Commits December 2013

commits@lists.einsteintoolkit.org

5 participants
79 discussions

[svn:einsteintoolkit] IDLinearWaves/ (Rev. 117)
by knarf＠cct.lsu.edu 02 Dec '13

02 Dec '13

User: knarf Date: 2013/12/02 06:35 PM Modified: /trunk/test/test_pw_ADM_leap/ gxx.dl, gxy.dl, gxz.dl, gyy.dl, gyz.dl, gzz.dl, kxx.dl, kxy.dl, kxz.dl, kyy.dl, kyz.dl, kzz.dl Log: set correct mime type Directory Changes: Directory: /svn:mime-type/ ========================== - video/dl + text/plain File Changes: Directory: /trunk/test/test_pw_ADM_leap/ ======================================== File [modified]: gxx.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/test_pw_ADM_leap/gxx.dl ___________________________________________________________________ File [modified]: gxy.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/test_pw_ADM_leap/gxy.dl ___________________________________________________________________ File [modified]: gxz.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/test_pw_ADM_leap/gxz.dl ___________________________________________________________________ File [modified]: gyy.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/test_pw_ADM_leap/gyy.dl ___________________________________________________________________ File [modified]: gyz.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/test_pw_ADM_leap/gyz.dl ___________________________________________________________________ File [modified]: gzz.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/test_pw_ADM_leap/gzz.dl ___________________________________________________________________ File [modified]: kxx.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/test_pw_ADM_leap/kxx.dl ___________________________________________________________________ File [modified]: kxy.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/test_pw_ADM_leap/kxy.dl ___________________________________________________________________ File [modified]: kxz.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/test_pw_ADM_leap/kxz.dl ___________________________________________________________________ File [modified]: kyy.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/test_pw_ADM_leap/kyy.dl ___________________________________________________________________ File [modified]: kyz.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/test_pw_ADM_leap/kyz.dl ___________________________________________________________________ File [modified]: kzz.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/test_pw_ADM_leap/kzz.dl ___________________________________________________________________

1 0

[svn:einsteintoolkit] IDLinearWaves/trunk/test/ (Rev. 116)
by knarf＠cct.lsu.edu 02 Dec '13

02 Dec '13

