Add test case for extracting multiple radii

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

41 files changed, 1159 insertions, 0 deletions

diff --git a/test/test_rads.par b/test/test_rads.par
new file mode 100644
index 0000000..726e81b
--- /dev/null
+++ b/test/test_rads.par
@@ -0,0 +1,92 @@
+
+ActiveThorns = "CoordBase SymBase Boundary CartGrid3d IOUtil Carpet CarpetLib  CarpetInterp AEILocalInterp InitBase Multipole"
+
+#############################################################
+# 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                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::poison_new_timelevels           = "yes"
+Carpet::check_for_poison                = "no"
+Carpet::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 = 3
+Multipole::radius[0] = 8.0
+Multipole::radius[1] = 6.0
+Multipole::radius[2] = 4.0
+Multipole::variables = "Multipole::harmonic_re{sw=-2 cmplx='Multipole::harmonic_im' name='harmonic'}"
+Multipole::l_max = 2
+
+Multipole::enable_test = "yes"
+Multipole::test_mode_proportional_to_r = "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
diff --git a/test/test_rads/mp_harmonic_im_r4.00.ph.asc b/test/test_rads/mp_harmonic_im_r4.00.ph.asc
new file mode 100644
index 0000000..373a9d6
--- /dev/null
+++ b/test/test_rads/mp_harmonic_im_r4.00.ph.asc
@@ -0,0 +1,104 @@
+"Time = 0
+0.000000 -7.401262412693975668e-17
+0.062832 0.2363324706322694757
+0.125664 0.4689329301629398539
+0.188496 0.6941366327500196887
+0.251327 0.908397208004283252
+0.314159 1.108337011388746962
+0.376991 1.290797316344458157
+0.439823 1.45289168829120241
+0.502655 1.592078014559861643
+0.565487 1.706165034348026133
+0.628319 1.793331857477059144
+0.691150 1.852225826015087229
+0.753982 1.881903335672080901
+0.816814 1.881903335672081345
+0.879646 1.852225826015088339
+0.942478 1.793331857477059588
+1.005310 1.706165034348027465
+1.068142 1.592078014559864529
+1.130973 1.452891688291203964
+1.193805 1.290797316344458157
+1.256637 1.108337011388747184
+1.319469 0.9083972080042843622
+1.382301 0.6941366327500200217
+1.445133 0.4689329301629400204
+1.507964 0.2363324706322695867
+1.570796 6.405154908074140523e-16
+1.633628 -0.2363324706322705027
+1.696460 -0.4689329301629376889
+1.759292 -0.6941366327500188005
+1.822124 -0.9083972080042839181
+1.884956 -1.108337011388750071
+1.947787 -1.290797316344459489
+2.010619 -1.452891688291201744
+2.073451 -1.592078014559864751
+2.136283 -1.706165034348028131
+2.199115 -1.793331857477057811
+2.261947 -1.852225826015088561
+2.324779 -1.88190333567208401
+2.387610 -1.881903335672080901
+2.450442 -1.852225826015087673
+2.513274 -1.793331857477061586
+2.576106 -1.706165034348026577
+2.638938 -1.592078014559861865
+2.701770 -1.452891688291202632
+2.764602 -1.29079731634445749
+2.827433 -1.108337011388747184
+2.890265 -0.908397208004283252
+2.953097 -0.6941366327500195776
+3.015929 -0.4689329301629353575
+3.078761 -0.2363324706322678659
+3.141593 -1.198810537017405911e-16
+3.204425 0.2363324706322707247
+3.267256 0.4689329301629380775
+3.330088 0.6941366327500188005
+3.392920 0.9083972080042829189
+3.455752 1.108337011388749183
+3.518584 1.290797316344457935
+3.581416 1.452891688291202188
+3.644247 1.592078014559862753
+3.707079 1.706165034348028131
+3.769911 1.793331857477061364
+3.832743 1.852225826015088783
+3.895575 1.881903335672082678
+3.958407 1.881903335672080235
+4.021239 1.852225826015087007
+4.084070 1.793331857477061364
+4.146902 1.706165034348026799
+4.209734 1.592078014559862309
+4.272566 1.45289168829120241
+4.335398 1.290797316344458379
+4.398230 1.108337011388751625
+4.461062 0.9083972080042822528
+4.523893 0.6941366327500189115
+4.586725 0.4689329301629360236
+4.649557 0.2363324706322711688
+4.712389 -9.147249194367437519e-16
+4.775221 -0.2363324706322693369
+4.838053 -0.4689329301629397428
+4.900885 -0.6941366327500239075
+4.963716 -0.908397208004283585
+5.026548 -1.108337011388748738
+5.089380 -1.290797316344457046
+5.152212 -1.4528916882912013
+5.215044 -1.592078014559863863
+5.277876 -1.706165034348028797
+5.340708 -1.793331857477062474
+5.403539 -1.852225826015087673
+5.466371 -1.881903335672080679
+5.529203 -1.881903335672078681
+5.592035 -1.852225826015087895
+5.654867 -1.793331857477059366
+5.717699 -1.706165034348028353
+5.780530 -1.592078014559863641
+5.843362 -1.452891688291201744
+5.906194 -1.290797316344458157
+5.969026 -1.108337011388746296
+6.031858 -0.9083972080042839181
+6.094690 -0.6941366327500184674
+6.157522 -0.4689329301629361346
+6.220353 -0.23633247063227103
+6.283185 2.446529326830952138e-15
+
+
diff --git a/test/test_rads/mp_harmonic_im_r4.00.th.asc b/test/test_rads/mp_harmonic_im_r4.00.th.asc
new file mode 100644
index 0000000..61531ba
--- /dev/null
+++ b/test/test_rads/mp_harmonic_im_r4.00.th.asc
@@ -0,0 +1,54 @@
+"Time = 0
+0.000000 2.523132522020167023
+0.062832 -0.2480546203736004318
+0.125664 0.02572324522335723077
+0.188496 1.300412972816408185e-15
+0.251327 -8.290482298185427189e-16
+0.314159 1.422905590373744042e-15
+0.376991 9.627294281313936891e-16
+0.439823 1.346828507873481372e-15
+0.502655 9.215938591053693621e-16
+0.565487 4.716979192346364955e-17
+0.628319 5.371537648620685327e-16
+0.691150 6.644518394028439732e-16
+0.753982 -7.401262412693975668e-17
+0.816814 7.306700155701662743e-16
+0.879646 1.96819601879617785e-17
+0.942478 6.4369229079478495e-16
+1.005310 7.048140315923532797e-16
+1.068142 2.298918245670697738e-16
+1.130973 3.951743798364942992e-16
+1.193805 4.253622622013715869e-16
+1.256637 3.019031988576083198e-16
+1.319469 2.176841464672257058e-16
+1.382301 -1.196040100677527523e-16
+1.445133 1.30652474834588205e-16
+1.507964 1.288455465056744061e-16
+1.570796 1.11826556979572894e-16
+1.633628 -1.487574209416853533e-17
+1.696460 8.632343551358583468e-17
+1.759292 9.693782735106190257e-17
+1.822124 3.833750546884154189e-17
+1.884956 3.194027942863962038e-18
+1.947787 5.188203904292530496e-17
+2.010619 2.951166179570749264e-17
+2.073451 3.581355224780174593e-17
+2.136283 1.661976823479674027e-17
+2.199115 1.9608615389292225e-17
+2.261947 3.259204477469943187e-17
+2.324779 2.125012219464427539e-17
+2.387610 1.229276775444637264e-17
+2.450442 2.740191622875059393e-18
+2.513274 1.061506189990507119e-17
+2.576106 4.12574668807967695e-18
+2.638938 2.121301176700787072e-18
+2.701770 5.196466700977786436e-19
+2.764602 1.479411516049145903e-18
+2.827433 9.68975177808027074e-19
+2.890265 -2.097518135615728137e-19
+2.953097 -9.541411057598607838e-20
+3.015929 3.419382822504113259e-20
+3.078761 6.141183624723773964e-21
+3.141593 -5.501569712789865975e-20
+
+
diff --git a/test/test_rads/mp_harmonic_im_r6.00.ph.asc b/test/test_rads/mp_harmonic_im_r6.00.ph.asc
new file mode 100644
index 0000000..85b13f5
--- /dev/null
+++ b/test/test_rads/mp_harmonic_im_r6.00.ph.asc
@@ -0,0 +1,104 @@
+"Time = 0
+0.000000 8.845964161748810866e-16
+0.062832 0.3544952786343023021
+0.125664 0.7034009541834183032
+0.188496 1.041214574194973652
+0.251327 1.36260656874174102
+0.314159 1.662508367570813217
+0.376991 1.936192150544441271
+0.439823 2.179342528910793941
+0.502655 2.388122556580707823
+0.565487 2.559241177622988861
+0.628319 2.689998927592172961
+0.691150 2.778332346497143224
+0.753982 2.822849635654231637
+0.816814 2.822849635654230749
+0.879646 2.778332346497144112
+0.942478 2.689998927592171185
+1.005310 2.559241177622988861
+1.068142 2.388122556580675848
+1.130973 2.179342528910796162
+1.193805 1.936192150544441715
+1.256637 1.662508367570813661
+1.319469 1.362606568741739022
+1.382301 1.041214574194972764
+1.445133 0.7034009541834186363
+1.507964 0.3544952786343024687
+1.570796 9.065659625476588288e-16
+1.633628 -0.3544952786343045781
+1.696460 -0.7034009541834195245
+1.759292 -1.041214574194974318
+1.822124 -1.362606568741739022
+1.884956 -1.662508367570814549
+1.947787 -1.936192150544441715
+2.010619 -2.179342528910794829
+2.073451 -2.38812255658067718
+2.136283 -2.559241177622988417
+2.199115 -2.689998927592174294
+2.261947 -2.77833234649715255
+2.324779 -2.822849635654211209
+2.387610 -2.822849635654230749
+2.450442 -2.778332346497146332
+2.513274 -2.689998927592172517
+2.576106 -2.559241177622987529
+2.638938 -2.388122556580681177
+2.701770 -2.179342528910788612
+2.764602 -1.93619215054443683
+2.827433 -1.662508367570811885
+2.890265 -1.362606568741741242
+2.953097 -1.041214574194969433
+3.015929 -0.7034009541834185253
+3.078761 -0.3544952786343018025
+3.141593 5.646914740910456145e-16
+3.204425 0.3544952786343023021
+3.267256 0.7034009541834180812
+3.330088 1.041214574194972542
+3.392920 1.362606568741740354
+3.455752 1.662508367570812773
+3.518584 1.936192150544445045
+3.581416 2.179342528910793941
+3.644247 2.388122556580678957
+3.707079 2.55924117762299419
+3.769911 2.689998927592168076
+3.832743 2.778332346497148553
+3.895575 2.822849635654221867
+3.958407 2.822849635654224976
+4.021239 2.778332346497146332
+4.084070 2.689998927592176514
+4.146902 2.559241177622981311
+4.209734 2.388122556580681621
+4.272566 2.179342528910790389
+4.335398 1.936192150544435941
+4.398230 1.662508367570813883
+4.461062 1.362606568741736801
+4.523893 1.041214574194973652
+4.586725 0.7034009541834155277
+4.649557 0.3544952786343045781
+4.712389 -1.156438725869027658e-15
+4.775221 -0.3544952786343038564
+4.838053 -0.7034009541834197465
+4.900885 -1.041214574194974096
+4.963716 -1.362606568741741464
+5.026548 -1.662508367570810774
+5.089380 -1.936192150544445045
+5.152212 -2.179342528910795718
+5.215044 -2.38812255658068473
+5.277876 -2.559241177622987973
+5.340708 -2.689998927592173406
+5.403539 -2.778332346497148109
+5.466371 -2.822849635654231637
+5.529203 -2.822849635654230305
+5.592035 -2.778332346497151217
+5.654867 -2.689998927592168076
+5.717699 -2.55924117762298442
+5.780530 -2.388122556580678957
+5.843362 -2.179342528910794829
+5.906194 -1.936192150544441493
+5.969026 -1.662508367570812107
+6.031858 -1.362606568741737467
+6.094690 -1.041214574194973208
+6.157522 -0.7034009541834156387
+6.220353 -0.3544952786343039119
+6.283185 2.446532740057408798e-15
+
+
diff --git a/test/test_rads/mp_harmonic_im_r6.00.th.asc b/test/test_rads/mp_harmonic_im_r6.00.th.asc
new file mode 100644
index 0000000..985ceee
--- /dev/null
+++ b/test/test_rads/mp_harmonic_im_r6.00.th.asc
@@ -0,0 +1,54 @@
+"Time = 0
+0.000000 3.784698783030249647
+0.062832 -0.05867611814867000619
+0.125664 9.980367879620083493e-16
+0.188496 4.836350568642748723e-15
+0.251327 2.345867367225603486e-15
+0.314159 2.164474559169993504e-15
+0.376991 2.351912289018611647e-15
+0.439823 -1.389321621009536033e-16
+0.502655 9.779789910089601505e-16
+0.565487 1.166733578558169415e-15
+0.628319 3.337346579419405878e-16
+0.691150 5.592788098842174113e-16
+0.753982 8.845964161748810866e-16
+0.816814 7.054036462425777245e-16
+0.879646 3.902081629518370634e-16
+0.942478 3.42701115138936474e-16
+1.005310 4.78945246718037378e-16
+1.068142 2.583707567632333068e-16
+1.130973 -2.025534730464371624e-17
+1.193805 1.895718701686133983e-16
+1.256637 2.949708720990815899e-16
+1.319469 2.459404920232426121e-16
+1.382301 2.113812588681842964e-16
+1.445133 -7.762333238577484763e-17
+1.507964 6.355535500769766477e-16
+1.570796 9.007234476487628386e-17
+1.633628 -4.313837884850812957e-17
+1.696460 6.424144629611896116e-17
+1.759292 1.300430028219536752e-16
+1.822124 9.469609087768473701e-17
+1.884956 6.925300632092745011e-17
+1.947787 1.62589964865188628e-16
+2.010619 5.538940923328626015e-17
+2.073451 3.06477084594907695e-17
+2.136283 1.703723395243381606e-17
+2.199115 2.294866843239971664e-17
+2.261947 -4.050214362245610545e-18
+2.324779 1.67706721814622975e-17
+2.387610 8.065815125815240363e-18
+2.450442 -5.114677309030565691e-18
+2.513274 -1.532689298504537983e-19
+2.576106 4.24746229250417635e-18
+2.638938 5.531742865870244717e-18
+2.701770 2.336701641655587213e-18
+2.764602 1.387624927722005422e-18
+2.827433 9.248123155517364725e-19
+2.890265 6.947150699519010796e-19
+2.953097 2.528196414160398378e-19
+3.015929 4.034572764754988926e-20
+3.078761 -1.121249133514266495e-21
+3.141593 -1.72734299268512402e-20
+
+
diff --git a/test/test_rads/mp_harmonic_im_r8.00.ph.asc b/test/test_rads/mp_harmonic_im_r8.00.ph.asc
new file mode 100644
index 0000000..005d5b7
--- /dev/null
+++ b/test/test_rads/mp_harmonic_im_r8.00.ph.asc
@@ -0,0 +1,104 @@
+"Time = 0
+0.000000 3.472115407824081654e-16
+0.062832 0.4726618019463603715
+0.125664 0.9378689567934566096
+0.188496 1.388286184243222587
+0.251327 1.816808870154915434
+0.314159 2.216679070880364311
+0.376991 2.581592212405786402
+0.439823 2.905790726107930322
+0.502655 3.184164147459177752
+0.565487 3.412320634332897562
+0.628319 3.586663500855614917
+0.691150 3.704442263462494278
+0.753982 3.763799576318072493
+0.816814 3.763799576318068496
+0.879646 3.704442263462494722
+0.942478 3.586663500855615805
+1.005310 3.41232063433289845
+1.068142 3.184164147459176419
+1.130973 2.905790726107930766
+1.193805 2.58159221240578507
+1.256637 2.216679070880365199
+1.319469 1.816808870154916988
+1.382301 1.388286184243221921
+1.445133 0.9378689567934563875
+1.507964 0.4726618019463570408
+1.570796 -2.608493726006150761e-17
+1.633628 -0.472661801946357818
+1.696460 -0.9378689567934566096
+1.759292 -1.388286184243223031
+1.822124 -1.816808870154911881
+1.884956 -2.216679070880366531
+1.947787 -2.581592212405784625
+2.010619 -2.905790726107928545
+2.073451 -3.184164147459178196
+2.136283 -3.412320634332898894
+2.199115 -3.586663500855617137
+2.261947 -3.704442263462492946
+2.324779 -3.763799576318064055
+2.387610 -3.763799576318069828
+2.450442 -3.704442263462496499
+2.513274 -3.586663500855617137
+2.576106 -3.412320634332899338
+2.638938 -3.184164147459175531
+2.701770 -2.905790726107929878
+2.764602 -2.581592212405784181
+2.827433 -2.216679070880361646
+2.890265 -1.816808870154914768
+2.953097 -1.388286184243222809
+3.015929 -0.9378689567934580529
+3.078761 -0.4726618019463565967
+3.141593 2.828760169978631935e-15
+3.204425 0.4726618019463585951
+3.267256 0.9378689567934570537
+3.330088 1.388286184243219257
+3.392920 1.816808870154918765
+3.455752 2.216679070880366975
+3.518584 2.581592212405790399
+3.581416 2.905790726107925437
+3.644247 3.184164147459177752
+3.707079 3.412320634332901559
+3.769911 3.586663500855625131
+3.832743 3.704442263462496943
+3.895575 3.763799576318069384
+3.958407 3.76379957631806894
+4.021239 3.704442263462497831
+4.084070 3.58666350085562291
+4.146902 3.412320634332899782
+4.209734 3.184164147459180416
+4.272566 2.905790726107924993
+4.335398 2.58159221240578729
+4.398230 2.216679070880366087
+4.461062 1.816808870154911437
+4.523893 1.388286184243221033
+4.586725 0.9378689567934513915
+4.649557 0.4726618019463559861
+4.712389 1.541953057049712655e-15
+4.775221 -0.4726618019463632581
+4.838053 -0.9378689567934559435
+4.900885 -1.38828618424322503
+4.963716 -1.816808870154915878
+5.026548 -2.21667907088036209
+5.089380 -2.581592212405787734
+5.152212 -2.905790726107930766
+5.215044 -3.18416414745918086
+5.277876 -3.412320634332902891
+5.340708 -3.586663500855618913
+5.403539 -3.704442263462496499
+5.466371 -3.763799576318070716
+5.529203 -3.763799576318068052
+5.592035 -3.70444226346249561
+5.654867 -3.586663500855617137
+5.717699 -3.412320634332899338
+5.780530 -3.184164147459186189
+5.843362 -2.905790726107925437
+5.906194 -2.581592212405785514
+5.969026 -2.216679070880362978
+6.031858 -1.816808870154911881
+6.094690 -1.388286184243219701
+6.157522 -0.9378689567934528348
+6.220353 -0.4726618019463564857
+6.283185 6.360959921110405284e-15
+
+
diff --git a/test/test_rads/mp_harmonic_im_r8.00.th.asc b/test/test_rads/mp_harmonic_im_r8.00.th.asc
new file mode 100644
index 0000000..ce3b81f
--- /dev/null
+++ b/test/test_rads/mp_harmonic_im_r8.00.th.asc
@@ -0,0 +1,54 @@
+"Time = 0
+0.000000 5.04626504404033227
+0.062832 0.05121456776014572665
+0.125664 3.725280847919165093e-15
+0.188496 2.05134023716644306e-15
+0.251327 -7.369714291519664029e-16
+0.314159 1.36932023392800103e-15
+0.376991 8.701382433310238401e-16
+0.439823 1.229339936299604151e-15
+0.502655 8.164745875824769803e-16
+0.565487 -2.118643647607437957e-16
+0.628319 5.036756239298998587e-16
+0.691150 7.623988744447549555e-16
+0.753982 3.472115407824081654e-16
+0.816814 4.405029649936108356e-16
+0.879646 4.776832573215842252e-16
+0.942478 3.700689835978649222e-16
+1.005310 -1.318568378411203561e-16
+1.068142 3.929033768675659405e-16
+1.130973 3.218447565677357695e-16
+1.193805 2.715013020910437942e-16
+1.256637 2.424638812734975442e-16
+1.319469 3.005975144276718949e-16
+1.382301 1.952325340163290623e-16
+1.445133 1.269808009268706646e-16
+1.507964 1.083378524380504138e-16
+1.570796 8.910988680428007735e-17
+1.633628 1.547636890587835784e-16
+1.696460 2.438993712677536739e-18
+1.759292 5.013229567452844107e-17
+1.822124 8.969833483855360034e-17
+1.884956 5.410928171474172311e-17
+1.947787 8.018181683091965715e-17
+2.010619 1.508813990596288782e-17
+2.073451 2.124607007949434441e-17
+2.136283 2.324979050360401287e-17
+2.199115 1.988133437726511766e-17
+2.261947 2.725534200886170096e-17
+2.324779 1.478838924484028583e-17
+2.387610 1.409594681164944945e-17
+2.450442 1.139628579327599753e-17
+2.513274 5.693442069264047822e-18
+2.576106 4.016347149417074138e-18
+2.638938 4.344093593422070616e-18
+2.701770 1.929694084034253861e-18
+2.764602 -2.261447074092380337e-18
+2.827433 1.525503799472493377e-18
+2.890265 -2.553102147588907258e-19
+2.953097 3.686876459396992572e-20
+3.015929 4.603861186098072422e-20
+3.078761 9.124113053196154538e-21
+3.141593 -7.504662048703381121e-21
+
+
diff --git a/test/test_rads/mp_harmonic_l0_m0_r4.00.asc b/test/test_rads/mp_harmonic_l0_m0_r4.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l0_m0_r4.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_rads/mp_harmonic_l0_m0_r6.00.asc b/test/test_rads/mp_harmonic_l0_m0_r6.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l0_m0_r6.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_rads/mp_harmonic_l0_m0_r8.00.asc b/test/test_rads/mp_harmonic_l0_m0_r8.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l0_m0_r8.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_rads/mp_harmonic_l1_m-1_r4.00.asc b/test/test_rads/mp_harmonic_l1_m-1_r4.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l1_m-1_r4.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_rads/mp_harmonic_l1_m-1_r6.00.asc b/test/test_rads/mp_harmonic_l1_m-1_r6.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l1_m-1_r6.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_rads/mp_harmonic_l1_m-1_r8.00.asc b/test/test_rads/mp_harmonic_l1_m-1_r8.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l1_m-1_r8.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_rads/mp_harmonic_l1_m0_r4.00.asc b/test/test_rads/mp_harmonic_l1_m0_r4.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l1_m0_r4.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_rads/mp_harmonic_l1_m0_r6.00.asc b/test/test_rads/mp_harmonic_l1_m0_r6.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l1_m0_r6.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_rads/mp_harmonic_l1_m0_r8.00.asc b/test/test_rads/mp_harmonic_l1_m0_r8.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l1_m0_r8.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_rads/mp_harmonic_l1_m1_r4.00.asc b/test/test_rads/mp_harmonic_l1_m1_r4.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l1_m1_r4.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_rads/mp_harmonic_l1_m1_r6.00.asc b/test/test_rads/mp_harmonic_l1_m1_r6.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l1_m1_r6.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_rads/mp_harmonic_l1_m1_r8.00.asc b/test/test_rads/mp_harmonic_l1_m1_r8.00.asc
new file mode 100644
index 0000000..23eaf79
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l1_m1_r8.00.asc
@@ -0,0 +1 @@
+0.000000 0 0
diff --git a/test/test_rads/mp_harmonic_l2_m-1_r4.00.asc b/test/test_rads/mp_harmonic_l2_m-1_r4.00.asc
new file mode 100644
index 0000000..1788393
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m-1_r4.00.asc
@@ -0,0 +1 @@
+0.000000 0.01472621395168234246 1.564338938176084565e-09
diff --git a/test/test_rads/mp_harmonic_l2_m-1_r6.00.asc b/test/test_rads/mp_harmonic_l2_m-1_r6.00.asc
new file mode 100644
index 0000000..c21deda
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m-1_r6.00.asc
@@ -0,0 +1 @@
+0.000000 0.02208932395513868788 -5.683896356897662323e-10
diff --git a/test/test_rads/mp_harmonic_l2_m-1_r8.00.asc b/test/test_rads/mp_harmonic_l2_m-1_r8.00.asc
new file mode 100644
index 0000000..e414a3c
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m-1_r8.00.asc
@@ -0,0 +1 @@
+0.000000 0.02945243063564317518 4.961101318244261637e-10
diff --git a/test/test_rads/mp_harmonic_l2_m-2_r4.00.asc b/test/test_rads/mp_harmonic_l2_m-2_r4.00.asc
new file mode 100644
index 0000000..82b84c2
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m-2_r4.00.asc
@@ -0,0 +1 @@
+0.000000 0.006666694481119555572 8.704132452431115113e-11
diff --git a/test/test_rads/mp_harmonic_l2_m-2_r6.00.asc b/test/test_rads/mp_harmonic_l2_m-2_r6.00.asc
new file mode 100644
index 0000000..23d230f
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m-2_r6.00.asc
@@ -0,0 +1 @@
+0.000000 0.01000000117781982625 -8.931185644313283081e-12
diff --git a/test/test_rads/mp_harmonic_l2_m-2_r8.00.asc b/test/test_rads/mp_harmonic_l2_m-2_r8.00.asc
new file mode 100644
index 0000000..037dc22
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m-2_r8.00.asc
@@ -0,0 +1 @@
+0.000000 0.01333333435361747782 7.795539574208940158e-12
diff --git a/test/test_rads/mp_harmonic_l2_m0_r4.00.asc b/test/test_rads/mp_harmonic_l2_m0_r4.00.asc
new file mode 100644
index 0000000..39f9102
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m0_r4.00.asc
@@ -0,0 +1 @@
+0.000000 0.02449492188562732778 -2.064761642782179971e-06
diff --git a/test/test_rads/mp_harmonic_l2_m0_r6.00.asc b/test/test_rads/mp_harmonic_l2_m0_r6.00.asc
new file mode 100644
index 0000000..645f68a
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m0_r6.00.asc
@@ -0,0 +1 @@
+0.000000 0.03674238032471683485 -1.07237621656233756e-06
diff --git a/test/test_rads/mp_harmonic_l2_m0_r8.00.asc b/test/test_rads/mp_harmonic_l2_m0_r8.00.asc
new file mode 100644
index 0000000..3b487b4
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m0_r8.00.asc
@@ -0,0 +1 @@
+0.000000 0.04898979147934824641 5.938470720142443481e-07
diff --git a/test/test_rads/mp_harmonic_l2_m1_r4.00.asc b/test/test_rads/mp_harmonic_l2_m1_r4.00.asc
new file mode 100644
index 0000000..5a94e53
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m1_r4.00.asc
@@ -0,0 +1 @@
+0.000000 0.0343626003613806269 -1.430760424557310765e-06
diff --git a/test/test_rads/mp_harmonic_l2_m1_r6.00.asc b/test/test_rads/mp_harmonic_l2_m1_r6.00.asc
new file mode 100644
index 0000000..86cfa08
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m1_r6.00.asc
@@ -0,0 +1 @@
+0.000000 0.0515423300344349869 -5.755202806191689259e-07
diff --git a/test/test_rads/mp_harmonic_l2_m1_r8.00.asc b/test/test_rads/mp_harmonic_l2_m1_r8.00.asc
new file mode 100644
index 0000000..5867ad5
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m1_r8.00.asc
@@ -0,0 +1 @@
+0.000000 0.06872183696740095937 5.023342253613046556e-07
diff --git a/test/test_rads/mp_harmonic_l2_m2_r4.00.asc b/test/test_rads/mp_harmonic_l2_m2_r4.00.asc
new file mode 100644
index 0000000..28f8829
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m2_r4.00.asc
@@ -0,0 +1 @@
+0.000000 4.035524347289729086 -3.074478609784238764e-05
diff --git a/test/test_rads/mp_harmonic_l2_m2_r6.00.asc b/test/test_rads/mp_harmonic_l2_m2_r6.00.asc
new file mode 100644
index 0000000..0fe95f0
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m2_r6.00.asc
@@ -0,0 +1 @@
+0.000000 6.056567657173909325 -9.156676133234360109e-06
diff --git a/test/test_rads/mp_harmonic_l2_m2_r8.00.asc b/test/test_rads/mp_harmonic_l2_m2_r8.00.asc
new file mode 100644
index 0000000..54c1f0f
--- /dev/null
+++ b/test/test_rads/mp_harmonic_l2_m2_r8.00.asc
@@ -0,0 +1 @@
+0.000000 8.073152867502308538 7.992267128060954068e-06
diff --git a/test/test_rads/mp_harmonic_re_r4.00.ph.asc b/test/test_rads/mp_harmonic_re_r4.00.ph.asc
new file mode 100644
index 0000000..cb91193
--- /dev/null
+++ b/test/test_rads/mp_harmonic_re_r4.00.ph.asc
@@ -0,0 +1,104 @@
+"Time = 0
+0.000000 1.885620462359344884
+0.062832 1.870749141423040163
+0.125664 1.826374560850948958
+0.188496 1.753201105915120683
+0.251327 1.652381650637045007
+0.314159 1.525502080469085975
+0.376991 1.374561142997398511
+0.439823 1.201941618805164591
+0.502655 1.010371314242054508
+0.565487 0.8028641818766975025
+0.628319 0.5826890978521712672
+0.691150 0.3533343149107346925
+0.753982 0.1184000260980698632
+0.816814 -0.1184000260980697106
+0.879646 -0.3533343149107350256
+0.942478 -0.5826890978521728215
+1.005310 -0.8028641818766973914
+1.068142 -1.010371314242056284
+1.130973 -1.201941618805165035
+1.193805 -1.374561142997398511
+1.256637 -1.525502080469087307
+1.319469 -1.652381650637046562
+1.382301 -1.753201105915121349
+1.445133 -1.826374560850949402
+1.507964 -1.870749141423043271
+1.570796 -1.88562046235934444
+1.633628 -1.870749141423040607
+1.696460 -1.826374560850951401
+1.759292 -1.753201105915120683
+1.822124 -1.652381650637046118
+1.884956 -1.525502080469084198
+1.947787 -1.374561142997399621
+2.010619 -1.201941618805164147
+2.073451 -1.010371314242055396
+2.136283 -0.8028641818766962812
+2.199115 -0.5826890978521711562
+2.261947 -0.3533343149107335268
+2.324779 -0.1184000260980704877
+2.387610 0.1184000260980687114
+2.450442 0.3533343149107354697
+2.513274 0.5826890978521720443
+2.576106 0.8028641818766971694
+2.638938 1.010371314242054508
+2.701770 1.201941618805162149
+2.764602 1.374561142997397845
+2.827433 1.525502080469085309
+2.890265 1.652381650637046562
+2.953097 1.753201105915121127
+3.015929 1.826374560850950735
+3.078761 1.870749141423038386
+3.141593 1.88562046235934222
+3.204425 1.870749141423040607
+3.267256 1.826374560850949402
+3.330088 1.753201105915121571
+3.392920 1.65238165063704634
+3.455752 1.525502080469085087
+3.518584 1.374561142997398067
+3.581416 1.201941618805163481
+3.644247 1.01037131424205362
+3.707079 0.8028641818766973914
+3.769911 0.5826890978521690467
+3.832743 0.3533343149107325276
+3.895575 0.1184000260980697106
+3.958407 -0.1184000260980703906
+4.021239 -0.3533343149107334713
+4.084070 -0.5826890978521710451
+4.146902 -0.8028641818766979465
+4.209734 -1.010371314242054286
+4.272566 -1.201941618805163925
+4.335398 -1.374561142997398511
+4.398230 -1.525502080469087973
+4.461062 -1.652381650637048782
+4.523893 -1.753201105915123126
+4.586725 -1.826374560850950957
+4.649557 -1.870749141423039275
+4.712389 -1.885620462359343108
+4.775221 -1.870749141423039497
+4.838053 -1.826374560850950513
+4.900885 -1.753201105915120905
+4.963716 -1.652381650637045452
+5.026548 -1.525502080469088417
+5.089380 -1.374561142997396956
+5.152212 -1.201941618805162149
+5.215044 -1.010371314242052732
+5.277876 -0.8028641818766969473
+5.340708 -0.5826890978521664932
+5.403539 -0.3533343149107341374
+5.466371 -0.1184000260980691416
+5.529203 0.11840002609807114
+5.592035 0.3533343149107336933
+5.654867 0.5826890978521720443
+5.717699 0.8028641818766985017
+5.780530 1.010371314242055618
+5.843362 1.201941618805164813
+5.906194 1.374561142997399399
+5.969026 1.525502080469087307
+6.031858 1.652381650637045896
+6.094690 1.753201105915121794
+6.157522 1.826374560850950513
+6.220353 1.870749141423040607
+6.283185 1.885620462359341998
+
+
diff --git a/test/test_rads/mp_harmonic_re_r4.00.th.asc b/test/test_rads/mp_harmonic_re_r4.00.th.asc
new file mode 100644
index 0000000..1c772ba
--- /dev/null
+++ b/test/test_rads/mp_harmonic_re_r4.00.th.asc
@@ -0,0 +1,54 @@
+"Time = 0
+0.000000 5.589453454777939113e-15
+0.062832 2.766210742652663424
+0.125664 2.477552789726411131
+0.188496 2.478638771237914984
+0.251327 2.444486262107310903
+0.314159 2.401152665929550345
+0.376991 2.349059914969390039
+0.439823 2.288712171368338932
+0.502655 2.220688897786289306
+0.565487 2.145637535535099172
+0.628319 2.064264629325662348
+0.691150 1.977326668124248732
+0.753982 1.885620462359344884
+0.816814 1.789972562259663169
+0.879646 1.691228737165055396
+0.942478 1.590243318593252742
+1.005310 1.487868452698897892
+1.068142 1.384943554366310448
+1.130973 1.282285299243998145
+1.193805 1.180677601977725422
+1.256637 1.080862950184381832
+1.319469 0.9835337302469323584
+1.382301 0.8893250652390954913
+1.445133 0.7988080548373885925
+1.507964 0.7124845678073858934
+1.570796 0.6307831305050356496
+1.633628 0.554055765908431308
+1.696460 0.4825756971413410201
+1.759292 0.4165371763889615941
+1.822124 0.3560561880437642102
+1.884956 0.3011721218200319461
+1.947787 0.2518505709448911056
+2.010619 0.2079877214934599261
+2.073451 0.1694151916375707945
+2.136283 0.1359064401175641146
+2.199115 0.107183232570289777
+2.261947 0.08292353831428847522
+2.324779 0.06276949643620742536
+2.387610 0.0463360109799808087
+2.450442 0.0332196466294391729
+2.513274 0.02300753995123193396
+2.576106 0.01528628745114547251
+2.638938 0.009651014104101776972
+2.701770 0.005713608487684786012
+2.764602 0.003110625215124666887
+2.827433 0.001511039386052391489
+2.890265 0.0006225879825632539251
+2.953097 0.0001978629585884935326
+3.015929 3.916648887955317704e-05
+3.078761 2.425936286859116353e-06
+3.141593 4.564233006487166158e-20
+
+
diff --git a/test/test_rads/mp_harmonic_re_r6.00.ph.asc b/test/test_rads/mp_harmonic_re_r6.00.ph.asc
new file mode 100644
index 0000000..b9461c4
--- /dev/null
+++ b/test/test_rads/mp_harmonic_re_r6.00.ph.asc
@@ -0,0 +1,104 @@
+"Time = 0
+0.000000 2.828430783489839229
+0.062832 2.806126641164295776
+0.125664 2.739568731384816491
+0.188496 2.629807592015437123
+0.251327 2.478573341847115419
+0.314159 2.288249669502208317
+0.376991 2.06183782017832451
+0.439823 1.802910771645378718
+0.502655 1.515549169477164115
+0.565487 1.204287651286701122
+0.628319 0.8740331218936553404
+0.691150 0.5299945458910427787
+0.753982 0.1775983529518873338
+0.816814 -0.1775983529518899429
+0.879646 -0.5299945458910432228
+0.942478 -0.8740331218936545632
+1.005310 -1.204287651286699568
+1.068142 -1.515549169477143909
+1.130973 -1.802910771645380494
+1.193805 -2.061837820178326286
+1.256637 -2.288249669502206984
+1.319469 -2.478573341847111866
+1.382301 -2.629807592015436679
+1.445133 -2.739568731384819156
+1.507964 -2.806126641164297997
+1.570796 -2.828430783489844558
+1.633628 -2.806126641164301994
+1.696460 -2.739568731384820488
+1.759292 -2.629807592015433126
+1.822124 -2.478573341847113642
+1.884956 -2.288249669502206984
+1.947787 -2.061837820178323177
+2.010619 -1.80291077164538005
+2.073451 -1.515549169477145686
+2.136283 -1.204287651286697125
+2.199115 -0.8740331218936561175
+2.261947 -0.5299945458910438889
+2.324779 -0.1775983529518876947
+2.387610 0.1775983529518926074
+2.450442 0.5299945458910428897
+2.513274 0.8740331218936550073
+2.576106 1.204287651286701122
+2.638938 1.515549169477149682
+2.701770 1.802910771645380272
+2.764602 2.061837820178322289
+2.827433 2.288249669502207428
+2.890265 2.478573341847114975
+2.953097 2.62980759201543135
+3.015929 2.739568731384818712
+3.078761 2.806126641164299329
+3.141593 2.828430783489841005
+3.204425 2.806126641164295776
+3.267256 2.739568731384822708
+3.330088 2.62980759201543135
+3.392920 2.478573341847117639
+3.455752 2.288249669502201211
+3.518584 2.061837820178325842
+3.581416 1.802910771645380716
+3.644247 1.515549169477149238
+3.707079 1.204287651286701566
+3.769911 0.8740331218936546742
+3.832743 0.5299945458910432228
+3.895575 0.1775983529518881665
+3.958407 -0.1775983529518892212
+4.021239 -0.5299945458910433338
+4.084070 -0.8740331218936561175
+4.146902 -1.204287651286699345
+4.209734 -1.515549169477148128
+4.272566 -1.802910771645384491
+4.335398 -2.061837820178320957
+4.398230 -2.288249669502205208
+4.461062 -2.478573341847117639
+4.523893 -2.629807592015432682
+4.586725 -2.739568731384818712
+4.649557 -2.806126641164301549
+4.712389 -2.8284307834898339
+4.775221 -2.806126641164297553
+4.838053 -2.739568731384822708
+4.900885 -2.629807592015433571
+4.963716 -2.478573341847116307
+5.026548 -2.288249669502203876
+5.089380 -2.061837820178327174
+5.152212 -1.802910771645381827
+5.215044 -1.51554916947714613
+5.277876 -1.204287651286699123
+5.340708 -0.8740331218936538971
+5.403539 -0.5299945458910401141
+5.466371 -0.1775983529518918302
+5.529203 0.1775983529518887494
+5.592035 0.5299945458910434448
+5.654867 0.8740331218936602253
+5.717699 1.20428765128669979
+5.780530 1.515549169477146574
+5.843362 1.802910771645381383
+5.906194 2.061837820178325842
+5.969026 2.288249669502211425
+6.031858 2.478573341847118972
+6.094690 2.629807592015434015
+6.157522 2.739568731384820488
+6.220353 2.806126641164301105
+6.283185 2.828430783489841449
+
+
diff --git a/test/test_rads/mp_harmonic_re_r6.00.th.asc b/test/test_rads/mp_harmonic_re_r6.00.th.asc
new file mode 100644
index 0000000..960d44a
--- /dev/null
+++ b/test/test_rads/mp_harmonic_re_r6.00.th.asc
@@ -0,0 +1,54 @@
+"Time = 0
+0.000000 8.630578307977105267e-15
+0.062832 3.835910342632805392
+0.125664 3.75491412303918537
+0.188496 3.717958204858506033
+0.251327 3.666729398496961778
+0.314159 3.601728969865258811
+0.376991 3.523589845435281198
+0.439823 3.433068194208238388
+0.502655 3.331033343393394741
+0.565487 3.218456267823690986
+0.628319 3.09639689712424282
+0.691150 2.96599003877175571
+0.753982 2.828430783489839229
+0.816814 2.68495891615730331
+0.879646 2.53684320835061472
+0.942478 2.38536505012181399
+1.005310 2.231802749767326954
+1.068142 2.077415414779697933
+1.130973 1.923428001793415776
+1.193805 1.771016477654032695
+1.256637 1.62129445373965142
+1.319469 1.475300618257998186
+1.382301 1.333987498395791782
+1.445133 1.19821189831058672
+1.507964 1.068726741928595292
+1.570796 0.9461746957575534189
+1.633628 0.8310835390801751821
+1.696460 0.7238633617665178877
+1.759292 0.6248056651205932122
+1.822124 0.5340843049532448816
+1.884956 0.4517582111931253142
+1.947787 0.3777759311047821922
+2.010619 0.3119816351676065747
+2.073451 0.25412287068658862
+2.136283 0.2038597308953270515
+2.199115 0.1607749210873694878
+2.261947 0.1243854100744692659
+2.324779 0.09415431742212695965
+2.387610 0.06950410642079878076
+2.450442 0.04982950652954279414
+2.513274 0.03451126306259526266
+2.576106 0.02292939569776776543
+2.638938 0.01447651787011384854
+2.701770 0.008570349887252423984
+2.764602 0.00466591080389366572
+2.827433 0.00226653005001414224
+2.890265 0.0009338873098514001783
+2.953097 0.000296842439520530446
+3.015929 5.882034786612743327e-05
+3.078761 3.679970000328322189e-06
+3.141593 1.420239752188108614e-20
+
+
diff --git a/test/test_rads/mp_harmonic_re_r8.00.ph.asc b/test/test_rads/mp_harmonic_re_r8.00.ph.asc
new file mode 100644
index 0000000..a7b3444
--- /dev/null
+++ b/test/test_rads/mp_harmonic_re_r8.00.ph.asc
@@ -0,0 +1,104 @@
+"Time = 0
+0.000000 3.771241044431569822
+0.062832 3.741503090232144668
+0.125664 3.65275996216742449
+0.188496 3.506411335140886543
+0.251327 3.304763630936279117
+0.314159 3.050998145393059691
+0.376991 2.749116782820001159
+0.439823 2.403880014988298619
+0.502655 2.020732258839770257
+0.565487 1.605716954602965707
+0.628319 1.165377656963972441
+0.691150 0.7066607121772099376
+0.753982 0.236798535923516934
+0.816814 -0.2367985359235170451
+0.879646 -0.7066607121772101596
+0.942478 -1.165377656963972441
+1.005310 -1.605716954602966373
+1.068142 -2.020732258839770701
+1.130973 -2.403880014988300839
+1.193805 -2.749116782820001603
+1.256637 -3.050998145393060135
+1.319469 -3.304763630936283114
+1.382301 -3.506411335140886987
+1.445133 -3.652759962167423602
+1.507964 -3.741503090232144224
+1.570796 -3.771241044431569822
+1.633628 -3.741503090232146889
+1.696460 -3.652759962167422714
+1.759292 -3.506411335140887875
+1.822124 -3.304763630936283114
+1.884956 -3.05099814539305747
+1.947787 -2.749116782820002935
+2.010619 -2.403880014988297731
+2.073451 -2.020732258839767592
+2.136283 -1.605716954602964153
+2.199115 -1.165377656963972441
+2.261947 -0.7066607121772116029
+2.324779 -0.2367985359235149634
+2.387610 0.2367985359235169063
+2.450442 0.7066607121772154887
+2.513274 1.165377656963968445
+2.576106 1.605716954602966595
+2.638938 2.020732258839769369
+2.701770 2.403880014988298619
+2.764602 2.749116782820003824
+2.827433 3.050998145393061023
+2.890265 3.304763630936281338
+2.953097 3.506411335140884322
+3.015929 3.652759962167421826
+3.078761 3.741503090232148221
+3.141593 3.771241044431564937
+3.204425 3.741503090232146889
+3.267256 3.652759962167424046
+3.330088 3.506411335140888763
+3.392920 3.304763630936281338
+3.455752 3.050998145393055694
+3.518584 2.749116782820001603
+3.581416 2.403880014988297287
+3.644247 2.02073225883976848
+3.707079 1.605716954602963487
+3.769911 1.165377656963969333
+3.832743 0.7066607121772109368
+3.895575 0.2367985359235172116
+3.958407 -0.2367985359235150467
+4.021239 -0.7066607121772155997
+4.084070 -1.165377656963969777
+4.146902 -1.605716954602965707
+4.209734 -2.020732258839769813
+4.272566 -2.403880014988298619
+4.335398 -2.749116782820007376
+4.398230 -3.050998145393059247
+4.461062 -3.304763630936284891
+4.523893 -3.506411335140886099
+4.586725 -3.652759962167423158
+4.649557 -3.741503090232146
+4.712389 -3.771241044431562273
+4.775221 -3.741503090232146
+4.838053 -3.652759962167422714
+4.900885 -3.506411335140884766
+4.963716 -3.304763630936278673
+5.026548 -3.050998145393062355
+5.089380 -2.749116782819999383
+5.152212 -2.403880014988298619
+5.215044 -2.020732258839765372
+5.277876 -1.605716954602963042
+5.340708 -1.165377656963966224
+5.403539 -0.7066607121772103817
+5.466371 -0.2367985359235177667
+5.529203 0.2367985359235197373
+5.592035 0.7066607121772152666
+5.654867 1.16537765696397666
+5.717699 1.605716954602965929
+5.780530 2.020732258839772477
+5.843362 2.40388001498830306
+5.906194 2.74911678282000338
+5.969026 3.050998145393060579
+6.031858 3.30476363093628267
+6.094690 3.506411335140886099
+6.157522 3.65275996216742449
+6.220353 3.741503090232145112
+6.283185 3.771241044431567602
+
+
diff --git a/test/test_rads/mp_harmonic_re_r8.00.th.asc b/test/test_rads/mp_harmonic_re_r8.00.th.asc
new file mode 100644
index 0000000..f4bd4cf
--- /dev/null
+++ b/test/test_rads/mp_harmonic_re_r8.00.th.asc
@@ -0,0 +1,54 @@
+"Time = 0
+0.000000 1.167173116941090532e-14
+0.062832 4.985097729690861179
+0.125664 5.006552178572228229
+0.188496 4.957277617810031067
+0.251327 4.888972535616995074
+0.314159 4.802305291558570666
+0.376991 4.698119794614672884
+0.439823 4.577424253796267273
+0.502655 4.441377782345509395
+0.565487 4.291275036053113112
+0.628319 4.128529188273529549
+0.691150 3.954653385908060592
+0.753982 3.771241044431569822
+0.816814 3.579945246363574807
+0.879646 3.382457606387820004
+0.942478 3.180486757449974
+1.005310 2.975736998698742308
+1.068142 2.769887245295235445
+1.130973 2.564570667936973525
+1.193805 2.361355311176972283
+1.256637 2.161725950254292261
+1.319469 1.967067462642798414
+1.382301 1.778649970797490321
+1.445133 1.597615829945008503
+1.507964 1.424968995236297831
+1.570796 1.261566261010071521
+1.633628 1.108111391438389992
+1.696460 0.9651511145529183544
+1.759292 0.8330741930972249687
+1.822124 0.7121123782364628951
+1.884956 0.6023442935255907127
+1.947787 0.5037012491113050938
+2.010619 0.4159755124358952005
+2.073451 0.3388305198377580263
+2.136283 0.2718129735360658161
+2.199115 0.2143665854040466279
+2.261947 0.1658472086862932959
+2.324779 0.125539114716664707
+2.387610 0.09267214167284612558
+2.450442 0.066439342918431607
+2.513274 0.04601500952466608396
+2.576106 0.03057253988520207064
+2.638938 0.0193020149811337291
+2.701770 0.01142712803495176017
+2.764602 0.006221215106142334063
+2.827433 0.003022038471575898885
+2.890265 0.001245187367502235487
+2.953097 0.0003958012513789150666
+3.015929 7.844165047749454404e-05
+3.078761 4.904765484158289893e-06
+3.141593 6.131879098763379389e-21
+
+
diff --git a/test/test_rads/test_rads.par b/test/test_rads/test_rads.par
new file mode 100644
index 0000000..726e81b
--- /dev/null
+++ b/test/test_rads/test_rads.par
@@ -0,0 +1,92 @@
+
+ActiveThorns = "CoordBase SymBase Boundary CartGrid3d IOUtil Carpet CarpetLib  CarpetInterp AEILocalInterp InitBase Multipole"
+
+#############################################################
+# 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                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::poison_new_timelevels           = "yes"
+Carpet::check_for_poison                = "no"
+Carpet::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 = 3
+Multipole::radius[0] = 8.0
+Multipole::radius[1] = 6.0
+Multipole::radius[2] = 4.0
+Multipole::variables = "Multipole::harmonic_re{sw=-2 cmplx='Multipole::harmonic_im' name='harmonic'}"
+Multipole::l_max = 2
+
+Multipole::enable_test = "yes"
+Multipole::test_mode_proportional_to_r = "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