diff options
Diffstat (limited to 'src/Initialisation.c')
| -rw-r--r-- | src/Initialisation.c | 186 |
1 files changed, 186 insertions, 0 deletions
diff --git a/src/Initialisation.c b/src/Initialisation.c index ba0bba8..f24c672 100644 --- a/src/Initialisation.c +++ b/src/Initialisation.c @@ -33,9 +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_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); /******************************************************************** ********************* Other Routine Prototypes ********************* @@ -249,6 +255,98 @@ void ADMBase_ShiftZero(CCTK_ARGUMENTS) } /*@@ + @routine ADMBase_DtlapseZero + @date Sep 06 2006 + @author Erik Schnetter + @desc + Scheduled routine to initialise 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 Sep 06 2006 + @author Erik Schnetter + @desc + Scheduled routine to initialise 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_SetShiftStateOn @date Thu Apr 25 23:12:18 2002 @author Tom Goodale @@ -292,6 +390,94 @@ void ADMBase_SetShiftStateOff(CCTK_ARGUMENTS) *shift_state = 0; } +/*@@ + @routine ADMBase_SetDtlapseStateOn + @date Sep 06 2006 + @author Erik Schnetter + @desc + Scheduled routine to set the value of the dtlapse state to on. + @enddesc + @calls + @calledby + @history + + @endhistory + + @@*/ +void ADMBase_SetDtlapseStateOn(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS; + DECLARE_CCTK_PARAMETERS; + + *dtlapse_state = 1; +} + +/*@@ + @routine ADMBase_SetDtlapseStateOff + @date Sep 06 2006 + @author Erik Schnetter + @desc + Scheduled routine to set the value of the dtlapse state to off. + @enddesc + @calls + @calledby + @history + + @endhistory + + @@*/ +void ADMBase_SetDtlapseStateOff(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS; + DECLARE_CCTK_PARAMETERS; + + *dtlapse_state = 0; +} + +/*@@ + @routine ADMBase_SetDtshiftStateOn + @date Sep 06 2006 + @author Erik Schnetter + @desc + Scheduled routine to set the value of the dtshift state to on. + @enddesc + @calls + @calledby + @history + + @endhistory + + @@*/ +void ADMBase_SetDtshiftStateOn(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS; + DECLARE_CCTK_PARAMETERS; + + *dtshift_state = 1; +} + +/*@@ + @routine ADMBase_SetDtshiftStateOff + @date Sep 06 2006 + @author Erik Schnetter + @desc + Scheduled routine to set the value of the dtshift state to off. + @enddesc + @calls + @calledby + @history + + @endhistory + + @@*/ +void ADMBase_SetDtshiftStateOff(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS; + DECLARE_CCTK_PARAMETERS; + + *dtshift_state = 0; +} + /******************************************************************** |