User: knarf Date: 2013/12/02 06:32 PM Removed: /trunk/test/ test_pw_ADM_icn_par, test_spw_ADM_icn_par /trunk/test/test_pw_ADM_icn/ Log: remove test that wasn't enabled for the last 14 years File Changes: Directory: /trunk/test/ ======================= File [removed]: test_pw_ADM_icn_par Delta lines: +0 -59 =================================================================== --- trunk/test/test_pw_ADM_icn_par 2013-11-25 15:39:38 UTC (rev 115) +++ trunk/test/test_pw_ADM_icn_par 2013-12-03 00:32:08 UTC (rev 116) @@ -1,59 +0,0 @@ -######################################################### -!DESC "Linearized plane wave along y-line, ADM, ICN" -######################################################### - -# Required thorns -ActiveThorns = "time ADMconstraints pugh cartgrid3d einstein ADM IDLinearWaves - ioascii ioutil iobasic" - -# GENERAL - -driver::global_nx = 3 -driver::global_ny = 200 -driver::global_nz = 3 -time::dtfac = 0.25 - -grid::type = "byspacing" -grid::dxyz = 0.01 -grid::domain = "full" - -cactus::cctk_initial_time = 0. -cactus::cctk_itlast = 100 - - -# Einstein - -einstein::initial_data = "planewaves" -einstein::evolution_system = "ADM" - -adm::method = "icn" -adm::bound = "flat" -einstein::slicing = "geodesic" - -# IDLinearwaves -idlinearwaves::amplitude = 0.001 -idlinearwaves::wavecenter = 0.0 -idlinearwaves::wavelength = 2.0 -idlinearwaves::wavepulse = 1.0 -idlinearwaves::wavetheta = 90.0 -idlinearwaves::wavephi = 90.0 -idlinearwaves::wavesgoing = "both" - -# Output - -IO::outdir = "test_pw_ADM_icn" -IOBasic::outScalar_every = 10 -IOBasic::outScalar_vars = "einstein::gzz einstein::kzz - admconstraints::ADMconstraints - einstein::grr einstein::alp" -IOASCII::out1D_every = 10 -IOASCII::out1D_vars = "einstein::gxx einstein::gxy einstein::gxz - einstein::gyy einstein::gyz einstein::gzz - einstein::kxx einstein::kxy einstein::kxz - einstein::kyy einstein::kyz einstein::kzz" - - -IOBasic::outInfo_every = 5 -IOBasic::outInfo_vars = "admconstraints::ham einstein::gzz" - -################################################################## File [removed]: test_spw_ADM_icn_par Delta lines: +0 -54 =================================================================== --- trunk/test/test_spw_ADM_icn_par 2013-11-25 15:39:38 UTC (rev 115) +++ trunk/test/test_spw_ADM_icn_par 2013-12-03 00:32:08 UTC (rev 116) @@ -1,54 +0,0 @@ -######################################################### -!DESC "Standing linearised plane wave along y-line, ADM, ICN" -######################################################### - -# Required thorns -ActiveThorns = "time ADMconstraints pugh cartgrid3d admbase staticconformal admcoupling admmacros boundary coordgauge pughslab spacemask ADM IDLinearWaves - ioascii ioutil iobasic pughreduce" - -# GENERAL - -driver::global_nx = 3 -driver::global_ny = 3 -driver::global_nz = 200 -time::dtfac = 0.25 - -grid::type = "byspacing" -grid::dxyz = 0.01 -grid::domain = "full" -grid::avoid_origin = no - -cactus::cctk_initial_time = 0. -cactus::cctk_itlast = 100 - - -# Einstein - -admbase::initial_data = "standing_planewaves" -admbase::evolution_method = "ADM" - -adm::method = "icn" -adm::bound = "static" -admbase::lapse_evolution_method = "geodesic" - -# IDLinearwaves -idlinearwaves::amplitude = 0.001 - -# Output - -IO::outdir = "test_spw_ADM_icn" -IOBasic::outScalar_every = 10 -IOBasic::outScalar_vars = "admbase::gzz admbase::kzz - admconstraints::ADMconstraints - admbase::alp" -IOASCII::out1D_every = 10 -IOASCII::out1D_vars = "admbase::gxx admbase::gxy admbase::gxz - admbase::gyy admbase::gyz admbase::gzz - admbase::kxx admbase::kxy admbase::kxz - admbase::kyy admbase::kyz admbase::kzz" - - -IOBasic::outInfo_every = 5 -IOBasic::outInfo_vars = "admconstraints::ham admbase::gzz" - -##################################################################

1 0

[svn:einsteintoolkit] IDConstraintViolate/trunk/test/ (Rev. 29)
by knarf＠cct.lsu.edu 02 Dec '13

02 Dec '13

User: knarf Date: 2013/12/02 05:37 PM Modified: /trunk/test/ ConstraintViolate.par /trunk/test/ConstraintViolate/ alp.dl, alp.xl, alp.yl, alp.zl, gxx.dl, gxx.xl, gxx.yl, gxx.zl, ham.dl, ham.xl, ham.yl, ham.zl, momx.dl, momx.xl, momx.yl, momx.zl, momy.dl, momy.xl, momy.yl, momy.zl, momz.dl, momz.xl, momz.yl, momz.zl Log: don't use evolution (ADM) in ID testsuite File Changes: Directory: /trunk/test/ConstraintViolate/ ========================================= File [modified]: alp.dl Delta lines: <Unavailable> File [modified]: alp.xl Delta lines: <Unavailable> File [modified]: alp.yl Delta lines: <Unavailable> File [modified]: alp.zl Delta lines: <Unavailable> File [modified]: gxx.dl Delta lines: <Unavailable> File [modified]: gxx.xl Delta lines: <Unavailable> File [modified]: gxx.yl Delta lines: <Unavailable> File [modified]: gxx.zl Delta lines: <Unavailable> File [modified]: ham.dl Delta lines: <Unavailable> File [modified]: ham.xl Delta lines: <Unavailable> File [modified]: ham.yl Delta lines: <Unavailable> File [modified]: ham.zl Delta lines: <Unavailable> File [modified]: momx.dl Delta lines: <Unavailable> File [modified]: momx.xl Delta lines: <Unavailable> File [modified]: momx.yl Delta lines: <Unavailable> File [modified]: momx.zl Delta lines: <Unavailable> File [modified]: momy.dl Delta lines: <Unavailable> File [modified]: momy.xl Delta lines: <Unavailable> File [modified]: momy.yl Delta lines: <Unavailable> File [modified]: momy.zl Delta lines: <Unavailable> File [modified]: momz.dl Delta lines: <Unavailable> File [modified]: momz.xl Delta lines: <Unavailable> File [modified]: momz.yl Delta lines: <Unavailable> File [modified]: momz.zl Delta lines: <Unavailable> Directory: /trunk/test/ ======================= File [modified]: ConstraintViolate.par Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] IDConstraintViolate/ (Rev. 28)
by knarf＠cct.lsu.edu 02 Dec '13

