From 46b6864c520fed2ff3d5a8709c7fbaab38d6644b Mon Sep 17 00:00:00 2001 From: hawke Date: Thu, 4 Sep 2003 08:35:11 +0000 Subject: Initialize all active timelevels. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinBase/ADMBase/trunk@23 d576a68a-b34a-40ae-82fc-004fa1a9d16f --- src/Initialisation.c | 101 +++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 101 insertions(+) diff --git a/src/Initialisation.c b/src/Initialisation.c index 271ef87..ba0bba8 100644 --- a/src/Initialisation.c +++ b/src/Initialisation.c @@ -88,6 +88,70 @@ void ADMBase_CartesianMinkowski(CCTK_ARGUMENTS) kxz[i] = 0.0; kyz[i] = 0.0; } + + if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::metric") > 1) + { + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + + gxx_p[i] = 1.0; + gyy_p[i] = 1.0; + gzz_p[i] = 1.0; + + gxy_p[i] = 0.0; + gxz_p[i] = 0.0; + gyz_p[i] = 0.0; + } + } + + if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::metric") > 2) + { + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + + gxx_p_p[i] = 1.0; + gyy_p_p[i] = 1.0; + gzz_p_p[i] = 1.0; + + gxy_p_p[i] = 0.0; + gxz_p_p[i] = 0.0; + gyz_p_p[i] = 0.0; + + } + + } + + if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::curv") > 1) + { + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + + kxx_p[i] = 0.0; + kyy_p[i] = 0.0; + kzz_p[i] = 0.0; + + kxy_p[i] = 0.0; + kxz_p[i] = 0.0; + kyz_p[i] = 0.0; + } + } + + if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::curv") > 2) + { + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + + kxx_p_p[i] = 0.0; + kyy_p_p[i] = 0.0; + kzz_p_p[i] = 0.0; + + kxy_p_p[i] = 0.0; + kxz_p_p[i] = 0.0; + kyz_p_p[i] = 0.0; + + } + + } } @@ -116,6 +180,23 @@ void ADMBase_LapseOne(CCTK_ARGUMENTS) { alp[i] = 1.0; } + + if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::lapse") > 1) + { + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + alp_p[i] = 1.0; + } + } + + if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::lapse") > 2) + { + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + alp_p_p[i] = 1.0; + } + } + } /*@@ @@ -144,6 +225,26 @@ void ADMBase_ShiftZero(CCTK_ARGUMENTS) betay[i] = 0.0; betaz[i] = 0.0; } + + if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::shift") > 1) + { + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + betax_p[i] = 0.0; + betay_p[i] = 0.0; + betaz_p[i] = 0.0; + } + } + + if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::shift") > 2) + { + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + betax_p_p[i] = 0.0; + betay_p_p[i] = 0.0; + betaz_p_p[i] = 0.0; + } + } } -- cgit v1.2.3