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, 17 insertions, 15 deletions

diff --git a/src/BrillLindquist.c b/src/BrillLindquist.c
index f0ecacd..74b56e6 100644
--- a/src/BrillLindquist.c
+++ b/src/BrillLindquist.c
@@ -16,6 +16,8 @@
 #include "cctk_Arguments.h"
 #include "cctk_Parameters.h"
 
+#include "IDAnalyticBH.h"
+
 static const char *rcsid = "$Header$";
 CCTK_FILEVERSION(CactusEinstein_IDAnalyticBH_BrillLindquist_c)
 
@@ -102,18 +104,18 @@ void BrillLindquist(CCTK_ARGUMENTS)
   
   if (make_conformal_derivs == 1)
   {
-    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);
     }
   }
 
@@ -231,10 +233,10 @@ void BrillLindquist(CCTK_ARGUMENTS)
   /*     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);
 }