Feature/mod adjoint

include/quda.h

@@ -863,6 +863,17 @@ extern "C" {
                                  with the underlying issue documented https://github.com/lattice/quda/issues/1315 */
   } QudaGaugeObservableParam;
 
+  typedef struct QudaFermMeasurements_s {
+    size_t struct_size; /**< Size of this struct in bytes.  Used to ensure that the host application and QUDA see the same struct*/
+    void **ppb;
+    void *ppb_t;
+    unsigned int meas_int;
+    void *meas_int_vec;
+    void *meas_list;
+    QudaBoolean take_meas;
+    QudaBoolean take_fwd_gflow;
+  } QudaFermMeasurements;
+
   typedef struct QudaGaugeSmearParam_s {
     size_t struct_size; /**< Size of this struct in bytes.  Used to ensure that the host application and QUDA see the same struct*/
     unsigned int n_steps; /**< The total number of smearing steps to perform. */
@@ -1094,6 +1105,14 @@ extern "C" {
    */
   QudaGaugeSmearParam newQudaGaugeSmearParam(void);
 
+  /**
+   * A new QudaFermMeasurements should always be initialized
+   * immediately after it's defined (and prior to explicitly setting
+   * its members) using this function.  Typical usage is as follows:
+   *
+   *   QudaFermMeasurements ferm_meas = newQudaFermMeasurements();
+   */
+  QudaFermMeasurements newQudaFermMeasurements(void);
   /**
    * A new QudaBLASParam should always be initialized immediately
    * after it's defined (and prior to explicitly setting its members)
@@ -1381,6 +1400,8 @@ extern "C" {
 
   void computeKSLinkQuda(void* fatlink, void* longlink, void* ulink, void* inlink,
                          double *path_coeff, QudaGaugeParam *param);
+
+  void set_act_path(double *act_path, int i);
 
   /**
    * Compute two-link field
@@ -1740,13 +1761,11 @@ extern "C" {
    * @param[in] h_in Input fermion field set
    * @param[in] inv_param Dirac/Laplacian and solver meta data
    * @param[in] smear_param Parameter struct that defines the computation parameters
-   * @param[in,out] obs_param Parameter struct that defines which
-   * observables we are making and the resulting observables.
+   * @param[in,out] ferm_meas Parameter struct that conducts measurements of fermionic observables
    * @param[in] nSpinors Number of spinors in the input and output fields
    */
-  void performAdjGFlowHier(void **h_out, void **h_in, QudaInvertParam *inv_param, QudaGaugeSmearParam *smear_param,
+  void performAdjGFlowHier(void **h_out, void **h_in, QudaInvertParam *inv_param, QudaGaugeSmearParam *smear_param, QudaFermMeasurements *ferm_meas,
                            size_t nSpinors);
-
   /**
    * @brief Calculates a variety of gauge-field observables.  If a
    * smeared gauge field is presently loaded (in gaugeSmeared) the

lib/check_params.h

@@ -1207,6 +1207,49 @@ void printQudaGaugeSmearParam(QudaGaugeSmearParam *param)
 #endif
 }
 
+#if defined INIT_PARAM
+QudaFermMeasurements newQudaFermMeasurements(void)
+{
+  QudaFermMeasurements ret;
+#elif defined CHECK_PARAM
+static void checkQudaFermMeasurements(QudaFermMeasurements *param)
+{
+#else
+void printQudaFermMeasurements(QudaFermMeasurements *param)
+{
+  printfQuda("QUDA Gauge Smear Parameters:\n");
+#endif
+
+#if defined CHECK_PARAM
+  if (param->struct_size != (size_t)INVALID_INT && param->struct_size != sizeof(*param))
+    errorQuda("Unexpected QudaGaugeSmearParam struct size %lu, expected %lu", param->struct_size, sizeof(*param));
+
+#else
+  P(struct_size, (size_t)INVALID_INT);
+#endif
+
+#ifdef INIT_PARAM
+  P(meas_int, 5);
+  P(meas_int_vec, nullptr);
+  P(take_meas, QUDA_BOOLEAN_FALSE);
+  P(take_fwd_gflow, QUDA_BOOLEAN_FALSE);
+  P(ppb, nullptr);
+  P(ppb_t, nullptr);
+  P(meas_list, nullptr);  
+#else
+  P(meas_int, (unsigned int)INVALID_INT);
+  P(take_meas, QUDA_BOOLEAN_INVALID);
+  P(take_fwd_gflow, QUDA_BOOLEAN_INVALID);
+#endif
+
+#ifdef INIT_PARAM
+  return ret;
+#endif
+}
+
+
+
+
 #if defined INIT_PARAM
 QudaBLASParam newQudaBLASParam(void)
 {