Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 481)
by rhaas＠tapir.caltech.edu 11 Feb '13

11 Feb '13

User: rhaas Date: 2013/02/11 10:30 AM Modified: /trunk/doc/ documentation.tex /trunk/src/ GRHydro_Con2PrimM.F90 Log: GRHydro: remove commented out code, update docs From: Roland Haas <rhaas(a)tapir.caltech.edu> File Changes: Directory: /trunk/doc/ ====================== File [modified]: documentation.tex Delta lines: +1 -1 =================================================================== --- trunk/doc/documentation.tex 2013-02-11 16:30:22 UTC (rev 480) +++ trunk/doc/documentation.tex 2013-02-11 16:30:23 UTC (rev 481) @@ -1991,7 +1991,7 @@ encouragement that it was possible! Thirdly, the support of the Cactus team, especially Tom Goodale, -Gabrielle Allen and Thomas Radke made (and continues to make) life a +Gabrielle Allen and Thomas Radke made life a lot easier. Finally, for their work in coding, ideas and suggestions, or just Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2PrimM.F90 Delta lines: +0 -12 =================================================================== --- trunk/src/GRHydro_Con2PrimM.F90 2013-02-11 16:30:22 UTC (rev 480) +++ trunk/src/GRHydro_Con2PrimM.F90 2013-02-11 16:30:23 UTC (rev 481) @@ -660,18 +660,6 @@ end if -! ! Again, reset 3-metric only for con2prim inversion. Restoring -! ! the current values for the 3-metric: -! if(.not.posdef)then -! g11(i,j,k) = g11c -! g12(i,j,k) = g12c -! g13(i,j,k) = g13c -! g22(i,j,k) = g22c -! g23(i,j,k) = g23c -! g33(i,j,k) = g33c -! posdef = .true. -! endif -! end do end do end do

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 480)
by rhaas＠tapir.caltech.edu 11 Feb '13

11 Feb '13

User: rhaas Date: 2013/02/11 10:30 AM Modified: /trunk/ schedule.ccl Log: GRHydro: synchronize press for initial guess of Con2Prim since the Polytype (!) con2prim routine uses it for its initial guess then continues to iterate over rho Do not sync variables not used for initial guess. This affects only runs with sync_only_conservatives=yes. From: Roland Haas <rhaas(a)tapir.caltech.edu> File Changes: Directory: /trunk/ ================== File [modified]: schedule.ccl Delta lines: +13 -13 =================================================================== --- trunk/schedule.ccl 2013-02-11 16:30:19 UTC (rev 479) +++ trunk/schedule.ccl 2013-02-11 16:30:22 UTC (rev 480) @@ -1268,9 +1268,9 @@ LANG: Fortran OPTIONS: LEVEL # SYNC: HydroBase::w_lorentz # not used as initial guess and hence not needed! -# SYNC: HydroBase::press # not used as initial guess and hence not needed! -# SYNC: HydroBase::vel -# SYNC: HydroBase::Bvec + SYNC: HydroBase::press # con2prim_ptPolytype uses this for its initial guess +# SYNC: HydroBase::vel # not needed; lvel is used as initial guess instead +# SYNC: HydroBase::Bvec # not needed; lBvec is used as initial guess instead # SYNC: hydrobase::entropy # not used as initial guess and hence not needed! # SYNC: hydrobase::Y_e # not used as initial guess and hence not needed! # SYNC: GRHydro_tracers # not used as initial guess and hence not needed! @@ -1287,10 +1287,10 @@ { LANG: Fortran OPTIONS: LEVEL - SYNC: HydroBase::w_lorentz - SYNC: HydroBase::press -# SYNC: HydroBase::vel -# SYNC: HydroBase::Bvec +# SYNC: HydroBase::w_lorentz # don't need that since computed in Con2Prim which runs in postregrid! + SYNC: HydroBase::press # con2prim_ptPolytype uses this for its initial guess +# SYNC: HydroBase::vel # not needed; obtained from lvel everywhere in post regrid after Con2Prim +# SYNC: HydroBase::Bvec # not needed; obtained from lBvec everywhere in post regrid after Con2Prim SYNC: hydrobase::entropy SYNC: hydrobase::Y_e SYNC: GRHydro_tracers @@ -1310,7 +1310,7 @@ LANG: Fortran OPTIONS: LEVEL # SYNC: HydroBase::w_lorentz # not used as initial guess and hence not needed! -# SYNC: HydroBase::press # not used as initial guess and hence not needed! + SYNC: HydroBase::press # con2prim_ptPolytype uses this for its initial guess # SYNC: hydrobase::entropy # not used as initial guess and hence not needed! # SYNC: hydrobase::Y_e # not used as initial guess and hence not needed! # SYNC: lvel @@ -1329,11 +1329,11 @@ { LANG: Fortran OPTIONS: LEVEL - SYNC: HydroBase::w_lorentz # not used as initial guess and hence not needed! - SYNC: HydroBase::press # not used as initial guess and hence not needed! - SYNC: hydrobase::entropy # not used as initial guess and hence not needed! - SYNC: hydrobase::Y_e # not used as initial guess and hence not needed! - SYNC: GRHydro_tracers # not used as initial guess and hence not needed! +# SYNC: HydroBase::w_lorentz # don't need that since computed in Con2Prim which runs in postregrid! + SYNC: HydroBase::press # con2prim_ptPolytype uses this for its initial guess + SYNC: hydrobase::entropy + SYNC: hydrobase::Y_e + SYNC: GRHydro_tracers SYNC: HydroBase::rho SYNC: HydroBase::eps SYNC: HydroBase::vel

