diff options
-rw-r--r-- | src/Norm1Functions.c | 4 | ||||
-rw-r--r-- | src/NormInfFunctions.c | 4 |
2 files changed, 4 insertions, 4 deletions
diff --git a/src/Norm1Functions.c b/src/Norm1Functions.c index 5e4713f..8f2082d 100644 --- a/src/Norm1Functions.c +++ b/src/Norm1Functions.c @@ -3306,7 +3306,7 @@ int LocalReduce_Norm1_COMPLEX16(int i, int weight_on, const void * const weight, #ifdef CCTK_REAL_PRECISION_4 #undef REDUCTION_PREOP_CAST #define REDUCTION_PREOP_CAST(inval, typed_vdata,sum_indices, out_type) \ - inval = (CCTK_REAL4) (CCTK_Cmplx8Abs(typed_vdata[sum_indices])); + inval = (CCTK_REAL4) (CCTK_Cmplx16Abs(typed_vdata[sum_indices])); #elif CCTK_REAL_PRECISION_8 #undef REDUCTION_PREOP_CAST #define REDUCTION_PREOP_CAST(inval, typed_vdata,sum_indices, out_type) \ @@ -3314,7 +3314,7 @@ int LocalReduce_Norm1_COMPLEX16(int i, int weight_on, const void * const weight, #elif CCTK_REAL_PRECISION_16 #undef REDUCTION_PREOP_CAST #define REDUCTION_PREOP_CAST(inval, typed_vdata,sum_indices, out_type) \ - inval = (CCTK_REAL16) (CCTK_Cmplx32Abs(typed_vdata[sum_indices])); + inval = (CCTK_REAL16) (CCTK_Cmplx16Abs(typed_vdata[sum_indices])); #endif #define REDUCTION_OPERATION(Norm1, scalar) Norm1 = Norm1 + scalar; diff --git a/src/NormInfFunctions.c b/src/NormInfFunctions.c index dee219d..0c61247 100644 --- a/src/NormInfFunctions.c +++ b/src/NormInfFunctions.c @@ -3716,7 +3716,7 @@ int LocalReduce_NormInf_COMPLEX16(int i, int weight_on, const void * const weigh #ifdef CCTK_REAL_PRECISION_4 #undef REDUCTION_PREOP_CAST #define REDUCTION_PREOP_CAST(inval, typed_vdata,sum_indices, out_type) \ - inval = (CCTK_REAL4) (CCTK_Cmplx8Abs(typed_vdata[sum_indices])); + inval = (CCTK_REAL4) (CCTK_Cmplx16Abs(typed_vdata[sum_indices])); #elif CCTK_REAL_PRECISION_8 #undef REDUCTION_PREOP_CAST #define REDUCTION_PREOP_CAST(inval, typed_vdata,sum_indices, out_type) \ @@ -3724,7 +3724,7 @@ int LocalReduce_NormInf_COMPLEX16(int i, int weight_on, const void * const weigh #elif CCTK_REAL_PRECISION_16 #undef REDUCTION_PREOP_CAST #define REDUCTION_PREOP_CAST(inval, typed_vdata,sum_indices, out_type) \ - inval = (CCTK_REAL16) (CCTK_Cmplx32Abs(typed_vdata[sum_indices])); + inval = (CCTK_REAL16) (CCTK_Cmplx16Abs(typed_vdata[sum_indices])); #endif #define REDUCTION_OPERATION(NormInf, scalar) NormInf = MAX(NormInf , ABS(scalar)); |