diff options
-rw-r--r-- | src/Initialisation.c | 33 |
1 files changed, 18 insertions, 15 deletions
diff --git a/src/Initialisation.c b/src/Initialisation.c index 1149897..4ed861f 100644 --- a/src/Initialisation.c +++ b/src/Initialisation.c @@ -13,11 +13,12 @@ void HydroBase_Zero (CCTK_ARGUMENTS) #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 @@ void HydroBase_Zero (CCTK_ARGUMENTS) 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 @@ void HydroBase_Zero (CCTK_ARGUMENTS) 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; } } |