diff options
| -rw-r--r-- | src/setup.cc | 41 |
1 files changed, 39 insertions, 2 deletions
diff --git a/src/setup.cc b/src/setup.cc index dc24e82..897efb3 100644 --- a/src/setup.cc +++ b/src/setup.cc @@ -101,8 +101,26 @@ extern "C" void SphericalSurface_Setup (CCTK_ARGUMENTS) sf_nphi[n] += 4 - (sf_nphi[n]-2*nghostsphi[n]) % 4 + 2*nghostsphi[n]; // consistency check - assert (sf_ntheta[n] >= 3*nghoststheta[n] && sf_ntheta[n] <= maxntheta); - assert (sf_nphi[n] >= 3*nghostsphi[n] && sf_nphi[n] <= maxnphi); + if (sf_ntheta[n] < 3*nghoststheta[n]) + { + CCTK_WARN(1, "AutoRes: Determined sf_ntheta[n] < 3*nghoststheta[n]! Setting sf_ntheta[n] = 3*nghoststheta[n]!"); + sf_ntheta[n] = 3*nghoststheta[n]; + } + if (sf_ntheta[n] > maxntheta) + { + CCTK_WARN(1, "AutoRes: Determined sf_ntheta[n] > maxntheta! You should increase maxntheta or auto_res_ratio[n]! Setting sf_ntheta[n] = maxntheta!"); + sf_ntheta[n] = maxntheta; + } + if (sf_nphi[n] < 3*nghostsphi[n]) + { + CCTK_WARN(1, "AutoRes: Determined sf_nphi[n] < 3*nghostsphi[n]! Setting sf_nphi[n] = 3*nghostsphi[n]!"); + sf_nphi[n] = 3*nghostsphi[n]; + } + if (sf_nphi[n] > maxnphi) + { + CCTK_WARN(1, "AutoRes: Determined sf_nphi[n] > maxnphi! You should increase maxnphi or auto_res_ratio[n]! Setting sf_nphi[n] = maxnphi!"); + sf_nphi[n] = maxnphi; + } // check, if there are enough maxntheta,maxnphi gridpoints to remove symmetries // (some thorns require that) @@ -338,6 +356,25 @@ extern "C" void SphericalSurface_SetupRes (CCTK_ARGUMENTS) dummy = 0; MPI_Allreduce(&sf_delta_phi_estimate[n], &dummy, 1, MPI_DOUBLE, MPI_MIN, MPI_COMM_WORLD); sf_delta_phi_estimate[n] = dummy; + + /* + // input arrays and output values + const CCTK_INT input_array_variable_indices[2] + = { CCTK_VarIndex("SphericalSurface::sf_delta_theta_estimate"), + CCTK_VarIndex("SphericalSurface::sf_delta_phi_estimate") }; + const CCTK_INT output_value_type_codes[2] + = { CCTK_VARIABLE_REAL, CCTK_VARIABLE_REAL }; + void *const output_numbers[2] + = { (void *) output_for_real_values, + (void *) output_for_real_values }; + + const int status + = CCTK_ReduceGridArrays(GH, + 0, + 2, input_array_variable_indices, + 2, output_value_type_codes, + output_values); + */ } // find maximum value of "sf_minreflevel" over all processors |