02 Dec '13

User: knarf Date: 2013/12/02 05:36 PM Modified: /trunk/test/ConstraintViolate/ alp.dl, gxx.dl, ham.dl, momx.dl, momy.dl, momz.dl Log: set correct mime type Directory Changes: Directory: /svn:mime-type/ ========================== - video/dl + text/plain File Changes: Directory: /trunk/test/ConstraintViolate/ ========================================= File [modified]: alp.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/ConstraintViolate/alp.dl ___________________________________________________________________ File [modified]: gxx.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/ConstraintViolate/gxx.dl ___________________________________________________________________ File [modified]: ham.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/ConstraintViolate/ham.dl ___________________________________________________________________ File [modified]: momx.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/ConstraintViolate/momx.dl ___________________________________________________________________ File [modified]: momy.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/ConstraintViolate/momy.dl ___________________________________________________________________ File [modified]: momz.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/ConstraintViolate/momz.dl ___________________________________________________________________

1 0

[svn:einsteintoolkit] IDAnalyticBH/trunk/test/ (Rev. 191)
by knarf＠cct.lsu.edu 02 Dec '13

02 Dec '13

User: knarf Date: 2013/12/02 05:34 PM Modified: /trunk/test/ kerr.par, test_bl.par, test_misner.par /trunk/test/kerr/ alp.dl, detg.dl, gpp.dl, gqp.dl, gqq.dl, grp.dl, grq.dl, grr.dl, grr_max.tl, grr_min.tl, grr_nm1.tl, grr_nm2.tl, gxx.dl, gxy.dl, gxz.dl, gyy.dl, gyz.dl, gzz.dl, kpp.dl, kqp.dl, kqq.dl, krp.dl, krq.dl, krr.dl, kxx.dl, kxy.dl, kxz.dl, kyy.dl, kyz.dl, kzz.dl, trK.dl, trK_max.tl, trK_min.tl, trK_nm1.tl, trK_nm2.tl /trunk/test/test_bl/ grr_max.tl, grr_min.tl, grr_nm1.tl, grr_nm2.tl /trunk/test/test_misner/ grr_max.tl, grr_min.tl, grr_nm1.tl, grr_nm2.tl, trK_max.tl, trK_min.tl, trK_nm1.tl, trK_nm2.tl Log: don't use evolution (ADM) in ID testsuite File Changes: Directory: /trunk/test/kerr/ ============================ File [modified]: alp.dl Delta lines: <Unavailable> File [modified]: detg.dl Delta lines: <Unavailable> File [modified]: gpp.dl Delta lines: <Unavailable> File [modified]: gqp.dl Delta lines: <Unavailable> File [modified]: gqq.dl Delta lines: <Unavailable> File [modified]: grp.dl Delta lines: <Unavailable> File [modified]: grq.dl Delta lines: <Unavailable> File [modified]: grr.dl Delta lines: <Unavailable> File [modified]: grr_max.tl Delta lines: <Unavailable> File [modified]: grr_min.tl Delta lines: <Unavailable> File [modified]: grr_nm1.tl Delta lines: <Unavailable> File [modified]: grr_nm2.tl Delta lines: <Unavailable> File [modified]: gxx.dl Delta lines: <Unavailable> File [modified]: gxy.dl Delta lines: <Unavailable> File [modified]: gxz.dl Delta lines: <Unavailable> File [modified]: gyy.dl Delta lines: <Unavailable> File [modified]: gyz.dl Delta lines: <Unavailable> File [modified]: gzz.dl Delta lines: <Unavailable> File [modified]: kpp.dl Delta lines: <Unavailable> File [modified]: kqp.dl Delta lines: <Unavailable> File [modified]: kqq.dl Delta lines: <Unavailable> File [modified]: krp.dl Delta lines: <Unavailable> File [modified]: krq.dl Delta lines: <Unavailable> File [modified]: krr.dl Delta lines: <Unavailable> File [modified]: kxx.dl Delta lines: <Unavailable> File [modified]: kxy.dl Delta lines: <Unavailable> File [modified]: kxz.dl Delta lines: <Unavailable> File [modified]: kyy.dl Delta lines: <Unavailable> File [modified]: kyz.dl Delta lines: <Unavailable> File [modified]: kzz.dl Delta lines: <Unavailable> File [modified]: trK.dl Delta lines: <Unavailable> File [modified]: trK_max.tl Delta lines: <Unavailable> File [modified]: trK_min.tl Delta lines: <Unavailable> File [modified]: trK_nm1.tl Delta lines: <Unavailable> File [modified]: trK_nm2.tl Delta lines: <Unavailable> Directory: /trunk/test/ ======================= File [modified]: kerr.par Delta lines: <Unavailable> File [modified]: test_bl.par Delta lines: <Unavailable> File [modified]: test_misner.par Delta lines: <Unavailable> Directory: /trunk/test/test_bl/ =============================== File [modified]: grr_max.tl Delta lines: <Unavailable> File [modified]: grr_min.tl Delta lines: <Unavailable> File [modified]: grr_nm1.tl Delta lines: <Unavailable> File [modified]: grr_nm2.tl Delta lines: <Unavailable> Directory: /trunk/test/test_misner/ =================================== File [modified]: grr_max.tl Delta lines: <Unavailable> File [modified]: grr_min.tl Delta lines: <Unavailable> File [modified]: grr_nm1.tl Delta lines: <Unavailable> File [modified]: grr_nm2.tl Delta lines: <Unavailable> File [modified]: trK_max.tl Delta lines: <Unavailable> File [modified]: trK_min.tl Delta lines: <Unavailable> File [modified]: trK_nm1.tl Delta lines: <Unavailable> File [modified]: trK_nm2.tl Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] IDAnalyticBH/ (Rev. 190)
by knarf＠cct.lsu.edu 02 Dec '13

