Change from using memset() to using a new function

IDAnalyticBH_zero_CCTK_REAL_array()  when setting a grid array to all zeros.
This is necessary because  memset() gives all bits 0, but we want
floating point (CCTK_REAL) 0.0, which might be a different bit pattern.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@126 6a3ddf76-46e1-4315-99d9-bc56cac1ef84

1 files changed, 20 insertions, 18 deletions

diff --git a/src/Misner_standard.c b/src/Misner_standard.c
index adde16b..cb14618 100644
--- a/src/Misner_standard.c
+++ b/src/Misner_standard.c
@@ -21,6 +21,8 @@
 #include "cctk_Arguments.h"
 #include "cctk_Parameters.h"
 
+#include "IDAnalyticBH.h"
+
 static const char *rcsid = "$Header$";
 
 CCTK_FILEVERSION(CactusEinstein_IDAnalyticBH_Misner_standard_c)
@@ -111,18 +113,18 @@ void Misner_standard(CCTK_ARGUMENTS)
    */
   if (make_conformal_derivs)
   {
-    memset (psix, 0, npoints * sizeof (CCTK_REAL));
-    memset (psiy, 0, npoints * sizeof (CCTK_REAL));
-    memset (psiz, 0, npoints * sizeof (CCTK_REAL));
+    IDAnalyticBH_zero_CCTK_REAL_array(npoints, psix);
+    IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiy);
+    IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiz);
 
     if(*conformal_state > 2)
     {
-      memset (psixx, 0, npoints * sizeof (CCTK_REAL));
-      memset (psixy, 0, npoints * sizeof (CCTK_REAL));
-      memset (psixz, 0, npoints * sizeof (CCTK_REAL));
-      memset (psiyy, 0, npoints * sizeof (CCTK_REAL));
-      memset (psiyz, 0, npoints * sizeof (CCTK_REAL));
-      memset (psizz, 0, npoints * sizeof (CCTK_REAL));
+      IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixx);
+      IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixy);
+      IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixz);
+      IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiyy);
+      IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiyz);
+      IDAnalyticBH_zero_CCTK_REAL_array(npoints, psizz);
     }
   }
 
@@ -259,19 +261,19 @@ void Misner_standard(CCTK_ARGUMENTS)
       gzz[i] = gxx[i];
     }
   }
-  memset (gxy, 0, npoints * sizeof (CCTK_REAL));
-  memset (gxz, 0, npoints * sizeof (CCTK_REAL));
-  memset (gyz, 0, npoints * sizeof (CCTK_REAL));
+  IDAnalyticBH_zero_CCTK_REAL_array(npoints, gxy);
+  IDAnalyticBH_zero_CCTK_REAL_array(npoints, gxz);
+  IDAnalyticBH_zero_CCTK_REAL_array(npoints, gyz);
 
   /*     Time-symmetric data
    *     -------------------
    */
-  memset (kxx, 0, npoints * sizeof (CCTK_REAL));
-  memset (kyy, 0, npoints * sizeof (CCTK_REAL));
-  memset (kzz, 0, npoints * sizeof (CCTK_REAL));
-  memset (kxy, 0, npoints * sizeof (CCTK_REAL));
-  memset (kxz, 0, npoints * sizeof (CCTK_REAL));
-  memset (kyz, 0, npoints * sizeof (CCTK_REAL));
+  IDAnalyticBH_zero_CCTK_REAL_array(npoints, kxx);
+  IDAnalyticBH_zero_CCTK_REAL_array(npoints, kxy);
+  IDAnalyticBH_zero_CCTK_REAL_array(npoints, kxz);
+  IDAnalyticBH_zero_CCTK_REAL_array(npoints, kyy);
+  IDAnalyticBH_zero_CCTK_REAL_array(npoints, kyz);
+  IDAnalyticBH_zero_CCTK_REAL_array(npoints, kzz);
 
   if (csch)
   {