Add test for trapezoidal integration method

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

15 files changed, 421 insertions, 0 deletions

diff --git a/test/test_trapezoidal.par b/test/test_trapezoidal.par
new file mode 100644
index 0000000..7c4f782
--- /dev/null
+++ b/test/test_trapezoidal.par
@@ -0,0 +1,92 @@
+
+ActiveThorns = "CoordBase SymBase Boundary CartGrid3d IOUtil Carpet CarpetLib  CarpetInterp AEILocalInterp InitBase Multipole LoopControl CarpetIOASCII"
+
+#############################################################
+# Grid
+#############################################################
+
+CartGrid3D::type                        = "coordbase"
+CartGrid3D::domain                      = "full"
+CartGrid3D::avoid_origin                = "no"
+
+CoordBase::domainsize                   = minmax
+CoordBase::xmin                         = -10
+CoordBase::ymin                         = -10
+CoordBase::zmin                         = -10
+CoordBase::xmax                         = 10
+CoordBase::ymax                         = 10
+CoordBase::zmax                         = 10
+CoordBase::dx                           = 0.2
+CoordBase::dy                           = 0.2
+CoordBase::dz                           = 0.2
+CoordBase::boundary_size_x_lower        = 2
+CoordBase::boundary_size_y_lower        = 2
+CoordBase::boundary_size_z_lower        = 2
+CoordBase::boundary_shiftout_x_lower    = 0
+CoordBase::boundary_shiftout_y_lower    = 0
+CoordBase::boundary_shiftout_z_lower    = 0
+CoordBase::boundary_size_x_upper        = 2
+CoordBase::boundary_size_y_upper        = 2
+CoordBase::boundary_size_z_upper        = 2
+CoordBase::boundary_shiftout_x_upper    = 0
+CoordBase::boundary_shiftout_y_upper    = 0
+CoordBase::boundary_shiftout_z_upper    = 0
+
+#############################################################
+# Carpet
+#############################################################
+
+Carpet::ghost_size                      = 3
+Carpet::domain_from_coordbase           = "yes"
+Carpet::poison_new_timelevels           = "yes"
+Carpet::check_for_poison                = "no"
+CarpetLib::poison_value                    = 113
+Carpet::init_fill_timelevels            = "yes"
+
+#############################################################
+# CarpetLib
+#############################################################
+
+CarpetLib::poison_new_memory            = "yes"
+
+#############################################################
+# Cactus
+#############################################################
+
+Cactus::terminate                     = "iteration"
+Cactus::cctk_itlast                   = 0
+
+#############################################################
+# Multipole
+#############################################################
+
+Multipole::nradii = 1
+Multipole::radius[0] = 8.0
+Multipole::variables = "Multipole::harmonic_re{sw=-2 cmplx='Multipole::harmonic_im' name='harmonic'}"
+Multipole::integration_method = "trapezoidal"
+
+Multipole::enable_test = "yes"
+Multipole::test_l = 2
+Multipole::test_m = 2
+Multipole::out_1d_every = 1
+
+#############################################################
+# Output
+#############################################################
+
+IO::out_dir                          = $parfile
+IO::out_fileinfo                     = "none"
+
+# Enabling 1D output for the test grid functions would be helpful for
+# localising any failures but it makes the tests dependent on the
+# number of processors, as CarpetIOASCII's output is dependent on
+# this.
+
+# CarpetIOASCII::out1d_vars = "Multipole::harmonics"
+# CarpetIOASCII::out1d_every = 1
+# CarpetIOASCII::out1d_x = yes
+# CarpetIOASCII::out1d_y = yes
+# CarpetIOASCII::out1d_z = yes
+
+CarpetIOASCII::out0d_vars = "Multipole::test_trapezoidal_convergence_order"
+CarpetIOASCII::out0d_every = 1
diff --git a/test/test_trapezoidal/mp_harmonic_im_r8.00.ph.asc b/test/test_trapezoidal/mp_harmonic_im_r8.00.ph.asc
new file mode 100644
index 0000000..4943e48
--- /dev/null
+++ b/test/test_trapezoidal/mp_harmonic_im_r8.00.ph.asc
@@ -0,0 +1,104 @@
+"Time = 0
+0.000000 5.341621717946576991e-17
+0.062832 0.05908272357625558163
+0.125664 0.1172335976390736084
+0.188496 0.1735357640795202605
+0.251327 0.2271010934402109549
+0.314159 0.2770848541179522506
+0.376991 0.3226990304560694023
+0.439823 0.3632238274429471225
+0.502655 0.3980205224812608522
+0.565487 0.4265400723078277556
+0.628319 0.4483329443869139563
+0.691150 0.4630553007318252146
+0.753982 0.4704749615530876006
+0.816814 0.4704749615530870455
+0.879646 0.4630553007318246594
+0.942478 0.4483329443869139563
+1.005310 0.4265400723078278111
+1.068142 0.3980205224812605191
+1.130973 0.3632238274429472891
+1.193805 0.3226990304560692913
+1.256637 0.2770848541179523616
+1.319469 0.2271010934402110937
+1.382301 0.1735357640795201495
+1.445133 0.1172335976390738305
+1.507964 0.05908272357625565102
+1.570796 -1.080937077952825899e-18
+1.633628 -0.05908272357625524857
+1.696460 -0.1172335976390735945
+1.759292 -0.1735357640795195111
+1.822124 -0.2271010934402106496
+1.884956 -0.2770848541179521396
+1.947787 -0.3226990304560689027
+2.010619 -0.3632238274429472336
+2.073451 -0.3980205224812608522
+2.136283 -0.4265400723078279777
+2.199115 -0.4483329443869139008
+2.261947 -0.4630553007318246594
+2.324779 -0.4704749615530868234
+2.387610 -0.4704749615530867679
+2.450442 -0.4630553007318249925
+2.513274 -0.4483329443869140674
+2.576106 -0.4265400723078281997
+2.638938 -0.3980205224812605191
+2.701770 -0.363223827442946956
+2.764602 -0.3226990304560695133
+2.827433 -0.2770848541179521951
+2.890265 -0.227101093440211288
+2.953097 -0.1735357640795202327
+3.015929 -0.1172335976390740525
+3.078761 -0.05908272357625561633
+3.141593 1.128851862588426379e-16
+3.204425 0.05908272357625510285
+3.267256 0.1172335976390735252
+3.330088 0.1735357640795197331
+3.392920 0.2271010934402107606
+3.455752 0.2770848541179523061
+3.518584 0.3226990304560693468
+3.581416 0.3632238274429465119
+3.644247 0.3980205224812606302
+3.707079 0.4265400723078277001
+3.769911 0.4483329443869141784
+3.832743 0.4630553007318256031
+3.895575 0.4704749615530871565
+3.958407 0.470474961553087101
+4.021239 0.4630553007318254366
+4.084070 0.4483329443869144004
+4.146902 0.4265400723078277556
+4.209734 0.3980205224812608522
+4.272566 0.3632238274429464564
+4.335398 0.3226990304560687917
+4.398230 0.2770848541179524727
+4.461062 0.2271010934402105108
+4.523893 0.1735357640795198719
+4.586725 0.1172335976390739554
+4.649557 0.05908272357625504734
+4.712389 2.284473287243558895e-16
+4.775221 -0.05908272357625498489
+4.838053 -0.1172335976390735807
+4.900885 -0.1735357640795197609
+4.963716 -0.2271010934402110382
+5.026548 -0.2770848541179522506
+5.089380 -0.3226990304560685141
+5.152212 -0.3632238274429464009
+5.215044 -0.3980205224812609632
+5.277876 -0.4265400723078279221
+5.340708 -0.4483329443869139008
+5.403539 -0.4630553007318252146
+5.466371 -0.470474961553087101
+5.529203 -0.4704749615530869344
+5.592035 -0.463055300731824937
+5.654867 -0.4483329443869135678
+5.717699 -0.4265400723078281997
+5.780530 -0.3980205224812614628
+5.843362 -0.363223827442946845
+5.906194 -0.3226990304560690692
+5.969026 -0.2770848541179523616
+6.031858 -0.22710109344021151
+6.094690 -0.1735357640795197609
+6.157522 -0.1172335976390737888
+6.220353 -0.05908272357625506122
+6.283185 -1.976266587454604297e-16
+
+
diff --git a/test/test_trapezoidal/mp_harmonic_im_r8.00.th.asc b/test/test_trapezoidal/mp_harmonic_im_r8.00.th.asc
new file mode 100644
index 0000000..7a0ece3
--- /dev/null
+++ b/test/test_trapezoidal/mp_harmonic_im_r8.00.th.asc
@@ -0,0 +1,54 @@
+"Time = 0
+0.000000 0.6307831305050434212
+0.062832 0.006414478478045328516
+0.125664 2.926186227276146829e-16
+0.188496 2.386143598765005219e-16
+0.251327 -1.420084232264331917e-16
+0.314159 1.777240527155209573e-16
+0.376991 1.638154941593066095e-16
+0.439823 1.580236523128021832e-16
+0.502655 9.684566415599272793e-17
+0.565487 -3.469487216420912885e-17
+0.628319 1.011041002333707997e-16
+0.691150 7.999484602692621821e-17
+0.753982 5.341621717946576991e-17
+0.816814 4.255698807620960161e-17
+0.879646 -2.755703360697257134e-17
+0.942478 3.312157109350618348e-17
+1.005310 9.164156058357568613e-18
+1.068142 4.136430182637461761e-17
+1.130973 3.93568926162343637e-17
+1.193805 2.959892296258195764e-17
+1.256637 2.515951604708110849e-17
+1.319469 3.916167830229968497e-17
+1.382301 1.059617470841231802e-17
+1.445133 6.098417063879056444e-18
+1.507964 1.027594448571300023e-17
+1.570796 -3.523247651629027143e-17
+1.633628 -1.456220359444229049e-18
+1.696460 -1.410063837397148702e-17
+1.759292 5.613015692432480741e-18
+1.822124 6.635796753069970284e-18
+1.884956 6.568542235803732263e-18
+1.947787 9.855341525364593679e-18
+2.010619 -5.508866192607108089e-18
+2.073451 4.622356468968703956e-18
+2.136283 6.611848247478812561e-18
+2.199115 2.858520867375647277e-18
+2.261947 2.590500913413465802e-19
+2.324779 1.520091510550632901e-18
+2.387610 -1.713262760772048037e-20
+2.450442 1.718382086966605172e-18
+2.513274 9.414189862246460786e-19
+2.576106 5.036919495096394639e-19
+2.638938 2.847075820349658673e-19
+2.701770 3.309989484742896642e-19
+2.764602 7.09384461930837558e-20
+2.827433 5.000736099493162766e-20
+2.890265 4.760554507445765016e-20
+2.953097 -1.805583750081199166e-20
+3.015929 5.343533318497654634e-21
+3.078761 6.3364941679916884e-22
+3.141593 -7.22781544709519842e-37
+
+
diff --git a/test/test_trapezoidal/mp_harmonic_l0_m0_r8.00.asc b/test/test_trapezoidal/mp_harmonic_l0_m0_r8.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_trapezoidal/mp_harmonic_l0_m0_r8.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_trapezoidal/mp_harmonic_l1_m-1_r8.00.asc b/test/test_trapezoidal/mp_harmonic_l1_m-1_r8.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_trapezoidal/mp_harmonic_l1_m-1_r8.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_trapezoidal/mp_harmonic_l1_m0_r8.00.asc b/test/test_trapezoidal/mp_harmonic_l1_m0_r8.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_trapezoidal/mp_harmonic_l1_m0_r8.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_trapezoidal/mp_harmonic_l1_m1_r8.00.asc b/test/test_trapezoidal/mp_harmonic_l1_m1_r8.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_trapezoidal/mp_harmonic_l1_m1_r8.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_trapezoidal/mp_harmonic_l2_m-1_r8.00.asc b/test/test_trapezoidal/mp_harmonic_l2_m-1_r8.00.asc
new file mode 100644
index 0000000..7041720
--- /dev/null
+++ b/test/test_trapezoidal/mp_harmonic_l2_m-1_r8.00.asc
@@ -0,0 +1 @@
+0.000000 -3.084588393705386578e-17 1.720378076789101466e-17
diff --git a/test/test_trapezoidal/mp_harmonic_l2_m-2_r8.00.asc b/test/test_trapezoidal/mp_harmonic_l2_m-2_r8.00.asc
new file mode 100644
index 0000000..9ccaf00
--- /dev/null
+++ b/test/test_trapezoidal/mp_harmonic_l2_m-2_r8.00.asc
@@ -0,0 +1 @@
+0.000000 1.138927479537339823e-10 1.118247678211310438e-17
diff --git a/test/test_trapezoidal/mp_harmonic_l2_m0_r8.00.asc b/test/test_trapezoidal/mp_harmonic_l2_m0_r8.00.asc
new file mode 100644
index 0000000..28ecad8
--- /dev/null
+++ b/test/test_trapezoidal/mp_harmonic_l2_m0_r8.00.asc
@@ -0,0 +1 @@
+0.000000 -1.972076544340019278e-17 7.195607036871429069e-08
diff --git a/test/test_trapezoidal/mp_harmonic_l2_m1_r8.00.asc b/test/test_trapezoidal/mp_harmonic_l2_m1_r8.00.asc
new file mode 100644
index 0000000..1f86d6f
--- /dev/null
+++ b/test/test_trapezoidal/mp_harmonic_l2_m1_r8.00.asc
@@ -0,0 +1 @@
+0.000000 9.233534510396670281e-17 3.678958802688812425e-16
diff --git a/test/test_trapezoidal/mp_harmonic_l2_m2_r8.00.asc b/test/test_trapezoidal/mp_harmonic_l2_m2_r8.00.asc
new file mode 100644
index 0000000..8e64289
--- /dev/null
+++ b/test/test_trapezoidal/mp_harmonic_l2_m2_r8.00.asc
@@ -0,0 +1 @@
+0.000000 0.9991533008868566101 7.11347051830939968e-18
diff --git a/test/test_trapezoidal/mp_harmonic_re_r8.00.ph.asc b/test/test_trapezoidal/mp_harmonic_re_r8.00.ph.asc
new file mode 100644
index 0000000..2c4f799
--- /dev/null
+++ b/test/test_trapezoidal/mp_harmonic_re_r8.00.ph.asc
@@ -0,0 +1,104 @@
+"Time = 0
+0.000000 0.4714051344795946452
+0.062832 0.467687868496530601
+0.125664 0.4565949764204048078
+0.188496 0.4383014236215818649
+0.251327 0.4130954585389419731
+0.314159 0.3813747754991172467
+0.376991 0.3436396201486748292
+0.439823 0.3004850309575128331
+0.502655 0.2525915434611012422
+0.565487 0.2007146483930815628
+0.628319 0.1456722117101460767
+0.691150 0.08833261385994624348
+0.753982 0.0295998309385308786
+0.816814 -0.02959983093853089942
+0.879646 -0.08833261385994582715
+0.942478 -0.1456722117101461045
+1.005310 -0.2007146483930816183
+1.068142 -0.2525915434611012977
+1.130973 -0.300485030957512611
+1.193805 -0.3436396201486749402
+1.256637 -0.3813747754991173022
+1.319469 -0.4130954585389421396
+1.382301 -0.4383014236215819759
+1.445133 -0.4565949764204046968
+1.507964 -0.4676878684965307675
+1.570796 -0.4714051344795952003
+1.633628 -0.4676878684965312116
+1.696460 -0.4565949764204050298
+1.759292 -0.4383014236215813653
+1.822124 -0.4130954585389423062
+1.884956 -0.3813747754991171912
+1.947787 -0.3436396201486755508
+2.010619 -0.3004850309575126666
+2.073451 -0.2525915434611009647
+2.136283 -0.2007146483930813963
+2.199115 -0.145672211710146271
+2.261947 -0.08833261385994620185
+2.324779 -0.02959983093853062186
+2.387610 0.02959983093853071207
+2.450442 0.08833261385994579939
+2.513274 0.1456722117101459657
+2.576106 0.2007146483930811465
+2.638938 0.2525915434611014088
+2.701770 0.3004850309575123335
+2.764602 0.3436396201486751623
+2.827433 0.3813747754991176908
+2.890265 0.4130954585389419731
+2.953097 0.4383014236215818649
+3.015929 0.4565949764204046968
+3.078761 0.4676878684965314892
+3.141593 0.4714051344795943677
+3.204425 0.4676878684965311006
+3.267256 0.4565949764204046413
+3.330088 0.4383014236215818649
+3.392920 0.4130954585389420841
+3.455752 0.3813747754991171912
+3.518584 0.3436396201486757729
+3.581416 0.3004850309575129441
+3.644247 0.2525915434611011867
+3.707079 0.2007146483930816183
+3.769911 0.145672211710146271
+3.832743 0.08833261385994643777
+3.895575 0.02959983093853109717
+3.958407 -0.02959983093853056288
+4.021239 -0.0883326138599464239
+4.084070 -0.1456722117101458269
+4.146902 -0.2007146483930817016
+4.209734 -0.2525915434611011312
+4.272566 -0.3004850309575127221
+4.335398 -0.3436396201486753843
+4.398230 -0.3813747754991170802
+4.461062 -0.4130954585389422506
+4.523893 -0.4383014236215819204
+4.586725 -0.4565949764204048633
+4.649557 -0.4676878684965311561
+4.712389 -0.4714051344795950338
+4.775221 -0.4676878684965316002
+4.838053 -0.4565949764204047523
+4.900885 -0.4383014236215819759
+4.963716 -0.4130954585389418066
+5.026548 -0.3813747754991175798
+5.089380 -0.3436396201486756063
+5.152212 -0.3004850309575128331
+5.215044 -0.2525915434611007981
+5.277876 -0.200714648393081424
+5.340708 -0.1456722117101463265
+5.403539 -0.08833261385994690962
+5.466371 -0.02959983093853164535
+5.529203 0.02959983093853053859
+5.592035 0.088332613859946299
+5.654867 0.1456722117101458269
+5.717699 0.2007146483930807301
+5.780530 0.2525915434611006316
+5.843362 0.3004850309575120559
+5.906194 0.3436396201486754398
+5.969026 0.3813747754991170247
+6.031858 0.4130954585389420286
+6.094690 0.4383014236215816428
+6.157522 0.4565949764204046968
+6.220353 0.4676878684965309896
+6.283185 0.4714051344795947562
+
+
diff --git a/test/test_trapezoidal/mp_harmonic_re_r8.00.th.asc b/test/test_trapezoidal/mp_harmonic_re_r8.00.th.asc
new file mode 100644
index 0000000..5a48b56
--- /dev/null
+++ b/test/test_trapezoidal/mp_harmonic_re_r8.00.th.asc
@@ -0,0 +1,54 @@
+"Time = 0
+0.000000 -8.80289179457185799e-16
+0.062832 0.6231245530661584375
+0.125664 0.6258190228933131483
+0.188496 0.6196596992337820575
+0.251327 0.6111215679566899261
+0.314159 0.6002881616211543925
+0.376991 0.5872649752392626876
+0.439823 0.5721780328572182439
+0.502655 0.5551722253737308987
+0.565487 0.5364093860695589022
+0.628319 0.5160661561144331833
+0.691150 0.4943316785539432767
+0.753982 0.4714051344795946452
+0.816814 0.447493156296473682
+0.879646 0.4228072006721167453
+0.942478 0.3975608424363436533
+1.005310 0.37196712173057489
+1.068142 0.3462359033535465924
+1.130973 0.3205713329354391572
+1.193805 0.2951694104757784931
+1.256637 0.2702157403654934509
+1.319469 0.2458834314661353659
+1.382301 0.2223312435844497503
+1.445133 0.1997019777222116266
+1.507964 0.1781211248076795539
+1.570796 0.1576957826262607443
+1.633628 0.1385139262975075558
+1.696460 0.1206438899337402004
+1.759292 0.1041342775866162806
+1.822124 0.08901404797365959631
+1.884956 0.07529303859782114638
+1.947787 0.06296265788326334267
+2.010619 0.05199693903018270252
+2.073451 0.04235381504000218511
+2.136283 0.03397662104341544337
+2.199115 0.02679582193746246274
+2.261947 0.02073089969599720139
+2.324779 0.0156923880855150584
+2.387610 0.01158401689597711795
+2.450442 0.008304916835505299355
+2.513274 0.005751875269321123606
+2.576106 0.003821566914013105852
+2.638938 0.002412751990280704532
+2.701770 0.001428390934720448913
+2.764602 0.0007776521592554753309
+2.827433 0.000377755173091246974
+2.890265 0.0001556483610775961923
+2.953097 4.947535271276211051e-05
+3.015929 9.805208296189968497e-06
+3.078761 6.131081706940525447e-07
+3.141593 -1.077203552463963929e-22
+
+
diff --git a/test/test_trapezoidal/test_trapezoidal_convergence_order..asc b/test/test_trapezoidal/test_trapezoidal_convergence_order..asc
new file mode 100644
index 0000000..32939d6
--- /dev/null
+++ b/test/test_trapezoidal/test_trapezoidal_convergence_order..asc
@@ -0,0 +1,4 @@
+# 0D ASCII output created by CarpetIOASCII
+#
+0	0	0 0 0	0 0 0	0	0 0 0	2.00042739440297
+