02 Dec '13

User: knarf Date: 2013/12/02 05:05 PM Modified: /trunk/test/kerr/ alp.dl, detg.dl, gpp.dl, gqp.dl, gqq.dl, grp.dl, grq.dl, grr.dl, gxx.dl, gxy.dl, gxz.dl, gyy.dl, gyz.dl, gzz.dl, kpp.dl, kqp.dl, kqq.dl, krp.dl, krq.dl, krr.dl, kxx.dl, kxy.dl, kxz.dl, kyy.dl, kyz.dl, kzz.dl, trK.dl Log: set correct mime type Directory Changes: Directory: /svn:mime-type/ ========================== - video/dl + text/plain File Changes: Directory: /trunk/test/kerr/ ============================ File [modified]: alp.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/alp.dl ___________________________________________________________________ File [modified]: detg.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/detg.dl ___________________________________________________________________ File [modified]: gpp.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/gpp.dl ___________________________________________________________________ File [modified]: gqp.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/gqp.dl ___________________________________________________________________ File [modified]: gqq.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/gqq.dl ___________________________________________________________________ File [modified]: grp.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/grp.dl ___________________________________________________________________ File [modified]: grq.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/grq.dl ___________________________________________________________________ File [modified]: grr.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/grr.dl ___________________________________________________________________ File [modified]: gxx.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/gxx.dl ___________________________________________________________________ File [modified]: gxy.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/gxy.dl ___________________________________________________________________ File [modified]: gxz.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/gxz.dl ___________________________________________________________________ File [modified]: gyy.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/gyy.dl ___________________________________________________________________ File [modified]: gyz.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/gyz.dl ___________________________________________________________________ File [modified]: gzz.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/gzz.dl ___________________________________________________________________ File [modified]: kpp.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/kpp.dl ___________________________________________________________________ File [modified]: kqp.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/kqp.dl ___________________________________________________________________ File [modified]: kqq.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/kqq.dl ___________________________________________________________________ File [modified]: krp.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/krp.dl ___________________________________________________________________ File [modified]: krq.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/krq.dl ___________________________________________________________________ File [modified]: krr.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/krr.dl ___________________________________________________________________ File [modified]: kxx.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/kxx.dl ___________________________________________________________________ File [modified]: kxy.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/kxy.dl ___________________________________________________________________ File [modified]: kxz.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/kxz.dl ___________________________________________________________________ File [modified]: kyy.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/kyy.dl ___________________________________________________________________ File [modified]: kyz.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/kyz.dl ___________________________________________________________________ File [modified]: kzz.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/kzz.dl ___________________________________________________________________ File [modified]: trK.dl Delta lines: +0 -0 =================================================================== Property changes on: trunk/test/kerr/trK.dl ___________________________________________________________________

