Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 533)
by rhaas＠tapir.caltech.edu 29 May '13

29 May '13

User: rhaas Date: 2013/05/29 11:35 AM Modified: /trunk/src/ GRHydro_ParamCheck.F90 Log: GRHydro: add PARAMCHECK for grhydro_stencil for WENO5 WENO5 needs at least a stencil width of 3 From: Roland Haas <rhaas(a)tapir.caltech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_ParamCheck.F90 Delta lines: +4 -0 =================================================================== --- trunk/src/GRHydro_ParamCheck.F90 2013-05-29 16:35:16 UTC (rev 532) +++ trunk/src/GRHydro_ParamCheck.F90 2013-05-29 16:35:18 UTC (rev 533) @@ -53,6 +53,10 @@ end if end if + if (CCTK_EQUALS(recon_method,"weno").and.(GRHydro_stencil < 3)) then + call CCTK_PARAMWARN("The stencil size must be at least 3 to use WENO reconstruction.") + end if + if (CCTK_EQUALS(recon_method,"eno").and.(GRHydro_stencil < eno_order)) then call CCTK_PARAMWARN("The stencil size must be at least the order of the reconstruction to use ENO reconstruction") end if

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 532)
by rhaas＠tapir.caltech.edu 29 May '13

29 May '13

User: rhaas Date: 2013/05/29 11:35 AM Modified: /trunk/src/ GRHydro_PPM.F90 Log: GRHydro: fixed critical bug in temperature reconstruction in ePPM From: Christian David Ott <cott(a)zwicky.(none)> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_PPM.F90 Delta lines: +7 -6 =================================================================== --- trunk/src/GRHydro_PPM.F90 2013-05-29 16:35:12 UTC (rev 531) +++ trunk/src/GRHydro_PPM.F90 2013-05-29 16:35:16 UTC (rev 532) @@ -1158,13 +1158,14 @@ tempminus(i+1) = tempplus(i) enddo - - !! Finally compute pressure gradient needed for flattening and shock detection - do i = 2, nx-1 - dpress(i) = press(i+1) - press(i-1) - end do - + endif + + !! Finally compute pressure gradient needed for flattening and shock detection + do i = 2, nx-1 + dpress(i) = press(i+1) - press(i-1) + end do + endif !!$ Zone flattening. See appendix of C&W, p. 197-8.

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 531)
by rhaas＠tapir.caltech.edu 29 May '13

29 May '13

