From a6877ae3a34dd99c8ae749ceda73cafa984e27d6 Mon Sep 17 00:00:00 2001 From: yye00 Date: Thu, 23 Feb 2006 21:48:01 +0000 Subject: 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 --- src/MaxFunctions.c | 252 ++++++++++++++++++++++++++--------------------------- 1 file 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 -- cgit v1.2.3