Handle all complex variable datatypes (not just the standard CCTK_COMPLEX).

This closes PR CactusPUGH/1475.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusPUGH/PUGHReduce/trunk@36 d60812e6-3970-4df4-986e-c251b06effeb

1 files changed, 28 insertions, 1 deletions

diff --git a/src/ReductionNormInf.c b/src/ReductionNormInf.c
index c304477..363ca5b 100644
--- a/src/ReductionNormInf.c
+++ b/src/ReductionNormInf.c
@@ -217,7 +217,7 @@ static int ReductionNormInf (const cGH *GH,
 {
   int i, total_outvals;
 #ifdef CCTK_MPI
-  pGH *pughGH;
+  const pGH *pughGH;
   CCTK_REAL *local_outvals;
 #endif
 
@@ -322,6 +322,33 @@ static int ReductionNormInf (const cGH *GH,
                        CCTK_REAL, outvals, num_outvals, total_outvals);
         break;
 
+#ifdef CCTK_REAL4
+#undef  REDUCTION_OPERATION
+#define REDUCTION_OPERATION(norm_inf, scalar)   norm_inf = MAX(norm_inf, CCTK_Cmplx8Abs(scalar))
+        ITERATE_ARRAY (CCTK_COMPLEX8, num_dims, inarrays[i],
+                       from, to, iterator, points_per_dim,
+                       CCTK_REAL4, outvals, num_outvals, total_outvals);
+        break;
+#endif
+
+#ifdef CCTK_REAL8
+#undef  REDUCTION_OPERATION
+#define REDUCTION_OPERATION(norm_inf, scalar)   norm_inf = MAX(norm_inf, CCTK_Cmplx16Abs(scalar))
+        ITERATE_ARRAY (CCTK_COMPLEX16, num_dims, inarrays[i],
+                       from, to, iterator, points_per_dim,
+                       CCTK_REAL8, outvals, num_outvals, total_outvals);
+        break;
+#endif
+
+#ifdef CCTK_REAL16
+#undef  REDUCTION_OPERATION
+#define REDUCTION_OPERATION(norm_inf, scalar)   norm_inf = MAX(norm_inf, CCTK_Cmplx32Abs(scalar))
+        ITERATE_ARRAY (CCTK_COMPLEX32, num_dims, inarrays[i],
+                       from, to, iterator, points_per_dim,
+                       CCTK_REAL16, outvals, num_outvals, total_outvals);
+        break;
+#endif
+
       default:
         CCTK_WARN (1, "PUGH_ReductionNormInf: Unknown variable type");
         return (-1);