[Commits] [svn:einsteintoolkit] HydroBase/trunk/src/ (Rev. 50)

[knarf at cct.lsu.edu]

User: knarf
Date: 2012/05/08 12:30 PM

Modified:
 /trunk/src/
  Initialisation.c

Log:
 initialize w_lorentz as well as the other declared quantities - especially as the velocities

File Changes:

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

File [modified]: Initialisation.c
Delta lines: +18 -15
===================================================================
--- trunk/src/Initialisation.c	2012-05-07 05:53:46 UTC (rev 49)
+++ trunk/src/Initialisation.c	2012-05-08 17:30:37 UTC (rev 50)
@@ -13,11 +13,12 @@
   
 #pragma omp parallel for
   for (int i=0; i<np; ++i) {
-    rho[i     ] = 0.0;
-    vel[i     ] = 0.0;
-    vel[i+  np] = 0.0;
-    vel[i+2*np] = 0.0;
-    eps[i     ] = 0.0;
+    rho[i     ]  = 0.0;
+    vel[i     ]  = 0.0;
+    vel[i+  np]  = 0.0;
+    vel[i+2*np]  = 0.0;
+    w_lorentz[i] = 1.0;
+    eps[i     ]  = 0.0;
   }
   
   if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::temperature") >= 1) {
@@ -45,11 +46,12 @@
     if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::rho") >= 2) {
 #pragma omp parallel for
       for (int i=0; i<np; ++i) {
-        rho_p[i     ] = 0.0;
-        vel_p[i     ] = 0.0;
-        vel_p[i+  np] = 0.0;
-        vel_p[i+2*np] = 0.0;
-        eps_p[i     ] = 0.0;
+        rho_p[i      ] = 0.0;
+        vel_p[i      ] = 0.0;
+        vel_p[i+  np ] = 0.0;
+        vel_p[i+2*np ] = 0.0;
+        w_lorentz_p[i] = 1.0;
+        eps_p[i      ] = 0.0;
       }
     }
 
@@ -70,11 +72,12 @@
     if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::rho") >= 3) {
 #pragma omp parallel for
       for (int i=0; i<np; ++i) {
-        rho_p_p[i     ] = 0.0;
-        vel_p_p[i     ] = 0.0;
-        vel_p_p[i+  np] = 0.0;
-        vel_p_p[i+2*np] = 0.0;
-        eps_p_p[i     ] = 0.0;
+        rho_p_p[i      ] = 0.0;
+        vel_p_p[i      ] = 0.0;
+        vel_p_p[i+  np ] = 0.0;
+        vel_p_p[i+2*np ] = 0.0;
+        w_lorentz_p_p[i] = 1.0;
+        eps_p_p[i      ] = 0.0;
       }
     }