CCTK_REAL dtlapse TYPE=gf TAGS='index={0 0 0}' { dtalp } "Time derivative of ADM lapse function"

# TODO: This is temporarily deactivated since BaikalX doesn't handle this
# variable correctly. Other respective parts of ADMBase have also been
# deactivated. We need to update BaikalX to set this variable as well.
# CCTK_REAL dtlapse TYPE=gf TAGS='index={0 0 0}' { dtalp } "Time derivative of ADM lapse function"

KEYWORD initial_dtlapse "Initial dtlapse value"
{
  "zero" :: "Dtlapse is zero"
} "zero"

# KEYWORD initial_dtlapse "Initial dtlapse value"
# {
#   "zero" :: "Dtlapse is zero"
# } "zero"

if (CCTK_EQUALS(initial_dtlapse, "zero")) {
  SCHEDULE ADMBase_initial_dtlapse IN ADMBase_InitialGauge
  {
    LANG: C
    WRITES: dtlapse(everywhere)
  } "Set dtlapse to zero"
}

# if (CCTK_EQUALS(initial_dtlapse, "zero")) {
#   SCHEDULE ADMBase_initial_dtlapse IN ADMBase_InitialGauge
#   {
#     LANG: C
#     WRITES: dtlapse(everywhere)
#   } "Set dtlapse to zero"
# }

extern "C" void ADMBase_initial_dtlapse(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_ADMBase_initial_dtlapse;
  DECLARE_CCTK_PARAMETERS;
  const GF3D<CCTK_REAL, 0, 0, 0> dtalp_(cctkGH, dtalp);

  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { dtalp_(p.I) = 0; });
}

// extern "C" void ADMBase_initial_dtlapse(CCTK_ARGUMENTS) {
//   DECLARE_CCTK_ARGUMENTS_ADMBase_initial_dtlapse;
//   DECLARE_CCTK_PARAMETERS;
// 
//   const GF3D<CCTK_REAL, 0, 0, 0> dtalp_(cctkGH, dtalp);
// 
//   loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { dtalp_(p.I) = 0; });
// }