* major changes to driver routines to find multiple horizons in parallel

  across multiple processors -- see src/driver/README.parallel for details
* drop convergence checks on ||Delta_h|| in param.ccl because they don't
  fit well with parallelization changes

==> With this changes, AHFinderDirect is now (I think)
    multiprocessor-ready!!


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

1 files changed, 19 insertions, 29 deletions

diff --git a/param.ccl b/param.ccl
index f5ce1e5..e54df30 100644
--- a/param.ccl
+++ b/param.ccl
@@ -43,14 +43,21 @@ boolean find_AHs_at_poststep \
 keyword method "what should this thorn do for each apparent horizon?"
 {
 # these options are mostly for testing/debugging
-"evaluate expansion"	  :: "evaluate the LHS function Theta(h)"
-"test expansion Jacobian" :: "compute/print the J[Theta(h)] Jacobian matrix \
-			      (possibly in several ways, depending on \
-			      the  test_all_Jacobian_methods  parameter"
+# ... in a multiprocessor Cactus run, the horizons are done sequentually
+#     on processor #0; the other processors do dummy computations
+"evaluate expansions"	  :: "evaluate the LHS function Theta(h)"
+
+# ... in a multiprocessor Cactus run, the Jacobian is computed on
+#     processor #0; the other processors do dummy computations
+"test expansion Jacobians" :: \
+  "compute/print horizon 1's J[Theta(h)] Jacobian matrix (possibly in \
+   multiple ways, depending on the  test_all_Jacobian_methods  parameter)"
 
 # this is for normal apparent horizon finding
-"find horizon"		  :: "find the apparent horizon"
-} "find horizon"
+# ... in a multiprocessor Cactus run, the horizons are done in parallel
+#     across processors; see src/driver/README.parallel for details
+"find horizons"		  :: "find the apparent horizon"
+} "find horizons"
 
 #
 # We support searching for up to N_horizons distinct apparent horizons
@@ -72,10 +79,8 @@ int N_horizons "number of apparent horizons to search for"
 } 1
 
 #
-# This parameter controls which (how many) informational (non-error)
-# messages are printed describing the operation of the apparent horizon
-# finder.  This is analogous to the -W Cactus command-line option, except
-# that this parameter controls (only) *non*-error messages.
+# This parameter controls how verbose this thorn is in printing
+# informational (non-error) messages describing what it's doing.
 #
 keyword verbose_level \
   "controls which (how many) messages to print describing AH finding"
@@ -200,35 +205,20 @@ int max_Newton_iterations__subsequent \
 # move in any single Newton iteration (i.e. the infinity-norm of Delta_h)
 # to <= this fraction of the mean horizon radius
 #
-real max_Delta_h_over_h \
+real max_allowable_Delta_h_over_h \
   "don't let horizon move > this fraction of mean radius in a Newton iteration"
 {
 (0.0:* :: "any positive real number"
 } 0.1
 
 #
-# we declare convergence if *either* of the following two criteria are met
+# convergence tolerance for the Newton iteration
 #
-real Theta_norm_for_convergence "declare convergence if ||Theta||_inf <= this"
+real Theta_norm_for_convergence \
+  "we declare the horizon to be found if ||Theta||_infinity <= this"
 {
 (0.0:* :: "any positive real number"
 } 1.0e-8
-real Delta_h_norm_for_convergence \
-  "declare convergence after any Newton step with ||Delta_h||_inf <= this"
-{
-(0.0:* :: "any positive real number"
-} 1.0e-8
-
-#
-# This only needs to be set to true for very careful convergence studies
-# etc.  On the other hand, setting it to true probably only slows down
-# the apparent horizon finder by a few percent.
-#
-boolean final_Theta_update_if_Delta_h_converged \
-  "should we do a final Theta(h) update if we terminate the \
-   Newton iteration by the small-||Delta_h|| convergence criterion?"
-{
-} "false"
 
 ################################################################################