1 0

[svn:einsteintoolkit] GRHydro/trunk/test/ (Rev. 479)
by rhaas＠tapir.caltech.edu 11 Feb '13

11 Feb '13

User: rhaas Date: 2013/02/11 10:30 AM Modified: /trunk/test/ test_tov_carpet_refined_nosync.par, tov_slowsector.par /trunk/test/test_tov_carpet_refined_nosync/ rho.x.asc, rho.y.asc, rho.z.asc /trunk/test/tov_slowsector/ gxx_x_[3][3].xg, gxx_y_[3][3].xg, gxx_z_[3][3].xg, rho_x_[3][3].xg, rho_y_[3][3].xg, rho_z_[3][3].xg Log: GRHydro: regenerate test data * after rev 72 of EOS_Omni changed constants of nature * increase con2prim accuracy to avoid O(1) errors in primitives due to O(eps) changes in conservatives * remove changes in K and Gamma to avoid setting points to atmosphere From: Roland Haas <rhaas(a)tapir.caltech.edu> File Changes: Directory: /trunk/test/ ======================= File [modified]: test_tov_carpet_refined_nosync.par Delta lines: <Unavailable> File [modified]: tov_slowsector.par Delta lines: <Unavailable> Directory: /trunk/test/test_tov_carpet_refined_nosync/ ====================================================== File [modified]: rho.x.asc Delta lines: <Unavailable> File [modified]: rho.y.asc Delta lines: <Unavailable> File [modified]: rho.z.asc Delta lines: <Unavailable> Directory: /trunk/test/tov_slowsector/ ====================================== File [modified]: gxx_x_[3][3].xg Delta lines: <Unavailable> File [modified]: gxx_y_[3][3].xg Delta lines: <Unavailable> File [modified]: gxx_z_[3][3].xg Delta lines: <Unavailable> File [modified]: rho_x_[3][3].xg Delta lines: <Unavailable> File [modified]: rho_y_[3][3].xg Delta lines: <Unavailable> File [modified]: rho_z_[3][3].xg Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] GRHydro/trunk/test/ (Rev. 478)
by rhaas＠tapir.caltech.edu 11 Feb '13

11 Feb '13

User: rhaas Date: 2013/02/11 10:30 AM Modified: /trunk/test/ test_tov_carpet_refined_nosync.par, tov_slowsector.par /trunk/test/test_tov_carpet_refined_nosync/ rho.x.asc, rho.y.asc, rho.z.asc /trunk/test/tov_slowsector/ rho_x_[3][3].xg, rho_y_[3][3].xg, rho_z_[3][3].xg Log: GRHydro: update tests where poly_gamma_ini was not set properly this affect the test_tov_carpet_refined_nosync and tov_slowsector tests. From: Roland Haas <roland.haas(a)physics.gatech.edu> File Changes: Directory: /trunk/test/ ======================= File [modified]: test_tov_carpet_refined_nosync.par Delta lines: <Unavailable> File [modified]: tov_slowsector.par Delta lines: <Unavailable> Directory: /trunk/test/test_tov_carpet_refined_nosync/ ====================================================== File [modified]: rho.x.asc Delta lines: <Unavailable> File [modified]: rho.y.asc Delta lines: <Unavailable> File [modified]: rho.z.asc Delta lines: <Unavailable> Directory: /trunk/test/tov_slowsector/ ====================================== File [modified]: rho_x_[3][3].xg Delta lines: <Unavailable> File [modified]: rho_y_[3][3].xg Delta lines: <Unavailable> File [modified]: rho_z_[3][3].xg Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 477)
by rhaas＠tapir.caltech.edu 11 Feb '13

11 Feb '13

