}

  DECLARE_CCTK_PARAMETERS;
  if (!poison_undefined_values)
    return;
  const valid_t &valid = scalargroupdata.valid.at(tl).at(vi);
  if (!valid.valid_int && !valid.valid_bnd)
    return;
  // scalar have no boundary so we expect it to alway be valid
  assert(valid.valid_bnd);
  atomic<bool> found_nan{false};
  if (valid.valid_int) {
    CCTK_REAL *restrict const ptr = scalargroupdata.data.at(tl).at(vi);
    if (CCTK_BUILTIN_EXPECT(isnan(*ptr), false)) {
      found_nan = true;
    }
  }
  if (CCTK_BUILTIN_EXPECT(found_nan, false)) {
    const char *where = valid.valid_int && valid.valid_bnd
                            ? "interior and boundary"
                            : valid.valid_int ? "interior" : "boundary";

                tl, where);
  }

}

    const auto &restrict globaldata = ghext->globaldata;

    const auto &restrict globaldata = ghext->globaldata;

    auto &restrict globaldata = ghext->globaldata;

    auto &restrict globaldata = ghext->globaldata;