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; } }