From 776418a82ced3c2ed3da4b9db71a57c97295e4db Mon Sep 17 00:00:00 2001 From: eschnett Date: Wed, 24 Nov 2010 01:57:07 +0000 Subject: 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 --- param.ccl | 17 +++++++++++++++++ 1 file changed, 17 insertions(+) (limited to 'param.ccl') diff --git a/param.ccl b/param.ccl index 1c03d3f..9164f7d 100644 --- a/param.ccl +++ b/param.ccl @@ -102,6 +102,23 @@ KEYWORD overwrite_boundary "Overwrite g and K on the boundary" "exact" :: "Use boundary data from an exact solution on a trivial slice" } "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 ############### ################################################################################ -- cgit v1.2.3