modify and regenerate test suites, they now pass (trivially)

* changed spatial prolongation order to 3 (from 1), mostly because the git
  version of Carpet seems to have a bug in the order 1 code. This was what
  prevented the testsuite from running with more than one processor

* increased resolution a bit to disentangle prolongation effects from outer
  boundary effects, also needed due to higher order spatial prolongation

* the ID testsuite produces almost identical data to the previous curie release
  version when the same parfile is used. The differences are only in the outer
  boundary and go away if ADMBase is reverted to revision 60 (which changed the
  outer boundary condition for "static" evolutions)

* the evolving tests do not agree as well. Agreement gets better if McLachlan
  is reverted to the Curie release version, but even so there are grid points
  at the edge of the prolongated region that do not match. This is very likey
  due to changes in when Psi4 is computed (Curie computes in PostRestrict,
  Maxwell does not). Since the differences creep in from these points, and
  these differences increase with time I decided to regenerate the test rather
  than finding out which precise change in McLachlan/Carpet/MoL/ADMBase/Kranc
  causes these to show up.

* this requires (as of 2011-09-29 a patch to Carpet/Initialise.cc/CallRegrid)



git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/WeylScal4/trunk@82 4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843

1 files changed, 27 insertions, 23 deletions

diff --git a/test/teukolskyID.par b/test/teukolskyID.par
index 9d0230a..36aefaf 100644
--- a/test/teukolskyID.par
+++ b/test/teukolskyID.par
@@ -7,13 +7,13 @@
 ActiveThorns = "CartGrid3D CoordBase SymBase Boundary Time Carpet CarpetLib CarpetSlab CarpetRegrid2 CarpetReduce InitBase CarpetIOAscii IOUtil IOBasic ADMBase ADMMacros ADMConstraints GenericFD TmunuBase ADMCoupling StaticConformal CoordGauge SpaceMask IDLinearWaves LoopControl GSL WeylScal4 ReflectionSymmetry MoL"
 
 # Grid
-driver::ghost_size   = 1
+driver::ghost_size   = 2
 
 CartGrid3D::avoid_origin = no
 
 Carpet::max_refinement_levels = 3
-Carpet::prolongation_order_space = 1
-Carpet::prolongation_order_time = 2
+Carpet::prolongation_order_space = 3 # 1 seems to be buggy in Carpet/git
+Carpet::prolongation_order_time = 0
 Carpet::domain_from_coordbase = yes
 
 Carpet::poison_new_timelevels = "yes"
@@ -29,30 +29,30 @@ CarpetRegrid2::ensure_proper_nesting = "yes"
 CarpetRegrid2::regrid_every          = 0
 CarpetRegrid2::num_centres           = 1
 CarpetRegrid2::num_levels_1          = 3
-CarpetRegrid2::position_x_1          = -0.0625
-CarpetRegrid2::position_y_1          = -0.0625
+CarpetRegrid2::position_x_1          = -0.062
+CarpetRegrid2::position_y_1          = -0.062
 CarpetRegrid2::position_z_1          = 0
 CarpetRegrid2::radius_1[1]           = 1.0
 CarpetRegrid2::radius_1[2]           = 0.5
 
 Grid::type = "coordbase"
 
-CoordBase::xmin = -2.0625
-CoordBase::ymin = -2.0625
+CoordBase::xmin = -2.062
+CoordBase::ymin = -2.062
 CoordBase::zmin = 0.0
-CoordBase::xmax = 1.9375
-CoordBase::ymax = 1.9375
+CoordBase::xmax = 1.938
+CoordBase::ymax = 1.938
 CoordBase::zmax = 2.0
-CoordBase::dx   =  0.5
-CoordBase::dy   =  0.5
-CoordBase::dz   =  0.5
+CoordBase::dx   =  0.25
+CoordBase::dy   =  0.25
+CoordBase::dz   =  0.25
 
-CoordBase::boundary_size_x_lower = 1
-CoordBase::boundary_size_y_lower = 1
-CoordBase::boundary_size_z_lower = 1
-CoordBase::boundary_size_x_upper = 1
-CoordBase::boundary_size_y_upper = 1
-CoordBase::boundary_size_z_upper = 1
+CoordBase::boundary_size_x_lower = 2
+CoordBase::boundary_size_y_lower = 2
+CoordBase::boundary_size_z_lower = 2
+CoordBase::boundary_size_x_upper = 2
+CoordBase::boundary_size_y_upper = 2
+CoordBase::boundary_size_z_upper = 2
 
 CoordBase::boundary_shiftout_z_lower = 1
 CoordBase::boundary_shiftout_z_upper = 0
@@ -77,9 +77,9 @@ ADMBase::dtshift_evolution_method = "static"
 # Evolution
 ADMBase::evolution_method = "static"
 
-ADMBase::lapse_timelevels = 3
-ADMBase::shift_timelevels = 3
-ADMBase::metric_timelevels = 3
+ADMBase::lapse_timelevels = 1
+ADMBase::shift_timelevels = 1
+ADMBase::metric_timelevels = 1
 
 # Output
 IO::out_dir = $parfile
@@ -87,6 +87,10 @@ IO::out_fileinfo = "none"
 IO::parfile_write = "no"
 
 IOASCII::out1D_every =  1
-IOASCII::out1d_zline_x = -0.0625
-IOASCII::out1d_zline_y = -0.0625
+IOASCII::out1d_yline_x = -0.062
+IOASCII::out1d_xline_y = -0.062
+IOASCII::out1d_zline_x = -0.062
+IOASCII::out1d_zline_y = -0.062
 IOASCII::out1D_vars  = "weylscal4::psi4r weylscal4::psi4i"
+
+WeylScal4::timelevels = 1