Added link updating. Modified code to work with QUDA master branch

Author: Justin Foley

milc_interface/Makefile

@@ -11,7 +11,7 @@ TEST_OBJS = ../tests/test_util.o ../tests/staggered_dslash_reference.o ../tests/
 LOCAL_HDRS = include/milc_utilities.h include/milc_timer.h include/milc_inverter_utilities.h
 EXTERNAL_HDR = external_headers/quda_milc_interface.h
 
-INTERFACE_OBJS = milc_init_finalize_interface.o milc_utilities.o milc_timer.o
+INTERFACE_OBJS = milc_init_finalize_interface.o milc_utilities.o milc_timer.o milc_update_interface.o
 
 ifeq ($(strip $(BUILD_STAGGERED_DIRAC)),yes)
 	INTERFACE_OBJS += milc_staggered_inverter_interface.o milc_inverter_utilities.o
@@ -53,5 +53,5 @@ clean:
 	-rm -f *.o
 
 %.o: %.cpp
-	$(CXX) $(CXXFLAGS) ${INC} $< -c -o $@
+	$(CXX) -g $(CXXFLAGS) ${INC} $< -c -o $@
 #	$(CXX) $(CXXFLAGS) -DTIME_INTERFACE ${INC} $< -c -o $@

milc_interface/external_headers/quda_milc_interface.h

@@ -2,6 +2,7 @@
 #define _QUDA_MILC_INTERFACE_H
 
 #include <enum_quda.h>
+#include <quda.h>
 
 #ifdef __cplusplus
 extern "C" {
@@ -53,7 +54,7 @@ extern "C" {
   void qudaHisqParamsInit(QudaHisqParams_t hisq_params);
 
 
-  void qudaLoadFatLink(int precision, QudaFatLinkArgs_t fatlink_args, const double act_path_coeff[6], void* inlink, void* outlink);
+  void qudaLoadKSLink(int precision, QudaFatLinkArgs_t fatlink_args, const double act_path_coeff[6], void* inlink, void* fatlink, void* longlink);
 
 
   void qudaLoadUnitarizedLink(int precision, QudaFatLinkArgs_t fatlink_args, const double path_coeff[6], void* inlink, void* fatlink, void* ulink);
@@ -67,6 +68,7 @@ extern "C" {
       double target_relresid,
       const void* const milc_fatlink,
       const void* const milc_longlink,
+      const double tadpole,
       void* source,
       void* solution,
       double* const final_resid,
@@ -103,6 +105,7 @@ extern "C" {
       const double* target_relative_residual,
       const void* const milc_fatlink,
       const void* const milc_longlink,
+      const double tadpole,
       void* source,
       void** solutionArray, 
       double* const final_residual,
@@ -128,6 +131,22 @@ extern "C" {
       );
 
 
+  void qudaCloverMultishiftInvert(int external_precision, 
+      int quda_precision,
+      int num_offsets,
+      double* const offset,
+      double kappa,
+      QudaInvertArgs_t inv_args,
+      const double* target_residual,
+      const void* milc_link,
+      void* milc_clover, 
+      void* milc_clover_inv,
+      void* source,
+      void** solutionArray,
+      double* const final_residual, 
+      int* num_iters
+      );
+
   void qudaHisqForce(
       int precision,
       const double level2_coeff[6],
@@ -141,6 +160,15 @@ extern "C" {
       void* const milc_momentum);
 
 
+  void qudaAsqtadForce(
+      int precision,
+      const double act_path_coeff[6],
+      const void* const one_link_src[4],
+      const void* const naik_src[4],
+      const void* const link,
+      void* const milc_momentum);
+
+
   void qudaGaugeForce(int precision,
       int num_loop_types,
       double milc_loop_coeff[3],
@@ -157,7 +185,26 @@ extern "C" {
       void* const one_link_src[4], 
       void* const three_link_src[4]);
 
+  void qudaComputeOprod(int precision,
+                        int num_terms,
+                        double** coeff,
+                        void** quark_field,
+                        void* oprod[2]);
+
+
+  void qudaUpdateU(int precision, 
+                   double eps,
+                   void* momentum, 
+                   void* link);
+
+
+  void qudaCloverDerivative(void* out, void* gauge, void* oprod, 
+                            int mu, int nu, int precision, int parity, int conjugate);
+
+
+  void* qudaCreateExtendedGaugeField(void* gauge, int geometry, int precision);
 
+  void qudaDestroyGaugeField(void* gauge);
 
 #ifdef __cplusplus
 }

milc_interface/include/milc_inverter_utilities.h

@@ -44,6 +44,7 @@ namespace milc_interface {
       QudaPrecision cuda_prec,
       QudaPrecision cuda_prec_sloppy,
       QudaPrecision cuda_prec_precondition,
+      const double tadpole,
       QudaGaugeParam *gaugeParam);