From c2f0a62812f40c9aef8ad75c28689776b6344da9 Mon Sep 17 00:00:00 2001 From: jthorn Date: Thu, 8 Jan 2004 12:51:03 +0000 Subject: Raise the limit for the maximum allowable number of AHs you can search for, from 10 all the way up to 100. Hopefully that will be enough for now... git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1234 f88db872-0e4f-0410-b76b-b9085cfa78c5 --- param.ccl | 68 +++++++++++++++++++++++++++++++-------------------------------- 1 file changed, 34 insertions(+), 34 deletions(-) (limited to 'param.ccl') diff --git a/param.ccl b/param.ccl index f22c574..dc9170d 100644 --- a/param.ccl +++ b/param.ccl @@ -75,8 +75,8 @@ keyword method "what should this thorn do for each apparent horizon?" # int N_horizons "number of apparent horizons to search for" { -0 :: "turn this thorn into a fancy no-op :)" -1:10 :: "search for this many apparent horizons" +0 :: "turn this thorn into a fancy no-op :)" +1:100 :: "search for this many apparent horizons" } 1 # @@ -89,8 +89,8 @@ int N_horizons "number of apparent horizons to search for" int which_horizon_to_announce_centroid \ "for which horizon should we announce the centroid?" { -0 :: "don't announce any centroid(s)" -1:10 :: "announce this horizon's centroid each time we find it" +0 :: "don't announce any centroid(s)" +1:100 :: "announce this horizon's centroid each time we find it" } 0 # @@ -130,7 +130,7 @@ boolean print_timing_stats \ # Theta(h) = 0 # This parameter allows the RHS to be set to any specified constant, # to find a surface of constant expansion. -real surface_expansion[11] "search for a surface with this (constant) expansion" +real surface_expansion[101] "search for a surface with this (constant) expansion" { *:* :: "any real number" } 0.0 @@ -147,7 +147,7 @@ real surface_expansion[11] "search for a surface with this (constant) expansion" # next attempting to find this horizon.) # -keyword initial_guess_method[11] \ +keyword initial_guess_method[101] \ "method used to set up initial guess for apparent horizon shape" { "read from named file" :: "read from explicitly-named input file" @@ -163,105 +163,105 @@ keyword initial_guess_method[11] \ } "coordinate sphere" # parameters for initial_guess_method = "read from named file" -string initial_guess__read_from_named_file__file_name[11] \ +string initial_guess__read_from_named_file__file_name[101] \ "file name to read initial guess from" { ".+" :: "file name to read initial guess from" } "h.gp" # parameters for initial_guess_method = "Kerr/Kerr" -real initial_guess__Kerr_Kerr__x_posn[11] "x coordinate of Kerr BH" +real initial_guess__Kerr_Kerr__x_posn[101] "x coordinate of Kerr BH" { *:* :: "any real number" } 0.0 -real initial_guess__Kerr_Kerr__y_posn[11] "y coordinate of Kerr BH" +real initial_guess__Kerr_Kerr__y_posn[101] "y coordinate of Kerr BH" { *:* :: "any real number" } 0.0 -real initial_guess__Kerr_Kerr__z_posn[11] "z coordinate of Kerr BH" +real initial_guess__Kerr_Kerr__z_posn[101] "z coordinate of Kerr BH" { *:* :: "any real number" } 0.0 -real initial_guess__Kerr_Kerr__mass[11] "mass of Kerr BH" +real initial_guess__Kerr_Kerr__mass[101] "mass of Kerr BH" { (0.0:* :: "BH mass = any real number > 0" } 1.0 # n.b. my convention is that a=J/m^2 is dimensionless, # while MTW take a=J/m=m * (my a) -real initial_guess__Kerr_Kerr__spin[11] "dimensionless spin a=J/m^2 of Kerr BH" +real initial_guess__Kerr_Kerr__spin[101] "dimensionless spin a=J/m^2 of Kerr BH" { (-1.0:1.0) :: \ "dimensionless BH spin = J/m^2 = any real number with absolute value < 1" } 0.6 # parameters for initial_guess_method = "Kerr/Kerr-Schild" -real initial_guess__Kerr_KerrSchild__x_posn[11] "x coordinate of Kerr BH" +real initial_guess__Kerr_KerrSchild__x_posn[101] "x coordinate of Kerr BH" { *:* :: "any real number" } 0.0 -real initial_guess__Kerr_KerrSchild__y_posn[11] "y coordinate of Kerr BH" +real initial_guess__Kerr_KerrSchild__y_posn[101] "y coordinate of Kerr BH" { *:* :: "any real number" } 0.0 -real initial_guess__Kerr_KerrSchild__z_posn[11] "z coordinate of Kerr BH" +real initial_guess__Kerr_KerrSchild__z_posn[101] "z coordinate of Kerr BH" { *:* :: "any real number" } 0.0 -real initial_guess__Kerr_KerrSchild__mass[11] "mass of Kerr BH" +real initial_guess__Kerr_KerrSchild__mass[101] "mass of Kerr BH" { (0.0:* :: "BH mass = any real number > 0" } 1.0 # n.b. my convention is that a=J/m^2 is dimensionless, # while MTW take a=J/m=m * (my a) -real initial_guess__Kerr_KerrSchild__spin[11] "dimensionless spin a=J/m^2 of Kerr BH" +real initial_guess__Kerr_KerrSchild__spin[101] "dimensionless spin a=J/m^2 of Kerr BH" { (-1.0:1.0) :: \ "dimensionless BH spin = J/m^2 = any real number with absolute value < 1" } 0.6 # parameters for initial_guess_method = "sphere" -real initial_guess__coord_sphere__x_center[11] "x coordinate of sphere center" +real initial_guess__coord_sphere__x_center[101] "x coordinate of sphere center" { *:* :: "any real number" } 0.0 -real initial_guess__coord_sphere__y_center[11] "y coordinate of sphere center" +real initial_guess__coord_sphere__y_center[101] "y coordinate of sphere center" { *:* :: "any real number" } 0.0 -real initial_guess__coord_sphere__z_center[11] "z coordinate of sphere center" +real initial_guess__coord_sphere__z_center[101] "z coordinate of sphere center" { *:* :: "any real number" } 0.0 -real initial_guess__coord_sphere__radius[11] "radius of sphere" +real initial_guess__coord_sphere__radius[101] "radius of sphere" { (0.0:* :: "any real number > 0.0" } 2.0 # parameters for initial_guess_method = "ellipsoid" -real initial_guess__coord_ellipsoid__x_center[11] \ +real initial_guess__coord_ellipsoid__x_center[101] \ "x coordinate of ellipsoid center" { *:* :: "any real number" } 0.0 -real initial_guess__coord_ellipsoid__y_center[11] \ +real initial_guess__coord_ellipsoid__y_center[101] \ "y coordinate of ellipsoid center" { *:* :: "any real number" } 0.0 -real initial_guess__coord_ellipsoid__z_center[11] \ +real initial_guess__coord_ellipsoid__z_center[101] \ "z coordinate of ellipsoid center" { *:* :: "any real number" } 0.0 -real initial_guess__coord_ellipsoid__x_radius[11] "x radius of ellipsoid" +real initial_guess__coord_ellipsoid__x_radius[101] "x radius of ellipsoid" { (0.0:* :: "any real number > 0.0" } 2.0 -real initial_guess__coord_ellipsoid__y_radius[11] "y radius of ellipsoid" +real initial_guess__coord_ellipsoid__y_radius[101] "y radius of ellipsoid" { (0.0:* :: "any real number > 0.0" } 2.0 -real initial_guess__coord_ellipsoid__z_radius[11] "z radius of ellipsoid" +real initial_guess__coord_ellipsoid__z_radius[101] "z radius of ellipsoid" { (0.0:* :: "any real number > 0.0" } 2.0 @@ -477,7 +477,7 @@ boolean set_mask_for_all_horizons \ { } "false" -boolean set_mask_for_individual_horizon[11] \ +boolean set_mask_for_individual_horizon[101] \ "should we set a mask grid function (or functions) for *this* horizon?" { } "false" @@ -762,7 +762,7 @@ real max_allowable_Delta_h_over_h \ # sphere a bit larger than you expect the horizon to be (eg a sphere with # areal radius 4m or so, where m is the ADM mass of the slice). # -real max_allowable_horizon_radius[11] \ +real max_allowable_horizon_radius[101] \ "max mean-coordinate-radius allowed for any trial surface \ before we give up and say we can't find this horizon" { @@ -816,15 +816,15 @@ real max_allowable_Theta \ # of the horizon radius is no problem, and even 1/2 the horizon radius # only slows the convergence by an extra iteration or two. # -real origin_x[11] "global x coordinate of patch system origin" +real origin_x[101] "global x coordinate of patch system origin" { *:* :: "any real number" } 0.0 -real origin_y[11] "global y coordinate of patch system origin" +real origin_y[101] "global y coordinate of patch system origin" { *:* :: "any real number" } 0.0 -real origin_z[11] "global z coordinate of patch system origin" +real origin_z[101] "global z coordinate of patch system origin" { *:* :: "any real number" } 0.0 @@ -834,7 +834,7 @@ real origin_z[11] "global z coordinate of patch system origin" # but alas they don't work yet for apparent horizon finding # (the Jacobian computation doesn't yet grok the nonlocal rotation BCs). # -keyword patch_system_type[11] "what type of patch system should we use?" +keyword patch_system_type[101] "what type of patch system should we use?" { # choose this for normal use "match Cactus grid symmetry" :: \ @@ -887,7 +887,7 @@ keyword patch_system_type[11] "what type of patch system should we use?" # ... 3rd power of this parameter # if Jacobian_store_solve_method = "row-oriented sparse matrix/ILUCG" # -int N_zones_per_right_angle[11] "sets angular resolution of patch systems" +int N_zones_per_right_angle[101] "sets angular resolution of patch systems" { 1:* :: "any integer >= 1; must be even for patch systems other than full-sphere" } 18 -- cgit v1.2.3