Changed the scale factor used for the long links.

Only matters for reconstruct 13 and reconstruct 9.

Author: Justin Foley

milc_interface/milc_inverter_utilities.cpp

@@ -151,7 +151,8 @@ namespace milc_interface {
     gaugeParam->t_boundary = QUDA_PERIODIC_T; // anti-periodic boundary conditions are built into the gauge field
     gaugeParam->gauge_order = QUDA_MILC_GAUGE_ORDER; 
     gaugeParam->ga_pad = dim[0]*dim[1]*dim[2]/2;
-    gaugeParam->scale = -1.0/(24.0*gaugeParam->tadpole_coeff*gaugeParam->tadpole_coeff);
+    //gaugeParam->scale = -1.0/(24.0*gaugeParam->tadpole_coeff*gaugeParam->tadpole_coeff);
+    gaugeParam->scale = -1.0/(24.0);
 
 
     // preconditioning...

milc_interface/milc_staggered_inverter_interface.cpp

@@ -183,15 +183,16 @@ void qudaMultishiftInvert(int external_precision,
     const int long_pad = 3*fat_pad;
     gaugeParam.type = QUDA_THREE_LINKS;
     gaugeParam.ga_pad = long_pad; // don't know if this will work
+    gaugeParam.reconstruct = gaugeParam.reconstruct_sloppy = QUDA_RECONSTRUCT_13;
     loadGaugeQuda(const_cast<void*>(longlink), &gaugeParam);
 #else // single-gpu code
     gaugeParam.type = QUDA_GENERAL_LINKS;
     gaugeParam.reconstruct = gaugeParam.reconstruct_sloppy = QUDA_RECONSTRUCT_NO;
     loadGaugeQuda(const_cast<void*>(fatlink), &gaugeParam);
 
     gaugeParam.type = QUDA_ASQTAD_LONG_LINKS;
-    gaugeParam.reconstruct = QUDA_RECONSTRUCT_NO;
-    gaugeParam.reconstruct_sloppy = QUDA_RECONSTRUCT_NO;
+    gaugeParam.reconstruct = QUDA_RECONSTRUCT_13;
+    gaugeParam.reconstruct_sloppy = QUDA_RECONSTRUCT_13;
     loadGaugeQuda(const_cast<void*>(longlink), &gaugeParam);
 #endif
     milc_interface::invalidate_quda_gauge = false;
@@ -346,15 +347,16 @@ void qudaInvert(int external_precision,
 
     gaugeParam.type = QUDA_THREE_LINKS;
     gaugeParam.ga_pad = long_pad; 
+    gaugeParam.reconstruct = gaugeParam.reconstruct_sloppy = QUDA_RECONSTRUCT_13;
     loadGaugeQuda(const_cast<void*>(longlink), &gaugeParam);
 #else // single-gpu code
     gaugeParam.type = QUDA_GENERAL_LINKS;
     gaugeParam.reconstruct = gaugeParam.reconstruct_sloppy = QUDA_RECONSTRUCT_NO;
     loadGaugeQuda(const_cast<void*>(fatlink), &gaugeParam);
 
     gaugeParam.type = QUDA_THREE_LINKS;
-    gaugeParam.reconstruct = QUDA_RECONSTRUCT_NO;
-    gaugeParam.reconstruct_sloppy = QUDA_RECONSTRUCT_NO;
+    gaugeParam.reconstruct = QUDA_RECONSTRUCT_13;
+    gaugeParam.reconstruct_sloppy = QUDA_RECONSTRUCT_13;
     loadGaugeQuda(const_cast<void*>(longlink), &gaugeParam);
 #endif
     milc_interface::invalidate_quda_gauge = false;