# 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"

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"
}

    # READS: dtlapse(everywhere)

    READS: dtlapse(everywhere)

    # WRITES: dtlapse(everywhere)

    WRITES: dtlapse(everywhere)

// 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; });
}

  // const GF3D<CCTK_REAL, 0, 0, 0> dtalp_(cctkGH, dtalp);

  const GF3D<CCTK_REAL, 0, 0, 0> dtalp_(cctkGH, dtalp);

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

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