Changes in methodology for filling in the interior of Cook-Pfeiffer data as

used in the "turducken" paper. This includes some spline blending from
Erik and the most recent method of filling in the interior using an elliptic
equation with the "good" points as boundary values. This method uses a
conjugate gradient solver (built in), that should be robust (i.e. it should
always work) but may not be the fastest. However, since this is only done on
initial data this shouldn't be an issue. The elliptic method currently support
a second order, fourth order and sixth order equation, that will give C0, C1
and C2 solutions, respectively, across the boundary of the fill in region.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/NoExcision/trunk@7 4ec1db94-0e4f-0410-ada3-8bed251432c9

1 files changed, 84 insertions, 6 deletions

diff --git a/param.ccl b/param.ccl
index 1fe7342..6d93bd3 100644
--- a/param.ccl
+++ b/param.ccl
@@ -5,13 +5,18 @@ BOOLEAN verbose "Produce some screen output"
 {
 } "no"
 
+KEYWORD method "Method to use"
+{
+  "old" :: "Use old method"
+  "new" :: "Use new method"
+} "old"
+
+
 INT num_regions "Number of no-excision regions"
 {
   0:10 :: ""
 } 0
 
-
-
 BOOLEAN overwrite_geometry[10] "Set the geometry to Minkowski"
 {
 } "yes"
@@ -40,27 +45,100 @@ BOOLEAN overwrite_shift[10] "Set the shift to zero"
 
 
 
-REAL centre_x[10] "X-coordinate of the centre of the region"
+REAL centre_x[10] "x-coordinate of the centre of the region"
 {
   (*:*) :: ""
 } 0.0
 
-REAL centre_y[10] "Y-coordinate of the centre of the region"
+REAL centre_y[10] "y-coordinate of the centre of the region"
 {
   (*:*) :: ""
 } 0.0
 
-REAL centre_z[10] "Z-coordinate of the centre of the region"
+REAL centre_z[10] "z-coordinate of the centre of the region"
 {
   (*:*) :: ""
 } 0.0
 
+KEYWORD region_shape[10] "Shape of the region"
+{
+  "sphere"    :: "use radius"
+  "ellipsoid" :: "use radius_x, radius_y, and radius_z"
+  "surface"   :: "use a spherical surface shape"
+} "sphere"
+
 REAL radius[10] "Radius of the region"
 {
   0.0:*) :: ""
 } 1.0
 
-REAL smoothing_zone_width[10] "Width of smoothing zone inside the region"
+REAL radius_x[10] "x-radius of the region"
+{
+  0.0:*) :: ""
+} 1.0
+
+REAL radius_y[10] "y-radius of the region"
 {
   0.0:*) :: ""
+} 1.0
+
+REAL radius_z[10] "z-radius of the region"
+{
+  0.0:*) :: ""
+} 1.0
+
+INT surface_index[10] "Spherical surface index"
+{
+  0:* :: "must be an index of a spherical surface"
+} 0
+
+
+
+BOOLEAN reduce_rhs[10] "Reduce RHS"
+{
+} "no"
+
+REAL reduction_factor[10] "Reduction factor for RHS (0=complete, 1=no reduction)"
+{
+  *:* :: ""
+} 0.0
+
+
+
+KEYWORD smoothing_function[10] "Smoothing function"
+{
+  "linear" :: "linear ramp"
+  "spline" :: "cubic spline ramp"
+  "cosine" :: "cosine ramp"
+} "linear"
+
+REAL smoothing_zone_width[10] "Relative width of smoothing zone inside the region"
+{
+  0.0:1.0 :: ""
 } 0.0
+
+INT smoothing_iterations "Smoothing iterations"
+{
+  0:* :: ""
+} 10
+
+REAL smoothing_factor "Initial moothing factor"
+{
+  (0:2) :: ""
+} 1.2
+
+
+
+BOOLEAN smooth_regions "Smooth overwritten regions?"
+{
+} "no"
+
+INT smoothing_order "Order of the derivatives used for CG smoothing"
+{
+  2:6:2 :: ""
+} 6
+
+REAL smoothing_eps "CG smoothing stop criteria"
+{
+  (0.0:* :: ""
+} 1e-6