[Commits] [svn:einsteintoolkit] TOVSolver/trunk/src/ (Rev. 113)

knarf at cct.lsu.edu knarf at cct.lsu.edu
Sun Jul 11 22:31:00 CDT 2010


User: knarf
Date: 2010/07/11 10:31 PM

Modified:
 /trunk/src/
  tov.c

Log:
 Clear old matter initial data if not specified to be used (even if 'clear initial data' was set to false

File Changes:

Directory: /trunk/src/
======================

File [modified]: tov.c
Delta lines: +16 -16
===================================================================
--- trunk/src/tov.c	2010-06-17 19:59:37 UTC (rev 112)
+++ trunk/src/tov.c	2010-07-12 03:31:00 UTC (rev 113)
@@ -604,23 +604,23 @@
         }
       }
     }
-
-    if (!TOV_Use_Old_Matter_Initial_Data)
-    {
-      TOV_C_fill(rho,        LSH_MAX_I+1, 0.0);
-      TOV_C_fill(dens,       LSH_MAX_I+1, 0.0);
-      TOV_C_fill(eps,        LSH_MAX_I+1, 0.0);
-      TOV_C_fill(press,      LSH_MAX_I+1, 0.0);
-      TOV_C_fill(tau,        LSH_MAX_I+1, 0.0);
-      TOV_C_fill(w_lorentz,  LSH_MAX_I+1, 0.0);
-      TOV_C_fill(sx,         LSH_MAX_I+1, 0.0);
-      TOV_C_fill(sy,         LSH_MAX_I+1, 0.0);
-      TOV_C_fill(sz,         LSH_MAX_I+1, 0.0);
-      TOV_C_fill(velx,       LSH_MAX_I+1, 0.0);
-      TOV_C_fill(vely,       LSH_MAX_I+1, 0.0);
-      TOV_C_fill(velz,       LSH_MAX_I+1, 0.0);
-    }
   }
+  if (!TOV_Use_Old_Matter_Initial_Data)
+  {
+    CCTK_INFO("Not using old matter initial data");
+    TOV_C_fill(rho,        LSH_MAX_I+1, 0.0);
+    TOV_C_fill(dens,       LSH_MAX_I+1, 0.0);
+    TOV_C_fill(eps,        LSH_MAX_I+1, 0.0);
+    TOV_C_fill(press,      LSH_MAX_I+1, 0.0);
+    TOV_C_fill(tau,        LSH_MAX_I+1, 0.0);
+    TOV_C_fill(w_lorentz,  LSH_MAX_I+1, 0.0);
+    TOV_C_fill(sx,         LSH_MAX_I+1, 0.0);
+    TOV_C_fill(sy,         LSH_MAX_I+1, 0.0);
+    TOV_C_fill(sz,         LSH_MAX_I+1, 0.0);
+    TOV_C_fill(velx,       LSH_MAX_I+1, 0.0);
+    TOV_C_fill(vely,       LSH_MAX_I+1, 0.0);
+    TOV_C_fill(velz,       LSH_MAX_I+1, 0.0);
+  }
   /* use the fast interpolation? only useful for testing this */
   if (TOV_Fast_Interpolation == 0)
       CCTK_INFO("Interpolating the slow way.");



More information about the Commits mailing list