diff options
Diffstat (limited to 'src/Norm2Functions.c')
-rw-r--r-- | src/Norm2Functions.c | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/Norm2Functions.c b/src/Norm2Functions.c index deb863f..bcbd458 100644 --- a/src/Norm2Functions.c +++ b/src/Norm2Functions.c @@ -2814,15 +2814,15 @@ int LocalReduce_Norm2_COMPLEX(int i, int weight_on, const void * const weight, C #ifdef CCTK_REAL_PRECISION_4 #undef REDUCTION_PREOP_CAST #define REDUCTION_PREOP_CAST(inval, typed_vdata,sum_indices, out_type) \ - inval = (CCTK_REAL4) ((typed_vdata[sum_indices]).Re); + inval = (CCTK_REAL4) CCTK_CmplxAbs(typed_vdata[sum_indices]); #elif CCTK_REAL_PRECISION_8 #undef REDUCTION_PREOP_CAST #define REDUCTION_PREOP_CAST(inval, typed_vdata,sum_indices, out_type) \ - inval = (CCTK_REAL8) ((typed_vdata[sum_indices]).Re); + inval = (CCTK_REAL8) CCTK_CmplxAbs(typed_vdata[sum_indices]); #elif CCTK_REAL_PRECISION_16 #undef REDUCTION_PREOP_CAST #define REDUCTION_PREOP_CAST(inval, typed_vdata,sum_indices, out_type) \ - inval = (CCTK_REAL16) ((typed_vdata[sum_indices]).Re); + inval = (CCTK_REAL16) CCTK_CmplxAbs(typed_vdata[sum_indices]); #endif #define REDUCTION_OPERATION(Norm2, scalar) Norm2 = Norm2 + SQR(scalar); |