Patch to change the initial value of Maximum reduction from 0.0 to the first element in the input array

git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalReduce/trunk@74 7daa882c-dc44-4453-834e-278d26b18e6a

1 files changed, 126 insertions, 126 deletions

diff --git a/src/MaxFunctions.c b/src/MaxFunctions.c
index 7a8b4c0..9fd7d57 100644
--- a/src/MaxFunctions.c
+++ b/src/MaxFunctions.c
@@ -39,7 +39,7 @@ int LocalReduce_Max_BYTE(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -112,7 +112,7 @@ int LocalReduce_Max_BYTE(int i, int weight_on, const void * const weight, CCTK_I
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -125,7 +125,7 @@ int LocalReduce_Max_BYTE(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -152,7 +152,7 @@ int LocalReduce_Max_BYTE(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -165,7 +165,7 @@ int LocalReduce_Max_BYTE(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -193,7 +193,7 @@ int LocalReduce_Max_BYTE(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -207,7 +207,7 @@ int LocalReduce_Max_BYTE(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -235,7 +235,7 @@ int LocalReduce_Max_BYTE(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -249,7 +249,7 @@ int LocalReduce_Max_BYTE(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
     break;
     #endif
@@ -280,7 +280,7 @@ int LocalReduce_Max_INT(int i, int weight_on, const void * const weight, CCTK_IN
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -352,7 +352,7 @@ int LocalReduce_Max_INT(int i, int weight_on, const void * const weight, CCTK_IN
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -365,7 +365,7 @@ int LocalReduce_Max_INT(int i, int weight_on, const void * const weight, CCTK_IN
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -392,7 +392,7 @@ int LocalReduce_Max_INT(int i, int weight_on, const void * const weight, CCTK_IN
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -405,7 +405,7 @@ int LocalReduce_Max_INT(int i, int weight_on, const void * const weight, CCTK_IN
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -433,7 +433,7 @@ int LocalReduce_Max_INT(int i, int weight_on, const void * const weight, CCTK_IN
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -446,7 +446,7 @@ int LocalReduce_Max_INT(int i, int weight_on, const void * const weight, CCTK_IN
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -474,7 +474,7 @@ int LocalReduce_Max_INT(int i, int weight_on, const void * const weight, CCTK_IN
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -487,7 +487,7 @@ int LocalReduce_Max_INT(int i, int weight_on, const void * const weight, CCTK_IN
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
     break;
     #endif
@@ -518,7 +518,7 @@ int LocalReduce_Max_INT1(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -590,7 +590,7 @@ int LocalReduce_Max_INT1(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -603,7 +603,7 @@ int LocalReduce_Max_INT1(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -630,7 +630,7 @@ int LocalReduce_Max_INT1(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -643,7 +643,7 @@ int LocalReduce_Max_INT1(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -671,7 +671,7 @@ int LocalReduce_Max_INT1(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -685,7 +685,7 @@ int LocalReduce_Max_INT1(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -713,7 +713,7 @@ int LocalReduce_Max_INT1(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -727,7 +727,7 @@ int LocalReduce_Max_INT1(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
     break;
     #endif
@@ -759,7 +759,7 @@ int LocalReduce_Max_INT2(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -831,7 +831,7 @@ int LocalReduce_Max_INT2(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -844,7 +844,7 @@ int LocalReduce_Max_INT2(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -871,7 +871,7 @@ int LocalReduce_Max_INT2(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -884,7 +884,7 @@ int LocalReduce_Max_INT2(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -912,7 +912,7 @@ int LocalReduce_Max_INT2(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -926,7 +926,7 @@ int LocalReduce_Max_INT2(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -954,7 +954,7 @@ int LocalReduce_Max_INT2(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -968,7 +968,7 @@ int LocalReduce_Max_INT2(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
     break;
     #endif
@@ -1000,7 +1000,7 @@ int LocalReduce_Max_INT4(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -1072,7 +1072,7 @@ int LocalReduce_Max_INT4(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -1085,7 +1085,7 @@ int LocalReduce_Max_INT4(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -1112,7 +1112,7 @@ int LocalReduce_Max_INT4(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -1125,7 +1125,7 @@ int LocalReduce_Max_INT4(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -1153,7 +1153,7 @@ int LocalReduce_Max_INT4(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -1167,7 +1167,7 @@ int LocalReduce_Max_INT4(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -1195,7 +1195,7 @@ int LocalReduce_Max_INT4(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -1209,7 +1209,7 @@ int LocalReduce_Max_INT4(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
     break;
     #endif
@@ -1241,7 +1241,7 @@ int LocalReduce_Max_INT8(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -1313,7 +1313,7 @@ int LocalReduce_Max_INT8(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -1326,7 +1326,7 @@ int LocalReduce_Max_INT8(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -1353,7 +1353,7 @@ int LocalReduce_Max_INT8(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -1366,7 +1366,7 @@ int LocalReduce_Max_INT8(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -1394,7 +1394,7 @@ int LocalReduce_Max_INT8(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -1408,7 +1408,7 @@ int LocalReduce_Max_INT8(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -1436,7 +1436,7 @@ int LocalReduce_Max_INT8(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -1450,7 +1450,7 @@ int LocalReduce_Max_INT8(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
     break;
     #endif
@@ -1481,7 +1481,7 @@ int LocalReduce_Max_REAL(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -1553,7 +1553,7 @@ int LocalReduce_Max_REAL(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -1566,7 +1566,7 @@ int LocalReduce_Max_REAL(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -1593,7 +1593,7 @@ int LocalReduce_Max_REAL(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -1606,7 +1606,7 @@ int LocalReduce_Max_REAL(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -1634,7 +1634,7 @@ int LocalReduce_Max_REAL(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -1648,7 +1648,7 @@ int LocalReduce_Max_REAL(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -1676,7 +1676,7 @@ int LocalReduce_Max_REAL(int i, int weight_on, const void * const weight, CCTK_I
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -1690,7 +1690,7 @@ int LocalReduce_Max_REAL(int i, int weight_on, const void * const weight, CCTK_I
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
     break;
     #endif
@@ -1721,7 +1721,7 @@ int LocalReduce_Max_REAL4(int i, int weight_on, const void * const weight, CCTK_
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -1793,7 +1793,7 @@ int LocalReduce_Max_REAL4(int i, int weight_on, const void * const weight, CCTK_
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -1806,7 +1806,7 @@ int LocalReduce_Max_REAL4(int i, int weight_on, const void * const weight, CCTK_
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -1833,7 +1833,7 @@ int LocalReduce_Max_REAL4(int i, int weight_on, const void * const weight, CCTK_
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -1846,7 +1846,7 @@ int LocalReduce_Max_REAL4(int i, int weight_on, const void * const weight, CCTK_
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -1874,7 +1874,7 @@ int LocalReduce_Max_REAL4(int i, int weight_on, const void * const weight, CCTK_
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -1888,7 +1888,7 @@ int LocalReduce_Max_REAL4(int i, int weight_on, const void * const weight, CCTK_
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -1916,7 +1916,7 @@ int LocalReduce_Max_REAL4(int i, int weight_on, const void * const weight, CCTK_
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -1930,7 +1930,7 @@ int LocalReduce_Max_REAL4(int i, int weight_on, const void * const weight, CCTK_
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
     break;
     #endif
@@ -1962,7 +1962,7 @@ int LocalReduce_Max_REAL8(int i, int weight_on, const void * const weight, CCTK_
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -2034,7 +2034,7 @@ int LocalReduce_Max_REAL8(int i, int weight_on, const void * const weight, CCTK_
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -2047,7 +2047,7 @@ int LocalReduce_Max_REAL8(int i, int weight_on, const void * const weight, CCTK_
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -2074,7 +2074,7 @@ int LocalReduce_Max_REAL8(int i, int weight_on, const void * const weight, CCTK_
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -2087,7 +2087,7 @@ int LocalReduce_Max_REAL8(int i, int weight_on, const void * const weight, CCTK_
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -2115,7 +2115,7 @@ int LocalReduce_Max_REAL8(int i, int weight_on, const void * const weight, CCTK_
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -2129,7 +2129,7 @@ int LocalReduce_Max_REAL8(int i, int weight_on, const void * const weight, CCTK_
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -2157,7 +2157,7 @@ int LocalReduce_Max_REAL8(int i, int weight_on, const void * const weight, CCTK_
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -2172,7 +2172,7 @@ int LocalReduce_Max_REAL8(int i, int weight_on, const void * const weight, CCTK_
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
     break;
     #endif
@@ -2204,7 +2204,7 @@ int LocalReduce_Max_REAL16(int i, int weight_on, const void * const weight, CCTK
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -2276,7 +2276,7 @@ int LocalReduce_Max_REAL16(int i, int weight_on, const void * const weight, CCTK
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -2289,7 +2289,7 @@ int LocalReduce_Max_REAL16(int i, int weight_on, const void * const weight, CCTK
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -2316,7 +2316,7 @@ int LocalReduce_Max_REAL16(int i, int weight_on, const void * const weight, CCTK
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -2329,7 +2329,7 @@ int LocalReduce_Max_REAL16(int i, int weight_on, const void * const weight, CCTK
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -2357,7 +2357,7 @@ int LocalReduce_Max_REAL16(int i, int weight_on, const void * const weight, CCTK
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -2371,7 +2371,7 @@ int LocalReduce_Max_REAL16(int i, int weight_on, const void * const weight, CCTK
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
 
     break;
@@ -2399,7 +2399,7 @@ int LocalReduce_Max_REAL16(int i, int weight_on, const void * const weight, CCTK
 #endif
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = typed_vdata[0]; (num).Im = 0;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -2413,7 +2413,7 @@ int LocalReduce_Max_REAL16(int i, int weight_on, const void * const weight, CCTK
         inval = (out_type) typed_vdata[sum_indices];
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = typed_vdata[0];
 #define EXTRA_STEP(a, b)
     break;
     #endif
@@ -2455,7 +2455,7 @@ int LocalReduce_Max_COMPLEX(int i, int weight_on, const void * const weight, CCT
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -2517,7 +2517,7 @@ int LocalReduce_Max_COMPLEX(int i, int weight_on, const void * const weight, CCT
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -2541,7 +2541,7 @@ int LocalReduce_Max_COMPLEX(int i, int weight_on, const void * const weight, CCT
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
     break;
@@ -2558,7 +2558,7 @@ int LocalReduce_Max_COMPLEX(int i, int weight_on, const void * const weight, CCT
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -2582,7 +2582,7 @@ int LocalReduce_Max_COMPLEX(int i, int weight_on, const void * const weight, CCT
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
     break;
@@ -2600,7 +2600,7 @@ int LocalReduce_Max_COMPLEX(int i, int weight_on, const void * const weight, CCT
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -2625,7 +2625,7 @@ int LocalReduce_Max_COMPLEX(int i, int weight_on, const void * const weight, CCT
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
     break;
@@ -2643,7 +2643,7 @@ int LocalReduce_Max_COMPLEX(int i, int weight_on, const void * const weight, CCT
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -2668,7 +2668,7 @@ int LocalReduce_Max_COMPLEX(int i, int weight_on, const void * const weight, CCT
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
     break;
     #endif
@@ -2711,7 +2711,7 @@ int LocalReduce_Max_COMPLEX8(int i, int weight_on, const void * const weight, CC
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -2773,7 +2773,7 @@ int LocalReduce_Max_COMPLEX8(int i, int weight_on, const void * const weight, CC
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -2797,7 +2797,7 @@ int LocalReduce_Max_COMPLEX8(int i, int weight_on, const void * const weight, CC
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
     break;
@@ -2814,7 +2814,7 @@ int LocalReduce_Max_COMPLEX8(int i, int weight_on, const void * const weight, CC
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -2838,7 +2838,7 @@ int LocalReduce_Max_COMPLEX8(int i, int weight_on, const void * const weight, CC
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
     break;
@@ -2856,7 +2856,7 @@ int LocalReduce_Max_COMPLEX8(int i, int weight_on, const void * const weight, CC
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -2881,7 +2881,7 @@ int LocalReduce_Max_COMPLEX8(int i, int weight_on, const void * const weight, CC
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
     break;
@@ -2899,7 +2899,7 @@ int LocalReduce_Max_COMPLEX8(int i, int weight_on, const void * const weight, CC
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -2924,7 +2924,7 @@ int LocalReduce_Max_COMPLEX8(int i, int weight_on, const void * const weight, CC
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
     break;
     #endif
@@ -2967,7 +2967,7 @@ int LocalReduce_Max_COMPLEX16(int i, int weight_on, const void * const weight, C
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -3029,7 +3029,7 @@ int LocalReduce_Max_COMPLEX16(int i, int weight_on, const void * const weight, C
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -3053,7 +3053,7 @@ int LocalReduce_Max_COMPLEX16(int i, int weight_on, const void * const weight, C
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
     break;
@@ -3070,7 +3070,7 @@ int LocalReduce_Max_COMPLEX16(int i, int weight_on, const void * const weight, C
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -3094,7 +3094,7 @@ int LocalReduce_Max_COMPLEX16(int i, int weight_on, const void * const weight, C
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
     break;
@@ -3112,7 +3112,7 @@ int LocalReduce_Max_COMPLEX16(int i, int weight_on, const void * const weight, C
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -3137,7 +3137,7 @@ int LocalReduce_Max_COMPLEX16(int i, int weight_on, const void * const weight, C
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
     break;
@@ -3155,7 +3155,7 @@ int LocalReduce_Max_COMPLEX16(int i, int weight_on, const void * const weight, C
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -3180,7 +3180,7 @@ int LocalReduce_Max_COMPLEX16(int i, int weight_on, const void * const weight, C
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
     break;
     #endif
@@ -3223,7 +3223,7 @@ int LocalReduce_Max_COMPLEX32(int i, int weight_on, const void * const weight, C
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
   switch (output_number_type_codes[i])
@@ -3285,7 +3285,7 @@ int LocalReduce_Max_COMPLEX32(int i, int weight_on, const void * const weight, C
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_COMPLEX, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -3309,7 +3309,7 @@ int LocalReduce_Max_COMPLEX32(int i, int weight_on, const void * const weight, C
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
     break;
@@ -3326,7 +3326,7 @@ int LocalReduce_Max_COMPLEX32(int i, int weight_on, const void * const weight, C
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
 #undef REDUCTION_OPERATION
@@ -3350,7 +3350,7 @@ int LocalReduce_Max_COMPLEX32(int i, int weight_on, const void * const weight, C
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
     break;
@@ -3368,7 +3368,7 @@ int LocalReduce_Max_COMPLEX32(int i, int weight_on, const void * const weight, C
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -3393,7 +3393,7 @@ int LocalReduce_Max_COMPLEX32(int i, int weight_on, const void * const weight, C
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
 
     break;
@@ -3411,7 +3411,7 @@ int LocalReduce_Max_COMPLEX32(int i, int weight_on, const void * const weight, C
 
 #define REDUCTION_OPERATION(Max, scalar)   if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight) if(ABS(weight)>SMALL_NUMBER){if ((Max).Re < (scalar).Re) (Max).Re = (scalar).Re;};
-#define REDUCTION_INITIAL(num) (num).Re = 0.0; (num).Im = 0.0;
+#define REDUCTION_INITIAL(num) (num).Re = (typed_vdata[0]).Re; (num).Im = (typed_vdata[0]).Im;
 #define EXTRA_STEP(a, b)
 
       ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], weight_on, weight, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,input_array_dims,product)
@@ -3436,7 +3436,7 @@ int LocalReduce_Max_COMPLEX32(int i, int weight_on, const void * const weight, C
 
 #define  REDUCTION_OPERATION(Max, scalar)   if (Max < scalar) Max = scalar;
 #define WEIGHTED_REDUCTION_OPERATION(Max, scalar, weight)   if (ABS(weight)> SMALL_NUMBER) {if (Max < scalar) Max = scalar;} ;
-#define REDUCTION_INITIAL(num) num = 0;
+#define REDUCTION_INITIAL(num) num = (typed_vdata[0]).Re;
 #define EXTRA_STEP(a, b)
     break;
     #endif