if (dimension == 0) {
      *ptr = 0.0 / 0.0;
    } else {
      const int *gsh = arraygroupdata.gsh;
      int n_elems = 1;
      for (int i = 0; i < dimension; i++) n_elems *= gsh[i];
      for (int i = 0; i < n_elems; i++) ptr[i] = 0.0 / 0.0;
    }

    const int *gsh = arraygroupdata.gsh;
    int n_elems = 1;
    for (int i = 0; i < dimension; i++) n_elems *= gsh[i];
    for (int i = 0; i < n_elems; i++) ptr[i] = 0.0 / 0.0;

    if (dimension == 0) {
      if (CCTK_BUILTIN_EXPECT(!CCTK_isfinite(*ptr), false)) {

    const int *gsh = arraygroupdata.gsh;
    int n_elems = 1;
    for (int i = 0; i < dimension; i++) n_elems *= gsh[i];
    for (int i = 0; i < n_elems; i++) {
      if (CCTK_BUILTIN_EXPECT(!CCTK_isfinite(ptr[i]), false)) {

    } else {
      const int *gsh = arraygroupdata.gsh;
      int n_elems = 1;
      for (int i = 0; i < dimension; i++) {
        n_elems *= gsh[i];
      }
      for (int i = 0; i < n_elems; i++) {
        if (CCTK_BUILTIN_EXPECT(!CCTK_isfinite(ptr[i]), false)) {
          ++nan_count;
        }
      }