CarpetProlongateTest: Test interpolation as well

Add a routine that interpolates to a grid array, testing whether
interpolation (including time interpolation) works correctly.

Add a test case for 3rd order spatial and 2nd order temporal
interpolation.

1 files changed, 112 insertions, 0 deletions

diff --git a/CarpetExtra/CarpetProlongateTest/test/test_interp/test_interp.par b/CarpetExtra/CarpetProlongateTest/test/test_interp/test_interp.par
new file mode 100644
index 000000000..06e11e855
--- /dev/null
+++ b/CarpetExtra/CarpetProlongateTest/test/test_interp/test_interp.par
@@ -0,0 +1,112 @@
+Cactus::cctk_itlast = 512
+
+
+
+ActiveThorns = "IOUtil"
+
+IO::out_dir = $parfile
+
+
+
+ActiveThorns = "AEILocalInterp"
+
+ActiveThorns = "InitBase"
+InitBase::initial_data_setup_method = "init_single_level"
+
+ActiveThorns = "LoopControl"
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce"
+
+Carpet::domain_from_coordbase = yes
+Carpet::max_refinement_levels = 10
+
+driver::ghost_size       = 2
+Carpet::use_buffer_zones = yes
+
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time  = 2
+
+Carpet::init_each_timelevel = yes
+
+Carpet::poison_new_timelevels = yes
+CarpetLib::poison_new_memory  = yes
+
+
+
+ActiveThorns = "Boundary CartGrid3D CoordBase SymBase"
+
+CoordBase::domainsize = "minmax"
+
+CoordBase::xmin = -1.00
+CoordBase::ymin = -1.00
+CoordBase::zmin = -1.00
+CoordBase::xmax = +1.00
+CoordBase::ymax = +1.00
+CoordBase::zmax = +1.00
+CoordBase::dx   =  0.02
+CoordBase::dy   =  0.02
+CoordBase::dz   =  0.02
+
+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
+
+CartGrid3D::type = "coordbase"
+
+
+
+ActiveThorns = "SphericalSurface"
+
+
+
+ActiveThorns = "CarpetRegrid2 CarpetTracker"
+
+CarpetRegrid2::regrid_every = 0
+CarpetRegrid2::verbose      = yes
+
+CarpetRegrid2::num_centres = 1
+
+CarpetRegrid2::num_levels_1 = 2
+CarpetRegrid2::radius_1[1]  = 0.2
+
+
+
+ActiveThorns = "MoL"
+
+
+
+ActiveThorns = "CarpetProlongateTest"
+
+CarpetProlongateTest::power_x = 3
+CarpetProlongateTest::power_y = 3
+CarpetProlongateTest::power_z = 3
+CarpetProlongateTest::power_t = 2
+
+CarpetProlongateTest::interpolator_options = "order=4"
+
+
+
+ActiveThorns = "CarpetIOASCII"
+
+IOASCII::one_file_per_group = yes
+
+IOASCII::out0D_every = 1
+IOASCII::out0D_vars  = "
+        CarpetProlongateTest::interp_errornorm
+"
+
+IOASCII::out1D_every = 1
+IOASCII::out1D_vars  = "
+        CarpetProlongateTest::scalar
+        CarpetProlongateTest::interp_difference
+"
+
+IOASCII::out3D_every = 1
+IOASCII::out3D_vars  = "
+        CarpetProlongateTest::interp_difference
+"