diff options
author | eschnett <eschnett@4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843> | 2012-06-16 15:39:07 +0000 |
---|---|---|
committer | eschnett <eschnett@4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843> | 2012-06-16 15:39:07 +0000 |
commit | 167b7353013f12c35d935a88a91c7c961ae5a2db (patch) | |
tree | a4ef14090fdbcdf49ee6a72c43a228eda8f3689c /test/test_driscollhealy/test_driscollhealy.par | |
parent | 74ad038e293cc3dfb5a819d7888b3473d6896885 (diff) |
Implement Driscoll&Healy integration
Implements a more accurate integration over the sphere using an
algorithm by Driscoll & Healy. This algorithm uses Gaussian
integration weights, leading (almost) to exponential convergence.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/Multipole/trunk@83 4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843
Diffstat (limited to 'test/test_driscollhealy/test_driscollhealy.par')
-rw-r--r-- | test/test_driscollhealy/test_driscollhealy.par | 90 |
1 files changed, 90 insertions, 0 deletions
diff --git a/test/test_driscollhealy/test_driscollhealy.par b/test/test_driscollhealy/test_driscollhealy.par new file mode 100644 index 0000000..4edac91 --- /dev/null +++ b/test/test_driscollhealy/test_driscollhealy.par @@ -0,0 +1,90 @@ + +ActiveThorns = "CoordBase SymBase Boundary CartGrid3d IOUtil Carpet CarpetLib CarpetInterp AEILocalInterp InitBase Multipole LoopControl" + +############################################################# +# 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::ntheta = 51 +Multipole::radius[0] = 8.0 +Multipole::variables = "Multipole::harmonic_re{sw=-2 cmplx='Multipole::harmonic_im' name='harmonic'}" +Multipole::integration_method = "DriscollHealy" + +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 |