This thorn EinsteinInitial/Exact uses finite differencing to calculate

the extrinsic curvature. It uses currently only a 2nd order stencil.
This patch implements 4th order finite differencing, which reduces the
error significantly.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/Exact/trunk@266 e296648e-0e4f-0410-bd07-d597d9acff87

1 files changed, 17 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
index 1c03d3f..9164f7d 100644
--- a/param.ccl
+++ b/param.ccl
@@ -103,6 +103,23 @@ KEYWORD overwrite_boundary "Overwrite g and K on the boundary"
 } "no"
 
 ################################################################################
+##### parameters for calculating the extrinsic curvature #######################
+################################################################################
+
+private:
+
+INT exact_order "finite differencing order"
+{
+2 :: "2"
+4 :: "4"
+} 2
+
+REAL exact_eps "finite differencing stencil size (in terms of the grid spacing)"
+{
+(0.0:* :: "Positive please"
+} 1.0e-6
+
+################################################################################
 ##### parameters for boosting any non-stress-energy-tensor model ###############
 ################################################################################