parameter files and shell scripts to test surface-integral computations

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

4 files changed, 129 insertions, 0 deletions

diff --git a/run/test-ahfinderdirect/surface/README b/run/test-ahfinderdirect/surface/README
new file mode 100644
index 0000000..7c33af9
--- /dev/null
+++ b/run/test-ahfinderdirect/surface/README
@@ -0,0 +1,3 @@
+This directory contains parameter files to test the surface integral
+computations used to estimate black hole areas, masses, and centroid
+positions.
diff --git a/run/test-ahfinderdirect/surface/doit.grep b/run/test-ahfinderdirect/surface/doit.grep
new file mode 100755
index 0000000..b006e82
--- /dev/null
+++ b/run/test-ahfinderdirect/surface/doit.grep
@@ -0,0 +1,19 @@
+#!/bin/sh
+
+	for method in trapezoid Simpson Simpson-variant
+	do
+	for offset in '' '-offset'
+	do
+	name="$method$offset"
+	echo $name
+	fgrep 'AH found' Schw*$name.log \
+	   | sed -e 's/^Schw//' \
+		 -e 's/-.*AH found: A=/	/' \
+		 -e 's/ m=/	/' \
+		 -e 's/ at (/	/' \
+		 -e 's/,/	/g' \
+		 -e 's/)$//' \
+	   | pre1st '# delta	area	mass	posn_x	posn_y	posn_z' \
+	   > $name.dat
+	done
+	done
diff --git a/run/test-ahfinderdirect/surface/doit.run b/run/test-ahfinderdirect/surface/doit.run
new file mode 100755
index 0000000..367b397
--- /dev/null
+++ b/run/test-ahfinderdirect/surface/doit.run
@@ -0,0 +1,40 @@
+#!/bin/sh
+
+	for offset in 0 1
+	do
+	if [ $offset -eq 1 ]
+	   then offset_name='-offset'
+		origin_x=0.5
+		origin_y=0.7
+		origin_z=0.6
+	   else offset_name=''
+		origin_x=0.0
+		origin_y=0.0
+		origin_z=0.0
+	fi
+
+		for delta in 7.5 5.0 3.0 2.25 1.5 1.0
+		do
+		for method in 'trapezoid' 'Simpson' 'Simpson (variant)'
+		do
+		if [ "$method" = 'Simpson (variant)' ]
+		   then method_name='Simpson-variant'
+		   else method_name="$method"
+		fi
+		base=Schw${delta}-${method_name}${offset_name}
+		##echo ":${delta}:${method}:${origin_x}:${origin_y}:${origin_z}:"
+		##echo ":${delta}:${method_name}:${origin_x}:${origin_y}:${origin_z}:"
+		##echo $base
+
+		sed -e "s#\$delta_drho_dsigma#$delta#" \
+		    -e "s#\$surface_integral_method#$method#" \
+		    -e "s#\$origin_x#$origin_x#" \
+		    -e "s#\$origin_y#$origin_y#" \
+		    -e "s#\$origin_z#$origin_z#" \
+		    -e "s#\$base_name#$base#" \
+		    <template.par >$base.par
+		cactus_test-ahfinderdirect $base.par 2>&1 | tee $base.log
+		
+		done
+		done
+	done
diff --git a/run/test-ahfinderdirect/surface/template.par b/run/test-ahfinderdirect/surface/template.par
new file mode 100644
index 0000000..a44a644
--- /dev/null
+++ b/run/test-ahfinderdirect/surface/template.par
@@ -0,0 +1,67 @@
+# parameter file for patch system test
+# uses variables
+#	$delta_drho_dsigma
+#	$surface_integral_method
+#	$origin_x
+#	$origin_y
+#	$origin_z
+#	$base_name
+ActiveThorns = "CartGrid3D LocalInterp PUGH ADMBase ADMCoupling CoordGauge Exact AHFinderDirect"
+
+# flesh
+cactus::cctk_itlast = 0
+
+# PUGH
+driver::ghost_size = 2
+driver::global_nx = 31
+driver::global_ny = 31
+driver::global_nz = 31
+
+# CartGrid3D
+grid::domain = "full"
+grid::avoid_origin = "false"
+grid::type = "byspacing"
+grid::dxyz = 0.2
+
+# ADMBase
+ADMBase::initial_lapse = "exact"
+ADMBase::initial_shift = "exact"
+ADMBase::initial_data = "exact"
+ADMBase::lapse_evolution_method = "static"
+ADMBase::shift_evolution_method = "static"
+ADMBase::metric_type = "physical"
+
+# Exact
+Exact::exact_model = "Schwarzschild/EF"
+Exact::Schwarzschild_EF__mass = 1.0
+
+AHFinderDirect::find_AHs_at_postinitial = "true"
+AHFinderDirect::find_AHs_at_poststep    = "false"
+AHFinderDirect::verbose_level = "algorithm details"
+AHFinderDirect::print_timing_stats = "true"
+
+AHFinderDirect::method = "horizon function"
+AHFinderDirect::geometry_method = "Schwarzschild/EF"
+AHFinderDirect::final_H_update_if_exit_x_H_small = "true"
+AHFinderDirect::surface_integral_method = "$surface_integral_method"
+##AHFinderDirect::output_initial_guess = "false"
+##AHFinderDirect::output_h = "false"
+##AHFinderDirect::output_H_of_h = "false"
+AHFinderDirect::h_base_file_name      = "$base_name.h"
+AHFinderDirect::H_of_h_base_file_name = "$base_name.H"
+
+AHFinderDirect::N_horizons = 1
+AHFinderDirect::origin_x[1] = $origin_x
+AHFinderDirect::origin_y[1] = $origin_y
+AHFinderDirect::origin_z[1] = $origin_z
+AHFinderDirect::N_ghost_points = 2
+AHFinderDirect::N_overlap_points = 1
+AHFinderDirect::delta_drho_dsigma = $delta_drho_dsigma
+
+AHFinderDirect::geometry_interpolator_pars = "order=3   out_of_range_tolerance={-1.0 -1.0 -1.0 -1.0 -1.0 -1.0}"
+
+AHFinderDirect::initial_guess_method = "sphere"
+AHFinderDirect::initial_guess__sphere__x_center[1] = 0.0
+AHFinderDirect::initial_guess__sphere__y_center[1] = 0.0
+AHFinderDirect::initial_guess__sphere__z_center[1] = 0.0
+AHFinderDirect::initial_guess__sphere__radius[1] = 2.0