User: rhaas Date: 2013/02/11 10:30 AM Modified: /trunk/src/ GRHydro_HLLEM.F90 Log: GRHydro: declare charXXX_dc variables in HLLEM OMP private From: Philipp Moesta <pmoesta(a)tapir.caltech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_HLLEM.F90 Delta lines: +2 -1 =================================================================== --- trunk/src/GRHydro_HLLEM.F90 2013-02-11 16:30:12 UTC (rev 476) +++ trunk/src/GRHydro_HLLEM.F90 2013-02-11 16:30:14 UTC (rev 477) @@ -136,7 +136,8 @@ !$OMP bdotvp,bdotvm,b2p,b2m,v2p,v2m,wp,wm,& !$OMP bxlowp,bxlowm,bylowp,bylowm,bzlowp,bzlowm,& !$OMP rhoenth_p,rhoenth_m,ab0p,ab0m,vA2p,vA2m,pressstarp,pressstarm,& - !$OMP usendh,psidcdiff,psidcfp,psidcfm,charmin,charmax,chartop,charpm,m,xtemp,& + !$OMP usendh,psidcdiff,psidcfp,psidcfm,charmin,charmax,chartop,charpm,& + !$OMP charmin_dc,charmax_dc,charpm_dc,m,xtemp,& !$OMP entropyconsp,entropyconsm,entropyp,entropym,entropyf,entropydiff,entropyfp,entropyfm) do k = GRHydro_stencil, cctk_lsh(3) - GRHydro_stencil + transport_constraints*(1-zoffset) do j = GRHydro_stencil, cctk_lsh(2) - GRHydro_stencil + transport_constraints*(1-yoffset)

1 0

[svn:einsteintoolkit] GRHydro/trunk/test/ (Rev. 476)
by rhaas＠tapir.caltech.edu 11 Feb '13

11 Feb '13

User: rhaas Date: 2013/02/11 10:30 AM Added: /trunk/test/ tov_carpetevolutionmask2.par /trunk/test/tov_carpetevolutionmask2/ evolution_mask.x.asc, evolution_mask.y.asc, evolution_mask.z.asc, rho.x.asc, rho.y.asc, rho.z.asc Modified: /trunk/test/ test.ccl Log: GRHydro: add evolution mask test with more levels and gamma law From: Roland Haas <rhaas(a)tapir.caltech.edu> File Changes: Directory: /trunk/test/ ======================= File [modified]: test.ccl Delta lines: <Unavailable> File [added]: tov_carpetevolutionmask2.par Delta lines: <Unavailable> Directory: /trunk/test/tov_carpetevolutionmask2/ ================================================ File [added]: evolution_mask.x.asc Delta lines: <Unavailable> File [added]: evolution_mask.y.asc Delta lines: <Unavailable> File [added]: evolution_mask.z.asc Delta lines: <Unavailable> File [added]: rho.x.asc Delta lines: <Unavailable> File [added]: rho.y.asc Delta lines: <Unavailable> File [added]: rho.z.asc Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] GRHydro/trunk/test/ (Rev. 475)
by rhaas＠tapir.caltech.edu 11 Feb '13

11 Feb '13

User: rhaas Date: 2013/02/11 10:30 AM Modified: /trunk/test/ tov_carpetevolutionmask.par Log: GRHydro: explicitly set EOS parameters in evolution_mask test remove unused Thorn from ActiveThorns. From: Roland Haas <rhaas(a)tapir.caltech.edu> File Changes: Directory: /trunk/test/ ======================= File [modified]: tov_carpetevolutionmask.par Delta lines: +6 -2 =================================================================== --- trunk/test/tov_carpetevolutionmask.par 2013-02-11 16:30:08 UTC (rev 474) +++ trunk/test/tov_carpetevolutionmask.par 2013-02-11 16:30:10 UTC (rev 475) @@ -1,5 +1,5 @@ -!DESC "Carpet mesh refined, one TOV, synchronize only conservatives" -ActiveThorns = "time MoL coordbase symbase cartgrid3d ioutil boundary admbase staticconformal spacemask admcoupling coordgauge admmacros hydrobase grhydro grhydro_initdata tovsolver eos_omni Constants LocalReduce Carpet CarpetLib CarpetInterp CarpetReduce CarpetIOScalar CarpetIOBasic LoopControl TmunuBase InitBase CarpetRegrid2 ReflectionSymmetry CarpetIOASCII CarpetEvolutionMask" +!DESC "Carpet mesh refined, one TOV, polytrope, test evolution mask" +ActiveThorns = "time MoL coordbase symbase cartgrid3d ioutil boundary admbase staticconformal spacemask admcoupling coordgauge admmacros hydrobase grhydro grhydro_initdata tovsolver eos_omni Constants Carpet CarpetLib CarpetInterp CarpetReduce CarpetIOScalar CarpetIOBasic LoopControl TmunuBase InitBase CarpetRegrid2 ReflectionSymmetry CarpetIOASCII CarpetEvolutionMask" TmunuBase::stress_energy_storage = yes TmunuBase::stress_energy_at_RHS = yes @@ -90,6 +90,10 @@ tovsolver::TOV_K[0] = 100.0 tovsolver::TOV_Combine_Method = "maximum" +eos_omni::poly_gamma=2.0 +eos_omni::poly_gamma_ini=2.0 +eos_omni::poly_k=100. + Carpet::prolongation_order_space = 3 Carpet::prolongation_order_time = 2

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 474)
by rhaas＠tapir.caltech.edu 11 Feb '13

11 Feb '13

User: rhaas Date: 2013/02/11 10:30 AM Modified: /trunk/src/ GRHydro_Con2Prim.F90 Log: GRHydro: add fuzzy comparison for evolution_mask evolution_mask is every only set to 0d0 or 1d0 (explicitly) right now. Should this change, this tests needs to be amended. From: Roland Haas <rhaas(a)tapir.caltech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2013-02-11 16:30:06 UTC (rev 473) +++ trunk/src/GRHydro_Con2Prim.F90 2013-02-11 16:30:08 UTC (rev 474) @@ -2283,7 +2283,7 @@ ! do not collapse conditions since Fortran does not guarantee an order if (check_evolution_mask.ne.0) then - if (evolution_mask(i,j,k).eq.0d0) then + if (abs(evolution_mask(i,j,k)).le.0.5d0) then cycle end if end if

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 473)
by rhaas＠tapir.caltech.edu 11 Feb '13

11 Feb '13

User: rhaas Date: 2013/02/11 10:30 AM Modified: /trunk/ interface.ccl, param.ccl, schedule.ccl /trunk/src/ GRHydro_Con2Prim.F90, GRHydro_ENOReconstruct_drv.F90, GRHydro_Eigenproblem.F90, GRHydro_Eigenproblem_Marquina.F90, GRHydro_HLLE.F90, GRHydro_Marquina.F90, GRHydro_PPM.F90, GRHydro_PPMReconstruct_drv.F90, GRHydro_Prim2Con.F90, GRHydro_TVDReconstruct_drv.F90, GRHydro_UpdateMask.F90, GRHydro_WENOReconstruct_drv.F90 Log: * add option to reconstruct temperature to GRHydro This does not yet work with the MHD part of the code, but Philipp is going to change this. * Improve TOVSolverHot; in particular, give functions unique names to avoid duplicate symbols with TOVSolverC From: Christian Ott <cott(a)bethe.tapir.caltech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: <Unavailable> File [modified]: GRHydro_ENOReconstruct_drv.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Eigenproblem.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Eigenproblem_Marquina.F90 Delta lines: <Unavailable> File [modified]: GRHydro_HLLE.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Marquina.F90 Delta lines: <Unavailable> File [modified]: GRHydro_PPM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_PPMReconstruct_drv.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Prim2Con.F90 Delta lines: <Unavailable> File [modified]: GRHydro_TVDReconstruct_drv.F90 Delta lines: <Unavailable> File [modified]: GRHydro_UpdateMask.F90 Delta lines: <Unavailable> File [modified]: GRHydro_WENOReconstruct_drv.F90 Delta lines: <Unavailable> Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: <Unavailable> File [modified]: param.ccl Delta lines: <Unavailable> File [modified]: schedule.ccl Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 472)
by rhaas＠tapir.caltech.edu 11 Feb '13

11 Feb '13

User: rhaas Date: 2013/02/11 10:30 AM Modified: /trunk/src/ GRHydro_Con2PrimAM.F90 Log: GRHydro: disable in-AH polytrope call by default for Avec From: Roland Haas <rhaas(a)tapir.caltech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2PrimAM.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_Con2PrimAM.F90 2013-01-20 18:20:48 UTC (rev 471) +++ trunk/src/GRHydro_Con2PrimAM.F90 2013-02-11 16:30:02 UTC (rev 472) @@ -373,7 +373,7 @@ local_pgam=log(xpress(1)/local_K)/log(xrho(1)) sc = local_K*dens(i,j,k) - if(sdet.ge.sqrtdet_thr) then + if(sqrtdet_thr.gt.0d0 .and. sdet.ge.sqrtdet_thr) then GRHydro_C2P_failed(i,j,k) = 0 rho_tmp = rho(i,j,k)

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits