Register symmetries of Ricci.

Apply symmetry and flat boundary conditions to Ricci.
Added testsuite for Ricci output.

Fixes CactusEinstein/1165.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/ADMAnalysis/trunk@14 80bd93c7-81bc-453a-9e3f-619c0b4f6fe4

31 files changed, 1095 insertions, 27 deletions

diff --git a/interface.ccl b/interface.ccl
index ec01397..434b6f4 100644
--- a/interface.ccl
+++ b/interface.ccl
@@ -5,6 +5,8 @@ implements: ADMAnalysis
 
 inherits: ADMBase, StaticConformal, Grid, ADMMacros
 
+USES INCLUDE HEADER: Symmetry.h
+
 # For evaltrK
 REAL trace_of_K TYPE = GF
 {
diff --git a/schedule.ccl b/schedule.ccl
index 59a1cfc..9f119ce 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -3,9 +3,14 @@
 
 SCHEDULE ADMAnalysis_ParamCheck AT CCTK_PARAMCHECK
 {
-  LANG:C
+  LANG: C
 } "Check that the metric_type is recognised"
 
+SCHEDULE ADMAnalysis_RegisterSymmetry AT CCTK_BASEGRID
+{
+  LANG: C
+} "Register symmetry of Ricci tensor and scalar"
+
 
 SCHEDULE ADMAnalysis_EvaltrK AT CCTK_ANALYSIS
 {
@@ -31,6 +36,6 @@ SCHEDULE ADMAnalysis_CurvCartToSphere AT CCTK_ANALYSIS
 SCHEDULE ADMAnalysis_Ricci AT CCTK_ANALYSIS
 {
   STORAGE: ricci_tensor, ricci_scalar, detofg
-  LANG:C
+  LANG: C
   TRIGGERS: ricci_tensor, ricci_scalar
 } "Calculation of the Ricci tensor"
diff --git a/src/Ricci.c b/src/Ricci.c
index e311ce5..74b11a9 100644
--- a/src/Ricci.c
+++ b/src/Ricci.c
@@ -10,9 +10,7 @@
  @@*/
 
 #include "cctk.h"
-
 #include "cctk_Arguments.h"
-#include "cctk_Parameters.h"
 
 #include "ADMAnalysis.h"
 
@@ -68,17 +66,16 @@ void ADMAnalysis_Ricci(CCTK_ARGUMENTS)
 {
 
   DECLARE_CCTK_ARGUMENTS;
-  DECLARE_CCTK_PARAMETERS;
 
-  CCTK_INT i,j,k, ijk, di, dj, dk;
+  CCTK_INT i,j,k, ijk, di, dj, dk, err, sw[3];
 
 #include "CactusEinstein/ADMMacros/src/macro/RICCI_declare.h"
 
-  for (k = 0; k < cctk_lsh[2]; k++) 
+  for (k = 1; k < cctk_lsh[2]-1; k++) 
   {
-    for (j = 0; j < cctk_lsh[1]; j++) 
+    for (j = 1; j < cctk_lsh[1]-1; j++) 
     {
-      for (i = 0; i < cctk_lsh[0]; i++) 
+      for (i = 1; i < cctk_lsh[0]-1; i++) 
       {
     
         ijk = CCTK_GFINDEX3D(cctkGH,i,j,k);
@@ -101,6 +98,24 @@ void ADMAnalysis_Ricci(CCTK_ARGUMENTS)
 
 #include "CactusEinstein/ADMMacros/src/macro/RICCI_undefine.h"
 
+  /* Execute Symmetry Conditions */
+  CartSymGN(cctkGH,"ADMAnalysis::ricci_scalar");
+  CartSymGN(cctkGH,"ADMAnalysis::ricci_tensor");
+
+  /* Apply Flat Boundary Condition */
+  sw[0] = 1; sw[1] = 1; sw[2] = 1;
+  err = BndFlatGN(cctkGH,sw,"ADMAnalysis::ricci_scalar");
+  err = BndFlatGN(cctkGH,sw,"ADMAnalysis::ricci_tensor");
+  if (err < 0)
+  {
+    CCTK_WARN(2,"Error in applying flat boundary condition to Ricci tensor");
+  }
+
+  /* WARNING:  Only flat boundary conditions are used here.  If the
+   *  value of the Ricci tensor is going to be used for something other
+   *  than output, then the boundaries should be handled more properly.
+   */
+
   ADMAnalysis_Trace(cctk_lsh, gxx, gxy, gxz, gyy, gyz, gzz, 
                     Ricci11, Ricci12, Ricci13, Ricci22, Ricci23, Ricci33,
                     Ricci, detg);
@@ -112,19 +127,3 @@ void ADMAnalysis_Ricci(CCTK_ARGUMENTS)
 /********************************************************************
  *********************     Local Routines   *************************
  ********************************************************************/
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
diff --git a/src/Symmetry.c b/src/Symmetry.c
new file mode 100644
index 0000000..96fa1ad
--- /dev/null
+++ b/src/Symmetry.c
@@ -0,0 +1,81 @@
+ /*@@
+   @file      Symmetry.c
+   @date      Tue Jul 30 11:15:23 CEST 2002
+   @author    David Rideout
+   @desc 
+   Symmetry registration stuff for ADMAnalysis
+   @enddesc
+   @version $Header$
+ @@*/
+
+#include "cctk.h"
+
+#include "cctk_Arguments.h"
+#include "Symmetry.h"
+
+static const char *rcsid = "$Header$";
+
+CCTK_FILEVERSION(CactusEinstein_ADMAnalysis_Symmetry_c)
+
+/********************************************************************
+ *********************     Local Data Types   ***********************
+ ********************************************************************/
+
+/********************************************************************
+ ********************* Local Routine Prototypes *********************
+ ********************************************************************/
+
+/********************************************************************
+ ***************** Scheduled Routine Prototypes *********************
+ ********************************************************************/
+
+
+/********************************************************************
+ ********************* Other Routine Prototypes *********************
+ ********************************************************************/
+
+/********************************************************************
+ *********************     Local Data   *****************************
+ ********************************************************************/
+
+/********************************************************************
+ *********************     External Routines   **********************
+
+ ********************************************************************/
+
+ /*@@
+   @routine    ADMAnalysis_RegisterSymmetry
+   @date       Tue Jul 30 11:31:50 CEST 2002
+   @author     David Rideout
+   @desc 
+   Scheduled routine to register symmetry of Ricci.
+   @enddesc 
+   @calls     
+   @calledby   
+   @history 
+ 
+   @endhistory 
+ 
+ @@*/
+void ADMAnalysis_RegisterSymmetry(CCTK_ARGUMENTS)
+{
+  DECLARE_CCTK_ARGUMENTS;
+  int sym[3];
+
+  /*CCTK_INFO("Registering symmetries for ricci_scalar and ricci_tensor");*/
+  sym[0] = 1; sym[1] = 1; sym[2] = 1;
+  SetCartSymVN(cctkGH, sym, "ADMAnalysis::Ricci");
+  SetCartSymVN(cctkGH, sym, "ADMAnalysis::Ricci11");
+  SetCartSymVN(cctkGH, sym, "ADMAnalysis::Ricci22");
+  SetCartSymVN(cctkGH, sym, "ADMAnalysis::Ricci33");
+  sym[0] = -1; sym[1] = -1; sym[2] = 1;
+  SetCartSymVN(cctkGH, sym, "ADMAnalysis::Ricci12");
+  sym[0] = -1; sym[1] = 1; sym[2] = -1;
+  SetCartSymVN(cctkGH, sym, "ADMAnalysis::Ricci13");
+  sym[0] = 1; sym[1] = -1; sym[2] = -1;
+  SetCartSymVN(cctkGH, sym, "ADMAnalysis::Ricci23");
+}
+
+/********************************************************************
+ *********************     Local Routines   *************************
+ ********************************************************************/
diff --git a/src/make.code.defn b/src/make.code.defn
index ec1603b..c3ae6cc 100644
--- a/src/make.code.defn
+++ b/src/make.code.defn
@@ -2,8 +2,7 @@
 # $Header$
 
 # Source files in this directory
-SRCS = ParamCheck.c Analysis.c CartToSphere.c Ricci.c Trace.c
+SRCS = ParamCheck.c Symmetry.c Analysis.c CartToSphere.c Ricci.c Trace.c
 
 # Subdirectories containing source files
 SUBDIRS = 
-
diff --git a/test/4BL_holes.par b/test/4BL_holes.par
new file mode 100644
index 0000000..b02631a
--- /dev/null
+++ b/test/4BL_holes.par
@@ -0,0 +1,53 @@
+!DESC "1+log sliced 4 BL black holes, ADM, double leapfrog"
+
+# This testsuite is taken from IDAnalyticBH/test/test_bl.par
+# but with modified output to make more use of this thorn.
+
+ActiveThorns = "Boundary Time ADMBase ADM ADMAnalysis CoordGauge ADMMacros SpaceMask StaticConformal CartGrid3D pugh pughreduce  pughslab IDanalyticBH IOUtil iobasic IOASCII"
+
+driver::global_nx          = 20
+driver::global_ny          = 20
+driver::global_nz          = 20
+
+grid::type                 = "byspacing"
+grid::dxyz                 =  0.8
+grid::domain               = "full"
+
+time::dtfac                =  0.25
+
+admbase::metric_type      = "static conformal"
+admbase::evolution_method = "ADM"
+admbase::lapse_evolution_method = "1+log"
+admbase::initial_data     = "bl_bh"
+
+idanalyticbh::bl_nbh       = 4
+idanalyticbh::bl_x0_1      = 3.5
+idanalyticbh::bl_y0_1      = 3.5
+idanalyticbh::bl_z0_1      = 3.5
+idanalyticbh::bl_M_1       = 0.1
+idanalyticbh::bl_x0_2      = 3.5
+idanalyticbh::bl_y0_2      = -3.5
+idanalyticbh::bl_z0_2      = 3.1
+idanalyticbh::bl_M_2       = 0.2
+idanalyticbh::bl_x0_3      = -3.5
+idanalyticbh::bl_y0_3      = 3.2
+idanalyticbh::bl_z0_3      = -3.5
+idanalyticbh::bl_M_3       = 0.3
+idanalyticbh::bl_x0_4      = -3.5
+idanalyticbh::bl_y0_4      = -3.3
+idanalyticbh::bl_z0_4      = 3.5
+idanalyticbh::bl_M_4       = 0.4
+
+cactus::cctk_itlast        = 4
+adm::bound                 = "flat"
+adm::method                = "leapfrog"
+
+IO::out_dir              = "4BL_holes"
+IO::out_fileinfo         = "none"
+IOBasic::outInfo_every   = 1
+IOBasic::outInfo_vars    = "admbase::gxx admanalysis::Ricci"
+IOBasic::outScalar_every = 2
+IOBasic::outScalar_vars  = "admanalysis::grr admanalysis::ricci admanalysis::ricci22 admanalysis::ricci13"
+IOBasic::outScalar_vars  = "admanalysis::ricci admanalysis::ricci22 admanalysis::ricci13"
+IO::out_every            = 2
+IOASCII::out1D_vars     = "admanalysis::grr admanalysis::ricci admanalysis::ricci22"
diff --git a/test/4BL_holes/4BL_holes.par b/test/4BL_holes/4BL_holes.par
new file mode 100644
index 0000000..b02631a
--- /dev/null
+++ b/test/4BL_holes/4BL_holes.par
@@ -0,0 +1,53 @@
+!DESC "1+log sliced 4 BL black holes, ADM, double leapfrog"
+
+# This testsuite is taken from IDAnalyticBH/test/test_bl.par
+# but with modified output to make more use of this thorn.
+
+ActiveThorns = "Boundary Time ADMBase ADM ADMAnalysis CoordGauge ADMMacros SpaceMask StaticConformal CartGrid3D pugh pughreduce  pughslab IDanalyticBH IOUtil iobasic IOASCII"
+
+driver::global_nx          = 20
+driver::global_ny          = 20
+driver::global_nz          = 20
+
+grid::type                 = "byspacing"
+grid::dxyz                 =  0.8
+grid::domain               = "full"
+
+time::dtfac                =  0.25
+
+admbase::metric_type      = "static conformal"
+admbase::evolution_method = "ADM"
+admbase::lapse_evolution_method = "1+log"
+admbase::initial_data     = "bl_bh"
+
+idanalyticbh::bl_nbh       = 4
+idanalyticbh::bl_x0_1      = 3.5
+idanalyticbh::bl_y0_1      = 3.5
+idanalyticbh::bl_z0_1      = 3.5
+idanalyticbh::bl_M_1       = 0.1
+idanalyticbh::bl_x0_2      = 3.5
+idanalyticbh::bl_y0_2      = -3.5
+idanalyticbh::bl_z0_2      = 3.1
+idanalyticbh::bl_M_2       = 0.2
+idanalyticbh::bl_x0_3      = -3.5
+idanalyticbh::bl_y0_3      = 3.2
+idanalyticbh::bl_z0_3      = -3.5
+idanalyticbh::bl_M_3       = 0.3
+idanalyticbh::bl_x0_4      = -3.5
+idanalyticbh::bl_y0_4      = -3.3
+idanalyticbh::bl_z0_4      = 3.5
+idanalyticbh::bl_M_4       = 0.4
+
+cactus::cctk_itlast        = 4
+adm::bound                 = "flat"
+adm::method                = "leapfrog"
+
+IO::out_dir              = "4BL_holes"
+IO::out_fileinfo         = "none"
+IOBasic::outInfo_every   = 1
+IOBasic::outInfo_vars    = "admbase::gxx admanalysis::Ricci"
+IOBasic::outScalar_every = 2
+IOBasic::outScalar_vars  = "admanalysis::grr admanalysis::ricci admanalysis::ricci22 admanalysis::ricci13"
+IOBasic::outScalar_vars  = "admanalysis::ricci admanalysis::ricci22 admanalysis::ricci13"
+IO::out_every            = 2
+IOASCII::out1D_vars     = "admanalysis::grr admanalysis::ricci admanalysis::ricci22"
diff --git a/test/4BL_holes/Ricci13_maximum.xg b/test/4BL_holes/Ricci13_maximum.xg
new file mode 100644
index 0000000..364946e
--- /dev/null
+++ b/test/4BL_holes/Ricci13_maximum.xg
@@ -0,0 +1,4 @@
+"Ricci13 v time
+0.0000000000000	2.5690402281264
+0.4000000000000	19.8893828270625
+0.8000000000000	12.5001146663148
diff --git a/test/4BL_holes/Ricci13_minimum.xg b/test/4BL_holes/Ricci13_minimum.xg
new file mode 100644
index 0000000..3a8d566
--- /dev/null
+++ b/test/4BL_holes/Ricci13_minimum.xg
@@ -0,0 +1,4 @@
+"Ricci13 v time
+0.0000000000000	-11.2693247293089
+0.4000000000000	-0.4648984705784
+0.8000000000000	-17.7893101514920
diff --git a/test/4BL_holes/Ricci13_norm1.xg b/test/4BL_holes/Ricci13_norm1.xg
new file mode 100644
index 0000000..79d752e
--- /dev/null
+++ b/test/4BL_holes/Ricci13_norm1.xg
@@ -0,0 +1,4 @@
+"Ricci13 v time
+0.0000000000000	0.0054408630140
+0.4000000000000	0.0067513049648
+0.8000000000000	0.0080817389600
diff --git a/test/4BL_holes/Ricci13_norm2.xg b/test/4BL_holes/Ricci13_norm2.xg
new file mode 100644
index 0000000..427f02b
--- /dev/null
+++ b/test/4BL_holes/Ricci13_norm2.xg
@@ -0,0 +1,4 @@
+"Ricci13 v time
+0.0000000000000	0.1308908534750
+0.4000000000000	0.2245754661059
+0.8000000000000	0.2441820052355
diff --git a/test/4BL_holes/Ricci22_3D_diagonal.xg b/test/4BL_holes/Ricci22_3D_diagonal.xg
new file mode 100644
index 0000000..7e1a638
--- /dev/null
+++ b/test/4BL_holes/Ricci22_3D_diagonal.xg
@@ -0,0 +1,69 @@
+
+
+"Time = 0.0000000000000
+-13.1635861375235		-0.0000747251238
+-11.7779454914684		-0.0000747251238
+-10.3923048454133		-0.0001261945344
+-9.0066641993582		-0.0002140016621
+-7.6210235533031		-0.0003552682038
+-6.2353829072480		-0.0005458121333
+-4.8497422611929		-0.0006890535227
+-3.4641016151378		-0.0005470719845
+-2.0784609690827		-0.0000454067838
+-0.6928203230276		0.0002884896026
+0.6928203230276		-0.0000856513648
+2.0784609690827		-0.0008959271032
+3.4641016151378		-0.0014019357606
+4.8497422611929		-0.0014973908407
+6.2353829072480		0.0060478745896
+7.6210235533031		-0.0004995056562
+9.0066641993582		-0.0003618856178
+10.3923048454133		-0.0002372710513
+11.7779454914684		-0.0001555561428
+13.1635861375235		-0.0001555561428
+
+
+"Time = 0.4000000000000
+-13.1635861375235		-0.0000789158525
+-11.7779454914684		-0.0000789158525
+-10.3923048454133		-0.0001261802045
+-9.0066641993582		-0.0002139992717
+-7.6210235533031		-0.0003553923357
+-6.2353829072480		-0.0005463724005
+-4.8497422611929		-0.0006902631246
+-3.4641016151378		-0.0005488641353
+-2.0784609690827		-0.0000508217721
+-0.6928203230276		0.0002780004002
+0.6928203230276		-0.0000944967705
+2.0784609690827		-0.0009131302196
+3.4641016151378		-0.0016772213592
+4.8497422611929		-0.0099479221735
+6.2353829072480		53.0582473241441
+7.6210235533031		-0.0055514603281
+9.0066641993582		-0.0004942349640
+10.3923048454133		-0.0002475393594
+11.7779454914684		-0.0001994784055
+13.1635861375235		-0.0001994784055
+
+
+"Time = 0.8000000000000
+-13.1635861375235		-0.0000909861730
+-11.7779454914684		-0.0000909861730
+-10.3923048454133		-0.0001272696582
+-9.0066641993582		-0.0002139900188
+-7.6210235533031		-0.0003557637025
+-6.2353829072480		-0.0005480681498
+-4.8497422611929		-0.0006939732115
+-3.4641016151378		-0.0005544559643
+-2.0784609690827		-0.0000672710971
+-0.6928203230276		0.0002463721185
+0.6928203230276		-0.0001211660220
+2.0784609690827		-0.0009660532816
+3.4641016151378		-0.0023154136682
+4.8497422611929		0.1690835190105
+6.2353829072480		46.7974935824031
+7.6210235533031		0.2043091989584
+9.0066641993582		-0.0008734701396
+10.3923048454133		-0.0002900867765
+11.7779454914684		-0.0003329825999
+13.1635861375235		-0.0003329825999
diff --git a/test/4BL_holes/Ricci22_maximum.xg b/test/4BL_holes/Ricci22_maximum.xg
new file mode 100644
index 0000000..93d1aa4
--- /dev/null
+++ b/test/4BL_holes/Ricci22_maximum.xg
@@ -0,0 +1,4 @@
+"Ricci22 v time
+0.0000000000000	2.2944275501020
+0.4000000000000	53.0582473241441
+0.8000000000000	46.7974935824031
diff --git a/test/4BL_holes/Ricci22_minimum.xg b/test/4BL_holes/Ricci22_minimum.xg
new file mode 100644
index 0000000..f1bc61a
--- /dev/null
+++ b/test/4BL_holes/Ricci22_minimum.xg
@@ -0,0 +1,4 @@
+"Ricci22 v time
+0.0000000000000	-6.1484696959425
+0.4000000000000	-1.6462268592568
+0.8000000000000	-0.3445701352369
diff --git a/test/4BL_holes/Ricci22_norm1.xg b/test/4BL_holes/Ricci22_norm1.xg
new file mode 100644
index 0000000..793672c
--- /dev/null
+++ b/test/4BL_holes/Ricci22_norm1.xg
@@ -0,0 +1,4 @@
+"Ricci22 v time
+0.0000000000000	0.0065607245458
+0.4000000000000	0.0126297980747
+0.8000000000000	0.0180908154083
diff --git a/test/4BL_holes/Ricci22_norm2.xg b/test/4BL_holes/Ricci22_norm2.xg
new file mode 100644
index 0000000..2827175
--- /dev/null
+++ b/test/4BL_holes/Ricci22_norm2.xg
@@ -0,0 +1,4 @@
+"Ricci22 v time
+0.0000000000000	0.0986483090900
+0.4000000000000	0.5977173888673
+0.8000000000000	0.6160958971265
diff --git a/test/4BL_holes/Ricci22_x_[10][10].xg b/test/4BL_holes/Ricci22_x_[10][10].xg
new file mode 100644
index 0000000..a905788
--- /dev/null
+++ b/test/4BL_holes/Ricci22_x_[10][10].xg
@@ -0,0 +1,69 @@
+
+
+"Time = 0.0000000000000
+-7.6000000000000		0.0001204281063
+-6.8000000000000		0.0001204281063
+-6.0000000000000		-0.0004636437054
+-5.2000000000000		-0.0012341070942
+-4.4000000000000		-0.0019639045417
+-3.6000000000000		-0.0023057952984
+-2.8000000000000		-0.0020584777243
+-2.0000000000000		-0.0013774811881
+-1.2000000000000		-0.0006406800837
+-0.4000000000000		-0.0001636864179
+0.4000000000000		-0.0000856513648
+1.2000000000000		-0.0004035255836
+2.0000000000000		-0.0009899805703
+2.8000000000000		-0.0015736879119
+3.6000000000000		-0.0018253792442
+4.4000000000000		-0.0016007395187
+5.2000000000000		-0.0010613721593
+6.0000000000000		-0.0004874429830
+6.8000000000000		-0.0000530654040
+7.6000000000000		-0.0000530654040
+
+
+"Time = 0.4000000000000
+-7.6000000000000		0.0000603213882
+-6.8000000000000		0.0000603213882
+-6.0000000000000		-0.0004743684468
+-5.2000000000000		-0.0012405104358
+-4.4000000000000		-0.0019552389342
+-3.6000000000000		-0.0022853662671
+-2.8000000000000		-0.0020462976382
+-2.0000000000000		-0.0013817013965
+-1.2000000000000		-0.0006519839227
+-0.4000000000000		-0.0001741050349
+0.4000000000000		-0.0000944967705
+1.2000000000000		-0.0004117243622
+2.0000000000000		-0.0009932113417
+2.8000000000000		-0.0015640520363
+3.6000000000000		-0.0018084135809
+4.4000000000000		-0.0015932078610
+5.2000000000000		-0.0010650808334
+6.0000000000000		-0.0004935507820
+6.8000000000000		-0.0000981366034
+7.6000000000000		-0.0000981366034
+
+
+"Time = 0.8000000000000
+-7.6000000000000		-0.0001116428405
+-6.8000000000000		-0.0001116428405
+-6.0000000000000		-0.0005163742887
+-5.2000000000000		-0.0012593805526
+-4.4000000000000		-0.0019293898315
+-3.6000000000000		-0.0022261539568
+-2.8000000000000		-0.0020103185588
+-2.0000000000000		-0.0013938789209
+-1.2000000000000		-0.0006858997300
+-0.4000000000000		-0.0002054813378
+0.4000000000000		-0.0001211660220
+1.2000000000000		-0.0004364930829
+2.0000000000000		-0.0010025208473
+2.8000000000000		-0.0015353484323
+3.6000000000000		-0.0017596212330
+4.4000000000000		-0.0015705143991
+5.2000000000000		-0.0010760889575
+6.0000000000000		-0.0005201751214
+6.8000000000000		-0.0002254238787
+7.6000000000000		-0.0002254238787
diff --git a/test/4BL_holes/Ricci22_y_[10][10].xg b/test/4BL_holes/Ricci22_y_[10][10].xg
new file mode 100644
index 0000000..8bbd6d0
--- /dev/null
+++ b/test/4BL_holes/Ricci22_y_[10][10].xg
@@ -0,0 +1,69 @@
+
+
+"Time = 0.0000000000000
+-7.6000000000000		-0.0004361340935
+-6.8000000000000		-0.0004361340935
+-6.0000000000000		0.0006765824207
+-5.2000000000000		0.0022993792007
+-4.4000000000000		0.0040068030025
+-3.6000000000000		0.0049485143397
+-2.8000000000000		0.0045149330434
+-2.0000000000000		0.0030262358883
+-1.2000000000000		0.0013845234309
+-0.4000000000000		0.0002749014029
+0.4000000000000		-0.0000856513648
+1.2000000000000		0.0002088770424
+2.0000000000000		0.0008854002941
+2.8000000000000		0.0015381566151
+3.6000000000000		0.0017646431341
+4.4000000000000		0.0014289485251
+5.2000000000000		0.0007494439371
+6.0000000000000		0.0000570581498
+6.8000000000000		-0.0004521306149
+7.6000000000000		-0.0004521306149
+
+
+"Time = 0.4000000000000
+-7.6000000000000		-0.0005405078627
+-6.8000000000000		-0.0005405078627
+-6.0000000000000		0.0006646916906
+-5.2000000000000		0.0022906314551
+-4.4000000000000		0.0040174830730
+-3.6000000000000		0.0049780053892
+-2.8000000000000		0.0045334143041
+-2.0000000000000		0.0030219749923
+-1.2000000000000		0.0013715015051
+-0.4000000000000		0.0002632484394
+0.4000000000000		-0.0000944967705
+1.2000000000000		0.0002029048873
+2.0000000000000		0.0008855095339
+2.8000000000000		0.0015479404647
+3.6000000000000		0.0017785367833
+4.4000000000000		0.0014358301831
+5.2000000000000		0.0007479888866
+6.0000000000000		0.0000529734237
+6.8000000000000		-0.0005036670975
+7.6000000000000		-0.0005036670975
+
+
+"Time = 0.8000000000000
+-7.6000000000000		-0.0008345167944
+-6.8000000000000		-0.0008345167944
+-6.0000000000000		0.0006105496525
+-5.2000000000000		0.0022638233536
+-4.4000000000000		0.0040493603665
+-3.6000000000000		0.0050677322930
+-2.8000000000000		0.0045895638236
+-2.0000000000000		0.0030087929516
+-1.2000000000000		0.0013321627858
+-0.4000000000000		0.0002281649485
+0.4000000000000		-0.0001211660220
+1.2000000000000		0.0001847876793
+2.0000000000000		0.0008856332769
+2.8000000000000		0.0015774357327
+3.6000000000000		0.0018204569030
+4.4000000000000		0.0014562256683
+5.2000000000000		0.0007432498700
+6.0000000000000		0.0000310284683
+6.8000000000000		-0.0006487840191
+7.6000000000000		-0.0006487840191
diff --git a/test/4BL_holes/Ricci22_z_[10][10].xg b/test/4BL_holes/Ricci22_z_[10][10].xg
new file mode 100644
index 0000000..3f0169c
--- /dev/null
+++ b/test/4BL_holes/Ricci22_z_[10][10].xg
@@ -0,0 +1,69 @@
+
+
+"Time = 0.0000000000000
+-7.6000000000000		0.0006954128388
+-6.8000000000000		0.0006954128388
+-6.0000000000000		0.0006400941167
+-5.2000000000000		0.0005256084273
+-4.4000000000000		0.0004005648388
+-3.6000000000000		0.0003473675467
+-2.8000000000000		0.0004065256111
+-2.0000000000000		0.0005138052218
+-1.2000000000000		0.0005371608407
+-0.4000000000000		0.0003582472897
+0.4000000000000		-0.0000856513648
+1.2000000000000		-0.0007781174685
+2.0000000000000		-0.0015778895440
+2.8000000000000		-0.0022033680700
+3.6000000000000		-0.0023623800612
+4.4000000000000		-0.0019866763657
+5.2000000000000		-0.0012943561909
+6.0000000000000		-0.0005860260387
+6.8000000000000		-0.0000434708775
+7.6000000000000		-0.0000434708775
+
+
+"Time = 0.4000000000000
+-7.6000000000000		0.0006875995113
+-6.8000000000000		0.0006875995113
+-6.0000000000000		0.0006359045500
+-5.2000000000000		0.0005230607417
+-4.4000000000000		0.0004040347160
+-3.6000000000000		0.0003554264020
+-2.8000000000000		0.0004108945622
+-2.0000000000000		0.0005110644630
+-1.2000000000000		0.0005306778062
+-0.4000000000000		0.0003504936324
+0.4000000000000		-0.0000944967705
+1.2000000000000		-0.0007858320779
+2.0000000000000		-0.0015767152425
+2.8000000000000		-0.0021868551735
+3.6000000000000		-0.0023407027931
+4.4000000000000		-0.0019775595587
+5.2000000000000		-0.0012986472709
+6.0000000000000		-0.0005940614417
+6.8000000000000		-0.0000984854691
+7.6000000000000		-0.0000984854691
+
+
+"Time = 0.8000000000000
+-7.6000000000000		0.0006646858395
+-6.8000000000000		0.0006646858395
+-6.0000000000000		0.0006224349207
+-5.2000000000000		0.0005158748183
+-4.4000000000000		0.0004145266499
+-3.6000000000000		0.0003785685687
+-2.8000000000000		0.0004238664521
+-2.0000000000000		0.0005033422437
+-1.2000000000000		0.0005113482067
+-0.4000000000000		0.0003271343856
+0.4000000000000		-0.0001211660220
+1.2000000000000		-0.0008089166035
+2.0000000000000		-0.0015726427515
+2.8000000000000		-0.0021382032706
+3.6000000000000		-0.0022783314585
+4.4000000000000		-0.0019502315291
+5.2000000000000		-0.0013111835131
+6.0000000000000		-0.0006274797336
+6.8000000000000		-0.0002549512732
+7.6000000000000		-0.0002549512732
diff --git a/test/4BL_holes/Ricci_3D_diagonal.xg b/test/4BL_holes/Ricci_3D_diagonal.xg
new file mode 100644
index 0000000..43f7be0
--- /dev/null
+++ b/test/4BL_holes/Ricci_3D_diagonal.xg
@@ -0,0 +1,69 @@
+
+
+"Time = 0.0000000000000
+-13.1635861375235		-0.0000000000000
+-11.7779454914684		-0.0000000000000
+-10.3923048454133		-0.0000000000000
+-9.0066641993582		-0.0000000000000
+-7.6210235533031		0.0000000000000
+-6.2353829072480		-0.0000000000000
+-4.8497422611929		0.0000000000000
+-3.4641016151378		-0.0000000000000
+-2.0784609690827		-0.0000000000000
+-0.6928203230276		-0.0000000000000
+0.6928203230276		-0.0000000000000
+2.0784609690827		0.0000000000000
+3.4641016151378		0.0000000000000
+4.8497422611929		-0.0000000000000
+6.2353829072480		0.0000000000037
+7.6210235533031		-0.0000000000000
+9.0066641993582		0.0000000000000
+10.3923048454133		0.0000000000000
+11.7779454914684		0.0000000000000
+13.1635861375235		0.0000000000000
+
+
+"Time = 0.4000000000000
+-13.1635861375235		-0.0000136974241
+-11.7779454914684		-0.0000136974241
+-10.3923048454133		0.0000001622725
+-9.0066641993582		0.0000001170628
+-7.6210235533031		-0.0000003927530
+-6.2353829072480		-0.0000019790236
+-4.8497422611929		-0.0000035362317
+-3.4641016151378		-0.0000021467627
+-2.0784609690827		-0.0000093435941
+-0.6928203230276		-0.0000255051305
+0.6928203230276		-0.0000234683420
+2.0784609690827		-0.0000493245972
+3.4641016151378		-0.0008018558746
+4.8497422611929		-0.0230359411924
+6.2353829072480		195.4472755875492
+7.6210235533031		-0.0146868874395
+9.0066641993582		-0.0003853632865
+10.3923048454133		-0.0000291265795
+11.7779454914684		-0.0001294277030
+13.1635861375235		-0.0001294277030
+
+
+"Time = 0.8000000000000
+-13.1635861375235		-0.0000530259555
+-11.7779454914684		-0.0000530259555
+-10.3923048454133		-0.0000027434771
+-9.0066641993582		0.0000004692947
+-7.6210235533031		-0.0000015675100
+-6.2353829072480		-0.0000079056714
+-4.8497422611929		-0.0000141905352
+-3.4641016151378		-0.0000088881963
+-2.0784609690827		-0.0000374787087
+-0.6928203230276		-0.0001023002122
+0.6928203230276		-0.0000941312023
+2.0784609690827		-0.0002006834298
+3.4641016151378		-0.0026303028058
+4.8497422611929		0.5264759249645
+6.2353829072480		169.4352957090077
+7.6210235533031		0.6171820964260
+9.0066641993582		-0.0014855472926
+10.3923048454133		-0.0001508671958
+11.7779454914684		-0.0005230960516
+13.1635861375235		-0.0005230960516
diff --git a/test/4BL_holes/Ricci_maximum.xg b/test/4BL_holes/Ricci_maximum.xg
new file mode 100644
index 0000000..7168f8e
--- /dev/null
+++ b/test/4BL_holes/Ricci_maximum.xg
@@ -0,0 +1,4 @@
+"Ricci v time
+0.0000000000000	0.0000000000037
+0.4000000000000	195.4472755875492
+0.8000000000000	169.4352957090077
diff --git a/test/4BL_holes/Ricci_minimum.xg b/test/4BL_holes/Ricci_minimum.xg
new file mode 100644
index 0000000..c45d0e5
--- /dev/null
+++ b/test/4BL_holes/Ricci_minimum.xg
@@ -0,0 +1,4 @@
+"Ricci v time
+0.0000000000000	-0.0000000000000
+0.4000000000000	-0.4745827436382
+0.8000000000000	-0.3512824219088
diff --git a/test/4BL_holes/Ricci_norm1.xg b/test/4BL_holes/Ricci_norm1.xg
new file mode 100644
index 0000000..217c294
--- /dev/null
+++ b/test/4BL_holes/Ricci_norm1.xg
@@ -0,0 +1,4 @@
+"Ricci v time
+0.0000000000000	0.0000000000000
+0.4000000000000	0.0321851594951
+0.8000000000000	0.0521728973688
diff --git a/test/4BL_holes/Ricci_norm2.xg b/test/4BL_holes/Ricci_norm2.xg
new file mode 100644
index 0000000..ab3c20d
--- /dev/null
+++ b/test/4BL_holes/Ricci_norm2.xg
@@ -0,0 +1,4 @@
+"Ricci v time
+0.0000000000000	0.0000000000000
+0.4000000000000	2.1988854043208
+0.8000000000000	2.1957617696270
diff --git a/test/4BL_holes/Ricci_x_[10][10].xg b/test/4BL_holes/Ricci_x_[10][10].xg
new file mode 100644
index 0000000..556fc7e
--- /dev/null
+++ b/test/4BL_holes/Ricci_x_[10][10].xg
@@ -0,0 +1,69 @@
+
+
+"Time = 0.0000000000000
+-7.6000000000000		0.0000000000000
+-6.8000000000000		0.0000000000000
+-6.0000000000000		0.0000000000000
+-5.2000000000000		-0.0000000000000
+-4.4000000000000		-0.0000000000000
+-3.6000000000000		-0.0000000000000
+-2.8000000000000		0.0000000000000
+-2.0000000000000		0.0000000000000
+-1.2000000000000		0.0000000000000
+-0.4000000000000		-0.0000000000000
+0.4000000000000		-0.0000000000000
+1.2000000000000		-0.0000000000000
+2.0000000000000		-0.0000000000000
+2.8000000000000		-0.0000000000000
+3.6000000000000		0.0000000000000
+4.4000000000000		0.0000000000000
+5.2000000000000		-0.0000000000000
+6.0000000000000		-0.0000000000000
+6.8000000000000		0.0000000000000
+7.6000000000000		0.0000000000000
+
+
+"Time = 0.4000000000000
+-7.6000000000000		-0.0002487692674
+-6.8000000000000		-0.0002487692674
+-6.0000000000000		-0.0000222584394
+-5.2000000000000		-0.0000015562741
+-4.4000000000000		0.0000515907769
+-3.6000000000000		0.0000905250189
+-2.8000000000000		0.0000628220739
+-2.0000000000000		0.0000060483928
+-1.2000000000000		-0.0000230286644
+-0.4000000000000		-0.0000256527470
+0.4000000000000		-0.0000234683420
+1.2000000000000		-0.0000221742501
+2.0000000000000		-0.0000077257407
+2.8000000000000		0.0000287272119
+3.6000000000000		0.0000493598569
+4.4000000000000		0.0000231344622
+5.2000000000000		-0.0000091309694
+6.0000000000000		-0.0000169744901
+6.8000000000000		-0.0001134798776
+7.6000000000000		-0.0001134798776
+
+
+"Time = 0.8000000000000
+-7.6000000000000		-0.0009535743579
+-6.8000000000000		-0.0009535743579
+-6.0000000000000		-0.0001320952684
+-5.2000000000000		-0.0000062069193
+-4.4000000000000		0.0002054294899
+-3.6000000000000		0.0003579732687
+-2.8000000000000		0.0002499051394
+-2.0000000000000		0.0000245947128
+-1.2000000000000		-0.0000923174828
+-0.4000000000000		-0.0001028616805
+0.4000000000000		-0.0000941312023
+1.2000000000000		-0.0000888495787
+2.0000000000000		-0.0000297381433
+2.8000000000000		0.0001148086778
+3.6000000000000		0.0001932566230
+4.4000000000000		0.0000927752724
+5.2000000000000		-0.0000360333256
+6.0000000000000		-0.0000883886717
+6.8000000000000		-0.0004341538810
+7.6000000000000		-0.0004341538810
diff --git a/test/4BL_holes/Ricci_y_[10][10].xg b/test/4BL_holes/Ricci_y_[10][10].xg
new file mode 100644
index 0000000..2ac717e
--- /dev/null
+++ b/test/4BL_holes/Ricci_y_[10][10].xg
@@ -0,0 +1,69 @@
+
+
+"Time = 0.0000000000000
+-7.6000000000000		-0.0000000000000
+-6.8000000000000		-0.0000000000000
+-6.0000000000000		-0.0000000000000
+-5.2000000000000		-0.0000000000000
+-4.4000000000000		-0.0000000000000
+-3.6000000000000		-0.0000000000000
+-2.8000000000000		-0.0000000000000
+-2.0000000000000		-0.0000000000000
+-1.2000000000000		-0.0000000000000
+-0.4000000000000		-0.0000000000000
+0.4000000000000		-0.0000000000000
+1.2000000000000		-0.0000000000000
+2.0000000000000		-0.0000000000000
+2.8000000000000		-0.0000000000000
+3.6000000000000		-0.0000000000000
+4.4000000000000		-0.0000000000000
+5.2000000000000		0.0000000000000
+6.0000000000000		0.0000000000000
+6.8000000000000		0.0000000000000
+7.6000000000000		0.0000000000000
+
+
+"Time = 0.4000000000000
+-7.6000000000000		-0.0002136955096
+-6.8000000000000		-0.0002136955096
+-6.0000000000000		-0.0000343328599
+-5.2000000000000		-0.0000249078632
+-4.4000000000000		0.0000393900466
+-3.6000000000000		0.0001041384441
+-2.8000000000000		0.0000676273496
+-2.0000000000000		-0.0000077402192
+-1.2000000000000		-0.0000356622781
+-0.4000000000000		-0.0000312818285
+0.4000000000000		-0.0000234683420
+1.2000000000000		-0.0000161281859
+2.0000000000000		0.0000018263526
+2.8000000000000		0.0000326673550
+3.6000000000000		0.0000457105919
+4.4000000000000		0.0000218002680
+5.2000000000000		-0.0000053814552
+6.0000000000000		-0.0000131351673
+6.8000000000000		-0.0001061996365
+7.6000000000000		-0.0001061996365
+
+
+"Time = 0.8000000000000
+-7.6000000000000		-0.0008178519560
+-6.8000000000000		-0.0008178519560
+-6.0000000000000		-0.0001742791709
+-5.2000000000000		-0.0000987198052
+-4.4000000000000		0.0001588966905
+-3.6000000000000		0.0004039970678
+-2.8000000000000		0.0002691189937
+-2.0000000000000		-0.0000289707729
+-1.2000000000000		-0.0001431062198
+-0.4000000000000		-0.0001255290390
+0.4000000000000		-0.0000941312023
+1.2000000000000		-0.0000646094627
+2.0000000000000		0.0000079986095
+2.8000000000000		0.0001295834267
+3.6000000000000		0.0001783004767
+4.4000000000000		0.0000869633624
+5.2000000000000		-0.0000212272109
+6.0000000000000		-0.0000717866450
+6.8000000000000		-0.0004063994741
+7.6000000000000		-0.0004063994741
diff --git a/test/4BL_holes/Ricci_z_[10][10].xg b/test/4BL_holes/Ricci_z_[10][10].xg
new file mode 100644
index 0000000..372576c
--- /dev/null
+++ b/test/4BL_holes/Ricci_z_[10][10].xg
@@ -0,0 +1,69 @@
+
+
+"Time = 0.0000000000000
+-7.6000000000000		-0.0000000000000
+-6.8000000000000		-0.0000000000000
+-6.0000000000000		-0.0000000000000
+-5.2000000000000		-0.0000000000000
+-4.4000000000000		0.0000000000000
+-3.6000000000000		0.0000000000000
+-2.8000000000000		0.0000000000000
+-2.0000000000000		0.0000000000000
+-1.2000000000000		0.0000000000000
+-0.4000000000000		-0.0000000000000
+0.4000000000000		-0.0000000000000
+1.2000000000000		-0.0000000000000
+2.0000000000000		0.0000000000000
+2.8000000000000		0.0000000000000
+3.6000000000000		-0.0000000000000
+4.4000000000000		-0.0000000000000
+5.2000000000000		-0.0000000000000
+6.0000000000000		-0.0000000000000
+6.8000000000000		-0.0000000000000
+7.6000000000000		-0.0000000000000
+
+
+"Time = 0.4000000000000
+-7.6000000000000		-0.0001024889326
+-6.8000000000000		-0.0001024889326
+-6.0000000000000		-0.0000109090607
+-5.2000000000000		-0.0000043016427
+-4.4000000000000		0.0000164219789
+-3.6000000000000		0.0000318351781
+-2.8000000000000		0.0000194196986
+-2.0000000000000		-0.0000041584501
+-1.2000000000000		-0.0000163735930
+-0.4000000000000		-0.0000203228830
+0.4000000000000		-0.0000234683420
+1.2000000000000		-0.0000198265877
+2.0000000000000		0.0000066152676
+2.8000000000000		0.0000516002358
+3.6000000000000		0.0000672129042
+4.4000000000000		0.0000312323082
+5.2000000000000		-0.0000087024466
+6.0000000000000		-0.0000211766399
+6.8000000000000		-0.0002038616387
+7.6000000000000		-0.0002038616387
+
+
+"Time = 0.8000000000000
+-7.6000000000000		-0.0003916924448
+-6.8000000000000		-0.0003916924448
+-6.0000000000000		-0.0000624396121
+-5.2000000000000		-0.0000170612695
+-4.4000000000000		0.0000655760483
+-3.6000000000000		0.0001260839623
+-2.8000000000000		0.0000776383895
+-2.0000000000000		-0.0000161790538
+-1.2000000000000		-0.0000655228781
+-0.4000000000000		-0.0000814762334
+0.4000000000000		-0.0000941312023
+1.2000000000000		-0.0000793381054
+2.0000000000000		0.0000272861363
+2.8000000000000		0.0002046894260
+3.6000000000000		0.0002637994146
+4.4000000000000		0.0001247144071
+5.2000000000000		-0.0000344905910
+6.0000000000000		-0.0001204413745
+6.8000000000000		-0.0007817705956
+7.6000000000000		-0.0007817705956
diff --git a/test/4BL_holes/grr_3D_diagonal.xg b/test/4BL_holes/grr_3D_diagonal.xg
new file mode 100644
index 0000000..6a98116
--- /dev/null
+++ b/test/4BL_holes/grr_3D_diagonal.xg
@@ -0,0 +1,69 @@
+
+
+"Time = 0.0000000000000
+-13.1635861375235		1.0000000000000
+-11.7779454914684		1.0000000000000
+-10.3923048454133		1.0000000000000
+-9.0066641993582		1.0000000000000
+-7.6210235533031		1.0000000000000
+-6.2353829072480		1.0000000000000
+-4.8497422611929		1.0000000000000
+-3.4641016151378		1.0000000000000
+-2.0784609690827		1.0000000000000
+-0.6928203230276		1.0000000000000
+0.6928203230276		1.0000000000000
+2.0784609690827		1.0000000000000
+3.4641016151378		1.0000000000000
+4.8497422611929		1.0000000000000
+6.2353829072480		1.0000000000000
+7.6210235533031		1.0000000000000
+9.0066641993582		1.0000000000000
+10.3923048454133		1.0000000000000
+11.7779454914684		1.0000000000000
+13.1635861375235		1.0000000000000
+
+
+"Time = 0.4000000000000
+-13.1635861375235		1.0000951955360
+-11.7779454914684		1.0000951955360
+-10.3923048454133		1.0001111745761
+-9.0066641993582		1.0001208692995
+-7.6210235533031		1.0001092693451
+-6.2353829072480		1.0000480952362
+-4.8497422611929		0.9999072460689
+-3.4641016151378		0.9997140902667
+-2.0784609690827		0.9996085166659
+-0.6928203230276		0.9996959606253
+0.6928203230276		0.9999044778080
+2.0784609690827		1.0002334049468
+3.4641016151378		1.0012270163223
+4.8497422611929		1.0110676342094
+6.2353829072480		2.1253455302950
+7.6210235533031		1.0057140732643
+9.0066641993582		1.0010092043593
+10.3923048454133		1.0003795088067
+11.7779454914684		1.0001999637180
+13.1635861375235		1.0001999637180
+
+
+"Time = 0.8000000000000
+-13.1635861375235		1.0003846615105
+-11.7779454914684		1.0003846615105
+-10.3923048454133		1.0004447624379
+-9.0066641993582		1.0004835613639
+-7.6210235533031		1.0004372180569
+-6.2353829072480		1.0001927040665
+-4.8497422611929		0.9996294952948
+-3.4641016151378		0.9988562300676
+-2.0784609690827		0.9984336455803
+-0.6928203230276		0.9987854982447
+0.6928203230276		0.9996200522743
+2.0784609690827		1.0009413255772
+3.4641016151378		1.0050636001104
+4.8497422611929		1.0538497273229
+6.2353829072480		2.0976489647219
+7.6210235533031		1.0268501521106
+9.0066641993582		1.0041162448038
+10.3923048454133		1.0015242829450
+11.7779454914684		1.0008258747060
+13.1635861375235		1.0008258747060
diff --git a/test/4BL_holes/grr_x_[10][10].xg b/test/4BL_holes/grr_x_[10][10].xg
new file mode 100644
index 0000000..5f68ae1
--- /dev/null
+++ b/test/4BL_holes/grr_x_[10][10].xg
@@ -0,0 +1,69 @@
+
+
+"Time = 0.0000000000000
+-7.6000000000000		1.0000000000000
+-6.8000000000000		1.0000000000000
+-6.0000000000000		1.0000000000000
+-5.2000000000000		1.0000000000000
+-4.4000000000000		1.0000000000000
+-3.6000000000000		1.0000000000000
+-2.8000000000000		1.0000000000000
+-2.0000000000000		1.0000000000000
+-1.2000000000000		1.0000000000000
+-0.4000000000000		1.0000000000000
+0.4000000000000		1.0000000000000
+1.2000000000000		1.0000000000000
+2.0000000000000		1.0000000000000
+2.8000000000000		1.0000000000000
+3.6000000000000		1.0000000000000
+4.4000000000000		1.0000000000000
+5.2000000000000		1.0000000000000
+6.0000000000000		1.0000000000000
+6.8000000000000		1.0000000000000
+7.6000000000000		1.0000000000000
+
+
+"Time = 0.4000000000000
+-7.6000000000000		1.0000058051803
+-6.8000000000000		1.0000062481189
+-6.0000000000000		0.9998564319489
+-5.2000000000000		0.9996657007687
+-4.4000000000000		0.9994883480149
+-3.6000000000000		0.9994024235733
+-2.8000000000000		0.9994494646761
+-2.0000000000000		0.9995860136920
+-1.2000000000000		0.9996966089212
+-0.4000000000000		0.9995986225336
+0.4000000000000		0.9999044778080
+1.2000000000000		1.0000667312522
+2.0000000000000		0.9999834613050
+2.8000000000000		0.9998655093659
+3.6000000000000		0.9997961076429
+4.4000000000000		0.9998144422162
+5.2000000000000		0.9998974640197
+6.0000000000000		0.9999929926990
+6.8000000000000		1.0000656030072
+7.6000000000000		1.0000669365079
+
+
+"Time = 0.8000000000000
+-7.6000000000000		1.0000522963029
+-6.8000000000000		1.0000540099857
+-6.0000000000000		0.9994275497055
+-5.2000000000000		0.9986634544247
+-4.4000000000000		0.9979508225911
+-3.6000000000000		0.9976045734015
+-2.8000000000000		0.9977940364782
+-2.0000000000000		0.9983433890828
+-1.2000000000000		0.9987875705983
+-0.4000000000000		0.9983962648543
+0.4000000000000		0.9996200522743
+1.2000000000000		1.0002685838598
+2.0000000000000		0.9999341591540
+2.8000000000000		0.9994597372605
+3.6000000000000		0.9991810018444
+4.4000000000000		0.9992564652227
+5.2000000000000		0.9995909350540
+6.0000000000000		0.9999736158751
+6.8000000000000		1.0002766391962
+7.6000000000000		1.0002819890258
diff --git a/test/4BL_holes/grr_y_[10][10].xg b/test/4BL_holes/grr_y_[10][10].xg
new file mode 100644
index 0000000..72fff59
--- /dev/null
+++ b/test/4BL_holes/grr_y_[10][10].xg
@@ -0,0 +1,69 @@
+
+
+"Time = 0.0000000000000
+-7.6000000000000		1.0000000000000
+-6.8000000000000		1.0000000000000
+-6.0000000000000		1.0000000000000
+-5.2000000000000		1.0000000000000
+-4.4000000000000		1.0000000000000
+-3.6000000000000		1.0000000000000
+-2.8000000000000		1.0000000000000
+-2.0000000000000		1.0000000000000
+-1.2000000000000		1.0000000000000
+-0.4000000000000		1.0000000000000
+0.4000000000000		1.0000000000000
+1.2000000000000		1.0000000000000
+2.0000000000000		1.0000000000000
+2.8000000000000		1.0000000000000
+3.6000000000000		1.0000000000000
+4.4000000000000		1.0000000000000
+5.2000000000000		1.0000000000000
+6.0000000000000		1.0000000000000
+6.8000000000000		1.0000000000000
+7.6000000000000		1.0000000000000
+
+
+"Time = 0.4000000000000
+-7.6000000000000		1.0000279423591
+-6.8000000000000		1.0000249030833
+-6.0000000000000		0.9998830701549
+-5.2000000000000		0.9996940043668
+-4.4000000000000		0.9995203581323
+-3.6000000000000		0.9994656099493
+-2.8000000000000		0.9995887522390
+-2.0000000000000		0.9998383666284
+-1.2000000000000		1.0001150660583
+-0.4000000000000		1.0002990188162
+0.4000000000000		0.9999044778080
+1.2000000000000		0.9999116088889
+2.0000000000000		0.9998860221316
+2.8000000000000		0.9998288539825
+3.6000000000000		0.9998035939762
+4.4000000000000		0.9998394531732
+5.2000000000000		0.9999185646811
+6.0000000000000		1.0000028022178
+6.8000000000000		1.0000666448909
+7.6000000000000		1.0000656318332
+
+
+"Time = 0.8000000000000
+-7.6000000000000		1.0001375591995
+-6.8000000000000		1.0001253793368
+-6.0000000000000		0.9995355296810
+-5.2000000000000		0.9987788712195
+-4.4000000000000		0.9980803551009
+-3.6000000000000		0.9978570137219
+-2.8000000000000		0.9983514320699
+-2.0000000000000		0.9993549730460
+-1.2000000000000		1.0004641411073
+-0.4000000000000		1.0011996388779
+0.4000000000000		0.9996200522743
+1.2000000000000		0.9996481183925
+2.0000000000000		0.9995442433108
+2.8000000000000		0.9993131856329
+3.6000000000000		0.9992112341451
+4.4000000000000		0.9993563418128
+5.2000000000000		0.9996746893905
+6.0000000000000		1.0000122410166
+6.8000000000000		1.0002797997222
+7.6000000000000		1.0002757344691
diff --git a/test/4BL_holes/grr_z_[10][10].xg b/test/4BL_holes/grr_z_[10][10].xg
new file mode 100644
index 0000000..ad6ddbf
--- /dev/null
+++ b/test/4BL_holes/grr_z_[10][10].xg
@@ -0,0 +1,69 @@
+
+
+"Time = 0.0000000000000
+-7.6000000000000		1.0000000000000
+-6.8000000000000		1.0000000000000
+-6.0000000000000		1.0000000000000
+-5.2000000000000		1.0000000000000
+-4.4000000000000		1.0000000000000
+-3.6000000000000		1.0000000000000
+-2.8000000000000		1.0000000000000
+-2.0000000000000		1.0000000000000
+-1.2000000000000		1.0000000000000
+-0.4000000000000		1.0000000000000
+0.4000000000000		1.0000000000000
+1.2000000000000		1.0000000000000
+2.0000000000000		1.0000000000000
+2.8000000000000		1.0000000000000
+3.6000000000000		1.0000000000000
+4.4000000000000		1.0000000000000
+5.2000000000000		1.0000000000000
+6.0000000000000		1.0000000000000
+6.8000000000000		1.0000000000000
+7.6000000000000		1.0000000000000
+
+
+"Time = 0.4000000000000
+-7.6000000000000		1.0000738477984
+-6.8000000000000		1.0000746738590
+-6.0000000000000		1.0000139294269
+-5.2000000000000		0.9999335668421
+-4.4000000000000		0.9998597654908
+-3.6000000000000		0.9998308035446
+-2.8000000000000		0.9998664602084
+-2.0000000000000		0.9999453985495
+-1.2000000000000		1.0000218531751
+-0.4000000000000		1.0000421488283
+0.4000000000000		0.9999044778080
+1.2000000000000		0.9998063593872
+2.0000000000000		0.9996727055706
+2.8000000000000		0.9995573668941
+3.6000000000000		0.9995275804985
+4.4000000000000		0.9996058075689
+5.2000000000000		0.9997528001224
+6.0000000000000		0.9999064286963
+6.8000000000000		1.0000260142473
+7.6000000000000		1.0000232047158
+
+
+"Time = 0.8000000000000
+-7.6000000000000		1.0003084927887
+-6.8000000000000		1.0003117716704
+-6.0000000000000		1.0000570336601
+-5.2000000000000		0.9997349437959
+-4.4000000000000		0.9994378872238
+-3.6000000000000		0.9993206181563
+-2.8000000000000		0.9994641000274
+-2.0000000000000		0.9997817333039
+-1.2000000000000		1.0000885881204
+-0.4000000000000		1.0001699135404
+0.4000000000000		0.9996200522743
+1.2000000000000		0.9992275404278
+2.0000000000000		0.9986912126658
+2.8000000000000		0.9982271205426
+3.6000000000000		0.9981071361077
+4.4000000000000		0.9984221876412
+5.2000000000000		0.9990124490126
+6.0000000000000		0.9996276248519
+6.8000000000000		1.0001285038909
+7.6000000000000		1.0001172624462