diff options
author | knarf <knarf@57fe0bb3-ccba-405f-9b23-de0201f165b7> | 2010-03-08 16:12:49 +0000 |
---|---|---|
committer | knarf <knarf@57fe0bb3-ccba-405f-9b23-de0201f165b7> | 2010-03-08 16:12:49 +0000 |
commit | 9565258e85778f68e55cf346399ba52a63687022 (patch) | |
tree | 3dc5be82c548d94b8c1bfc0136634a64d67b35f9 /src/Initialisation.c | |
parent | 284f26cffee4e51ec42b9b47b1bf615232b269ed (diff) |
patch from Christian: add entropy and temperature
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinBase/HydroBase/trunk@23 57fe0bb3-ccba-405f-9b23-de0201f165b7
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 |