diff options
Diffstat (limited to 'src/Initialisation.c')
-rw-r--r-- | src/Initialisation.c | 128 |
1 files changed, 110 insertions, 18 deletions
diff --git a/src/Initialisation.c b/src/Initialisation.c index f24c672..10fae5e 100644 --- a/src/Initialisation.c +++ b/src/Initialisation.c @@ -33,15 +33,15 @@ CCTK_FILEVERSION(CactusEinstein_ADMBase_Initialisation_c) void ADMBase_CartesianMinkowski(CCTK_ARGUMENTS); void ADMBase_LapseOne(CCTK_ARGUMENTS); void ADMBase_ShiftZero(CCTK_ARGUMENTS); -void ADMBase_DtlapseZero(CCTK_ARGUMENTS); -void ADMBase_DtshiftZero(CCTK_ARGUMENTS); +void ADMBase_DtLapseZero(CCTK_ARGUMENTS); +void ADMBase_DtShiftZero(CCTK_ARGUMENTS); void ADMBase_SetShiftStateOn(CCTK_ARGUMENTS); void ADMBase_SetShiftStateOff(CCTK_ARGUMENTS); -void ADMBase_SetDtlapseStateOn(CCTK_ARGUMENTS); -void ADMBase_SetDtlapseStateOff(CCTK_ARGUMENTS); -void ADMBase_SetDtshiftStateOn(CCTK_ARGUMENTS); -void ADMBase_SetDtshiftStateOff(CCTK_ARGUMENTS); +void ADMBase_SetDtLapseStateOn(CCTK_ARGUMENTS); +void ADMBase_SetDtLapseStateOff(CCTK_ARGUMENTS); +void ADMBase_SetDtShiftStateOn(CCTK_ARGUMENTS); +void ADMBase_SetDtShiftStateOff(CCTK_ARGUMENTS); /******************************************************************** ********************* Other Routine Prototypes ********************* @@ -255,6 +255,98 @@ void ADMBase_ShiftZero(CCTK_ARGUMENTS) } /*@@ + @routine ADMBase_DtLapseZero + @date Oct 18 2007 + @author Erik Schnetter + @desc + Scheduled routine to initialise the dtlapse to zero. + @enddesc + @calls + @calledby + @history + + @endhistory + + @@*/ +void ADMBase_DtLapseZero(CCTK_ARGUMENTS) +{ + int i; + DECLARE_CCTK_ARGUMENTS; + DECLARE_CCTK_PARAMETERS; + + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + dtalp[i] = 0.0; + } + + if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::dtlapse") > 1) + { + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + dtalp_p[i] = 0.0; + } + } + + if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::dtlapse") > 2) + { + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + dtalp_p_p[i] = 0.0; + } + } + +} + +/*@@ + @routine ADMBase_DtShiftZero + @date Oct 18 2007 + @author Erik Schnetter + @desc + Scheduled routine to initialise the dtshift to zero. + @enddesc + @calls + @calledby + @history + + @endhistory + + @@*/ +void ADMBase_DtShiftZero(CCTK_ARGUMENTS) +{ + int i; + DECLARE_CCTK_ARGUMENTS; + DECLARE_CCTK_PARAMETERS; + + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + dtbetax[i] = 0.0; + dtbetay[i] = 0.0; + dtbetaz[i] = 0.0; + } + + if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::dtshift") > 1) + { + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + dtbetax_p[i] = 0.0; + dtbetay_p[i] = 0.0; + dtbetaz_p[i] = 0.0; + } + } + + if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::dtshift") > 2) + { + for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++) + { + dtbetax_p_p[i] = 0.0; + dtbetay_p_p[i] = 0.0; + dtbetaz_p_p[i] = 0.0; + } + } + +} + +/*@@ @routine ADMBase_DtlapseZero @date Sep 06 2006 @author Erik Schnetter @@ -391,8 +483,8 @@ void ADMBase_SetShiftStateOff(CCTK_ARGUMENTS) } /*@@ - @routine ADMBase_SetDtlapseStateOn - @date Sep 06 2006 + @routine ADMBase_SetDtLapseStateOn + @date Oct 18 2007 @author Erik Schnetter @desc Scheduled routine to set the value of the dtlapse state to on. @@ -404,7 +496,7 @@ void ADMBase_SetShiftStateOff(CCTK_ARGUMENTS) @endhistory @@*/ -void ADMBase_SetDtlapseStateOn(CCTK_ARGUMENTS) +void ADMBase_SetDtLapseStateOn(CCTK_ARGUMENTS) { DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; @@ -413,8 +505,8 @@ void ADMBase_SetDtlapseStateOn(CCTK_ARGUMENTS) } /*@@ - @routine ADMBase_SetDtlapseStateOff - @date Sep 06 2006 + @routine ADMBase_SetDtLapseStateOff + @date Oct 18 2007 @author Erik Schnetter @desc Scheduled routine to set the value of the dtlapse state to off. @@ -426,7 +518,7 @@ void ADMBase_SetDtlapseStateOn(CCTK_ARGUMENTS) @endhistory @@*/ -void ADMBase_SetDtlapseStateOff(CCTK_ARGUMENTS) +void ADMBase_SetDtLapseStateOff(CCTK_ARGUMENTS) { DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; @@ -435,8 +527,8 @@ void ADMBase_SetDtlapseStateOff(CCTK_ARGUMENTS) } /*@@ - @routine ADMBase_SetDtshiftStateOn - @date Sep 06 2006 + @routine ADMBase_SetDtShiftStateOn + @date Oct 18 2007 @author Erik Schnetter @desc Scheduled routine to set the value of the dtshift state to on. @@ -448,7 +540,7 @@ void ADMBase_SetDtlapseStateOff(CCTK_ARGUMENTS) @endhistory @@*/ -void ADMBase_SetDtshiftStateOn(CCTK_ARGUMENTS) +void ADMBase_SetDtShiftStateOn(CCTK_ARGUMENTS) { DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; @@ -457,8 +549,8 @@ void ADMBase_SetDtshiftStateOn(CCTK_ARGUMENTS) } /*@@ - @routine ADMBase_SetDtshiftStateOff - @date Sep 06 2006 + @routine ADMBase_SetDtShiftStateOff + @date Oct 18 2007 @author Erik Schnetter @desc Scheduled routine to set the value of the dtshift state to off. @@ -470,7 +562,7 @@ void ADMBase_SetDtshiftStateOn(CCTK_ARGUMENTS) @endhistory @@*/ -void ADMBase_SetDtshiftStateOff(CCTK_ARGUMENTS) +void ADMBase_SetDtShiftStateOff(CCTK_ARGUMENTS) { DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; |