diff options
97 files changed, 2003 insertions, 120 deletions
@@ -33,6 +33,11 @@ Copyright As noted above, this code has been hacked on by many different people, so its' copyright status is somewhat unclear. +Therefore here is a (incomplete) list of authors and their wishes: + +Jonathan Thornburg: +------------------- + However, the new code added in the June 2002 revisions is copyright (C) 2001-2002 by Jonathan Thornburg <jthorn@aei.mpg.de>. @@ -50,3 +55,19 @@ You should have received a copy of the GNU General Public License along with this thorn (see the file COPYING in this directory); if not, write to the Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA + +Frank Loeffler +-------------- + +Code added by myself (which is mainly Minkowski_conf_wave and the +code using the conformal factor added in January 2003) is free +software, you can redistribute it and/or modify it under the terms +of the GNU General Public License as published bythe Free Software +Foundation; either version 2 of the License, or(at your option) +any later version. + +This new code is distributed in the hope that it will be useful, but +WITHOUT ANY WARRANTY; without even the implied warranty of +MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +GNU General Public License for more details. + diff --git a/doc/how_to_add_a_new_model b/doc/how_to_add_a_new_model index 86e8eed..2a49c63 100644 --- a/doc/how_to_add_a_new_model +++ b/doc/how_to_add_a_new_model @@ -1,7 +1,7 @@ ************************************************ ***** How to Add a New Model to this Thorn ***** ************************************************ -$Header: /usr/local/svn/cvs-repositories/numrelcvs/AEIThorns/Exact/doc/how_to_add_a_new_model,v 1.2 2002-11-17 15:44:32 jthorn Exp $ +$Header: /usr/local/svn/cvs-repositories/numrelcvs/AEIThorns/Exact/doc/how_to_add_a_new_model,v 1.3 2003-01-28 14:49:40 knarf Exp $ A general note: Please keep the sections for various models in the same order in all the files. This makes it _much_ easier for people @@ -51,10 +51,7 @@ Then you need to modify... src/metrics/make.code.defn: Add the new file name to the list of files to be compiled. -If your model needs to set up a stress-energy tensor (n.b. as described -in doc/documentation.tex, this also includes the case of a vacuum solution -with a cosmological constant), then you need to modify some additional -files: +You need to modify some additional files: interface.ccl: Add copies of all your model's CCTK_INT and/or CCTK_REAL parameters to the Exact_pars_int and/or Exact_pars_real variable groups as diff --git a/interface.ccl b/interface.ccl index 7ab7310..d7db872 100644 --- a/interface.ccl +++ b/interface.ccl @@ -2,7 +2,7 @@ # $Header$ implements: exact -inherits: admbase grid coordgauge +inherits: admbase grid coordgauge staticconformal friend: ADMCoupling # evolution thorns need to include these files @@ -39,6 +39,10 @@ EXTENDS KEYWORD evolution_method "slice" :: "Evolve arbitrary slice and extract Cauchy data" } +shares:StaticConformal + +USES KEYWORD conformal_storage + ################################################################################ ##### parameters for the blended boundary conditions ########################### ################################################################################ @@ -162,6 +166,7 @@ KEYWORD exact_model "The exact solution/coordinates used in thorn exact" "Minkowski/shift" :: "Minkowski spacetime with time-dependent shift vector" "Minkowski/funny" :: "Minkowski spacetime in non-trivial spatial coordinates" "Minkowski/gauge wave" :: "Minkowski spacetime in gauge-wave coordinates" +"Minkowski/conf wave" :: "Minkowski spacetime with 'waves' in conformal factor" # # black hole spacetimes # @@ -285,6 +290,30 @@ REAL Minkowski_gauge_wave__phase \ } 0.0 ################################################################################ + +# +# parameters for Minkowski spacetime in gauge-wave coordinates +# + +REAL Minkowski_conf_wave__amplitude \ + "Minkowski/conf wave: amplitude of the variation of the conformal factor" +{ + 0:* :: "any positive real number" +} 0.5 + +REAL Minkowski_conf_wave__wavelength \ + "Minkowski/conf wave: wave length in cactus units" +{ + 0:* :: "any positive real number" +} 1.0 + +INT Minkowski_conf_wave__direction \ + "Minkowski/conf wave: direction of 'wave' 0,1,2 : x,y,z" +{ + 0:2 :: "0, 1 or 2 for x, y or z" +} 0 + +################################################################################ ##### black hole spacetimes #################################################### ################################################################################ diff --git a/src/Bona_Masso_data.F77 b/src/Bona_Masso_data.F77 index 9805e0a..8318c01 100644 --- a/src/Bona_Masso_data.F77 +++ b/src/Bona_Masso_data.F77 @@ -13,6 +13,8 @@ c $Header$ $ x, y, z, t, $ gxx, gyy, gzz, gxy, gyz, gxz, $ hxx, hyy, hzz, hxy, hyz, hxz, + $ psi, psix, psiy, psiz, + $ psixx, psiyy, psizz, psixy, psiyz, psixz, $ dxgxx, dxgyy, dxgzz, dxgxy, dxgyz, dxgxz, $ dygxx, dygyy, dygzz, dygxy, dygyz, dygxz, $ dzgxx, dzgyy, dzgzz, dzgxy, dzgyz, dzgxz, @@ -21,9 +23,12 @@ c $Header$ implicit none CCTK_INT decoded_exact_model + CCTK_INT psi_on CCTK_REAL x, y, z, t, $ gxx, gyy, gzz, gxy, gyz, gxz, $ hxx, hyy, hzz, hxy, hyz, hxz, + $ psi, psix, psiy, psiz, + $ psixx, psiyy, psizz, psixy, psiyz, psixz, $ dxgxx, dxgyy, dxgzz, dxgxy, dxgyz, dxgxz, $ dygxx, dygyy, dygzz, dygxy, dygyz, dygxz, $ dzgxx, dzgyy, dzgzz, dzgxy, dzgyz, dzgxz, @@ -37,7 +42,7 @@ C alp is N, betax is N^x etc. C bxy is (/2) dN^y / dx (sic and sic!) C ax is dN / dx / N (sic!) - CCTK_REAL eps, + CCTK_REAL eps, four_eps2, $ gdtt, gdtx, gdty, gdtz, $ gutt, gutx, guty, gutz, $ guxx, guyy, guzz, guxy, guyz, guxz, @@ -45,13 +50,36 @@ C ax is dN / dx / N (sic!) $ gdxx_p, gdyy_p, gdzz_p, gdxy_p, gdyz_p, gdxz_p, $ gutt_p, gutx_p, guty_p, gutz_p, $ guxx_p, guyy_p, guzz_p, guxy_p, guyz_p, guxz_p, + $ psix_p, psiy_p, psiz_p, $ gdtt_m, gdtx_m, gdty_m, gdtz_m, $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, $ gutt_m, gutx_m, guty_m, gutz_m, - $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, + $ psix_m, psiy_m, psiz_m, + $ psix_px_p, psix_py_p, psix_pz_p, + $ psiy_py_p, psiy_pz_p, + $ psiz_pz_p, + $ + $ psix_mx_m, psix_my_m, psix_mz_m, + $ psiy_my_m, psiy_mz_m, + $ psiz_mz_m, + $ + $ psix_py_m, psix_pz_m, + $ psiy_pz_m, + $ + $ psix_my_p, psix_mz_p, + $ psiy_mz_p parameter (eps=1.d-6) +C Save, if we have to provide the conformal factor + + if (psi .gt. 0.0d0) then + psi_on=1 + else + psi_on=0 + end if + C Get the spacetime metric and its inverse at the base point. call Exact__metric( @@ -60,7 +88,7 @@ C Get the spacetime metric and its inverse at the base point. $ gdtt, gdtx, gdty, gdtz, $ gxx, gyy, gzz, gxy, gyz, gxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, psi) C Calculate lapse and shift from the upper metric. @@ -82,14 +110,14 @@ C Calculate x-derivatives. $ gdtt_p, gdtx_p, gdty_p, gdtz_p, $ gdxx_p, gdyy_p, gdzz_p, gdxy_p, gdyz_p, gdxz_p, $ gutt_p, gutx_p, guty_p, gutz_p, - $ guxx_p, guyy_p, guzz_p, guxy_p, guyz_p, guxz_p) + $ guxx_p, guyy_p, guzz_p, guxy_p, guyz_p, guxz_p, psix_p) call Exact__metric( $ decoded_exact_model, $ x-eps, y, z, t, $ gdtt_m, gdtx_m, gdty_m, gdtz_m, $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, $ gutt_m, gutx_m, guty_m, gutz_m, - $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m) + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psix_m) dxgxx = (gdxx_p - gdxx_m) / 4.d0 / eps dxgyy = (gdyy_p - gdyy_m) / 4.d0 / eps @@ -104,6 +132,7 @@ C Calculate x-derivatives. bxy = ((- guty_p / gutt_p) - (- guty_m / gutt_m)) / 4.d0 / eps bxz = ((- gutz_p / gutt_p) - (- gutz_m / gutt_m)) / 4.d0 / eps + psix = (psix_p - psix_m) / 2.0d0 / eps C Calculate y-derivatives. call Exact__metric( @@ -112,14 +141,14 @@ C Calculate y-derivatives. $ gdtt_p, gdtx_p, gdty_p, gdtz_p, $ gdxx_p, gdyy_p, gdzz_p, gdxy_p, gdyz_p, gdxz_p, $ gutt_p, gutx_p, guty_p, gutz_p, - $ guxx_p, guyy_p, guzz_p, guxy_p, guyz_p, guxz_p) + $ guxx_p, guyy_p, guzz_p, guxy_p, guyz_p, guxz_p, psiy_p) call Exact__metric( $ decoded_exact_model, $ x, y-eps, z, t, $ gdtt_m, gdtx_m, gdty_m, gdtz_m, $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, $ gutt_m, gutx_m, guty_m, gutz_m, - $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m) + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psiy_m) dygxx = (gdxx_p - gdxx_m) / 4.d0 / eps dygyy = (gdyy_p - gdyy_m) / 4.d0 / eps @@ -134,6 +163,7 @@ C Calculate y-derivatives. byy = ((- guty_p / gutt_p) - (- guty_m / gutt_m)) / 4.d0 / eps byz = ((- gutz_p / gutt_p) - (- gutz_m / gutt_m)) / 4.d0 / eps + psiy = (psiy_p - psiy_m) / 2.0d0 / eps C Calculate z-derivatives. call Exact__metric( @@ -142,14 +172,14 @@ C Calculate z-derivatives. $ gdtt_p, gdtx_p, gdty_p, gdtz_p, $ gdxx_p, gdyy_p, gdzz_p, gdxy_p, gdyz_p, gdxz_p, $ gutt_p, gutx_p, guty_p, gutz_p, - $ guxx_p, guyy_p, guzz_p, guxy_p, guyz_p, guxz_p) + $ guxx_p, guyy_p, guzz_p, guxy_p, guyz_p, guxz_p, psiz_p) call Exact__metric( $ decoded_exact_model, $ x, y, z-eps, t, $ gdtt_m, gdtx_m, gdty_m, gdtz_m, $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, $ gutt_m, gutx_m, guty_m, gutz_m, - $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m) + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psiz_m) dzgxx = (gdxx_p - gdxx_m) / 4.d0 / eps dzgyy = (gdyy_p - gdyy_m) / 4.d0 / eps @@ -164,6 +194,7 @@ C Calculate z-derivatives. bzy = ((- guty_p / gutt_p) - (- guty_m / gutt_m)) / 4.d0 / eps bzz = ((- gutz_p / gutt_p) - (- gutz_m / gutt_m)) / 4.d0 / eps + psiz = (psiz_p - psiz_m) / 2.0d0 / eps C Calculate t-derivatives, and extrinsic curvature. call Exact__metric( @@ -172,14 +203,14 @@ C Calculate t-derivatives, and extrinsic curvature. $ gdtt_p, gdtx_p, gdty_p, gdtz_p, $ gdxx_p, gdyy_p, gdzz_p, gdxy_p, gdyz_p, gdxz_p, $ gutt_p, gutx_p, guty_p, gutz_p, - $ guxx_p, guyy_p, guzz_p, guxy_p, guyz_p, guxz_p) + $ guxx_p, guyy_p, guzz_p, guxy_p, guyz_p, guxz_p, psi) call Exact__metric( $ decoded_exact_model, $ x, y, z, t-eps, $ gdtt_m, gdtx_m, gdty_m, gdtz_m, $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, $ gutt_m, gutx_m, guty_m, gutz_m, - $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m) + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psi) hxx = - (gdxx_p - gdxx_m) / 4.d0 / eps / alp $ + (dxgxx * betax + dygxx * betay + dzgxx * betaz @@ -208,6 +239,145 @@ C Calculate t-derivatives, and extrinsic curvature. $ + bxx * gxz + bxy * gyz + bxz * gzz $ + bzx * gxx + bzy * gxy + bzz * gxz) / alp +C Calculate second derivatives of the conformal factor + if (psi_on .eq. 1) then + call Exact__metric( + $ decoded_exact_model, + $ x+eps+eps, y, z, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psix_px_p) + call Exact__metric( + $ decoded_exact_model, + $ x, y+eps+eps, z, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psiy_py_p) + call Exact__metric( + $ decoded_exact_model, + $ x, y, z+eps+eps, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psiz_pz_p) + call Exact__metric( + $ decoded_exact_model, + $ x-eps-eps, y, z, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psix_mx_m) + call Exact__metric( + $ decoded_exact_model, + $ x, y-eps-eps, z, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psiy_my_m) + call Exact__metric( + $ decoded_exact_model, + $ x, y, z-eps-eps, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psiz_mz_m) + call Exact__metric( + $ decoded_exact_model, + $ x+eps, y+eps, z, t, + $ gdtt_p, gdtx_p, gdty_p, gdtz_p, + $ gdxx_p, gdyy_p, gdzz_p, gdxy_p, gdyz_p, gdxz_p, + $ gutt_p, gutx_p, guty_p, gutz_p, + $ guxx_p, guyy_p, guzz_p, guxy_p, guyz_p, guxz_p, psix_py_p) + call Exact__metric( + $ decoded_exact_model, + $ x+eps, y, z+eps, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psix_pz_p) + call Exact__metric( + $ decoded_exact_model, + $ x, y+eps, z+eps, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psiy_pz_p) + call Exact__metric( + $ decoded_exact_model, + $ x-eps, y-eps, z, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psix_my_m) + call Exact__metric( + $ decoded_exact_model, + $ x-eps, y, z-eps, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psix_mz_m) + call Exact__metric( + $ decoded_exact_model, + $ x, y-eps, z-eps, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psiy_mz_m) + call Exact__metric( + $ decoded_exact_model, + $ x+eps, y-eps, z, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psix_py_m) + call Exact__metric( + $ decoded_exact_model, + $ x+eps, y, z-eps, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psix_pz_m) + call Exact__metric( + $ decoded_exact_model, + $ x-eps, y+eps, z, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psix_my_p) + call Exact__metric( + $ decoded_exact_model, + $ x-eps, y, z+eps, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psix_mz_p) + call Exact__metric( + $ decoded_exact_model, + $ x, y+eps, z-eps, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psiy_pz_m) + call Exact__metric( + $ decoded_exact_model, + $ x, y-eps, z+eps, t, + $ gdtt_m, gdtx_m, gdty_m, gdtz_m, + $ gdxx_m, gdyy_m, gdzz_m, gdxy_m, gdyz_m, gdxz_m, + $ gutt_m, gutx_m, guty_m, gutz_m, + $ guxx_m, guyy_m, guzz_m, guxy_m, guyz_m, guxz_m, psiy_mz_p) + + four_eps2 = 4.0d0 * eps**2 + + psixx = (psix_mx_m- 2*psi+ psix_px_p) / four_eps2 + psiyy = (psiy_my_m- 2*psi+ psiy_py_p) / four_eps2 + psizz = (psiz_mz_m- 2*psi+ psiz_pz_p) / four_eps2 + psixy = (psix_my_m- psix_my_p- psix_py_m+ psix_py_p) / four_eps2 + psixz = (psix_mz_m- psix_mz_p- psix_pz_m+ psix_pz_p) / four_eps2 + psiyz = (psiy_mz_m- psiy_mz_p- psiy_pz_m+ psiy_pz_p) / four_eps2 + + end if return end diff --git a/src/ParamCheck.c b/src/ParamCheck.c index ed19b5b..c66125e 100644 --- a/src/ParamCheck.c +++ b/src/ParamCheck.c @@ -64,12 +64,12 @@ void Exact_ParamCheck(CCTK_ARGUMENTS); void Exact_ParamCheck(CCTK_ARGUMENTS) { DECLARE_CCTK_PARAMETERS; - +/* if(! CCTK_EQUALS(metric_type, "physical")) { CCTK_PARAMWARN("Unknown ADMBase::metric_type - known types are \"physical\""); } - +*/ if ( CCTK_EQUALS(shift_evolution_method, "exact") && CCTK_EQUALS(initial_shift, "none") ) { diff --git a/src/blended_boundary.F77 b/src/blended_boundary.F77 index 804f679..c51ac61 100644 --- a/src/blended_boundary.F77 +++ b/src/blended_boundary.F77 @@ -24,6 +24,12 @@ C $Header$ CCTK_REAL dygxxe, dygyye, dygzze, dygxye, dygyze, dygxze CCTK_REAL dzgxxe, dzgyye, dzgzze, dzgxye, dzgyze, dzgxze + CCTK_REAL + $ exact_psi, + $ exact_psix, exact_psiy, exact_psiz, + $ exact_psixx, exact_psiyy, exact_psizz, + $ exact_psixy, exact_psiyz, exact_psixz + CCTK_REAL alpe, axe, aye, aze CCTK_REAL betaxe,betaye,betaze CCTK_REAL bxxe,bxye,bxze,byxe,byye,byze,bzxe,bzye,bzze @@ -84,11 +90,32 @@ c there. if (r(i,j,k) .ge. rinner) then +C Initialize the psi of exact +C (also to tell the models about the conformal_state) + if (conformal_state .ne. 0) then + exact_psi = 1.0D0 + else + exact_psi = 0.0D0 + end if + exact_psix = 0.0D0 + exact_psiy = 0.0D0 + exact_psiz = 0.0D0 + exact_psixx = 0.0D0 + exact_psixx = 0.0D0 + exact_psizz = 0.0D0 + exact_psixy = 0.0D0 + exact_psiyz = 0.0D0 + exact_psixz = 0.0D0 + call Exact__Bona_Masso_data( $ decoded_exact_model, $ x(i,j,k), y(i,j,k), z(i,j,k), time, $ gxxe, gyye, gzze, gxye, gyze, gxze, $ kxxe, kyye, kzze, kxye, kyze, kxze, + $ exact_psi, + $ exact_psix, exact_psiy, exact_psiz, + $ exact_psixx, exact_psiyy, exact_psizz, + $ exact_psixy, exact_psiyz, exact_psixz, $ dxgxxe, dxgyye, dxgzze, dxgxye, dxgyze, dxgxze, $ dygxxe, dygyye, dygzze, dygxye, dygyze, dygxze, $ dzgxxe, dzgyye, dzgzze, dzgxye, dzgyze, dzgxze, diff --git a/src/boundary.F77 b/src/boundary.F77 index 1862e43..3e65f53 100644 --- a/src/boundary.F77 +++ b/src/boundary.F77 @@ -25,6 +25,11 @@ C $Header$ $ bxxjunk, bxyjunk, bxzjunk, $ byxjunk, byyjunk, byzjunk, $ bzxjunk, bzyjunk, bzzjunk + CCTK_REAL + $ exact_psi, + $ exact_psix, exact_psiy, exact_psiz, + $ exact_psixx, exact_psiyy, exact_psizz, + $ exact_psixy, exact_psiyz, exact_psixz C Grid parameters. @@ -32,6 +37,23 @@ C Grid parameters. ny = cctk_lsh(2) nz = cctk_lsh(3) +C Initialize the psi of exact +C (also to tell the models about the conformal_state) + if (conformal_state .ne. 0) then + exact_psi = 1.0D0 + else + exact_psi = 0.0D0 + end if + exact_psix = 0.0D0 + exact_psiy = 0.0D0 + exact_psiz = 0.0D0 + exact_psixx = 0.0D0 + exact_psixx = 0.0D0 + exact_psizz = 0.0D0 + exact_psixy = 0.0D0 + exact_psiyz = 0.0D0 + exact_psixz = 0.0D0 + C Set all initial data including dijk and vi on all points which C are on the boundary of the domain if it really is the boundary C of the complete grid. Treat all six sides of the grid cube this way. @@ -55,6 +77,10 @@ C overwritten by Exact__gauge. gxy(i,j,k), gyz(i,j,k), gxz(i,j,k), \ kxx(i,j,k), kyy(i,j,k), kzz(i,j,k), \ kxy(i,j,k), kyz(i,j,k), kxz(i,j,k), \ + exact_psi, \ + exact_psix, exact_psiy, exact_psiz, \ + exact_psixx, exact_psiyy, exact_psizz, \ + exact_psixy, exact_psiyz, exact_psixz, \ dxgxxjunk, dxgyyjunk, dxgzzjunk, \ dxgxyjunk, dxgyzjunk, dxgxzjunk, \ dygxxjunk, dygyyjunk, dygzzjunk, \ diff --git a/src/decode_pars.F77 b/src/decode_pars.F77 index 1368f11..7790171 100644 --- a/src/decode_pars.F77 +++ b/src/decode_pars.F77 @@ -57,6 +57,8 @@ c Minkowski spacetime decoded_exact_model = EXACT__Minkowski_funny elseif (CCTK_Equals(exact_model, "Minkowski/gauge wave") .ne. 0) then decoded_exact_model = EXACT__Minkowski_gauge_wave + elseif (CCTK_Equals(exact_model, "Minkowski/conf wave") .ne. 0) then + decoded_exact_model = EXACT__Minkowski_conf_wave c black hole spacetimes elseif (CCTK_Equals(exact_model, "Schwarzschild/EF") .ne. 0) then diff --git a/src/gauge.F77 b/src/gauge.F77 index f5c2e88..4df2cb7 100644 --- a/src/gauge.F77 +++ b/src/gauge.F77 @@ -34,6 +34,11 @@ C $Header$ $ bxxjunk, bxyjunk, bxzjunk, $ byxjunk, byyjunk, byzjunk, $ bzxjunk, bzyjunk, bzzjunk + CCTK_REAL + $ exact_psi, + $ exact_psix, exact_psiy, exact_psiz, + $ exact_psixx, exact_psiyy, exact_psizz, + $ exact_psixy, exact_psiyz, exact_psixz C Grid parameters. @@ -58,6 +63,23 @@ C Set both lapse and shift. do j=1,ny do i=1,nx +C Initialize the psi of exact +C (also to tell the models about the conformal_state) + if (conformal_state .ne. 0) then + exact_psi = 1.0D0 + else + exact_psi = 0.0D0 + end if + exact_psix = 0.0D0 + exact_psiy = 0.0D0 + exact_psiz = 0.0D0 + exact_psixx = 0.0D0 + exact_psixx = 0.0D0 + exact_psizz = 0.0D0 + exact_psixy = 0.0D0 + exact_psiyz = 0.0D0 + exact_psixz = 0.0D0 + call Exact__Bona_Masso_data( $ decoded_exact_model, $ x(i,j,k), y(i,j,k), z(i,j,k), tplushalf, @@ -65,6 +87,10 @@ C Set both lapse and shift. $ gxyjunk, gyzjunk, gxzjunk, $ hxxjunk, hyyjunk, hzzjunk, $ hxyjunk, hyzjunk, hxzjunk, + $ exact_psi, + $ exact_psix, exact_psiy, exact_psiz, + $ exact_psixx, exact_psiyy, exact_psizz, + $ exact_psixy, exact_psiyz, exact_psixz, $ dxgxxjunk, dxgyyjunk, dxgzzjunk, $ dxgxyjunk, dxgyzjunk, dxgxzjunk, $ dygxxjunk, dygyyjunk, dygzzjunk, @@ -89,6 +115,23 @@ C Set lapse only. do j=1,ny do i=1,nx +C Initialize the psi of exact +C (also to tell the models about the conformal_state) + if (conformal_state .ne. 0) then + exact_psi = 1.0D0 + else + exact_psi = 0.0D0 + end if + exact_psix = 0.0D0 + exact_psiy = 0.0D0 + exact_psiz = 0.0D0 + exact_psixx = 0.0D0 + exact_psixx = 0.0D0 + exact_psizz = 0.0D0 + exact_psixy = 0.0D0 + exact_psiyz = 0.0D0 + exact_psixz = 0.0D0 + call Exact__Bona_Masso_data( $ decoded_exact_model, $ x(i,j,k), y(i,j,k), z(i,j,k), tplushalf, @@ -96,6 +139,10 @@ C Set lapse only. $ gxyjunk, gyzjunk, gxzjunk, $ hxxjunk, hyyjunk, hzzjunk, $ hxyjunk, hyzjunk, hxzjunk, + $ exact_psi, + $ exact_psix, exact_psiy, exact_psiz, + $ exact_psixx, exact_psiyy, exact_psizz, + $ exact_psixy, exact_psiyz, exact_psixz, $ dxgxxjunk, dxgyyjunk, dxgzzjunk, $ dxgxyjunk, dxgyzjunk, dxgxzjunk, $ dygxxjunk, dygyyjunk, dygzzjunk, @@ -121,6 +168,23 @@ C Set shift only. do j=1,ny do i=1,nx +C Initialize the psi of exact +C (also to tell the models about the conformal_state) + if (conformal_state .ne. 0) then + exact_psi = 1.0D0 + else + exact_psi = 0.0D0 + end if + exact_psix = 0.0D0 + exact_psiy = 0.0D0 + exact_psiz = 0.0D0 + exact_psixx = 0.0D0 + exact_psixx = 0.0D0 + exact_psizz = 0.0D0 + exact_psixy = 0.0D0 + exact_psiyz = 0.0D0 + exact_psixz = 0.0D0 + call Exact__Bona_Masso_data( $ decoded_exact_model, $ x(i,j,k), y(i,j,k), z(i,j,k), tplushalf, @@ -128,6 +192,10 @@ C Set shift only. $ gxyjunk, gyzjunk, gxzjunk, $ hxxjunk, hyyjunk, hzzjunk, $ hxyjunk, hyzjunk, hxzjunk, + $ exact_psi, + $ exact_psix, exact_psiy, exact_psiz, + $ exact_psixx, exact_psiyy, exact_psizz, + $ exact_psixy, exact_psiyz, exact_psixz, $ dxgxxjunk, dxgyyjunk, dxgzzjunk, $ dxgxyjunk, dxgyzjunk, dxgxzjunk, $ dygxxjunk, dygyyjunk, dygzzjunk, diff --git a/src/include/Scalar_CalcTmunu.inc b/src/include/Scalar_CalcTmunu.inc index 4f5f240..50020cd 100644 --- a/src/include/Scalar_CalcTmunu.inc +++ b/src/include/Scalar_CalcTmunu.inc @@ -40,6 +40,8 @@ c no stress-energy tensor in this model c no stress-energy tensor in this model elseif (decoded_exact_model .eq. EXACT__Minkowski_gauge_wave) then c no stress-energy tensor in this model + elseif (decoded_exact_model .eq. EXACT__Minkowski_conf_wave) then +c no stress-energy tensor in this model cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc ccccc black hole spacetimes cccccccccccccccccccccccccccccccccccccccccccccccccccc diff --git a/src/include/param_defs.inc b/src/include/param_defs.inc index df25a06..45642ba 100644 --- a/src/include/param_defs.inc +++ b/src/include/param_defs.inc @@ -34,6 +34,7 @@ c Minkowski spacetime #define EXACT__Minkowski_shift 2 #define EXACT__Minkowski_funny 3 #define EXACT__Minkowski_gauge_wave 4 +#define EXACT__Minkowski_conf_wave 5 c black hole spacetimes #define EXACT__Schwarzschild_EF 10 diff --git a/src/initialize.F77 b/src/initialize.F77 index 00fca60..74d47f6 100644 --- a/src/initialize.F77 +++ b/src/initialize.F77 @@ -30,9 +30,25 @@ C $Header$ $ dygxyjunk, dygyzjunk, dygxzjunk, $ dzgxxjunk, dzgyyjunk, dzgzzjunk, $ dzgxyjunk, dzgyzjunk, dzgxzjunk + CCTK_REAL + $ exact_psi, + $ exact_psix, exact_psiy, exact_psiz, + $ exact_psixx, exact_psiyy, exact_psizz, + $ exact_psixy, exact_psiyz, exact_psixz call CCTK_INFO('setting exact data on slice') +C Set conformal state + if (CCTK_EQUALS(metric_type, "static conformal")) then + conformal_state=1 + if (CCTK_EQUALS(conformal_storage,"factor+derivs")) then + conformal_state = 2 + else if (CCTK_EQUALS(conformal_storage, + $ "factor+derivs+2nd derivs")) then + conformal_state = 3 + end if + end if + C Note I assume time has been initialized to physical time. C Set data pointwise. @@ -43,7 +59,24 @@ C Set data pointwise. do k=1,nz do j=1,ny do i=1,nx - + +C Initialize the psi of exact +C (also to tell the models about the conformal_state) + if (conformal_state .ne. 0) then + exact_psi = 1.0D0 + else + exact_psi = 0.0D0 + end if + exact_psix = 0.0D0 + exact_psiy = 0.0D0 + exact_psiz = 0.0D0 + exact_psixx = 0.0D0 + exact_psixx = 0.0D0 + exact_psizz = 0.0D0 + exact_psixy = 0.0D0 + exact_psiyz = 0.0D0 + exact_psixz = 0.0D0 + call Exact__Bona_Masso_data( $ decoded_exact_model, $ x(i,j,k), y(i,j,k), z(i,j,k), cctk_time, @@ -51,6 +84,10 @@ C Set data pointwise. $ gxy(i,j,k), gyz(i,j,k), gxz(i,j,k), $ kxx(i,j,k), kyy(i,j,k), kzz(i,j,k), $ kxy(i,j,k), kyz(i,j,k), kxz(i,j,k), + $ exact_psi, + $ exact_psix, exact_psiy, exact_psiz, + $ exact_psixx, exact_psiyy, exact_psizz, + $ exact_psixy, exact_psiyz, exact_psixz, $ dxgxxjunk, dxgyyjunk, dxgzzjunk, $ dxgxyjunk, dxgyzjunk, dxgxzjunk, $ dygxxjunk, dygyyjunk, dygzzjunk, @@ -63,6 +100,24 @@ C Set data pointwise. $ byxjunk, byyjunk, byzjunk, $ bzxjunk, bzyjunk, bzzjunk) +C Save the conformal factor if wanted + if (conformal_state .ne. 0) then + psi(i,j,k) = exact_psi + if (conformal_state .gt. 1) then + psix(i,j,k) = exact_psix + psiy(i,j,k) = exact_psiy + psiz(i,j,k) = exact_psiz + if (conformal_state .gt. 2) then + psixx(i,j,k) = exact_psixx + psiyy(i,j,k) = exact_psiyy + psizz(i,j,k) = exact_psizz + psixy(i,j,k) = exact_psixy + psiyz(i,j,k) = exact_psiyz + psixz(i,j,k) = exact_psixz + end if + end if + end if + end do end do end do @@ -71,8 +126,7 @@ C Tell the code there is no need to treat the conformal factor C as a separate field. That is, we have set the physical metric here. c Commented out in einstein revamp, now Exact doesnot inherit anything c about the conformal factor -c conformal_state = 0 -c psi = 1.0D0 +c Now ot does again (see above, knarf) return end diff --git a/src/metric.F77 b/src/metric.F77 index c547b8d..32cbf81 100644 --- a/src/metric.F77 +++ b/src/metric.F77 @@ -14,7 +14,7 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz,rama) + $ guxx, guyy, guzz, guxy, guyz, guxz, psi, rama) implicit none DECLARE_CCTK_FUNCTIONS @@ -25,7 +25,7 @@ c arguments $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz, rama + $ guxx, guyy, guzz, guxy, guyz, guxz, psi, rama c local variables character*100 warn_buffer @@ -42,7 +42,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Minkowski_shift) then call Exact__Minkowski_shift( @@ -50,7 +51,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Minkowski_funny) then call Exact__Minkowski_funny( @@ -58,7 +60,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Minkowski_gauge_wave) then call Exact__Minkowski_gauge_wave( @@ -66,7 +69,17 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) + + elseif (decoded_exact_model .eq. EXACT__Minkowski_conf_wave) then + call Exact__Minkowski_conf_wave( + $ x, y, z, t, + $ gdtt, gdtx, gdty, gdtz, + $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, + $ gutt, gutx, guty, gutz, + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc @@ -80,7 +93,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Schwarzschild_PG) then call Exact__Schwarzschild_PG( @@ -88,14 +102,16 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Schwarzschild_Novikov) then call Exact__Schwarzschild_Novikov(x,y,z,t, $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Kerr_BoyerLindquist) then call Exact__Kerr_BoyerLindquist( @@ -103,7 +119,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Kerr_KerrSchild) then call Exact__Kerr_KerrSchild( @@ -111,7 +128,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Schwarzschild_Lemaitre) then call Exact__Schwarzschild_Lemaitre( @@ -119,7 +137,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__multi_BH) then call Exact__multi_BH( @@ -127,7 +146,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) c c not fully implemented yet -- see Nina Jansen for details @@ -138,7 +158,8 @@ c $ x, y, z, t, c $ gdtt, gdtx, gdty, gdtz, c $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, c $ gutt, gutx, guty, gutz, -c $ guxx, guyy, guzz, guxy, guyz, guxz) +c $ guxx, guyy, guzz, guxy, guyz, guxz, +c $ psi) c elseif (decoded_exact_model .eq. EXACT__Thorne_fakebinary) then @@ -147,7 +168,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc @@ -161,7 +183,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Robertson_Walker) then call Exact__Robertson_Walker( @@ -169,7 +192,7 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz,rama) + $ guxx, guyy, guzz, guxy, guyz, guxz, psi, rama) elseif (decoded_exact_model .eq. EXACT__de_Sitter) then call Exact__de_Sitter( @@ -177,7 +200,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__de_Sitter_Lambda) then call Exact__de_Sitter_Lambda( @@ -185,7 +209,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__anti_de_Sitter_Lambda) then call Exact__anti_de_Sitter_Lambda( @@ -193,7 +218,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Bianchi_I) then call Exact__Bianchi_I( @@ -201,7 +227,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Goedel) then call Exact__Goedel( @@ -209,7 +236,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Bertotti) then call Exact__Bertotti( @@ -217,7 +245,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Kasner_like) then call Exact__Kasner_like( @@ -225,7 +254,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Kasner_axisymmetric) then call Exact__Kasner_axisymmetric( @@ -233,7 +263,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Kasner_generalized) then call Exact__Kasner_generalized( @@ -241,7 +272,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Gowdy_wave) then call Exact__Gowdy_wave( @@ -249,7 +281,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__Milne) then call Exact__Milne( @@ -257,7 +290,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc @@ -271,7 +305,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__bowl) then call Exact__bowl( @@ -279,7 +314,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) elseif (decoded_exact_model .eq. EXACT__constant_density_star) then call Exact__constant_density_star( @@ -287,7 +323,8 @@ c $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc diff --git a/src/metrics/Alvi.F77 b/src/metrics/Alvi.F77 index 81789e6..7ee21fb 100644 --- a/src/metrics/Alvi.F77 +++ b/src/metrics/Alvi.F77 @@ -11,7 +11,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -30,7 +31,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi CCTK_REAL rin1,rin2,rout,x1,x2, r1, r2, r3, results(10) diff --git a/src/metrics/Bertotti.F77 b/src/metrics/Bertotti.F77 index f9b32df..5533422 100644 --- a/src/metrics/Bertotti.F77 +++ b/src/metrics/Bertotti.F77 @@ -10,7 +10,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -20,7 +21,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/Bianchi_I.F77 b/src/metrics/Bianchi_I.F77 index 0af7b12..feece73 100644 --- a/src/metrics/Bianchi_I.F77 +++ b/src/metrics/Bianchi_I.F77 @@ -12,7 +12,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -22,7 +23,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/Goedel.F77 b/src/metrics/Goedel.F77 index fec2c73..513e3e1 100644 --- a/src/metrics/Goedel.F77 +++ b/src/metrics/Goedel.F77 @@ -12,7 +12,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -22,7 +23,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/Gowdy.F77 b/src/metrics/Gowdy.F77 index 06370bf..8573710 100644 --- a/src/metrics/Gowdy.F77 +++ b/src/metrics/Gowdy.F77 @@ -20,7 +20,8 @@ C @@ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -30,7 +31,8 @@ C @@ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/Kasner_axisymmetric.F77 b/src/metrics/Kasner_axisymmetric.F77 index 0dc98d1..d6480dd 100644 --- a/src/metrics/Kasner_axisymmetric.F77 +++ b/src/metrics/Kasner_axisymmetric.F77 @@ -10,7 +10,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -20,7 +21,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi gdtt = -t**(-0.5d0) gdtx = 0.d0 diff --git a/src/metrics/Kasner_generalized.F77 b/src/metrics/Kasner_generalized.F77 index 9ec71d5..c3e8a63 100644 --- a/src/metrics/Kasner_generalized.F77 +++ b/src/metrics/Kasner_generalized.F77 @@ -10,7 +10,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -20,7 +21,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/Kasner_like.F77 b/src/metrics/Kasner_like.F77 index 15e5cff..bd24380 100644 --- a/src/metrics/Kasner_like.F77 +++ b/src/metrics/Kasner_like.F77 @@ -14,7 +14,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -24,7 +25,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/Kerr_BoyerLindquist.F77 b/src/metrics/Kerr_BoyerLindquist.F77 index 3a16cc6..076d9f8 100644 --- a/src/metrics/Kerr_BoyerLindquist.F77 +++ b/src/metrics/Kerr_BoyerLindquist.F77 @@ -11,7 +11,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -21,7 +22,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/Kerr_KerrSchild.F77 b/src/metrics/Kerr_KerrSchild.F77 index d998fb4..ad9ac27 100644 --- a/src/metrics/Kerr_KerrSchild.F77 +++ b/src/metrics/Kerr_KerrSchild.F77 @@ -12,7 +12,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -23,7 +24,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi CCTK_REAL gamma, t0, z0, x0, y0, rho02, r02, r0, costheta0, $ lt0, lx0, ly0, lz0, hh, lt, lx, ly, lz diff --git a/src/metrics/Lemaitre.F77 b/src/metrics/Lemaitre.F77 index 2062f4b..9819c45 100644 --- a/src/metrics/Lemaitre.F77 +++ b/src/metrics/Lemaitre.F77 @@ -12,7 +12,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -24,7 +25,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/Milne.F77 b/src/metrics/Milne.F77 index 0167873..d441cb6 100644 --- a/src/metrics/Milne.F77 +++ b/src/metrics/Milne.F77 @@ -12,7 +12,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -26,6 +27,7 @@ c Output. $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi, $ coef, x2,y2,z2,t2 x2= x*x diff --git a/src/metrics/Minkowski.F77 b/src/metrics/Minkowski.F77 index 5318332..88031c2 100644 --- a/src/metrics/Minkowski.F77 +++ b/src/metrics/Minkowski.F77 @@ -8,7 +8,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -21,7 +22,8 @@ c Output. CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi gdtt = -1.d0 gdtx = 0.d0 diff --git a/src/metrics/Minkowski_conf_wave.F77 b/src/metrics/Minkowski_conf_wave.F77 new file mode 100644 index 0000000..301fb62 --- /dev/null +++ b/src/metrics/Minkowski_conf_wave.F77 @@ -0,0 +1,75 @@ +c $Header$ + +C Author: Frank Loeffler (frank.loeffler@aei.mpg.de) +C Licence: GPL 2 or later + +#include "cctk.h" +#include "cctk_Arguments.h" +#include "cctk_Parameters.h" + +#define Pi (4 * atan(1.d0)) + + + subroutine Exact__Minkowski_conf_wave( + $ x, y, z, t, + $ gdtt, gdtx, gdty, gdtz, + $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, + $ gutt, gutx, guty, gutz, + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) + + implicit none + DECLARE_CCTK_PARAMETERS + DECLARE_CCTK_FUNCTIONS + + CCTK_REAL x, y, z, t + CCTK_REAL gdtt, gdtx, gdty, gdtz, + $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, + $ gutt, gutx, guty, gutz, + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi + +C write conformal factor + if (Minkowski_conf_wave__direction .eq. 0) then + psi = Minkowski_conf_wave__amplitude * + . sin(2.0d0*Pi/Minkowski_conf_wave__wavelength* x ) + 1.0d0 + else if (Minkowski_conf_wave__direction .eq. 1) then + psi = Minkowski_conf_wave__amplitude * + . sin(2.0d0*Pi/Minkowski_conf_wave__wavelength* y ) + 1.0d0 + else if (Minkowski_conf_wave__direction .eq. 2) then + psi = Minkowski_conf_wave__amplitude * + . sin(2.0d0*Pi/Minkowski_conf_wave__wavelength* z ) + 1.0d0 + end if + +C write metric. + + gdxx = psi**(-4.0d0) + gdyy = gdxx + gdzz = gdxx + + gdxy = 0.0d0 + gdyz = 0.0d0 + gdzx = 0.0d0 + + gdtt = -1.0d0 + gdtx = 0.0d0 + gdty = 0.0d0 + gdtz = 0.0d0 + +C and upper metric. + + guxx = psi**4.0d0 + guyy = guxx + guzz = guxx + + guxy = 0.0d0 + guyz = 0.0d0 + guzx = 0.0d0 + + gutt = -1.0d0 + gutx = 0.0d0 + guty = 0.0d0 + gutz = 0.0d0 + + return + end diff --git a/src/metrics/Minkowski_funny.F77 b/src/metrics/Minkowski_funny.F77 index 5a9c8e2..a1b5af7 100644 --- a/src/metrics/Minkowski_funny.F77 +++ b/src/metrics/Minkowski_funny.F77 @@ -19,7 +19,8 @@ c $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -34,7 +35,8 @@ c Output. CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi c Internal. diff --git a/src/metrics/Minkowski_gauge_wave.F77 b/src/metrics/Minkowski_gauge_wave.F77 index 5e67a70..a4f8a6f 100644 --- a/src/metrics/Minkowski_gauge_wave.F77 +++ b/src/metrics/Minkowski_gauge_wave.F77 @@ -12,7 +12,8 @@ c $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none DECLARE_CCTK_PARAMETERS @@ -22,7 +23,8 @@ c $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/Minkowski_shift.F77 b/src/metrics/Minkowski_shift.F77 index f999d15..ab8bca4 100644 --- a/src/metrics/Minkowski_shift.F77 +++ b/src/metrics/Minkowski_shift.F77 @@ -21,7 +21,8 @@ c $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -36,7 +37,8 @@ c Output. CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi c Internal. diff --git a/src/metrics/Robertson_Walker.F77 b/src/metrics/Robertson_Walker.F77 index 85e0739..8335762 100644 --- a/src/metrics/Robertson_Walker.F77 +++ b/src/metrics/Robertson_Walker.F77 @@ -13,7 +13,7 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx,rama) + $ guxx, guyy, guzz, guxy, guyz, guzx, psi, rama) implicit none @@ -25,7 +25,7 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx,rama + $ guxx, guyy, guzz, guxy, guyz, guzx, psi, rama logical firstcall diff --git a/src/metrics/Schwarzschild_EF.F77 b/src/metrics/Schwarzschild_EF.F77 index a53858d..4f69802 100644 --- a/src/metrics/Schwarzschild_EF.F77 +++ b/src/metrics/Schwarzschild_EF.F77 @@ -10,7 +10,8 @@ c $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -20,7 +21,8 @@ c $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/Schwarzschild_Lemaitre.F77 b/src/metrics/Schwarzschild_Lemaitre.F77 index 9ff7835..6d9d033 100644 --- a/src/metrics/Schwarzschild_Lemaitre.F77 +++ b/src/metrics/Schwarzschild_Lemaitre.F77 @@ -12,7 +12,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -22,7 +23,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/Schwarzschild_Novikov.F77 b/src/metrics/Schwarzschild_Novikov.F77 index 61a2d1a..38dafc7 100644 --- a/src/metrics/Schwarzschild_Novikov.F77 +++ b/src/metrics/Schwarzschild_Novikov.F77 @@ -14,7 +14,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -29,7 +30,8 @@ c Output. CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi c Internal. logical firstcall diff --git a/src/metrics/Schwarzschild_PG.F77 b/src/metrics/Schwarzschild_PG.F77 index 37072fc..13db0f1 100644 --- a/src/metrics/Schwarzschild_PG.F77 +++ b/src/metrics/Schwarzschild_PG.F77 @@ -11,7 +11,8 @@ c $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -21,7 +22,8 @@ c $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/Thorne_fakebinary.F77 b/src/metrics/Thorne_fakebinary.F77 index 2cd954b..8c3f842 100644 --- a/src/metrics/Thorne_fakebinary.F77 +++ b/src/metrics/Thorne_fakebinary.F77 @@ -15,7 +15,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz) + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi) implicit none DECLARE_CCTK_PARAMETERS @@ -28,7 +29,8 @@ C output CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guxz + $ guxx, guyy, guzz, guxy, guyz, guxz, + $ psi C static diff --git a/src/metrics/anti_de_Sitter_Lambda.F77 b/src/metrics/anti_de_Sitter_Lambda.F77 index 5e2733c..483d847 100644 --- a/src/metrics/anti_de_Sitter_Lambda.F77 +++ b/src/metrics/anti_de_Sitter_Lambda.F77 @@ -10,7 +10,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -20,7 +21,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/boost_rotation_symmetric.F77 b/src/metrics/boost_rotation_symmetric.F77 index ec3dbf9..0a69f8d 100644 --- a/src/metrics/boost_rotation_symmetric.F77 +++ b/src/metrics/boost_rotation_symmetric.F77 @@ -10,7 +10,8 @@ c $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -20,7 +21,8 @@ c $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi CCTK_REAL a, b, mu0, mu1, lam1, mu2, lam2, $ lam3, mu4, lam4, mu5, lam5, num, div, f, $ elam, emu0, delta, gfunc, tmp diff --git a/src/metrics/bowl.F77 b/src/metrics/bowl.F77 index 192c4e0..36a1cf7 100644 --- a/src/metrics/bowl.F77 +++ b/src/metrics/bowl.F77 @@ -57,7 +57,8 @@ c $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -71,7 +72,8 @@ c Output. CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall,evolve diff --git a/src/metrics/constant_density_star.F77 b/src/metrics/constant_density_star.F77 index ab3aa30..aecb4ab 100644 --- a/src/metrics/constant_density_star.F77 +++ b/src/metrics/constant_density_star.F77 @@ -15,7 +15,8 @@ c $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) c c The metric is given as a conformally flat metric. @@ -39,7 +40,8 @@ c Output. CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi c Internal. diff --git a/src/metrics/de_Sitter.F77 b/src/metrics/de_Sitter.F77 index 2ced3e5..f073d73 100644 --- a/src/metrics/de_Sitter.F77 +++ b/src/metrics/de_Sitter.F77 @@ -16,7 +16,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -26,7 +27,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/de_Sitter_Lambda.F77 b/src/metrics/de_Sitter_Lambda.F77 index 381159b..868557a 100644 --- a/src/metrics/de_Sitter_Lambda.F77 +++ b/src/metrics/de_Sitter_Lambda.F77 @@ -10,7 +10,8 @@ C $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -20,7 +21,8 @@ C $Header$ CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/metrics/make.code.defn b/src/metrics/make.code.defn index 9c33332..ccd943d 100644 --- a/src/metrics/make.code.defn +++ b/src/metrics/make.code.defn @@ -10,6 +10,7 @@ SRCS = Minkowski.F77 \ Minkowski_funny.F77 \ Minkowski_shift.F77 \ Minkowski_gauge_wave.F77 \ + Minkowski_conf_wave.F77 \ \ Schwarzschild_EF.F77 \ Schwarzschild_PG.F77 \ diff --git a/src/metrics/multi_BH.F77 b/src/metrics/multi_BH.F77 index a942c67..0068db0 100644 --- a/src/metrics/multi_BH.F77 +++ b/src/metrics/multi_BH.F77 @@ -24,7 +24,8 @@ c $Header$ $ gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx) + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi) implicit none @@ -37,7 +38,8 @@ c Output. CCTK_REAL gdtt, gdtx, gdty, gdtz, $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, $ gutt, gutx, guty, gutz, - $ guxx, guyy, guzz, guxy, guyz, guzx + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi logical firstcall diff --git a/src/slice_data.F b/src/slice_data.F index f7edb5e..088f635 100644 --- a/src/slice_data.F +++ b/src/slice_data.F @@ -18,7 +18,7 @@ c #define-ing the symbol EXACT_NO_F90 will turn this subroutine into a no-op CCTK_REAL s1d(4,3), nd(4), nu(4), norm, gd(4,4), gu(4,4), g3(3,3), $ gd_p(4,4), gd_m(4,4), gd1d(4,4,4), s2d(4,3,3), k3(3,3), - $ ex_eps, d3(3,3,3), dx, dy, dz + $ ex_eps, d3(3,3,3), dx, dy, dz, exact_psi parameter (ex_eps=1.d-6) C Grid parameters. @@ -55,6 +55,13 @@ C Calculate first derivatives of slice coordinates. s1d(4,2) = 0.5d0 * (slicet(i,j+1,k) - slicet(i,j-1,k))/dy s1d(4,3) = 0.5d0 * (slicet(i,j,k+1) - slicet(i,j,k-1))/dz +C Preset the conformal factor + if (conformal_state .ne. 0) then + exact_psi=1.0d0 + else + exact_psi=0.0d0 + end if + C Now we need the exact solution metric in the preferred coordinates x^A. call Exact__metric( @@ -66,7 +73,7 @@ C Now we need the exact solution metric in the preferred coordinates x^A. $ gd(1,2), gd(2,3), gd(1,3), $ gu(4,4), gu(1,4), gu(2,4), gu(3,4), $ gu(1,1), gu(2,2), gu(3,3), - $ gu(1,2), gu(2,3), gu(1,3)) + $ gu(1,2), gu(2,3), gu(1,3), exact_psi) C Calculate n^A and dx^A/dt @@ -165,6 +172,12 @@ C with an editor (hint for proofreading). C Calculate g_AB,C. Need to sum explicitly over C. Do this with C the editor. +C Preset the conformal factor + if (conformal_state .ne. 0) then + exact_psi=1.0d0 + else + exact_psi=0.0d0 + end if call Exact__metric( $ decoded_exact_model, $ slicex(i,j,k)+ex_eps, slicey(i,j,k), slicez(i,j,k), @@ -174,7 +187,14 @@ C the editor. $ gd_p(1,2), gd_p(2,3), gd_p(1,3), $ gu(4,4), gu(1,4), gu(2,4), gu(3,4), $ gu(1,1), gu(2,2), gu(3,3), - $ gu(1,2), gu(2,3), gu(1,3)) + $ gu(1,2), gu(2,3), gu(1,3), exact_psi) + +C Preset the conformal factor + if (conformal_state .ne. 0) then + exact_psi=1.0d0 + else + exact_psi=0.0d0 + end if call Exact__metric( $ decoded_exact_model, $ slicex(i,j,k)-ex_eps, slicey(i,j,k), slicez(i,j,k), @@ -184,13 +204,19 @@ C the editor. $ gd_m(1,2), gd_m(2,3), gd_m(1,3), $ gu(4,4), gu(1,4), gu(2,4), gu(3,4), $ gu(1,1), gu(2,2), gu(3,3), - $ gu(1,2), gu(2,3), gu(1,3)) + $ gu(1,2), gu(2,3), gu(1,3), exact_psi) do m=1,4 do n=m,4 gd1d(m,n,1) = (gd_p(m,n) - gd_m(m,n)) / (2.*ex_eps) end do end do +C Preset the conformal factor + if (conformal_state .ne. 0) then + exact_psi=1.0d0 + else + exact_psi=0.0d0 + end if call Exact__metric( $ decoded_exact_model, $ slicex(i,j,k), slicey(i,j,k)+ex_eps, slicez(i,j,k), @@ -200,7 +226,14 @@ C the editor. $ gd_p(1,2), gd_p(2,3), gd_p(1,3), $ gu(4,4), gu(1,4), gu(2,4), gu(3,4), $ gu(1,1), gu(2,2), gu(3,3), - $ gu(1,2), gu(2,3), gu(1,3)) + $ gu(1,2), gu(2,3), gu(1,3), exact_psi) + +C Preset the conformal factor + if (conformal_state .ne. 0) then + exact_psi=1.0d0 + else + exact_psi=0.0d0 + end if call Exact__metric( $ decoded_exact_model, $ slicex(i,j,k), slicey(i,j,k)-ex_eps, slicez(i,j,k), @@ -210,13 +243,20 @@ C the editor. $ gd_m(1,2), gd_m(2,3), gd_m(1,3), $ gu(4,4), gu(1,4), gu(2,4), gu(3,4), $ gu(1,1), gu(2,2), gu(3,3), - $ gu(1,2), gu(2,3), gu(1,3)) + $ gu(1,2), gu(2,3), gu(1,3), exact_psi) do m=1,4 do n=m,4 gd1d(m,n,2) = (gd_p(m,n) - gd_m(m,n)) / (2.*ex_eps) end do end do + +C Preset the conformal factor + if (conformal_state .ne. 0) then + exact_psi=1.0d0 + else + exact_psi=0.0d0 + end if call Exact__metric( $ decoded_exact_model, $ slicex(i,j,k), slicey(i,j,k), slicez(i,j,k)+ex_eps, @@ -226,7 +266,14 @@ C the editor. $ gd_p(1,2), gd_p(2,3), gd_p(1,3), $ gu(4,4), gu(1,4), gu(2,4), gu(3,4), $ gu(1,1), gu(2,2), gu(3,3), - $ gu(1,2), gu(2,3), gu(1,3)) + $ gu(1,2), gu(2,3), gu(1,3), exact_psi) + +C Preset the conformal factor + if (conformal_state .ne. 0) then + exact_psi=1.0d0 + else + exact_psi=0.0d0 + end if call Exact__metric( $ decoded_exact_model, $ slicex(i,j,k), slicey(i,j,k), slicez(i,j,k)-ex_eps, @@ -236,13 +283,20 @@ C the editor. $ gd_m(1,2), gd_m(2,3), gd_m(1,3), $ gu(4,4), gu(1,4), gu(2,4), gu(3,4), $ gu(1,1), gu(2,2), gu(3,3), - $ gu(1,2), gu(2,3), gu(1,3)) + $ gu(1,2), gu(2,3), gu(1,3), exact_psi) do m=1,4 do n=m,4 gd1d(m,n,3) = (gd_p(m,n) - gd_m(m,n)) / (2.*ex_eps) end do end do + +C Preset the conformal factor + if (conformal_state .ne. 0) then + exact_psi=1.0d0 + else + exact_psi=0.0d0 + end if call Exact__metric( $ decoded_exact_model, $ slicex(i,j,k), slicey(i,j,k), slicez(i,j,k), @@ -252,7 +306,14 @@ C the editor. $ gd_p(1,2), gd_p(2,3), gd_p(1,3), $ gu(4,4), gu(1,4), gu(2,4), gu(3,4), $ gu(1,1), gu(2,2), gu(3,3), - $ gu(1,2), gu(2,3), gu(1,3)) + $ gu(1,2), gu(2,3), gu(1,3), exact_psi) + +C Preset the conformal factor + if (conformal_state .ne. 0) then + exact_psi=1.0d0 + else + exact_psi=0.0d0 + end if call Exact__metric( $ decoded_exact_model, $ slicex(i,j,k), slicey(i,j,k), slicez(i,j,k), @@ -262,7 +323,7 @@ C the editor. $ gd_m(1,2), gd_m(2,3), gd_m(1,3), $ gu(4,4), gu(1,4), gu(2,4), gu(3,4), $ gu(1,1), gu(2,2), gu(3,3), - $ gu(1,2), gu(2,3), gu(1,3)) + $ gu(1,2), gu(2,3), gu(1,3), exact_psi) do m=1,4 do n=m,4 gd1d(m,n,4) = (gd_p(m,n) - gd_m(m,n)) / (2.*ex_eps) diff --git a/src/xyz_blended_boundary.F77 b/src/xyz_blended_boundary.F77 index 9c6ef44..ee0cd41 100644 --- a/src/xyz_blended_boundary.F77 +++ b/src/xyz_blended_boundary.F77 @@ -33,6 +33,11 @@ C $Header$ CCTK_REAL bxxe,bxye,bxze,byxe,byye,byze,bzxe,bzye,bzze CCTK_REAL det, uxx, uxy, uxz, uyy, uyz, uzz CCTK_REAL vxe,vye,vze,sav + CCTK_REAL + $ exact_psi, + $ exact_psix, exact_psiy, exact_psiz, + $ exact_psixx, exact_psiyy, exact_psizz, + $ exact_psixy, exact_psiyz, exact_psixz CCTK_REAL dx,dy,dz,time integer ierr @@ -86,11 +91,32 @@ c We only do anything if in the blending region $ z(i,j,k) .ge. zmax - zblend .or. $ z(i,j,k) .le. zmin + zblend) then +C Initialize the psi of exact +C (also to tell the models about the conformal_state) + if (conformal_state .ne. 0) then + exact_psi = 1.0D0 + else + exact_psi = 0.0D0 + end if + exact_psix = 0.0D0 + exact_psiy = 0.0D0 + exact_psiz = 0.0D0 + exact_psixx = 0.0D0 + exact_psixx = 0.0D0 + exact_psizz = 0.0D0 + exact_psixy = 0.0D0 + exact_psiyz = 0.0D0 + exact_psixz = 0.0D0 + call Exact__Bona_Masso_data( $ decoded_exact_model, $ x(i,j,k), y(i,j,k), z(i,j,k), time, $ gxxe, gyye, gzze, gxye, gyze, gxze, $ kxxe, kyye, kzze, kxye, kyze, kxze, + $ exact_psi, + $ exact_psix, exact_psiy, exact_psiz, + $ exact_psixx, exact_psiyy, exact_psizz, + $ exact_psixy, exact_psiyz, exact_psixz, $ dxgxxe, dxgyye, dxgzze, dxgxye, dxgyze, dxgxze, $ dygxxe, dygyye, dygzze, dygxye, dygyze, dygxze, $ dzgxxe, dzgyye, dzgzze, dzgxye, dzgyze, dzgxze, diff --git a/test/schwarzschild_ef.par b/test/schwarzschild_ef.par new file mode 100644 index 0000000..bd8fc16 --- /dev/null +++ b/test/schwarzschild_ef.par @@ -0,0 +1,54 @@ +# Initial data and short evolution (ADM) for DeSitter spacetime +# Author: Mitica Vulcanov <vulcan@aei.mpg.de> +# $Header$ + +# Author: Frank Loeffler (frank.loeffler@aei.mpg.de) + +# Required thorns +ActiveThorns = "CartGrid3D Boundary Time PUGH PUGHReduce PUGHSlab IOAscii IOUtil IOBasic ADMBase ADMMacros ADMConstraints MoL BSSN_MoL ADMCoupling StaticConformal CoordGauge SpaceMask Exact NaNChecker" + +# Grid +driver::global_nx = 33 +driver::global_ny = 33 +driver::global_nz = 33 +grid::type = "byrange" +grid::domain = "full" + +grid::xmin = -5.0 +grid::xmax = 5.1 +grid::ymin = -5.0 +grid::ymax = 5.1 +grid::zmin = -5.0 +grid::zmax = 5.1 + +cactus::cctk_itlast = 0 +time::dtfac = 0.25 + +# Initial data +ADMBase::initial_data = "exact" +ADMBase::initial_shift = "exact" +Exact::exact_model = "Schwarzschild/EF" + +# Gauge +ADMBase::lapse_evolution_method = "exact" +ADMBase::shift_evolution_method = "exact" + +# Evolution +ADMBase::evolution_method = "adm_bssn" +ADM_BSSN::bound = "radiative" +boundary::radpower = 2 + +methodoflines::ode_method = "icn" + +# Output +IO::out_dir = "schwarzschild_ef" + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = "admconstraints::ham" + +IOBasic::outScalar_every = 1 +IOBasic::outScalar_vars = "admbase::gxx admbase::gyy admbase::gzz admbase::gxy admconstraints::ham admbase::lapse" + +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = "admbase::gxx admbase::gyy admbase::gzz admbase::gxy admconstraints::ham admbase::lapse" + diff --git a/test/schwarzschild_ef/alp_3D_diagonal.xg b/test/schwarzschild_ef/alp_3D_diagonal.xg new file mode 100644 index 0000000..1b50c87 --- /dev/null +++ b/test/schwarzschild_ef/alp_3D_diagonal.xg @@ -0,0 +1,40 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label diagonal +"y-label ADMBASE::alp + + +"Time = 0.0000000000000 +-8.6602540378444 0.9013252521712 +-8.1135755017055 0.8956818637024 +-7.5668969655665 0.8893513361568 +-7.0202184294276 0.8821994554954 +-6.4735398932887 0.8740544236021 +-5.9268613571498 0.8646926211453 +-5.3801828210108 0.8538173125557 +-4.8335042848719 0.8410257993211 +-4.2868257487330 0.8257568854970 +-3.7401472125940 0.8072031220124 +-3.1934686764551 0.7841561845659 +-2.6467901403162 0.7547154219106 +-2.1001116041773 0.7156873618401 +-1.5534330680383 0.6611837984071 +-1.0067545318994 0.5786458057228 +-0.4600759957605 0.4324546046304 +0.0866025403784 0.2037255270177 +0.6332810765174 0.4903990794737 +1.1799596126563 0.6091479250571 +1.7266381487952 0.6806785316502 +2.2733166849342 0.7293692178168 +2.8199952210731 0.7648933450237 +3.3666737572120 0.7920414666707 +3.9133522933509 0.8135000735488 +4.4600308294899 0.8309055026211 +5.0067093656288 0.8453160969210 +5.5533879017677 0.8574486299279 +6.1000664379066 0.8678066864076 +6.6467449740456 0.8767548618666 +7.1934235101845 0.8845638390699 +7.7401020463234 0.8914389167433 +8.2867805824623 0.8975386961475 +8.8334591186013 0.9029876848104 diff --git a/test/schwarzschild_ef/alp_maximum.xg b/test/schwarzschild_ef/alp_maximum.xg new file mode 100644 index 0000000..a2414a5 --- /dev/null +++ b/test/schwarzschild_ef/alp_maximum.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::alp +"alp v time +0.0000000000000 0.9029876848104 diff --git a/test/schwarzschild_ef/alp_minimum.xg b/test/schwarzschild_ef/alp_minimum.xg new file mode 100644 index 0000000..ba7582b --- /dev/null +++ b/test/schwarzschild_ef/alp_minimum.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label time +"y-label ADMBASE::alp +"alp v time +0.0000000000000 0.2037255270177 diff --git a/test/schwarzschild_ef/alp_norm1.xg b/test/schwarzschild_ef/alp_norm1.xg new file mode 100644 index 0000000..828150e --- /dev/null +++ b/test/schwarzschild_ef/alp_norm1.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::alp +"alp v time +0.0000000000000 0.8345487599747 diff --git a/test/schwarzschild_ef/alp_norm2.xg b/test/schwarzschild_ef/alp_norm2.xg new file mode 100644 index 0000000..3a64bd8 --- /dev/null +++ b/test/schwarzschild_ef/alp_norm2.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::alp +"alp v time +0.0000000000000 0.8360803073831 diff --git a/test/schwarzschild_ef/alp_x_[16][16].xg b/test/schwarzschild_ef/alp_x_[16][16].xg new file mode 100644 index 0000000..8fdeec0 --- /dev/null +++ b/test/schwarzschild_ef/alp_x_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label x +"y-label ADMBASE::alp (y = 0.0500000000000, z = 0.0500000000000), (yi = 16, zi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 0.8451663268086 +-4.6843750000000 0.8371491202724 +-4.3687500000000 0.8282480583464 +-4.0531250000000 0.8183066754501 +-3.7375000000000 0.8071287088617 +-3.4218750000000 0.7944643560816 +-3.1062500000000 0.7799902886628 +-2.7906250000000 0.7632797531455 +-2.4750000000000 0.7437563197794 +-2.1593750000000 0.7206193824028 +-1.8437500000000 0.6927179917005 +-1.5281250000000 0.6583230648496 +-1.2125000000000 0.6146795018325 +-0.8968750000000 0.5570130331675 +-0.5812500000000 0.4758842147107 +-0.2656250000000 0.3476076502339 +0.0500000000000 0.2037255270177 +0.3656250000000 0.3961966505475 +0.6812500000000 0.5050703361853 +0.9968750000000 0.5772308535168 +1.3125000000000 0.6297405399634 +1.6281250000000 0.6700628103203 +1.9437500000000 0.7021636168282 +2.2593750000000 0.7284019913708 +2.5750000000000 0.7502894531721 +2.8906250000000 0.7688475843882 +3.2062500000000 0.7847954234481 +3.5218750000000 0.7986555817630 +3.8375000000000 0.8108179872091 +4.1531250000000 0.8215800036862 +4.4687500000000 0.8311726601213 +4.7843750000000 0.8397783498435 +5.1000000000000 0.8475430993433 diff --git a/test/schwarzschild_ef/alp_y_[16][16].xg b/test/schwarzschild_ef/alp_y_[16][16].xg new file mode 100644 index 0000000..2e9a179 --- /dev/null +++ b/test/schwarzschild_ef/alp_y_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label y +"y-label ADMBASE::alp (x = 0.0500000000000, z = 0.0500000000000), (xi = 16, zi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 0.8451663268086 +-4.6843750000000 0.8371491202724 +-4.3687500000000 0.8282480583464 +-4.0531250000000 0.8183066754501 +-3.7375000000000 0.8071287088617 +-3.4218750000000 0.7944643560816 +-3.1062500000000 0.7799902886628 +-2.7906250000000 0.7632797531455 +-2.4750000000000 0.7437563197794 +-2.1593750000000 0.7206193824028 +-1.8437500000000 0.6927179917005 +-1.5281250000000 0.6583230648496 +-1.2125000000000 0.6146795018325 +-0.8968750000000 0.5570130331675 +-0.5812500000000 0.4758842147107 +-0.2656250000000 0.3476076502339 +0.0500000000000 0.2037255270177 +0.3656250000000 0.3961966505475 +0.6812500000000 0.5050703361853 +0.9968750000000 0.5772308535168 +1.3125000000000 0.6297405399634 +1.6281250000000 0.6700628103203 +1.9437500000000 0.7021636168282 +2.2593750000000 0.7284019913708 +2.5750000000000 0.7502894531721 +2.8906250000000 0.7688475843882 +3.2062500000000 0.7847954234481 +3.5218750000000 0.7986555817630 +3.8375000000000 0.8108179872091 +4.1531250000000 0.8215800036862 +4.4687500000000 0.8311726601213 +4.7843750000000 0.8397783498435 +5.1000000000000 0.8475430993433 diff --git a/test/schwarzschild_ef/alp_z_[16][16].xg b/test/schwarzschild_ef/alp_z_[16][16].xg new file mode 100644 index 0000000..9f99b1f --- /dev/null +++ b/test/schwarzschild_ef/alp_z_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label z +"y-label ADMBASE::alp (x = 0.0500000000000, y = 0.0500000000000), (xi = 16, yi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 0.8451663268086 +-4.6843750000000 0.8371491202724 +-4.3687500000000 0.8282480583464 +-4.0531250000000 0.8183066754501 +-3.7375000000000 0.8071287088617 +-3.4218750000000 0.7944643560816 +-3.1062500000000 0.7799902886628 +-2.7906250000000 0.7632797531455 +-2.4750000000000 0.7437563197794 +-2.1593750000000 0.7206193824028 +-1.8437500000000 0.6927179917005 +-1.5281250000000 0.6583230648496 +-1.2125000000000 0.6146795018325 +-0.8968750000000 0.5570130331675 +-0.5812500000000 0.4758842147107 +-0.2656250000000 0.3476076502339 +0.0500000000000 0.2037255270177 +0.3656250000000 0.3961966505475 +0.6812500000000 0.5050703361853 +0.9968750000000 0.5772308535168 +1.3125000000000 0.6297405399634 +1.6281250000000 0.6700628103203 +1.9437500000000 0.7021636168282 +2.2593750000000 0.7284019913708 +2.5750000000000 0.7502894531721 +2.8906250000000 0.7688475843882 +3.2062500000000 0.7847954234481 +3.5218750000000 0.7986555817630 +3.8375000000000 0.8108179872091 +4.1531250000000 0.8215800036862 +4.4687500000000 0.8311726601213 +4.7843750000000 0.8397783498435 +5.1000000000000 0.8475430993433 diff --git a/test/schwarzschild_ef/gxx_3D_diagonal.xg b/test/schwarzschild_ef/gxx_3D_diagonal.xg new file mode 100644 index 0000000..95cef73 --- /dev/null +++ b/test/schwarzschild_ef/gxx_3D_diagonal.xg @@ -0,0 +1,40 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label diagonal +"y-label ADMBASE::gxx + + +"Time = 0.0000000000000 +-8.6602540378444 1.0769800358920 +-8.1135755017055 1.0821668161622 +-7.5668969655665 1.0881030453699 +-7.0202184294276 1.0949638068058 +-6.4735398932887 1.1029833256080 +-5.9268613571498 1.1124822442257 +-5.3801828210108 1.1239115265866 +-4.8335042848719 1.1379261561334 +-4.2868257487330 1.1555152240241 +-3.7401472125940 1.1782461033678 +-3.1934686764551 1.2087594193680 +-2.6467901403162 1.2518774180514 +-2.1001116041773 1.3174434469771 +-1.5534330680383 1.4291569945196 +-1.0067545318994 1.6621938571351 +-0.4600759957605 2.4490359697308 +0.0866025403784 8.6980035891950 +0.6332810765174 2.0527184395480 +1.1799596126563 1.5649910891152 +1.7266381487952 1.3861067630944 +2.2733166849342 1.2932572795884 +2.8199952210731 1.2364070200137 +3.3666737572120 1.1980193849311 +3.9133522933509 1.1703569259020 +4.4600308294899 1.1494757978484 +5.0067093656288 1.1331546566780 +5.5533879017677 1.1200468395976 +6.1000664379066 1.1092884271758 +6.6467449740456 1.1002997210319 +7.1934235101845 1.0926772441137 +7.7401020463234 1.0861315086903 +8.2867805824623 1.0804494170001 +8.8334591186013 1.0754706234235 diff --git a/test/schwarzschild_ef/gxx_maximum.xg b/test/schwarzschild_ef/gxx_maximum.xg new file mode 100644 index 0000000..f8f6e24 --- /dev/null +++ b/test/schwarzschild_ef/gxx_maximum.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gxx +"gxx v time +0.0000000000000 8.6980035891950 diff --git a/test/schwarzschild_ef/gxx_minimum.xg b/test/schwarzschild_ef/gxx_minimum.xg new file mode 100644 index 0000000..618d0de --- /dev/null +++ b/test/schwarzschild_ef/gxx_minimum.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gxx +"gxx v time +0.0000000000000 1.0000133254896 diff --git a/test/schwarzschild_ef/gxx_norm1.xg b/test/schwarzschild_ef/gxx_norm1.xg new file mode 100644 index 0000000..2dd06ae --- /dev/null +++ b/test/schwarzschild_ef/gxx_norm1.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gxx +"gxx v time +0.0000000000000 1.1524177407289 diff --git a/test/schwarzschild_ef/gxx_norm2.xg b/test/schwarzschild_ef/gxx_norm2.xg new file mode 100644 index 0000000..b32be46 --- /dev/null +++ b/test/schwarzschild_ef/gxx_norm2.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gxx +"gxx v time +0.0000000000000 1.1660598040677 diff --git a/test/schwarzschild_ef/gxx_x_[16][16].xg b/test/schwarzschild_ef/gxx_x_[16][16].xg new file mode 100644 index 0000000..09248e1 --- /dev/null +++ b/test/schwarzschild_ef/gxx_x_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label x +"y-label ADMBASE::gxx (y = 0.0500000000000, z = 0.0500000000000), (yi = 16, zi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 1.3998800299930 +-4.6843750000000 1.4268054152043 +-4.3687500000000 1.4576170162950 +-4.0531250000000 1.4932212211925 +-3.7375000000000 1.5348298773928 +-3.4218750000000 1.5841007174696 +-3.1062500000000 1.6433630284612 +-2.7906250000000 1.7159956650215 +-2.4750000000000 1.8070924308486 +-2.1593750000000 1.9247061878960 +-1.8437500000000 2.0823569212870 +-1.5281250000000 2.3046011441221 +-1.2125000000000 2.6411053252770 +-0.8968750000000 3.2093335816411 +-0.5812500000000 4.3658652699413 +-0.2656250000000 7.7945223276589 +0.0500000000000 8.6980035891950 +0.3656250000000 6.1769420829572 +0.6812500000000 3.8889678377421 +0.9968750000000 2.9912227568063 +1.3125000000000 2.5171992397353 +1.6281250000000 2.2249395006047 +1.9437500000000 2.0268997428902 +2.2593750000000 1.8839015960002 +2.5750000000000 1.7758213205098 +2.8906250000000 1.6912713212121 +3.2062500000000 1.6233268610343 +3.5218750000000 1.5675361229725 +3.8375000000000 1.5209073234332 +4.1531250000000 1.4813557674072 +4.4687500000000 1.4473844135480 +4.7843750000000 1.4178905035399 +5.1000000000000 1.3920438112308 diff --git a/test/schwarzschild_ef/gxx_y_[16][16].xg b/test/schwarzschild_ef/gxx_y_[16][16].xg new file mode 100644 index 0000000..7d06172 --- /dev/null +++ b/test/schwarzschild_ef/gxx_y_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label y +"y-label ADMBASE::gxx (x = 0.0500000000000, z = 0.0500000000000), (xi = 16, zi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 1.0000399880030 +-4.6843750000000 1.0000486257845 +-4.3687500000000 1.0000599415250 +-4.0531250000000 1.0000750588260 +-3.7375000000000 1.0000957179230 +-3.4218750000000 1.0001247094795 +-3.1062500000000 1.0001666952776 +-2.7906250000000 1.0002298515519 +-2.4750000000000 1.0003293918706 +-2.1593750000000 1.0004957784375 +-1.8437500000000 1.0007959878536 +-1.5281250000000 1.0013966899306 +-1.2125000000000 1.0027906988207 +-0.8968750000000 1.0068665322743 +-0.5812500000000 1.0249063911754 +-0.2656250000000 1.2407470887032 +0.0500000000000 8.6980035891950 +0.3656250000000 1.0968147544187 +0.6812500000000 1.0155621531534 +0.9968750000000 1.0050093161992 +1.3125000000000 1.0022018310962 +1.6281250000000 1.0011552584619 +1.9437500000000 1.0006794965266 +2.2593750000000 1.0004328797686 +2.5750000000000 1.0002925144012 +2.8906250000000 1.0002068259556 +3.2062500000000 1.0001515866956 +3.5218750000000 1.0001143893632 +3.8375000000000 1.0000884308287 +4.1531250000000 1.0000697679855 +4.4687500000000 1.0000560078292 +4.7843750000000 1.0000456406779 +5.1000000000000 1.0000376820272 diff --git a/test/schwarzschild_ef/gxx_z_[16][16].xg b/test/schwarzschild_ef/gxx_z_[16][16].xg new file mode 100644 index 0000000..643a4f9 --- /dev/null +++ b/test/schwarzschild_ef/gxx_z_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label z +"y-label ADMBASE::gxx (x = 0.0500000000000, y = 0.0500000000000), (xi = 16, yi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 1.0000399880030 +-4.6843750000000 1.0000486257845 +-4.3687500000000 1.0000599415250 +-4.0531250000000 1.0000750588260 +-3.7375000000000 1.0000957179230 +-3.4218750000000 1.0001247094795 +-3.1062500000000 1.0001666952776 +-2.7906250000000 1.0002298515519 +-2.4750000000000 1.0003293918706 +-2.1593750000000 1.0004957784375 +-1.8437500000000 1.0007959878536 +-1.5281250000000 1.0013966899306 +-1.2125000000000 1.0027906988207 +-0.8968750000000 1.0068665322743 +-0.5812500000000 1.0249063911754 +-0.2656250000000 1.2407470887032 +0.0500000000000 8.6980035891950 +0.3656250000000 1.0968147544187 +0.6812500000000 1.0155621531534 +0.9968750000000 1.0050093161992 +1.3125000000000 1.0022018310962 +1.6281250000000 1.0011552584619 +1.9437500000000 1.0006794965266 +2.2593750000000 1.0004328797686 +2.5750000000000 1.0002925144012 +2.8906250000000 1.0002068259556 +3.2062500000000 1.0001515866956 +3.5218750000000 1.0001143893632 +3.8375000000000 1.0000884308287 +4.1531250000000 1.0000697679855 +4.4687500000000 1.0000560078292 +4.7843750000000 1.0000456406779 +5.1000000000000 1.0000376820272 diff --git a/test/schwarzschild_ef/gxy_3D_diagonal.xg b/test/schwarzschild_ef/gxy_3D_diagonal.xg new file mode 100644 index 0000000..e43b863 --- /dev/null +++ b/test/schwarzschild_ef/gxy_3D_diagonal.xg @@ -0,0 +1,40 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label diagonal +"y-label ADMBASE::gxy + + +"Time = 0.0000000000000 +-8.6602540378444 0.0769800358920 +-8.1135755017055 0.0821668161622 +-7.5668969655665 0.0881030453699 +-7.0202184294276 0.0949638068058 +-6.4735398932887 0.1029833256080 +-5.9268613571498 0.1124822442257 +-5.3801828210108 0.1239115265866 +-4.8335042848719 0.1379261561334 +-4.2868257487330 0.1555152240241 +-3.7401472125940 0.1782461033678 +-3.1934686764551 0.2087594193680 +-2.6467901403162 0.2518774180514 +-2.1001116041773 0.3174434469771 +-1.5534330680383 0.4291569945196 +-1.0067545318994 0.6621938571351 +-0.4600759957605 1.4490359697308 +0.0866025403784 7.6980035891950 +0.6332810765174 1.0527184395480 +1.1799596126563 0.5649910891152 +1.7266381487952 0.3861067630944 +2.2733166849342 0.2932572795884 +2.8199952210731 0.2364070200137 +3.3666737572120 0.1980193849311 +3.9133522933509 0.1703569259020 +4.4600308294899 0.1494757978484 +5.0067093656288 0.1331546566780 +5.5533879017677 0.1200468395976 +6.1000664379066 0.1092884271758 +6.6467449740456 0.1002997210319 +7.1934235101845 0.0926772441137 +7.7401020463234 0.0861315086903 +8.2867805824623 0.0804494170001 +8.8334591186013 0.0754706234235 diff --git a/test/schwarzschild_ef/gxy_maximum.xg b/test/schwarzschild_ef/gxy_maximum.xg new file mode 100644 index 0000000..f889a63 --- /dev/null +++ b/test/schwarzschild_ef/gxy_maximum.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gxy +"gxy v time +0.0000000000000 7.6980035891950 diff --git a/test/schwarzschild_ef/gxy_minimum.xg b/test/schwarzschild_ef/gxy_minimum.xg new file mode 100644 index 0000000..1921fc1 --- /dev/null +++ b/test/schwarzschild_ef/gxy_minimum.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gxy +"gxy v time +0.0000000000000 -2.0663480485652 diff --git a/test/schwarzschild_ef/gxy_norm1.xg b/test/schwarzschild_ef/gxy_norm1.xg new file mode 100644 index 0000000..c44eda4 --- /dev/null +++ b/test/schwarzschild_ef/gxy_norm1.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gxy +"gxy v time +0.0000000000000 0.1040489502246 diff --git a/test/schwarzschild_ef/gxy_norm2.xg b/test/schwarzschild_ef/gxy_norm2.xg new file mode 100644 index 0000000..17a3787 --- /dev/null +++ b/test/schwarzschild_ef/gxy_norm2.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gxy +"gxy v time +0.0000000000000 0.1435988504864 diff --git a/test/schwarzschild_ef/gxy_x_[16][16].xg b/test/schwarzschild_ef/gxy_x_[16][16].xg new file mode 100644 index 0000000..e01e50b --- /dev/null +++ b/test/schwarzschild_ef/gxy_x_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label x +"y-label ADMBASE::gxy (y = 0.0500000000000, z = 0.0500000000000), (yi = 16, zi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 -0.0039988002999 +-4.6843750000000 -0.0045556281810 +-4.3687500000000 -0.0052373907444 +-4.0531250000000 -0.0060844560826 +-3.7375000000000 -0.0071549147477 +-3.4218750000000 -0.0085348050041 +-3.1062500000000 -0.0103559441201 +-2.7906250000000 -0.0128285897428 +-2.4750000000000 -0.0163048975929 +-2.1593750000000 -0.0214114312682 +-1.8437500000000 -0.0293520521027 +-1.5281250000000 -0.0426863360040 +-1.2125000000000 -0.0676744464032 +-0.8968750000000 -0.1231684226699 +-0.5812500000000 -0.2895367974143 +-0.2656250000000 -1.2789689087358 +0.0500000000000 7.6980035891950 +0.3656250000000 0.7079578916864 +0.6812500000000 0.2120343367150 +0.9968750000000 0.0998732417207 +1.3125000000000 0.0577980662756 +1.6281250000000 0.0376181036654 +1.9437500000000 0.0264154274698 +2.2593750000000 0.0195607545450 +2.5750000000000 0.0150644916604 +2.8906250000000 0.0119571255561 +3.2062500000000 0.0097204968582 +3.5218750000000 0.0080573007698 +3.8375000000000 0.0067870661034 +4.1531250000000 0.0057951032946 +4.4687500000000 0.0050056997320 +4.7843750000000 0.0043672423623 +5.1000000000000 0.0038435667768 diff --git a/test/schwarzschild_ef/gxy_y_[16][16].xg b/test/schwarzschild_ef/gxy_y_[16][16].xg new file mode 100644 index 0000000..d5103c1 --- /dev/null +++ b/test/schwarzschild_ef/gxy_y_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label y +"y-label ADMBASE::gxy (x = 0.0500000000000, z = 0.0500000000000), (xi = 16, zi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 -0.0039988002999 +-4.6843750000000 -0.0045556281810 +-4.3687500000000 -0.0052373907444 +-4.0531250000000 -0.0060844560826 +-3.7375000000000 -0.0071549147477 +-3.4218750000000 -0.0085348050041 +-3.1062500000000 -0.0103559441201 +-2.7906250000000 -0.0128285897428 +-2.4750000000000 -0.0163048975929 +-2.1593750000000 -0.0214114312682 +-1.8437500000000 -0.0293520521027 +-1.5281250000000 -0.0426863360040 +-1.2125000000000 -0.0676744464032 +-0.8968750000000 -0.1231684226699 +-0.5812500000000 -0.2895367974143 +-0.2656250000000 -1.2789689087358 +0.0500000000000 7.6980035891950 +0.3656250000000 0.7079578916864 +0.6812500000000 0.2120343367150 +0.9968750000000 0.0998732417207 +1.3125000000000 0.0577980662756 +1.6281250000000 0.0376181036654 +1.9437500000000 0.0264154274698 +2.2593750000000 0.0195607545450 +2.5750000000000 0.0150644916604 +2.8906250000000 0.0119571255561 +3.2062500000000 0.0097204968582 +3.5218750000000 0.0080573007698 +3.8375000000000 0.0067870661034 +4.1531250000000 0.0057951032946 +4.4687500000000 0.0050056997320 +4.7843750000000 0.0043672423623 +5.1000000000000 0.0038435667768 diff --git a/test/schwarzschild_ef/gxy_z_[16][16].xg b/test/schwarzschild_ef/gxy_z_[16][16].xg new file mode 100644 index 0000000..37c8113 --- /dev/null +++ b/test/schwarzschild_ef/gxy_z_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label z +"y-label ADMBASE::gxy (x = 0.0500000000000, y = 0.0500000000000), (xi = 16, yi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 0.0000399880030 +-4.6843750000000 0.0000486257845 +-4.3687500000000 0.0000599415250 +-4.0531250000000 0.0000750588260 +-3.7375000000000 0.0000957179230 +-3.4218750000000 0.0001247094795 +-3.1062500000000 0.0001666952776 +-2.7906250000000 0.0002298515519 +-2.4750000000000 0.0003293918706 +-2.1593750000000 0.0004957784375 +-1.8437500000000 0.0007959878536 +-1.5281250000000 0.0013966899306 +-1.2125000000000 0.0027906988207 +-0.8968750000000 0.0068665322743 +-0.5812500000000 0.0249063911754 +-0.2656250000000 0.2407470887032 +0.0500000000000 7.6980035891950 +0.3656250000000 0.0968147544187 +0.6812500000000 0.0155621531534 +0.9968750000000 0.0050093161992 +1.3125000000000 0.0022018310962 +1.6281250000000 0.0011552584619 +1.9437500000000 0.0006794965266 +2.2593750000000 0.0004328797686 +2.5750000000000 0.0002925144012 +2.8906250000000 0.0002068259556 +3.2062500000000 0.0001515866956 +3.5218750000000 0.0001143893632 +3.8375000000000 0.0000884308287 +4.1531250000000 0.0000697679855 +4.4687500000000 0.0000560078292 +4.7843750000000 0.0000456406779 +5.1000000000000 0.0000376820272 diff --git a/test/schwarzschild_ef/gyy_3D_diagonal.xg b/test/schwarzschild_ef/gyy_3D_diagonal.xg new file mode 100644 index 0000000..553cd77 --- /dev/null +++ b/test/schwarzschild_ef/gyy_3D_diagonal.xg @@ -0,0 +1,40 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label diagonal +"y-label ADMBASE::gyy + + +"Time = 0.0000000000000 +-8.6602540378444 1.0769800358920 +-8.1135755017055 1.0821668161622 +-7.5668969655665 1.0881030453699 +-7.0202184294276 1.0949638068058 +-6.4735398932887 1.1029833256080 +-5.9268613571498 1.1124822442257 +-5.3801828210108 1.1239115265866 +-4.8335042848719 1.1379261561334 +-4.2868257487330 1.1555152240241 +-3.7401472125940 1.1782461033678 +-3.1934686764551 1.2087594193680 +-2.6467901403162 1.2518774180514 +-2.1001116041773 1.3174434469771 +-1.5534330680383 1.4291569945196 +-1.0067545318994 1.6621938571351 +-0.4600759957605 2.4490359697308 +0.0866025403784 8.6980035891950 +0.6332810765174 2.0527184395480 +1.1799596126563 1.5649910891152 +1.7266381487952 1.3861067630944 +2.2733166849342 1.2932572795884 +2.8199952210731 1.2364070200137 +3.3666737572120 1.1980193849311 +3.9133522933509 1.1703569259020 +4.4600308294899 1.1494757978484 +5.0067093656288 1.1331546566780 +5.5533879017677 1.1200468395976 +6.1000664379066 1.1092884271758 +6.6467449740456 1.1002997210319 +7.1934235101845 1.0926772441137 +7.7401020463234 1.0861315086903 +8.2867805824623 1.0804494170001 +8.8334591186013 1.0754706234235 diff --git a/test/schwarzschild_ef/gyy_maximum.xg b/test/schwarzschild_ef/gyy_maximum.xg new file mode 100644 index 0000000..14825a7 --- /dev/null +++ b/test/schwarzschild_ef/gyy_maximum.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gyy +"gyy v time +0.0000000000000 8.6980035891950 diff --git a/test/schwarzschild_ef/gyy_minimum.xg b/test/schwarzschild_ef/gyy_minimum.xg new file mode 100644 index 0000000..5c26dcd --- /dev/null +++ b/test/schwarzschild_ef/gyy_minimum.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gyy +"gyy v time +0.0000000000000 1.0000133254896 diff --git a/test/schwarzschild_ef/gyy_norm1.xg b/test/schwarzschild_ef/gyy_norm1.xg new file mode 100644 index 0000000..b67f248 --- /dev/null +++ b/test/schwarzschild_ef/gyy_norm1.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gyy +"gyy v time +0.0000000000000 1.1524177407289 diff --git a/test/schwarzschild_ef/gyy_norm2.xg b/test/schwarzschild_ef/gyy_norm2.xg new file mode 100644 index 0000000..efffa68 --- /dev/null +++ b/test/schwarzschild_ef/gyy_norm2.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gyy +"gyy v time +0.0000000000000 1.1660598040677 diff --git a/test/schwarzschild_ef/gyy_x_[16][16].xg b/test/schwarzschild_ef/gyy_x_[16][16].xg new file mode 100644 index 0000000..18ce660 --- /dev/null +++ b/test/schwarzschild_ef/gyy_x_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label x +"y-label ADMBASE::gyy (y = 0.0500000000000, z = 0.0500000000000), (yi = 16, zi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 1.0000399880030 +-4.6843750000000 1.0000486257845 +-4.3687500000000 1.0000599415250 +-4.0531250000000 1.0000750588260 +-3.7375000000000 1.0000957179230 +-3.4218750000000 1.0001247094795 +-3.1062500000000 1.0001666952776 +-2.7906250000000 1.0002298515519 +-2.4750000000000 1.0003293918706 +-2.1593750000000 1.0004957784375 +-1.8437500000000 1.0007959878536 +-1.5281250000000 1.0013966899306 +-1.2125000000000 1.0027906988207 +-0.8968750000000 1.0068665322743 +-0.5812500000000 1.0249063911754 +-0.2656250000000 1.2407470887032 +0.0500000000000 8.6980035891950 +0.3656250000000 1.0968147544187 +0.6812500000000 1.0155621531534 +0.9968750000000 1.0050093161992 +1.3125000000000 1.0022018310962 +1.6281250000000 1.0011552584619 +1.9437500000000 1.0006794965266 +2.2593750000000 1.0004328797686 +2.5750000000000 1.0002925144012 +2.8906250000000 1.0002068259556 +3.2062500000000 1.0001515866956 +3.5218750000000 1.0001143893632 +3.8375000000000 1.0000884308287 +4.1531250000000 1.0000697679855 +4.4687500000000 1.0000560078292 +4.7843750000000 1.0000456406779 +5.1000000000000 1.0000376820272 diff --git a/test/schwarzschild_ef/gyy_y_[16][16].xg b/test/schwarzschild_ef/gyy_y_[16][16].xg new file mode 100644 index 0000000..2b5e9ff --- /dev/null +++ b/test/schwarzschild_ef/gyy_y_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label y +"y-label ADMBASE::gyy (x = 0.0500000000000, z = 0.0500000000000), (xi = 16, zi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 1.3998800299930 +-4.6843750000000 1.4268054152043 +-4.3687500000000 1.4576170162950 +-4.0531250000000 1.4932212211925 +-3.7375000000000 1.5348298773928 +-3.4218750000000 1.5841007174696 +-3.1062500000000 1.6433630284612 +-2.7906250000000 1.7159956650215 +-2.4750000000000 1.8070924308486 +-2.1593750000000 1.9247061878960 +-1.8437500000000 2.0823569212870 +-1.5281250000000 2.3046011441221 +-1.2125000000000 2.6411053252770 +-0.8968750000000 3.2093335816411 +-0.5812500000000 4.3658652699413 +-0.2656250000000 7.7945223276589 +0.0500000000000 8.6980035891950 +0.3656250000000 6.1769420829572 +0.6812500000000 3.8889678377421 +0.9968750000000 2.9912227568063 +1.3125000000000 2.5171992397353 +1.6281250000000 2.2249395006047 +1.9437500000000 2.0268997428902 +2.2593750000000 1.8839015960002 +2.5750000000000 1.7758213205098 +2.8906250000000 1.6912713212121 +3.2062500000000 1.6233268610343 +3.5218750000000 1.5675361229725 +3.8375000000000 1.5209073234332 +4.1531250000000 1.4813557674072 +4.4687500000000 1.4473844135480 +4.7843750000000 1.4178905035399 +5.1000000000000 1.3920438112308 diff --git a/test/schwarzschild_ef/gyy_z_[16][16].xg b/test/schwarzschild_ef/gyy_z_[16][16].xg new file mode 100644 index 0000000..790fd2c --- /dev/null +++ b/test/schwarzschild_ef/gyy_z_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label z +"y-label ADMBASE::gyy (x = 0.0500000000000, y = 0.0500000000000), (xi = 16, yi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 1.0000399880030 +-4.6843750000000 1.0000486257845 +-4.3687500000000 1.0000599415250 +-4.0531250000000 1.0000750588260 +-3.7375000000000 1.0000957179230 +-3.4218750000000 1.0001247094795 +-3.1062500000000 1.0001666952776 +-2.7906250000000 1.0002298515519 +-2.4750000000000 1.0003293918706 +-2.1593750000000 1.0004957784375 +-1.8437500000000 1.0007959878536 +-1.5281250000000 1.0013966899306 +-1.2125000000000 1.0027906988207 +-0.8968750000000 1.0068665322743 +-0.5812500000000 1.0249063911754 +-0.2656250000000 1.2407470887032 +0.0500000000000 8.6980035891950 +0.3656250000000 1.0968147544187 +0.6812500000000 1.0155621531534 +0.9968750000000 1.0050093161992 +1.3125000000000 1.0022018310962 +1.6281250000000 1.0011552584619 +1.9437500000000 1.0006794965266 +2.2593750000000 1.0004328797686 +2.5750000000000 1.0002925144012 +2.8906250000000 1.0002068259556 +3.2062500000000 1.0001515866956 +3.5218750000000 1.0001143893632 +3.8375000000000 1.0000884308287 +4.1531250000000 1.0000697679855 +4.4687500000000 1.0000560078292 +4.7843750000000 1.0000456406779 +5.1000000000000 1.0000376820272 diff --git a/test/schwarzschild_ef/gzz_3D_diagonal.xg b/test/schwarzschild_ef/gzz_3D_diagonal.xg new file mode 100644 index 0000000..392c146 --- /dev/null +++ b/test/schwarzschild_ef/gzz_3D_diagonal.xg @@ -0,0 +1,40 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label diagonal +"y-label ADMBASE::gzz + + +"Time = 0.0000000000000 +-8.6602540378444 1.0769800358920 +-8.1135755017055 1.0821668161622 +-7.5668969655665 1.0881030453699 +-7.0202184294276 1.0949638068058 +-6.4735398932887 1.1029833256080 +-5.9268613571498 1.1124822442257 +-5.3801828210108 1.1239115265866 +-4.8335042848719 1.1379261561334 +-4.2868257487330 1.1555152240241 +-3.7401472125940 1.1782461033678 +-3.1934686764551 1.2087594193680 +-2.6467901403162 1.2518774180514 +-2.1001116041773 1.3174434469771 +-1.5534330680383 1.4291569945196 +-1.0067545318994 1.6621938571351 +-0.4600759957605 2.4490359697308 +0.0866025403784 8.6980035891950 +0.6332810765174 2.0527184395480 +1.1799596126563 1.5649910891152 +1.7266381487952 1.3861067630944 +2.2733166849342 1.2932572795884 +2.8199952210731 1.2364070200137 +3.3666737572120 1.1980193849311 +3.9133522933509 1.1703569259020 +4.4600308294899 1.1494757978484 +5.0067093656288 1.1331546566780 +5.5533879017677 1.1200468395976 +6.1000664379066 1.1092884271758 +6.6467449740456 1.1002997210319 +7.1934235101845 1.0926772441137 +7.7401020463234 1.0861315086903 +8.2867805824623 1.0804494170001 +8.8334591186013 1.0754706234235 diff --git a/test/schwarzschild_ef/gzz_maximum.xg b/test/schwarzschild_ef/gzz_maximum.xg new file mode 100644 index 0000000..fc4407e --- /dev/null +++ b/test/schwarzschild_ef/gzz_maximum.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gzz +"gzz v time +0.0000000000000 8.6980035891950 diff --git a/test/schwarzschild_ef/gzz_minimum.xg b/test/schwarzschild_ef/gzz_minimum.xg new file mode 100644 index 0000000..3c42326 --- /dev/null +++ b/test/schwarzschild_ef/gzz_minimum.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gzz +"gzz v time +0.0000000000000 1.0000133254896 diff --git a/test/schwarzschild_ef/gzz_norm1.xg b/test/schwarzschild_ef/gzz_norm1.xg new file mode 100644 index 0000000..856c9c4 --- /dev/null +++ b/test/schwarzschild_ef/gzz_norm1.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gzz +"gzz v time +0.0000000000000 1.1524177407289 diff --git a/test/schwarzschild_ef/gzz_norm2.xg b/test/schwarzschild_ef/gzz_norm2.xg new file mode 100644 index 0000000..41f8f56 --- /dev/null +++ b/test/schwarzschild_ef/gzz_norm2.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:38 +"x-label time +"y-label ADMBASE::gzz +"gzz v time +0.0000000000000 1.1660598040678 diff --git a/test/schwarzschild_ef/gzz_x_[16][16].xg b/test/schwarzschild_ef/gzz_x_[16][16].xg new file mode 100644 index 0000000..c3f9448 --- /dev/null +++ b/test/schwarzschild_ef/gzz_x_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label x +"y-label ADMBASE::gzz (y = 0.0500000000000, z = 0.0500000000000), (yi = 16, zi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 1.0000399880030 +-4.6843750000000 1.0000486257845 +-4.3687500000000 1.0000599415250 +-4.0531250000000 1.0000750588260 +-3.7375000000000 1.0000957179230 +-3.4218750000000 1.0001247094795 +-3.1062500000000 1.0001666952776 +-2.7906250000000 1.0002298515519 +-2.4750000000000 1.0003293918706 +-2.1593750000000 1.0004957784375 +-1.8437500000000 1.0007959878536 +-1.5281250000000 1.0013966899306 +-1.2125000000000 1.0027906988207 +-0.8968750000000 1.0068665322743 +-0.5812500000000 1.0249063911754 +-0.2656250000000 1.2407470887032 +0.0500000000000 8.6980035891950 +0.3656250000000 1.0968147544187 +0.6812500000000 1.0155621531534 +0.9968750000000 1.0050093161992 +1.3125000000000 1.0022018310962 +1.6281250000000 1.0011552584619 +1.9437500000000 1.0006794965266 +2.2593750000000 1.0004328797686 +2.5750000000000 1.0002925144012 +2.8906250000000 1.0002068259556 +3.2062500000000 1.0001515866956 +3.5218750000000 1.0001143893632 +3.8375000000000 1.0000884308287 +4.1531250000000 1.0000697679855 +4.4687500000000 1.0000560078292 +4.7843750000000 1.0000456406779 +5.1000000000000 1.0000376820272 diff --git a/test/schwarzschild_ef/gzz_y_[16][16].xg b/test/schwarzschild_ef/gzz_y_[16][16].xg new file mode 100644 index 0000000..d7717c1 --- /dev/null +++ b/test/schwarzschild_ef/gzz_y_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label y +"y-label ADMBASE::gzz (x = 0.0500000000000, z = 0.0500000000000), (xi = 16, zi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 1.0000399880030 +-4.6843750000000 1.0000486257845 +-4.3687500000000 1.0000599415250 +-4.0531250000000 1.0000750588260 +-3.7375000000000 1.0000957179230 +-3.4218750000000 1.0001247094795 +-3.1062500000000 1.0001666952776 +-2.7906250000000 1.0002298515519 +-2.4750000000000 1.0003293918706 +-2.1593750000000 1.0004957784375 +-1.8437500000000 1.0007959878536 +-1.5281250000000 1.0013966899306 +-1.2125000000000 1.0027906988207 +-0.8968750000000 1.0068665322743 +-0.5812500000000 1.0249063911754 +-0.2656250000000 1.2407470887032 +0.0500000000000 8.6980035891950 +0.3656250000000 1.0968147544187 +0.6812500000000 1.0155621531534 +0.9968750000000 1.0050093161992 +1.3125000000000 1.0022018310962 +1.6281250000000 1.0011552584619 +1.9437500000000 1.0006794965266 +2.2593750000000 1.0004328797686 +2.5750000000000 1.0002925144012 +2.8906250000000 1.0002068259556 +3.2062500000000 1.0001515866956 +3.5218750000000 1.0001143893632 +3.8375000000000 1.0000884308287 +4.1531250000000 1.0000697679855 +4.4687500000000 1.0000560078292 +4.7843750000000 1.0000456406779 +5.1000000000000 1.0000376820272 diff --git a/test/schwarzschild_ef/gzz_z_[16][16].xg b/test/schwarzschild_ef/gzz_z_[16][16].xg new file mode 100644 index 0000000..a0160bc --- /dev/null +++ b/test/schwarzschild_ef/gzz_z_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label z +"y-label ADMBASE::gzz (x = 0.0500000000000, y = 0.0500000000000), (xi = 16, yi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 1.3998800299930 +-4.6843750000000 1.4268054152043 +-4.3687500000000 1.4576170162950 +-4.0531250000000 1.4932212211925 +-3.7375000000000 1.5348298773928 +-3.4218750000000 1.5841007174696 +-3.1062500000000 1.6433630284612 +-2.7906250000000 1.7159956650215 +-2.4750000000000 1.8070924308486 +-2.1593750000000 1.9247061878960 +-1.8437500000000 2.0823569212870 +-1.5281250000000 2.3046011441221 +-1.2125000000000 2.6411053252770 +-0.8968750000000 3.2093335816411 +-0.5812500000000 4.3658652699413 +-0.2656250000000 7.7945223276589 +0.0500000000000 8.6980035891950 +0.3656250000000 6.1769420829572 +0.6812500000000 3.8889678377421 +0.9968750000000 2.9912227568063 +1.3125000000000 2.5171992397353 +1.6281250000000 2.2249395006047 +1.9437500000000 2.0268997428902 +2.2593750000000 1.8839015960002 +2.5750000000000 1.7758213205098 +2.8906250000000 1.6912713212121 +3.2062500000000 1.6233268610343 +3.5218750000000 1.5675361229725 +3.8375000000000 1.5209073234332 +4.1531250000000 1.4813557674072 +4.4687500000000 1.4473844135480 +4.7843750000000 1.4178905035399 +5.1000000000000 1.3920438112308 diff --git a/test/schwarzschild_ef/ham_3D_diagonal.xg b/test/schwarzschild_ef/ham_3D_diagonal.xg new file mode 100644 index 0000000..499f248 --- /dev/null +++ b/test/schwarzschild_ef/ham_3D_diagonal.xg @@ -0,0 +1,40 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:35 +"x-label diagonal +"y-label ADMCONSTRAINTS::ham + + +"Time = 0.0000000000000 +-8.6602540378444 -0.0000528143163 +-8.1135755017055 -0.0000528143163 +-7.5668969655665 -0.0000738064947 +-7.0202184294276 -0.0001056436990 +-6.4735398932887 -0.0001554524304 +-5.9268613571498 -0.0002362390033 +-5.3801828210108 -0.0003729387474 +-4.8335042848719 -0.0006162125391 +-4.2868257487330 -0.0010763685394 +-3.7401472125940 -0.0020146182917 +-3.1934686764551 -0.0041167645698 +-2.6467901403162 -0.0094289875687 +-2.1001116041773 -0.0250666506171 +-1.5534330680383 -0.0788972982249 +-1.0067545318994 -0.1456530657695 +-0.4600759957605 23.1353161983048 +0.0866025403784 -143.1316515243370 +0.6332810765174 4.6775934471488 +1.1799596126563 -0.1632285923450 +1.7266381487952 -0.0540460590569 +2.2733166849342 -0.0180535878881 +2.8199952210731 -0.0071484962168 +3.3666737572120 -0.0032469897201 +3.9133522933509 -0.0016377810215 +4.4600308294899 -0.0008959497205 +5.0067093656288 -0.0005226976644 +5.5533879017677 -0.0003212505513 +6.1000664379066 -0.0002061151283 +6.6467449740456 -0.0001370997391 +7.1934235101845 -0.0000940331339 +7.7401020463234 -0.0000662195428 +8.2867805824623 -0.0000477152109 +8.8334591186013 -0.0000477152109 diff --git a/test/schwarzschild_ef/ham_maximum.xg b/test/schwarzschild_ef/ham_maximum.xg new file mode 100644 index 0000000..9981410 --- /dev/null +++ b/test/schwarzschild_ef/ham_maximum.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label time +"y-label ADMCONSTRAINTS::ham +"ham v time +0.0000000000000 44.3957814069886 diff --git a/test/schwarzschild_ef/ham_minimum.xg b/test/schwarzschild_ef/ham_minimum.xg new file mode 100644 index 0000000..efa3b16 --- /dev/null +++ b/test/schwarzschild_ef/ham_minimum.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label time +"y-label ADMCONSTRAINTS::ham +"ham v time +0.0000000000000 -143.1316515243370 diff --git a/test/schwarzschild_ef/ham_norm1.xg b/test/schwarzschild_ef/ham_norm1.xg new file mode 100644 index 0000000..1c8e050 --- /dev/null +++ b/test/schwarzschild_ef/ham_norm1.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label time +"y-label ADMCONSTRAINTS::ham +"ham v time +0.0000000000000 0.0189288891482 diff --git a/test/schwarzschild_ef/ham_norm2.xg b/test/schwarzschild_ef/ham_norm2.xg new file mode 100644 index 0000000..4daefa2 --- /dev/null +++ b/test/schwarzschild_ef/ham_norm2.xg @@ -0,0 +1,6 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:36 +"x-label time +"y-label ADMCONSTRAINTS::ham +"ham v time +0.0000000000000 0.9460364102600 diff --git a/test/schwarzschild_ef/ham_x_[16][16].xg b/test/schwarzschild_ef/ham_x_[16][16].xg new file mode 100644 index 0000000..b92c30d --- /dev/null +++ b/test/schwarzschild_ef/ham_x_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:35 +"x-label x +"y-label ADMCONSTRAINTS::ham (y = 0.0500000000000, z = 0.0500000000000), (yi = 16, zi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 -0.0003365562508 +-4.6843750000000 -0.0003365562508 +-4.3687500000000 -0.0004635290728 +-4.0531250000000 -0.0006520686426 +-3.7375000000000 -0.0009395404067 +-3.4218750000000 -0.0013911020344 +-3.1062500000000 -0.0021244945253 +-2.7906250000000 -0.0033604935452 +-2.4750000000000 -0.0055274430879 +-2.1593750000000 -0.0094712284664 +-1.8437500000000 -0.0167785572162 +-1.5281250000000 -0.0292520976357 +-1.2125000000000 -0.0352711711531 +-0.8968750000000 0.1726939295639 +-0.5812500000000 3.3613260068408 +-0.2656250000000 44.3957814069885 +0.0500000000000 -143.1316515243370 +0.3656250000000 7.5689433864156 +0.6812500000000 1.4688650877804 +0.9968750000000 0.0367687147962 +1.3125000000000 -0.0372626774353 +1.6281250000000 -0.0247622584878 +1.9437500000000 -0.0139674110549 +2.2593750000000 -0.0079512484929 +2.5750000000000 -0.0047003412163 +2.8906250000000 -0.0028943111817 +3.2062500000000 -0.0018510929707 +3.5218750000000 -0.0012245677433 +3.8375000000000 -0.0008345477230 +4.1531250000000 -0.0005838064260 +4.4687500000000 -0.0004179151147 +4.7843750000000 -0.0003053235832 +5.1000000000000 -0.0003053235832 diff --git a/test/schwarzschild_ef/ham_y_[16][16].xg b/test/schwarzschild_ef/ham_y_[16][16].xg new file mode 100644 index 0000000..3d244cf --- /dev/null +++ b/test/schwarzschild_ef/ham_y_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:35 +"x-label y +"y-label ADMCONSTRAINTS::ham (x = 0.0500000000000, z = 0.0500000000000), (xi = 16, zi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 -0.0003365562508 +-4.6843750000000 -0.0003365562508 +-4.3687500000000 -0.0004635290728 +-4.0531250000000 -0.0006520686426 +-3.7375000000000 -0.0009395404067 +-3.4218750000000 -0.0013911020343 +-3.1062500000000 -0.0021244945253 +-2.7906250000000 -0.0033604935452 +-2.4750000000000 -0.0055274430879 +-2.1593750000000 -0.0094712284664 +-1.8437500000000 -0.0167785572162 +-1.5281250000000 -0.0292520976357 +-1.2125000000000 -0.0352711711531 +-0.8968750000000 0.1726939295639 +-0.5812500000000 3.3613260068408 +-0.2656250000000 44.3957814069886 +0.0500000000000 -143.1316515243370 +0.3656250000000 7.5689433864156 +0.6812500000000 1.4688650877804 +0.9968750000000 0.0367687147962 +1.3125000000000 -0.0372626774353 +1.6281250000000 -0.0247622584878 +1.9437500000000 -0.0139674110549 +2.2593750000000 -0.0079512484929 +2.5750000000000 -0.0047003412163 +2.8906250000000 -0.0028943111817 +3.2062500000000 -0.0018510929707 +3.5218750000000 -0.0012245677433 +3.8375000000000 -0.0008345477230 +4.1531250000000 -0.0005838064260 +4.4687500000000 -0.0004179151147 +4.7843750000000 -0.0003053235832 +5.1000000000000 -0.0003053235832 diff --git a/test/schwarzschild_ef/ham_z_[16][16].xg b/test/schwarzschild_ef/ham_z_[16][16].xg new file mode 100644 index 0000000..aab285c --- /dev/null +++ b/test/schwarzschild_ef/ham_z_[16][16].xg @@ -0,0 +1,41 @@ +"Parameter file arrangements/AEIThorns/Exact_old/test/schwarzschild_ef.par +"Created Jan 28 2003 15:34:35 +"x-label z +"y-label ADMCONSTRAINTS::ham (x = 0.0500000000000, y = 0.0500000000000), (xi = 16, yi = 16) + + + +"Time = 0.0000000000000 +-5.0000000000000 -0.0003365562503 +-4.6843750000000 -0.0003365562503 +-4.3687500000000 -0.0004635290728 +-4.0531250000000 -0.0006520686744 +-3.7375000000000 -0.0009395404067 +-3.4218750000000 -0.0013911020066 +-3.1062500000000 -0.0021244944914 +-2.7906250000000 -0.0033604935452 +-2.4750000000000 -0.0055274431962 +-2.1593750000000 -0.0094712284664 +-1.8437500000000 -0.0167785572162 +-1.5281250000000 -0.0292520976357 +-1.2125000000000 -0.0352711711531 +-0.8968750000000 0.1726939295639 +-0.5812500000000 3.3613260112518 +-0.2656250000000 44.3957813759909 +0.0500000000000 -143.1316515243370 +0.3656250000000 7.5689433708177 +0.6812500000000 1.4688650901672 +0.9968750000000 0.0367687141579 +1.3125000000000 -0.0372626774353 +1.6281250000000 -0.0247622584878 +1.9437500000000 -0.0139674110549 +2.2593750000000 -0.0079512486725 +2.5750000000000 -0.0047003410830 +2.8906250000000 -0.0028943111552 +3.2062500000000 -0.0018510929173 +3.5218750000000 -0.0012245677516 +3.8375000000000 -0.0008345477448 +4.1531250000000 -0.0005838064188 +4.4687500000000 -0.0004179151147 +4.7843750000000 -0.0003053236164 +5.1000000000000 -0.0003053236164 diff --git a/test/schwarzschild_ef/schwarzschild_ef.par b/test/schwarzschild_ef/schwarzschild_ef.par new file mode 100644 index 0000000..bd8fc16 --- /dev/null +++ b/test/schwarzschild_ef/schwarzschild_ef.par @@ -0,0 +1,54 @@ +# Initial data and short evolution (ADM) for DeSitter spacetime +# Author: Mitica Vulcanov <vulcan@aei.mpg.de> +# $Header$ + +# Author: Frank Loeffler (frank.loeffler@aei.mpg.de) + +# Required thorns +ActiveThorns = "CartGrid3D Boundary Time PUGH PUGHReduce PUGHSlab IOAscii IOUtil IOBasic ADMBase ADMMacros ADMConstraints MoL BSSN_MoL ADMCoupling StaticConformal CoordGauge SpaceMask Exact NaNChecker" + +# Grid +driver::global_nx = 33 +driver::global_ny = 33 +driver::global_nz = 33 +grid::type = "byrange" +grid::domain = "full" + +grid::xmin = -5.0 +grid::xmax = 5.1 +grid::ymin = -5.0 +grid::ymax = 5.1 +grid::zmin = -5.0 +grid::zmax = 5.1 + +cactus::cctk_itlast = 0 +time::dtfac = 0.25 + +# Initial data +ADMBase::initial_data = "exact" +ADMBase::initial_shift = "exact" +Exact::exact_model = "Schwarzschild/EF" + +# Gauge +ADMBase::lapse_evolution_method = "exact" +ADMBase::shift_evolution_method = "exact" + +# Evolution +ADMBase::evolution_method = "adm_bssn" +ADM_BSSN::bound = "radiative" +boundary::radpower = 2 + +methodoflines::ode_method = "icn" + +# Output +IO::out_dir = "schwarzschild_ef" + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = "admconstraints::ham" + +IOBasic::outScalar_every = 1 +IOBasic::outScalar_vars = "admbase::gxx admbase::gyy admbase::gzz admbase::gxy admconstraints::ham admbase::lapse" + +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = "admbase::gxx admbase::gyy admbase::gzz admbase::gxy admconstraints::ham admbase::lapse" + |