reduction avg

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

1 files changed, 954 insertions, 274 deletions

diff --git a/src/ReductionAvg.c b/src/ReductionAvg.c
index fb466a2..b844bb7 100644
--- a/src/ReductionAvg.c
+++ b/src/ReductionAvg.c
@@ -18,6 +18,11 @@ static const char *rcsid = "$Id$";
 
 CCTK_FILEVERSION(CCTDevelopment_LocalReduce_ReductionAvg_c);
 
+
+/* Define the reduction operations */
+
+
+
 /* local function prototypes */
 static int ReductionAvg (int N_dims, int operator_handle, 
                         int param_table_handle,   int N_input_arrays,
@@ -132,10 +137,8 @@ static int ReductionAvg (int N_dims, int operator_handle,
   int i, j, k, flag, product, num_points;
   int ierr;
   int * iters_per_dim;
-  void * data_pointer;
-  void * output_buffer;
-
-  output_buffer = (CCTK_REAL *) malloc (M_output_numbers*sizeof(CCTK_REAL));
+  void * data_pointer[10];
+  void * output_buffer[10];
 
   /* indices to hold the temp indices of size N_dims  and iteration indices*/
   int * indices;
@@ -160,6 +163,9 @@ static int ReductionAvg (int N_dims, int operator_handle,
   CCTK_INT mask_offset;
   CCTK_INT mask_time_level;
 
+
+
+
   /* set the number of points */
   num_points = 0;
 
@@ -235,322 +241,996 @@ static int ReductionAvg (int N_dims, int operator_handle,
     iters_per_dim [j]        = (int) ((input_array_max_subscripts[j] - input_array_min_subscripts[j])/input_array_strides[j]);
     max_iter *= iters_per_dim [j];
   }
-                                                                                
+
+#undef REDUCTION_OPERATION
+#undef REDUCTION_INITIAL
+#undef EXTRA_STEP
+
+#define REDUCTION_OPERATION(avg, scalar)   avg = avg + scalar;
+#define REDUCTION_INITIAL(num) num = 1.0;
+#define EXTRA_STEP(a, b) a = a/b;
+
   for (i = 0; i < N_input_arrays; i++)
   { 
-/*
+    /* Do the type matching */
+    switch (input_array_type_codes[i])
+    {
+    /* in values type switches*/
+      case CCTK_VARIABLE_CHAR:
       switch (output_number_type_codes[i])
       {
+        /* out values type switches*/
         case CCTK_VARIABLE_CHAR:
-        output_buffer = (CCTK_BYTE *) output_buffer;
-        output_buffer = (CCTK_BYTE *) malloc (M_output_numbers*sizeof(CCTK_BYTE));
+        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
         break;
         case CCTK_VARIABLE_INT:
-        output_buffer = (CCTK_INT *) output_buffer;
-        output_buffer = (CCTK_INT *) malloc (M_output_numbers*sizeof(CCTK_INT));
+        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
         break;
         #ifdef CCTK_INT1
         case CCTK_VARIABLE_INT1:
-        output_buffer = (CCTK_INT1 *) output_buffer;
-        output_buffer = (CCTK_INT1 *) malloc (M_output_numbers*sizeof(CCTK_INT1));
+        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
         break;
         #endif                                                              
         #ifdef CCTK_INT2
         case CCTK_VARIABLE_INT2:
-        output_buffer = (CCTK_INT2 *) output_buffer;
-        output_buffer = (CCTK_INT2 *) malloc (M_output_numbers*sizeof(CCTK_INT2));
+        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
         break;
         #endif
         #ifdef CCTK_INT4
         case CCTK_VARIABLE_INT4:
-        output_buffer = (CCTK_INT4 *) output_buffer;
-        output_buffer = (CCTK_INT4 *) malloc (M_output_numbers*sizeof(CCTK_INT4));
+        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
         break;
         #endif
         #ifdef CCTK_INT8
         case CCTK_VARIABLE_INT8:
-        output_buffer = (CCTK_INT8 *) output_buffer;
-        output_buffer = (CCTK_INT8 *) malloc (M_output_numbers*sizeof(CCTK_INT8));
+        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
         break;                                                                
         #endif
         case CCTK_VARIABLE_REAL:
-        #define type  CCTK_REAL
-        output_buffer = (CCTK_REAL *) output_buffer;
-        output_buffer = (CCTK_REAL *) malloc (M_output_numbers*sizeof(CCTK_REAL));
+        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
         break;
         #ifdef CCTK_REAL4
         case CCTK_VARIABLE_REAL4:
-        output_buffer = (CCTK_REAL4 *) output_buffer;
-        output_buffer = (CCTK_REAL4 *) malloc (M_output_numbers*sizeof(CCTK_REAL4));
+        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
         break;
         #endif
         #ifdef CCTK_REAL8
         case CCTK_VARIABLE_REAL8:
-        output_buffer = (CCTK_REAL8 *) output_buffer;
-        output_buffer = (CCTK_REAL8 *) malloc (M_output_numbers*sizeof(CCTK_REAL8));
+        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
         break;
         #endif
         #ifdef CCTK_REAL6
         case CCTK_VARIABLE_REAL6:
-        output_buffer = (CCTK_REAL6 *) output_buffer;
-        output_buffer = (CCTK_REAL6 *) malloc (M_output_numbers*sizeof(CCTK_REAL6));
+        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
         break;
         #endif
-        default:   
-        CCTK_WARN (1, "Average reduction: Unknown variable type");
-        return (-1);
-      }*/
-    if ( mask_on == 1)
-    {
-      if ( input_array_offsets[i] == 0)
+        case CCTK_VARIABLE_COMPLEX:
+/*        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_BYTE, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+      }
+      break;
+      case CCTK_VARIABLE_INT:
+      switch (output_number_type_codes[i])
       {
-        output_buffer = 0;
-        iter = 0;
-        sum_indices = 0;
-
-        while (iter < max_iter)
-        {
-          
-          ((CCTK_REAL *)output_buffer)[i] += ((CCTK_REAL* )data_pointer)[sum_indices];
-          num_points++;
-
-          printf("\nvalue is:%f",((CCTK_REAL* )data_pointer)[sum_indices]);
-
-          flag = 0;
-          for (k=0;k<N_dims;k++)
-          {
-            if( indices[k] < iters_per_dim[k]-1)
-            {
-              if (flag==1)
-              {
-                actual_indices[k] += input_array_strides[k-1];
-                indices[k]++;
-                iter++;
-                flag = 0;
-                break;
-              }
-              indices[k]++;
-              actual_indices[k] += input_array_strides[k];
-              iter++;
-              break;
-            }
-            else if (indices[k] == iters_per_dim[k]-1)
-            { 
-              indices[k] = 0;
-              actual_indices[k] = input_array_min_subscripts[k];
-              flag = 1;
-              continue;
-            }
-            else
-            {
-              CCTK_WARN(1,"indices out of bounds, this should not happen");
-              return -1;
-            }
-          }
-          sum_indices = actual_indices[0];
-          for (k=N_dims-1;k>0;k--)
-          {
-            product = 1;           
-            for (j=k-1;j>=0;j--)
-            {
-              product *= actual_iters_per_dim[j];
-            }
-            sum_indices += actual_indices[k]*product;
-          }          
-        }
+        /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        case CCTK_VARIABLE_INT:
+        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL6
+        case CCTK_VARIABLE_REAL6:
+        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+/*        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_INT, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
       }
-      else
+      break;
+      #ifdef CCTK_INT1
+      case CCTK_VARIABLE_INT1:
+      switch (output_number_type_codes[i])
       {
-        output_buffer = 0;
-        iter = 0;
-        sum_indices = 0;
-
-        while (iter < max_iter)
-        { 
-          /* prevent offset from giving segfaults */
-          if (sum_indices >= max_iter)
-          {
-            CCTK_WARN(1,"offsets and strides access unallocated memory");
-            return -1;
-          }
-
-          ((CCTK_REAL *)output_buffer)[i] += ((CCTK_REAL* )data_pointer)[sum_indices];
-          num_points++;
-
-          printf("\nvalue is:%f",((CCTK_REAL* )data_pointer)[sum_indices]);
-
-          flag = 0;
-          for (k=0;k<N_dims;k++)
-          {
-            if( indices[k] < iters_per_dim[k]-1)
-            {
-              if (flag==1)
-              {
-                actual_indices[k] += input_array_strides[k-1];
-                indices[k]++;
-                iter++;
-                flag = 0;
-                break;
-              }
-              indices[k]++;
-              actual_indices[k] += input_array_strides[k];
-              iter++;
-              break;
-            }
-            else if (indices[k] == iters_per_dim[k]-1)
-            { 
-              indices[k] = 0;
-              actual_indices[k] = input_array_min_subscripts[k];
-              flag = 1;
-              continue;
-            }
-            else
-            {
-              CCTK_WARN(1,"indices out of bounds, this should not happen");
-              return -1;
-            }
-          }
-          sum_indices = actual_indices[0]+input_array_offsets[i];
-          for (k=N_dims-1;k>0;k--)
-          {
-            product = 1;           
-            for (j=k-1;j>=0;j--)
-            {
-              product *= actual_iters_per_dim[j];
-            }
-            sum_indices += actual_indices[k]*product;
-          }          
-        }
-      }             
-    }
-    else if (mask_on == 0)
-    {
-      if ( input_array_offsets[i] == 0)
+        /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        case CCTK_VARIABLE_INT:
+        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL6
+        case CCTK_VARIABLE_REAL6:
+        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+/*        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_INT1, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+      }
+      break;
+      #endif                                                              
+      #ifdef CCTK_INT2
+      case CCTK_VARIABLE_INT2:
+      switch (output_number_type_codes[i])
       {
-        output_buffer = 0;
-        iter = 0;
-        sum_indices = 0;
-
-        while (iter < max_iter)
-        {
-          
-          ((CCTK_REAL *)output_buffer)[i] += ((CCTK_REAL* )data_pointer)[sum_indices];
-          num_points++;
-
-          printf("\nvalue is:%f",((CCTK_REAL* )data_pointer)[sum_indices]);
-
-          flag = 0;
-          for (k=0;k<N_dims;k++)
-          {
-            if( indices[k] < iters_per_dim[k]-1)
-            {
-              if (flag==1)
-              {
-                actual_indices[k] += input_array_strides[k-1];
-                indices[k]++;
-                iter++;
-                flag = 0;
-                break;
-              }
-              indices[k]++;
-              actual_indices[k] += input_array_strides[k];
-              iter++;
-              break;
-            }
-            else if (indices[k] == iters_per_dim[k]-1)
-            { 
-              indices[k] = 0;
-              actual_indices[k] = input_array_min_subscripts[k];
-              flag = 1;
-              continue;
-            }
-            else
-            {
-              CCTK_WARN(1,"indices out of bounds, this should not happen");
-              return -1;
-            }
-          }
-          sum_indices = actual_indices[0];
-          for (k=N_dims-1;k>0;k--)
-          {
-            product = 1;           
-            for (j=k-1;j>=0;j--)
-            {
-              product *= actual_iters_per_dim[j];
-            }
-            sum_indices += actual_indices[k]*product;
-          }          
-        }
+        /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        case CCTK_VARIABLE_INT:
+        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL6
+        case CCTK_VARIABLE_REAL6:
+        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+/*        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_INT2, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
       }
-      else
+      break;
+      #endif
+      #ifdef CCTK_INT4
+      case CCTK_VARIABLE_INT4:
+      switch (output_number_type_codes[i])
       {
-        output_buffer = 0;
-        iter = 0;
-        sum_indices = 0;
-
-        while (iter < max_iter)
-        { 
-          /* prevent offset from giving segfaults */
-          if (sum_indices >= max_iter)
-          {
-            CCTK_WARN(1,"offsets and strides access unallocated memory");
-            return -1;
-          }
-
-          ((CCTK_REAL *)output_buffer)[i] += ((CCTK_REAL* )data_pointer)[sum_indices];
-          num_points++;
-
-          printf("\nvalue is:%f",((CCTK_REAL* )data_pointer)[sum_indices]);
-
-          flag = 0;
-          for (k=0;k<N_dims;k++)
-          {
-            if( indices[k] < iters_per_dim[k]-1)
-            {
-              if (flag==1)
-              {
-                actual_indices[k] += input_array_strides[k-1];
-                indices[k]++;
-                iter++;
-                flag = 0;
-                break;
-              }
-              indices[k]++;
-              actual_indices[k] += input_array_strides[k];
-              iter++;
-              break;
-            }
-            else if (indices[k] == iters_per_dim[k]-1)
-            { 
-              indices[k] = 0;
-              actual_indices[k] = input_array_min_subscripts[k];
-              flag = 1;
-              continue;
-            }
-            else
-            {
-              CCTK_WARN(1,"indices out of bounds, this should not happen");
-              return -1;
-            }
-          }
-          sum_indices = actual_indices[0]+input_array_offsets[i];
-          for (k=N_dims-1;k>0;k--)
-          {
-            product = 1;           
-            for (j=k-1;j>=0;j--)
-            {
-              product *= actual_iters_per_dim[j];
-            }
-            sum_indices += actual_indices[k]*product;
-          }          
-        }
-      }         
-    }  
-    else
-    {
-      CCTK_WARN(1, "mask_on is not set to a valid value");
+        /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        case CCTK_VARIABLE_INT:
+        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL6
+        case CCTK_VARIABLE_REAL6:
+        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+/*        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_INT4, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+      }
+      break;
+      #endif
+      #ifdef CCTK_INT8
+      case CCTK_VARIABLE_INT8:
+      switch (output_number_type_codes[i])
+      {
+        /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        case CCTK_VARIABLE_INT:
+        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL6
+        case CCTK_VARIABLE_REAL6:
+        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+/*        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_INT8, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+      }
+      break;                                                                
+      #endif
+      case CCTK_VARIABLE_REAL:
+      switch (output_number_type_codes[i])
+      {
+        /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        case CCTK_VARIABLE_INT:
+        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL6
+        case CCTK_VARIABLE_REAL6:
+        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+/*        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_REAL, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+      }
+      break;
+      #ifdef CCTK_REAL4
+      case CCTK_VARIABLE_REAL4:
+      switch (output_number_type_codes[i])
+      {
+        /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        case CCTK_VARIABLE_INT:
+        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL6
+        case CCTK_VARIABLE_REAL6:
+        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+/*        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_REAL4, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+      }
+      break;
+      #endif
+      #ifdef CCTK_REAL8
+      case CCTK_VARIABLE_REAL8:
+      switch (output_number_type_codes[i])
+      {
+        /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        case CCTK_VARIABLE_INT:
+        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;	
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL6
+        case CCTK_VARIABLE_REAL6:
+        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+/*        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_REAL8, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+      }
+      break;
+      #endif
+      #ifdef CCTK_REAL6
+      case CCTK_VARIABLE_REAL6:
+      switch (output_number_type_codes[i])
+      {
+        /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        case CCTK_VARIABLE_INT:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;	
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL6
+        case CCTK_VARIABLE_REAL6:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_REAL16, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+      }
+      break;
+      #endif  
+      case CCTK_VARIABLE_COMPLEX:
+      switch (output_number_type_codes[i])
+      {
+        /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+/*        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        case CCTK_VARIABLE_INT:
+/*        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+/*        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;	
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+/*        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+/*        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+/*        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+/*        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+/*        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+/*        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #endif
+        #ifdef CCTK_REAL6
+        case CCTK_VARIABLE_REAL6:
+/*        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+/*        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)*/
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+      }
+      break;
+      #ifdef CCTK_COMPLEX8
+      case CCTK_VARIABLE_COMPLEX8:
+      switch (output_number_type_codes[i])
+      {
+        /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        case CCTK_VARIABLE_INT:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;	
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL6
+        case CCTK_VARIABLE_REAL6:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX8, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+      }
+      break;
+      #endif
+      #ifdef CCTK_COMPLEX16
+      case CCTK_VARIABLE_COMPLEX16:
+      switch (output_number_type_codes[i])
+      {
+        /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        case CCTK_VARIABLE_INT:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;	
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL6
+        case CCTK_VARIABLE_REAL6:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX16, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+      }
+      break;
+      #endif
+      #ifdef CCTK_COMPLEX32
+      case CCTK_VARIABLE_COMPLEX32:
+      switch (output_number_type_codes[i])
+      {
+        /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_BYTE, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        case CCTK_VARIABLE_INT:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_INT, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_INT1, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;	
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_INT2, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_INT4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_INT8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_REAL, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_REAL4, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_REAL8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_REAL6
+        case CCTK_VARIABLE_REAL6:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_REAL6, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_COMPLEX, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_COMPLEX8, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_COMPLEX16, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+        ITERATE_ON_ARRAY(i,CCTK_COMPLEX32, input_arrays[i], CCTK_COMPLEX32, output_numbers[i], mask_on, input_array_offsets[i], indices, sum_indices, max_iter, iter, flag, actual_indices,input_array_strides, input_array_min_subscripts,product)
+        break;
+        #endif
+      }
+      break;
+      #endif
     }
   }
+                                                                                
   return (0);
 }