1 0

[svn:einsteintoolkit] RotatingDBHIVP/trunk/test/ (Rev. 61)
by knarf＠cct.lsu.edu 02 Dec '13

02 Dec '13

User: knarf Date: 2013/12/02 04:53 PM Modified: /trunk/test/ E2xeon_test_rdbh.par Log: remove dependency on ADM File Changes: Directory: /trunk/test/ ======================= File [modified]: E2xeon_test_rdbh.par Delta lines: +1 -2 =================================================================== --- trunk/test/E2xeon_test_rdbh.par 2013-12-02 19:58:30 UTC (rev 60) +++ trunk/test/E2xeon_test_rdbh.par 2013-12-02 22:53:56 UTC (rev 61) @@ -2,7 +2,7 @@ #!DESC "Initial data for axisymmetric distored black hole" ###################################################################### -ActiveThorns = "coordgauge admmacros admanalysis tmunubase admbase staticconformal time boundary ellbase pugh pughinterp pughreduce pughslab cartgrid3d ADM RotatingDBHIVP ioascii iobasic ioutil ADMCoupling CoordBase SpaceMask LocalReduce LocalInterp SymBase" +ActiveThorns = "coordgauge admmacros admanalysis tmunubase admbase staticconformal time boundary ellbase pugh pughinterp pughreduce pughslab cartgrid3d RotatingDBHIVP ioascii iobasic ioutil ADMCoupling CoordBase SpaceMask LocalReduce LocalInterp SymBase" # GENERAL @@ -34,7 +34,6 @@ rotatingdbhivp::interpolation_order = 1 # GAUGE -admbase::lapse_evolution_method = "geodesic" admbase::metric_type = "static conformal" staticconformal::conformal_storage = "factor+derivs+2nd derivs"

1 0

[svn:einsteintoolkit] DistortedBHIVP/trunk/test/ (Rev. 67)
by knarf＠cct.lsu.edu 02 Dec '13

02 Dec '13

User: knarf Date: 2013/12/02 04:49 PM Removed: /trunk/test/ xeon_test_dbh.par /trunk/test/xeon_test_dbh/ Log: remove outdated test File Changes: Directory: /trunk/test/ ======================= File [removed]: xeon_test_dbh.par Delta lines: +0 -55 =================================================================== --- trunk/test/xeon_test_dbh.par 2013-12-02 20:48:49 UTC (rev 66) +++ trunk/test/xeon_test_dbh.par 2013-12-02 22:49:21 UTC (rev 67) @@ -1,55 +0,0 @@ -###################################################################### -#!DESC "Initial data for axisymmetric distored black hole" -###################################################################### - -ActiveThorns = "time boundary iobasic ADMconstraints pugh pughreduce pughslab pughinterp cartgrid3d einstein ADM DistortedBHIVP ioascii ioutil" - -# GENERAL - -einstein::evolution_system = "ADM" - -driver::global_nx = 32 -driver::global_ny = 32 -driver::global_nz = 32 - -grid::type = "byspacing" -grid::dxyz = 0.4 -grid::domain = "full" - -time::dtfac = 0.25 - -cactus::cctk_initial_time = 0. -cactus::cctk_itlast = 1 - -adm::method = "stagleap" -adm::bound = "static" -admconstraints::bound = "static" - -# MODEL - -einstein::initial_data = "distortedbh" -distortedbhivp::amp = -0.1 -distortedbhivp::eta0 = 0.0 -distortedbhivp::c = 0.5 -distortedbhivp::sigma = 1.0 -distortedbhivp::etamax = 6.0 -distortedbhivp::n = 4 -distortedbhivp::neta = 102 -distortedbhivp::ntheta = 27 -distortedbhivp::nphi = 5 -distortedbhivp::interpolation_order = 1 - - -# GAUGE -einstein::slicing = "geodesic" - -# OUTPUT ######################################################## -IO::outdir = "test_dbh" -IOASCII::out1d_every = 1 -IOBasic::outScalar_vars = "einstein::metric einstein::curv admconstraints::ham" -IOASCII::out1D_vars = "einstein::metric einstein::curv einstein::lapse admconstraints::ham einstein::grr einstein::psi einstein::psix einstein::psiy einstein::psiz einstein::psixx einstein::psixy einstein::psixz einstein::psiyy einstein::psiyz einstein::psizz" -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = "admconstraints::ham einstein::gxx" - -################################################################## -

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 594)
by rhaas＠tapir.caltech.edu 02 Dec '13

