Introduce a parameter TP_epsilon that can be used to smooth out

initial data.  It has the same meaning as for the Brill-Lindquist
initial data in thorn AEIThorns/Exact, namely by

        r -> (r^4 + epsilon^4)^(1/4)

By default, epsilon=0, which introduces no change.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@51 b2a53a04-0f4f-0410-87ed-f9f25ced00cf

1 files changed, 5 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
index 8764a16..d578502 100644
--- a/param.ccl
+++ b/param.ccl
@@ -71,6 +71,11 @@ INT Newton_maxit "Maximum number of Newton iterations"
   0:* :: ""
 } 5
 
+REAL TP_epsilon "A small number to smooth out singularities at the puncture locations"
+{
+  0:* :: ""
+} 0.0
+
 REAL TP_Tiny "Tiny number to avoid nans near or at the pucture locations"
 {
   0:* :: "anything positive, usually very small"