User: rhaas Date: 2013/05/29 11:35 AM Modified: /trunk/src/ GRHydro_ENOReconstruct_drv.F90, GRHydro_PPMMReconstruct_drv.F90, GRHydro_PPMReconstruct_drv.F90, GRHydro_TVDReconstruct_drv.F90, GRHydro_WENOReconstruct_drv.F90 Log: GRHydro: initialy set y_e_plus/minus to cell center value in reconstruct This ensures that all EOS calls have some sensible y_e to work with. From: Christian Ott <cott(a)bethe.tapir.caltech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_ENOReconstruct_drv.F90 Delta lines: +5 -2 =================================================================== --- trunk/src/GRHydro_ENOReconstruct_drv.F90 2013-05-29 16:35:08 UTC (rev 530) +++ trunk/src/GRHydro_ENOReconstruct_drv.F90 2013-05-29 16:35:12 UTC (rev 531) @@ -157,8 +157,11 @@ endif if (evolve_Y_e .ne. 0) then - Y_e_plus(i,j,k) = 0.0d0 - Y_e_minus(i,j,k) = 0.0d0 + ! set this to the cell center values + ! to make sure we have good Y_e even in + ! the boundary region (for full GF EOS calls) + Y_e_plus(i,j,k) = Y_e(i,j,k) + Y_e_minus(i,j,k) = Y_e(i,j,k) endif if (evolve_temper .ne. 0) then File [modified]: GRHydro_PPMMReconstruct_drv.F90 Delta lines: +8 -2 =================================================================== --- trunk/src/GRHydro_PPMMReconstruct_drv.F90 2013-05-29 16:35:08 UTC (rev 530) +++ trunk/src/GRHydro_PPMMReconstruct_drv.F90 2013-05-29 16:35:12 UTC (rev 531) @@ -175,10 +175,16 @@ endif if (evolve_Y_e .ne. 0) then - Y_e_plus(i,j,k) = 0.0d0 - Y_e_minus(i,j,k) = 0.0d0 + Y_e_plus(i,j,k) = Y_e(i,j,k) + Y_e_minus(i,j,k) = Y_e(i,j,k) endif + if(evolve_temper .ne. 0) then + tempplus(i,j,k) = temperature(i,j,k) + tempminus(i,j,k) = temperature(i,j,k) + endif + + enddo enddo enddo File [modified]: GRHydro_PPMReconstruct_drv.F90 Delta lines: +5 -2 =================================================================== --- trunk/src/GRHydro_PPMReconstruct_drv.F90 2013-05-29 16:35:08 UTC (rev 530) +++ trunk/src/GRHydro_PPMReconstruct_drv.F90 2013-05-29 16:35:12 UTC (rev 531) @@ -142,8 +142,11 @@ endif if (evolve_Y_e .ne. 0) then - Y_e_plus(i,j,k) = 0.0d0 - Y_e_minus(i,j,k) = 0.0d0 + ! set this to the cell center values + ! to make sure we have good Y_e even in + ! the boundary region (for full GF EOS calls) + Y_e_plus(i,j,k) = Y_e(i,j,k) + Y_e_minus(i,j,k) = Y_e(i,j,k) endif if (evolve_temper .ne. 0) then File [modified]: GRHydro_TVDReconstruct_drv.F90 Delta lines: +5 -2 =================================================================== --- trunk/src/GRHydro_TVDReconstruct_drv.F90 2013-05-29 16:35:08 UTC (rev 530) +++ trunk/src/GRHydro_TVDReconstruct_drv.F90 2013-05-29 16:35:12 UTC (rev 531) @@ -189,8 +189,11 @@ endif if (evolve_Y_e .ne. 0) then - Y_e_plus(i,j,k) = 0.0d0 - Y_e_minus(i,j,k) = 0.0d0 + ! set this to the cell center values + ! to make sure we have good Y_e even in + ! the boundary region (for full GF EOS calls) + Y_e_plus(i,j,k) = Y_e(i,j,k) + Y_e_minus(i,j,k) = Y_e(i,j,k) endif if(evolve_temper .ne. 0) then File [modified]: GRHydro_WENOReconstruct_drv.F90 Delta lines: +5 -2 =================================================================== --- trunk/src/GRHydro_WENOReconstruct_drv.F90 2013-05-29 16:35:08 UTC (rev 530) +++ trunk/src/GRHydro_WENOReconstruct_drv.F90 2013-05-29 16:35:12 UTC (rev 531) @@ -234,8 +234,11 @@ endif if (evolve_Y_e .ne. 0) then - Y_e_plus(i,j,k) = 0.0d0 - Y_e_minus(i,j,k) = 0.0d0 + ! set this to the cell center values + ! to make sure we have good Y_e even in + ! the boundary region (for full GF EOS calls) + Y_e_plus(i,j,k) = Y_e(i,j,k) + Y_e_minus(i,j,k) = Y_e(i,j,k) endif if (evolve_temper .ne. 0) then

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 530)
by rhaas＠tapir.caltech.edu 29 May '13

29 May '13

User: rhaas Date: 2013/05/29 11:35 AM Modified: /trunk/ param.ccl Log: GRHydro: correct parameter description text From: Roland Haas <rhaas(a)tapir.caltech.edu> File Changes: Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +1 -1 =================================================================== --- trunk/param.ccl 2013-05-29 16:35:04 UTC (rev 529) +++ trunk/param.ccl 2013-05-29 16:35:08 UTC (rev 530) @@ -83,7 +83,7 @@ restricted: -KEYWORD use_evolution_mask "Set this to 'true' to skip validity tests in regions where CarpetEvolutionMask::evolution_mask vanishes." STEERABLE=ALWAYS +KEYWORD use_evolution_mask "Set this to 'always' to skip validity tests in regions where CarpetEvolutionMask::evolution_mask vanishes." STEERABLE=ALWAYS { "always" :: "use the mask" "auto" :: "check if CarpetEvolutionMask is active, then use the mask"

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 529)
by rhaas＠tapir.caltech.edu 29 May '13

29 May '13

