* general tidying of multiprocessor code

* add current iteration number to what's sent to processor #0 for
  printing at each iteration
* add a new verbosity level "algorithm debug" for what used to be
  #ifdef DEBUG debugging code
* add new parameters & logic to set CCTK_Warn() levels for various
  warnings (eg AH is outside the gid, metric isn't +ve definite, etc)
* restore the ability to output the horizon shape to a data file,
  which I mistakenly deleted a few days ago :(


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@947 f88db872-0e4f-0410-b76b-b9085cfa78c5

1 files changed, 83 insertions, 2 deletions

diff --git a/param.ccl b/param.ccl
index e54df30..f104150 100644
--- a/param.ccl
+++ b/param.ccl
@@ -99,6 +99,10 @@ keyword verbose_level \
 # lots of details tracing what the code is doing
 "algorithm details"    :: \
   "physics details + lots of messages about the AH-finding algorithm"
+
+# even more details tracing what the code is doing :)
+"algorithm debug"    :: \
+  "physics details + lots and lots of messages about the AH-finding algorithm"
 } "algorithm highlights"
 
 # n.b. printing timing stats is independent of  verbose_level
@@ -260,8 +264,11 @@ keyword horizon_file_format \
 "HDF5"			:: "HDF5 surface format (alas not implemented yet)"
 } "ASCII (gnuplot)"
 
+# n.b. this refers to the patch system (angular) interpatch ghost zones,
+#      *not* the Cactus interprocessor ghost zones
 boolean output_ghost_zones_for_h \
-  "should we include the ghost zones in h data files?"
+  "should we include the patch system (angular) interpatch ghost zones \
+   in h data files?"
 {
 } "false"
 
@@ -470,7 +477,7 @@ keyword patch_system_type[5] "what type of patch system should we use?"
 int N_zones_per_right_angle[5] "sets angular resolution of patch systems"
 {
 1:* :: "any integer >= 1; must be even for patch systems other than full-sphere"
-} 12
+} 18
 
 #
 # This parameter sets the width of the interpatch ghost zones in the
@@ -814,6 +821,80 @@ real initial_guess__coord_ellipsoid__z_radius[5] "z radius of ellipsoid"
 ################################################################################
 
 #
+# ***** parameters controlling how we handle errors/warnings *****
+#
+
+#
+# We currently "handle" the following errors/warnings by CCTK_VWarn();
+# these parameters give the warning levels (in the sense of the -W flag
+# to Cactus, and the first argument to CCTK_VWarn()) for various conditions:
+#
+
+# point outside (or too close to boundary of) Cactus grid,
+# i.e. geometry interpolator returns CCTK_ERROR_INTERP_POINT_OUTSIDE
+# ... warning level if error occurs on first Newton iteration,
+#     i.e. when evaluating expansion/Jacobian for initial guess
+# ... this may mean that  driver::ghost_size  is too small, especially
+#     if any horizon crosses a symmetry boundary
+#     ==> visible message by default
+int warn_level__point_outside__initial \
+  "warning level for point outside (or too close to boundary of) Cactus grid \
+   (error occurs on first Newton iteration)"
+{
+-1:* :: "any valid Cactus warning level"
+} 1
+# ... warning level if error occurs on a subsequent Newton iteration
+# ... this probably "just" means there's no horizon in the grid
+#     ==> no visible message by default
+int warn_level__point_outside__subsequent \
+  "warning level for point outside (or too close to boundary of) Cactus grid \
+   (error occurs on subsequent Newton iteration)"
+{
+-1:* :: "any valid Cactus warning level"
+} 2
+
+# the Cactus configure process didn't find a finite() function,
+# so we're skipping the "check that the geometry is finite" tests
+# even though the user set  check_that_geometry_is_finite = "true"
+int warn_level__skipping_finite_check \
+  "warning level if the user sets  check_that_geometry_is_finite \
+   but the Cactus configure process doesn't find a  finite()  function \
+   so we have to skip the finite-geometry check"
+{
+-1:* :: "any valid Cactus warning level"
+} 3
+
+# infinity and/or NaN in interpolated {g_ij, partial_k g_ij, K_ij}
+int warn_level__nonfinite_geometry \
+  "warning level if we find infinity and/or NaN in the interpolated geometry \
+   values {g_ij, partial_k g_ij, K_ij}"
+{
+-1:* :: "any valid Cactus warning level"
+} 1
+
+# interpolated g_{ij} isn't positive definite
+# (usually this means we're too close to a singularity)
+# ... warning level if error occurs on first Newton iteration,
+#     i.e. when evaluating expansion/Jacobian for initial guess
+int warn_level__gij_not_positive_definite__initial \
+  "warning level if the interpolated g_{ij} isn't positive definite \
+   (usually this means we're too close to a singularity) \
+   (error occurs on first Newton iteration)"
+{
+-1:* :: "any valid Cactus warning level"
+} 2
+# ... warning level if error occurs on a subsequent Newton iteration
+int warn_level__gij_not_positive_definite__subsequent \
+  "warning level if the interpolated g_{ij} isn't positive definite \
+   (usually this means we're too close to a singularity) \
+   (error occurs on subsequent Newton iteration)"
+{
+-1:* :: "any valid Cactus warning level"
+} 2
+
+################################################################################
+
+#
 # ***** parameters for the test driver "src/patch/test_patch_system.cc" *****
 #
 # By default this test driver isn't compiled into the cactus executable,