diff options
Diffstat (limited to 'src/Initialisation.c')
| -rw-r--r-- | src/Initialisation.c | 45 |
1 files changed, 44 insertions, 1 deletions
diff --git a/src/Initialisation.c b/src/Initialisation.c index 3391fde..7e79528 100644 --- a/src/Initialisation.c +++ b/src/Initialisation.c @@ -20,6 +20,20 @@ void HydroBase_Zero (CCTK_ARGUMENTS) eps[i ] = 0.0; } + if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::temperature") >= 1) { +#pragma omp parallel for + for (int i=0; i<np; ++i) { + temperature[i ] = 0.0; + } + } + + if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::entropy") >= 1) { +#pragma omp parallel for + for (int i=0; i<np; ++i) { + entropy[i ] = 0.0; + } + } + if (CCTK_EQUALS (initial_data_setup_method, "init_some_levels") || CCTK_EQUALS (initial_data_setup_method, "init_single_levels")) { @@ -38,7 +52,21 @@ void HydroBase_Zero (CCTK_ARGUMENTS) eps_p[i ] = 0.0; } } - + + if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::temperature") >= 2) { +#pragma omp parallel for + for (int i=0; i<np; ++i) { + temperature_p[i ] = 0.0; + } + } + + if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::entropy") >= 2) { +#pragma omp parallel for + for (int i=0; i<np; ++i) { + entropy_p[i ] = 0.0; + } + } + if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::rho") >= 3) { #pragma omp parallel for for (int i=0; i<np; ++i) { @@ -54,6 +82,21 @@ void HydroBase_Zero (CCTK_ARGUMENTS) CCTK_WARN (CCTK_WARN_ABORT, "Too many active time levels for HydroBase variables"); } + + if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::temperature") >= 3) { +#pragma omp parallel for + for (int i=0; i<np; ++i) { + temperature_p_p[i ] = 0.0; + } + } + + if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::entropy") >= 3) { +#pragma omp parallel for + for (int i=0; i<np; ++i) { + entropy_p_p[i ] = 0.0; + } + } + } else |