User: rhaas Date: 2013/05/29 11:35 AM Modified: /trunk/ param.ccl /trunk/src/ GRHydro_Con2Prim.F90 Log: GRHydro: add option to control action taken when con2prim failed From: Roland Haas <rhaas(a)caltech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: +8 -1 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2013-05-29 16:35:01 UTC (rev 528) +++ trunk/src/GRHydro_Con2Prim.F90 2013-05-29 16:35:04 UTC (rev 529) @@ -2288,6 +2288,7 @@ DECLARE_CCTK_PARAMETERS DECLARE_CCTK_ARGUMENTS + DECLARE_CCTK_FUNCTIONS integer :: i, j, k, nx, ny, nz character(len=300) warnline @@ -2401,7 +2402,13 @@ write(warnline,'(a32,4g16.7)') 'lapse, shift: ',& alp(i,j,k),betax(i,j,k),betay(i,j,k),betaz(i,j,k) call CCTK_WARN(1,warnline) - call CCTK_WARN(0,"Aborting.") + if (CCTK_EQUALS(GRHydro_c2p_failed_action, "terminate")) then + call CCTK_TerminateNext(cctkGH) + else if (CCTK_EQUALS(GRHydro_c2p_failed_action, "abort")) then + call CCTK_ERROR("Aborting.") + else + call CCTK_ERROR("Internal error, unknown action") + end if !$OMP END CRITICAL end if Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +6 -0 =================================================================== --- trunk/param.ccl 2013-05-29 16:35:01 UTC (rev 528) +++ trunk/param.ccl 2013-05-29 16:35:04 UTC (rev 529) @@ -397,6 +397,12 @@ 0:1 :: "Either 0 or 1" } 0 +keyword GRHydro_c2p_failed_action "what to do when we detect a c2p failure" +{ + "abort" :: "abort with error" + "terminate" :: "request termination" +} "abort" + REAL sqrtdet_thr "Threshold to apply cons rescalings deep inside the horizon" STEERABLE=ALWAYS { 1.0: :: "Larger values guarantees this sort of rescaling only deep inside the horizon"

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 528)
by rhaas＠tapir.caltech.edu 29 May '13

29 May '13

User: rhaas Date: 2013/05/29 11:35 AM Modified: /trunk/src/ GRHydro_MP5Reconstruct.F90 Log: GRHydro: MP5: Make epsilon parameter scale dependent on reconstructed field. From: Christian Reisswig <reisswig(a)scriwalker.(none)> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_MP5Reconstruct.F90 Delta lines: +5 -2 =================================================================== --- trunk/src/GRHydro_MP5Reconstruct.F90 2013-05-29 16:34:58 UTC (rev 527) +++ trunk/src/GRHydro_MP5Reconstruct.F90 2013-05-29 16:35:01 UTC (rev 528) @@ -48,7 +48,7 @@ logical, dimension(nx) :: excise logical :: normal_m5 - CCTK_REAL :: vl, vmp, djm1, dj, djp1, dm4jph, dm4jmh, vul, vav, vmd, vlc, vmin, vmax + CCTK_REAL :: vl, vmp, djm1, dj, djp1, dm4jph, dm4jmh, vul, vav, vmd, vlc, vmin, vmax, vnorm ! sign requires its arguments to be of identical KIND CCTK_REAL, parameter :: one = 1d0 @@ -89,6 +89,9 @@ end if if (normal_m5) then + + vnorm = sqrt(v(i-2)**2 + v(i-1)**2 + v(i)**2 + v(i+1)**2 + v(i+2)**2) + #define MINMOD(x,y) \ 0.5d0*(sign(one,x) + sign(one,y)) * min(abs(x), abs(y)) @@ -98,7 +101,7 @@ #define MP5(am2, am1, a, ap1, ap2, arecon) \ vl = (2.0d0*am2 - 13.0d0*am1 + 47.0d0*a + 27.0d0*ap1 - 3.0d0*ap2)/60.0d0 &&\ vmp = a + MINMOD( ap1-a, mp5_alpha*(a-am1) ) &&\ - if ((vl-a)*(vl-vmp) .le. mp5_eps) then &&\ + if ((vl-a)*(vl-vmp) .le. mp5_eps*vnorm) then &&\ arecon = vl &&\ else &&\ djm1 = am2 -2.0d0*am1 + a &&\

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 527)
by rhaas＠tapir.caltech.edu 29 May '13

29 May '13

