From 8913a21ca51af0664ae90e4520349749c9b023a1 Mon Sep 17 00:00:00 2001
From: yye00 <yye00@7daa882c-dc44-4453-834e-278d26b18e6a>
Date: Mon, 30 Aug 2004 15:12:41 +0000
Subject: more conditioning

git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalReduce/trunk@36 7daa882c-dc44-4453-834e-278d26b18e6a
---
 src/ReductionNorm2.c | 78 ++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 78 insertions(+)

(limited to 'src/ReductionNorm2.c')

diff --git a/src/ReductionNorm2.c b/src/ReductionNorm2.c
index 97e8f4c..6a75ea6 100644
--- a/src/ReductionNorm2.c
+++ b/src/ReductionNorm2.c
@@ -1234,6 +1234,8 @@ static int ReductionL2 (int N_dims, int operator_handle,
     }
   }
 
+  num_points--;
+
   /* store the number of points in the paramater table */
   if ( Util_TableQueryNKeys(param_table_handle) != 0)
   {
@@ -1241,6 +1243,82 @@ static int ReductionL2 (int N_dims, int operator_handle,
     ierr = Util_TableSetInt(param_table_handle, 3,"mpi_operation");
     ierr = Util_TableSetInt(param_table_handle, 0,"perform_division");
   }
+  else
+  {
+    for (i = 0; i< M_output_numbers; i++)
+    {
+      switch (output_number_type_codes[i])
+      {
+       /* out values type switches*/
+        case CCTK_VARIABLE_CHAR:
+          *( (CCTK_CHAR *) output_numbers[i]) = *( (CCTK_CHAR *) output_numbers[i]) / num_points; 
+        break;
+        case CCTK_VARIABLE_INT:
+          *( (CCTK_INT *) output_numbers[i]) = *( (CCTK_INT *) output_numbers[i]) / num_points; 
+        break;
+        #ifdef CCTK_INT1
+        case CCTK_VARIABLE_INT1:
+          *( (CCTK_INT1 *) output_numbers[i]) = *( (CCTK_INT1 *) output_numbers[i]) / num_points; 
+        break;
+        #endif                                                              
+        #ifdef CCTK_INT2
+        case CCTK_VARIABLE_INT2:
+          *( (CCTK_INT2 *) output_numbers[i]) = *( (CCTK_INT2 *) output_numbers[i]) / num_points; 
+        break;
+        #endif
+        #ifdef CCTK_INT4
+        case CCTK_VARIABLE_INT4:
+          *( (CCTK_INT4 *) output_numbers[i]) = *( (CCTK_INT4 *) output_numbers[i]) / num_points; 
+        break;
+        #endif
+        #ifdef CCTK_INT8
+        case CCTK_VARIABLE_INT8:
+          *( (CCTK_INT8 *) output_numbers[i]) = *( (CCTK_INT8 *) output_numbers[i]) / num_points; 
+        break;                                                                
+        #endif
+        case CCTK_VARIABLE_REAL:
+          *( (CCTK_REAL *) output_numbers[i]) = *( (CCTK_REAL *) output_numbers[i]) / num_points; 
+        break;
+        #ifdef CCTK_REAL4
+        case CCTK_VARIABLE_REAL4:
+          *( (CCTK_REAL4 *) output_numbers[i]) = *( (CCTK_REAL4 *) output_numbers[i]) / num_points; 
+        break;
+        #endif
+        #ifdef CCTK_REAL8
+        case CCTK_VARIABLE_REAL8:
+          *( (CCTK_REAL8 *) output_numbers[i]) = *( (CCTK_REAL8 *) output_numbers[i]) / num_points; 
+        break;
+        #endif
+        #ifdef CCTK_REAL16
+        case CCTK_VARIABLE_REAL6:
+          *( (CCTK_REAL16 *) output_numbers[i]) = *( (CCTK_REAL16 *) output_numbers[i]) / num_points; 
+        break;
+        #endif
+        case CCTK_VARIABLE_COMPLEX:
+          (*( (CCTK_COMPLEX *) output_numbers[i])).Re = (*( (CCTK_COMPLEX *) output_numbers[i])).Re / num_points; 
+          (*( (CCTK_COMPLEX *) output_numbers[i])).Im = (*( (CCTK_COMPLEX *) output_numbers[i])).Im / num_points; 
+        break;
+        #ifdef CCTK_COMPLEX8
+        case CCTK_VARIABLE_COMPLEX8:
+          (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re / num_points; 
+          (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im / num_points; 
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX16
+        case CCTK_VARIABLE_COMPLEX16:
+          (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re / num_points; 
+          (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im / num_points; 
+        break;
+        #endif
+        #ifdef CCTK_COMPLEX32
+        case CCTK_VARIABLE_COMPLEX32:
+          (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re / num_points; 
+          (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im / num_points; 
+        break;
+        #endif
+      }
+    }
+  }
 
   /* free memory */
   free (iters_per_dim);
-- 
cgit v1.2.3