02 Dec '13

User: rhaas Date: 2013/12/02 03:50 PM Modified: /trunk/src/ GRHydro_ENOReconstruct_drv.F90, GRHydro_MP5Reconstruct_drv.F90, GRHydro_PPMMReconstruct_drv.F90, GRHydro_PPMReconstruct_drv.F90, GRHydro_PPMReconstruct_drv_opt.F90, GRHydro_ReconstructPoly.F90, GRHydro_TVDReconstruct_drv.F90, GRHydro_WENOReconstruct_drv.F90 Log: GRHydro: use ash instead of lsh when allocating arrays File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_ENOReconstruct_drv.F90 Delta lines: +4 -4 =================================================================== --- trunk/src/GRHydro_ENOReconstruct_drv.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_ENOReconstruct_drv.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -61,7 +61,7 @@ CCTK_INT :: ierr - allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if @@ -69,9 +69,9 @@ !!$ The dum variables are used as dummies if MHD on but divergence cleaning off - allocate(dum(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dump(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dumm(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(dum(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dump(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dumm(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with dum dump dumm") File [modified]: GRHydro_MP5Reconstruct_drv.F90 Delta lines: +4 -4 =================================================================== --- trunk/src/GRHydro_MP5Reconstruct_drv.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_MP5Reconstruct_drv.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -65,15 +65,15 @@ CCTK_INT :: ierr - allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if !!$ The dum variables are used as dummies if MHD on but divergence cleaning off - allocate(dum(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dump(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dumm(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(dum(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dump(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dumm(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with dum dump dumm") File [modified]: GRHydro_PPMMReconstruct_drv.F90 Delta lines: +4 -4 =================================================================== --- trunk/src/GRHydro_PPMMReconstruct_drv.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_PPMMReconstruct_drv.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -61,15 +61,15 @@ CCTK_INT :: ierr - allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if !!$ The dum variables are used as dummies if MHD on but divergence cleaning off - allocate(dum(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dump(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dumm(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(dum(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dump(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dumm(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) call SpaceMask_GetTypeBits(type_bitsx, "Hydro_RiemannProblemX") call SpaceMask_GetStateBits(trivialx, "Hydro_RiemannProblemX", & File [modified]: GRHydro_PPMReconstruct_drv.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_PPMReconstruct_drv.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_PPMReconstruct_drv.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -57,7 +57,7 @@ logical :: apply_enhanced_ppm - allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if File [modified]: GRHydro_PPMReconstruct_drv_opt.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_PPMReconstruct_drv_opt.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_PPMReconstruct_drv_opt.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -90,7 +90,7 @@ pvup = loc(vel) end if - allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if File [modified]: GRHydro_ReconstructPoly.F90 Delta lines: +8 -8 =================================================================== --- trunk/src/GRHydro_ReconstructPoly.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_ReconstructPoly.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -74,18 +74,18 @@ CCTK_INT :: GRHydro_UseGeneralCoordinates - allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if - allocate(psi4(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - lbetax(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - lbetay(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - lbetaz(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dum(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dump(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dumm(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(psi4(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + lbetax(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + lbetay(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + lbetaz(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dum(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dump(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dumm(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with lbeta") File [modified]: GRHydro_TVDReconstruct_drv.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_TVDReconstruct_drv.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_TVDReconstruct_drv.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -84,7 +84,7 @@ CCTK_REAL, DIMENSION(:,:,:,:), ALLOCATABLE :: wvel ! vel*w_lorentz CCTK_REAL :: agxx,agxy,agxz,agyy,agyz,agzz ! metric components - allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if File [modified]: GRHydro_WENOReconstruct_drv.F90 Delta lines: +4 -4 =================================================================== --- trunk/src/GRHydro_WENOReconstruct_drv.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_WENOReconstruct_drv.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -125,15 +125,15 @@ ! CCTK_REAL, dimension(ny) :: vxY,vyY,vzY ! CCTK_REAL, dimension(nz) :: vxZ,vyZ,vzZ - allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if !!$ The dum variables are used as dummies if MHD on but divergence cleaning off - allocate(dum(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dump(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dumm(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(dum(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dump(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dumm(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with dum dump dumm")

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 593)
by rhaas＠tapir.caltech.edu 02 Dec '13

02 Dec '13

User: rhaas Date: 2013/12/02 03:50 PM Modified: /trunk/src/ GRHydro_Con2PrimM.F90, GRHydro_UpdateMask.F90, GRHydro_UpdateMaskM.F90 Log: GRHydro: Better handling of excision for MHD. File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2PrimM.F90 Delta lines: +8 -43 =================================================================== --- trunk/src/GRHydro_Con2PrimM.F90 2013-11-25 15:39:05 UTC (rev 592) +++ trunk/src/GRHydro_Con2PrimM.F90 2013-12-02 21:50:25 UTC (rev 593) @@ -164,10 +164,6 @@ !do not compute if in atmosphere region! if (atmosphere_mask(i,j,k) .gt. 0) cycle - !do not compute if in atmosphere or in excised region - !if ((atmosphere_mask(i,j,k) .ne. 0) .or. & - ! (hydro_excision_mask(i,j,k) .ne. 0)) cycle - epsnegative = 0 sdet = sdetg(i,j,k) @@ -192,43 +188,6 @@ endif - ! In excised region, set to atmosphere! - if (GRHydro_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .gt. 0)) then - SET_ATMO_MIN(dens(i,j,k), sdet*GRHydro_rho_min, r(i,j,k)) !sqrt(det)*GRHydro_rho_min !/(1.d0+GRHydro_atmo_tolerance) - SET_ATMO_MIN(rho(i,j,k), GRHydro_rho_min, r(i,j,k)) !GRHydro_rho_min - scon(i,j,k,:) = 0.d0 - vup(i,j,k,:) = 0.d0 - w_lorentz(i,j,k) = 1.d0 - - if(evolve_temper.ne.0) then - ! set hot atmosphere values - temperature(i,j,k) = grhydro_hot_atmo_temp - y_e(i,j,k) = grhydro_hot_atmo_Y_e - y_e_con(i,j,k) = y_e(i,j,k) * dens(i,j,k) - keytemp = 1 - call EOS_Omni_press(GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,n,& - rho(i,j,k),eps(i,j,k),temperature(i,j,k),y_e(i,j,k),& - press(i,j,k),keyerr,anyerr) - else - keytemp = 0 - call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& - rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),keyerr,anyerr) - call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& - rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),eps(i,j,k),keyerr,anyerr) - endif - - !!$ tau does need to take into account the existing B-field - !!$ with w_lorentz=1, we find tau = sqrtdet*(rho (1+eps+b^2/2)) - dens [Press drops out] - tau(i,j,k) = sdet * (rho(i,j,k)*(1.0+eps(i,j,k)+b2/2.0)) - dens(i,j,k) - - if(tau(i,j,k).le.sdet*b2*0.5d0)then - tau(i,j,k) = GRHydro_tau_min + sdet*b2*0.5d0 - endif - - cycle - - endif - if(evolve_Y_e.ne.0) then Y_e(i,j,k) = max(min(Y_e_con(i,j,k) / dens(i,j,k),GRHydro_Y_e_max),& GRHydro_Y_e_min) @@ -239,8 +198,8 @@ 2.0*(g12(i,j,k)*Bprim(i,j,k,1)*Bprim(i,j,k,2)+g13(i,j,k)*Bprim(i,j,k,1)*Bprim(i,j,k,3)+ & g23(i,j,k)*Bprim(i,j,k,2)*Bprim(i,j,k,3)) - - if ( dens(i,j,k) .le. sdet*GRHydro_rho_min*(1.d0+GRHydro_atmo_tolerance) ) then + + if ( (GRHydro_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .gt. 0)) .or. (dens(i,j,k) .le. sdet*GRHydro_rho_min*(1.d0+GRHydro_atmo_tolerance))) then !call CCTK_WARN(1,"Con2Prim: Resetting to atmosphere") !write(warnline,"(3i5,1P10E15.6)") i,j,k,x(i,j,k),y(i,j,k),z(i,j,k) @@ -256,6 +215,10 @@ vup(i,j,k,:) = 0.d0 w_lorentz(i,j,k) = 1.d0 + Bprim(i,j,k,:) = 0.0d0 + Bcons(i,j,k,:) = 0.0d0 + b2 = 0 + if(evolve_temper.ne.0) then ! set hot atmosphere values temperature(i,j,k) = grhydro_hot_atmo_temp @@ -294,6 +257,8 @@ end if + + if(evolve_temper.eq.0) then File [modified]: GRHydro_UpdateMask.F90 Delta lines: +25 -0 =================================================================== --- trunk/src/GRHydro_UpdateMask.F90 2013-11-25 15:39:05 UTC (rev 592) +++ trunk/src/GRHydro_UpdateMask.F90 2013-12-02 21:50:25 UTC (rev 593) @@ -10,6 +10,7 @@ #include "cctk.h" #include "cctk_Parameters.h" #include "cctk_Arguments.h" +#include "cctk_Functions.h" #include "GRHydro_Macros.h" #include "SpaceMask.h" @@ -30,12 +31,24 @@ DECLARE_CCTK_ARGUMENTS DECLARE_CCTK_PARAMETERS + DECLARE_CCTK_FUNCTIONS CCTK_INT :: i,j,k CCTK_REAL :: frac + logical :: evolve_Bvec + logical :: evolve_Avec frac = CCTK_DELTA_TIME + evolve_Bvec = .false. + evolve_Avec = .false. + + if (CCTK_EQUALS(Bvec_evolution_method,"GRHydro")) then + evolve_Bvec = .true. + elseif (CCTK_EQUALS(Bvec_evolution_method,"GRHydro_Avec") .or. CCTK_EQUALS(Bvec_evolution_method,"GRHydro_Avec_centered")) then + evolve_Avec = .true. + endif + if(evolve_temper.ne.1.and.evolve_Y_e.ne.1) then !$OMP PARALLEL DO PRIVATE(k,j,i) do k = 1+GRHydro_stencil, cctk_lsh(3)-GRHydro_stencil @@ -48,6 +61,12 @@ densrhs(i,j,k) = 0.0d0 srhs(i,j,k,:) = 0.0d0 taurhs(i,j,k) = 0.0d0 + + if (evolve_Bvec) then + Bconsrhs(i,j,k,:) = 0.0d0 + endif + + ! TODO: Need to set Avecrhs and Aphirhs to zero for centered / or staggered case if vector potential is evolved ! Set real-valued mask! This will be sync'ed and right after syncing translated to ! our standard integer based mask (so that atmosphere_mask is still valid!). @@ -70,7 +89,13 @@ densrhs(i,j,k) = 0.0d0 srhs(i,j,k,:) = 0.0d0 taurhs(i,j,k) = 0.0d0 + + if (evolve_Bvec) then + Bconsrhs(i,j,k,:) = 0.0d0 + endif + ! TODO: Need to set Avecrhs and Aphirhs to zero for centered / or staggered case if vector potential is evolved + ! Set real-valued mask! This will be sync'ed and right after syncing translated to ! our standard integer based mask (so that atmosphere_mask is still valid!). atmosphere_mask_real(i,j,k) = 1 File [modified]: GRHydro_UpdateMaskM.F90 Delta lines: +2 -0 =================================================================== --- trunk/src/GRHydro_UpdateMaskM.F90 2013-11-25 15:39:05 UTC (rev 592) +++ trunk/src/GRHydro_UpdateMaskM.F90 2013-12-02 21:50:25 UTC (rev 593) @@ -115,6 +115,8 @@ vup(i,j,k,3) = 0.0d0 sdet = sdetg(i,j,k) + Bprim(i,j,k,:) = 0.0d0 + if(evolve_temper.ne.0) then ! set the temperature to be relatively low

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits December 2013