User: rhaas Date: 2013/05/29 11:34 AM Modified: /trunk/src/ GRHydro_Eigenproblem.F90 Log: GRHydro: add error checking after call for cs2 in Eigenproblem From: Christian David Ott <cott(a)zwicky-b.(none)> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Eigenproblem.F90 Delta lines: +11 -0 =================================================================== --- trunk/src/GRHydro_Eigenproblem.F90 2013-05-29 16:34:50 UTC (rev 526) +++ trunk/src/GRHydro_Eigenproblem.F90 2013-05-29 16:34:58 UTC (rev 527) @@ -146,6 +146,8 @@ CCTK_REAL lam1,lam2,lam3,lamm,lamp,lamm_nobeta,lamp_nobeta CCTK_INT handle,ii,jj,kk + character(len=512) :: warnline + ! begin EOS Omni vars integer :: n,keytemp,anyerr,keyerr(1) real*8 :: xpress,xeps @@ -159,6 +161,15 @@ !!$ Set required fluid quantities call EOS_Omni_cs2(handle,keytemp,GRHydro_eos_rf_prec,n,& rho,eps,temp,ye,cs2,keyerr,anyerr) + if(anyerr.ne.0) then + !$OMP CRITICAL + call CCTK_WARN(1,"EOS ERROR in eigenvalues_hot") + write(warnline,"(A10,i5,A10,i5)") "keyerr: ", keyerr, "keytemp: ",keytemp + call CCTK_WARN(1,warnline) + write(warnline,"(1P10E15.6)") rho,eps,temp,ye,cs2 + call CCTK_WARN(0,warnline) + !$OMP END CRITICAL + endif vlowx = gxx*velx + gxy*vely + gxz*velz vlowy = gxy*velx + gyy*vely + gyz*velz

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 526)
by rhaas＠tapir.caltech.edu 29 May '13

29 May '13

User: rhaas Date: 2013/05/29 11:34 AM Modified: /trunk/src/ GRHydro_HLLE.F90 Log: GRHydro: Fixing index offset bug in tempplus(i,j,k) for x direction flux. From: Philipp Moesta <pmoesta(a)tapir.caltech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_HLLE.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_HLLE.F90 2013-05-27 22:32:27 UTC (rev 525) +++ trunk/src/GRHydro_HLLE.F90 2013-05-29 16:34:50 UTC (rev 526) @@ -377,7 +377,7 @@ rhoplus(i,j,k),& velxplus(i,j,k),velyplus(i,j,k),& velzplus(i,j,k),epsplus(i,j,k), & - tempplus(i+xoffset,j+yoffset,k+zoffset),& + tempplus(i,j,k),& y_e_plus(i,j,k),& w_lorentzplus(i,j,k),& lamplus,gxxh,gxyh,gxzh,gyyh,&

1 0

[barrywardell/EinsteinExact] 287fe7: Tidy up auto-generated documentation after generat...
by Barry Wardell 28 May '13

28 May '13

Branch: refs/heads/master Home: https://github.com/barrywardell/EinsteinExact Commit: 287fe7532aeb1d64008a838645e3cbb68da338f0 https://github.com/barrywardell/EinsteinExact/commit/287fe7532aeb1d64008a83… Author: Barry Wardell <barry.wardell(a)gmail.com> Date: 2013-05-28 (Tue, 28 May 2013) Changed paths: M m/Makefile A m/docfix.patch Log Message: ----------- Tidy up auto-generated documentation after generating it Some of the autogenerated documentation is ugly and doesn't fit on a standard page. We have to manually tidy this up. Previously this was completely manual, now instead we make it semi-automatic by applying the fixup changes as a patch. Commit: 27756b1393f8fd888b1ff14484d9891b0e903977 https://github.com/barrywardell/EinsteinExact/commit/27756b1393f8fd888b1ff1… Author: Barry Wardell <barry.wardell(a)gmail.com> Date: 2013-05-28 (Tue, 28 May 2013) Changed paths: M doc/spacetimes.tex Log Message: ----------- Regenerate code (documentation) Compare: https://github.com/barrywardell/EinsteinExact/compare/9651db00ff16...27756b…

1 0

[barrywardell/EinsteinExact] 9651db: Regenerate code (Kranc commit 519f2cd182c4828e7442...
by Ian Hinder 28 May '13

28 May '13

Branch: refs/heads/master Home: https://github.com/barrywardell/EinsteinExact Commit: 9651db00ff160b2bef996b92a878e966e510fb4e https://github.com/barrywardell/EinsteinExact/commit/9651db00ff160b2bef996b… Author: Ian Hinder <ian.hinder(a)aei.mpg.de> Date: 2013-05-28 (Tue, 28 May 2013) Changed paths: M GaugeWave/schedule.ccl M KerrSchild/schedule.ccl M Minkowski/schedule.ccl M ModifiedSchwarzschildBL/schedule.ccl M ShiftedGaugeWave/schedule.ccl M Vaidya2/schedule.ccl M doc/spacetimes.tex Log Message: ----------- Regenerate code (Kranc commit 519f2cd182c4828e74429e883139fbb8b1a11d3c)

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits