Commits May 2013

commits@lists.einsteintoolkit.org

9 participants
106 discussions

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 535)
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 /trunk/test/ GRHydro_test_shock_mp5.par /trunk/test/GRHydro_test_shock_mp5/ dens_3D_diagonal.xg, dens_x_[7][7].xg, dens_y_[7][7].xg, dens_z_[7][7].xg, eps_3D_diagonal.xg, eps_x_[7][7].xg, eps_y_[7][7].xg, eps_z_[7][7].xg, press_3D_diagonal.xg, press_x_[7][7].xg, press_y_[7][7].xg, press_z_[7][7].xg, rho_3D_diagonal.xg, rho_x_[7][7].xg, rho_y_[7][7].xg, rho_z_[7][7].xg, scon[0]_3D_diagonal.xg, scon[0]_x_[7][7].xg, scon[0]_y_[7][7].xg, scon[0]_z_[7][7].xg, scon[1]_3D_diagonal.xg, scon[1]_x_[7][7].xg, scon[1]_y_[7][7].xg, scon[1]_z_[7][7].xg, scon[2]_3D_diagonal.xg, scon[2]_x_[7][7].xg, scon[2]_y_[7][7].xg, scon[2]_z_[7][7].xg, tau_3D_diagonal.xg, tau_x_[7][7].xg, tau_y_[7][7].xg, tau_z_[7][7].xg, vel[0]_3D_diagonal.xg, vel[0]_x_[7][7].xg, vel[0]_y_[7][7].xg, vel[0]_z_[7][7].xg, vel[1]_3D_diagonal.xg, vel[1]_x_[7][7].xg, vel[1]_y_[7][7].xg, vel[1]_z_[7][7].xg, vel[2]_3D_diagonal.xg, vel[2]_x_[7][7].xg, vel[2]_y_[7][7].xg, vel[2]_z_[7][7].xg, w_lorentz_3D_diagonal.xg, w_lorentz_x_[7][7].xg, w_lorentz_y_[7][7].xg, w_lorentz_z_[7][7].xg Log: GRHydro: add check for stencil width for MP5 regenerate test that was using wrong stencil before From: Roland Haas <rhaas(a)tapir.caltech.edu> File Changes: Directory: /trunk/test/GRHydro_test_shock_mp5/ ============================================== File [modified]: dens_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: dens_x_[7][7].xg Delta lines: <Unavailable> File [modified]: dens_y_[7][7].xg Delta lines: <Unavailable> File [modified]: dens_z_[7][7].xg Delta lines: <Unavailable> File [modified]: eps_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: eps_x_[7][7].xg Delta lines: <Unavailable> File [modified]: eps_y_[7][7].xg Delta lines: <Unavailable> File [modified]: eps_z_[7][7].xg Delta lines: <Unavailable> File [modified]: press_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: press_x_[7][7].xg Delta lines: <Unavailable> File [modified]: press_y_[7][7].xg Delta lines: <Unavailable> File [modified]: press_z_[7][7].xg Delta lines: <Unavailable> File [modified]: rho_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: rho_x_[7][7].xg Delta lines: <Unavailable> File [modified]: rho_y_[7][7].xg Delta lines: <Unavailable> File [modified]: rho_z_[7][7].xg Delta lines: <Unavailable> File [modified]: scon[0]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: scon[0]_x_[7][7].xg Delta lines: <Unavailable> File [modified]: scon[0]_y_[7][7].xg Delta lines: <Unavailable> File [modified]: scon[0]_z_[7][7].xg Delta lines: <Unavailable> File [modified]: scon[1]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: scon[1]_x_[7][7].xg Delta lines: <Unavailable> File [modified]: scon[1]_y_[7][7].xg Delta lines: <Unavailable> File [modified]: scon[1]_z_[7][7].xg Delta lines: <Unavailable> File [modified]: scon[2]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: scon[2]_x_[7][7].xg Delta lines: <Unavailable> File [modified]: scon[2]_y_[7][7].xg Delta lines: <Unavailable> File [modified]: scon[2]_z_[7][7].xg Delta lines: <Unavailable> File [modified]: tau_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: tau_x_[7][7].xg Delta lines: <Unavailable> File [modified]: tau_y_[7][7].xg Delta lines: <Unavailable> File [modified]: tau_z_[7][7].xg Delta lines: <Unavailable> File [modified]: vel[0]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: vel[0]_x_[7][7].xg Delta lines: <Unavailable> File [modified]: vel[0]_y_[7][7].xg Delta lines: <Unavailable> File [modified]: vel[0]_z_[7][7].xg Delta lines: <Unavailable> File [modified]: vel[1]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: vel[1]_x_[7][7].xg Delta lines: <Unavailable> File [modified]: vel[1]_y_[7][7].xg Delta lines: <Unavailable> File [modified]: vel[1]_z_[7][7].xg Delta lines: <Unavailable> File [modified]: vel[2]_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: vel[2]_x_[7][7].xg Delta lines: <Unavailable> File [modified]: vel[2]_y_[7][7].xg Delta lines: <Unavailable> File [modified]: vel[2]_z_[7][7].xg Delta lines: <Unavailable> File [modified]: w_lorentz_3D_diagonal.xg Delta lines: <Unavailable> File [modified]: w_lorentz_x_[7][7].xg Delta lines: <Unavailable> File [modified]: w_lorentz_y_[7][7].xg Delta lines: <Unavailable> File [modified]: w_lorentz_z_[7][7].xg Delta lines: <Unavailable> Directory: /trunk/src/ ====================== File [modified]: GRHydro_ParamCheck.F90 Delta lines: <Unavailable> Directory: /trunk/test/ ======================= File [modified]: GRHydro_test_shock_mp5.par Delta lines: <Unavailable>

1 0

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

29 May '13

User: rhaas Date: 2013/05/29 11:35 AM Modified: /trunk/src/ GRHydro_Boundaries.F90 Log: GRHydro: avoid implict save for ierr in boundary registration From: Roland Haas <rhaas(a)tapir.caltech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Boundaries.F90 Delta lines: +11 -5 =================================================================== --- trunk/src/GRHydro_Boundaries.F90 2013-05-29 16:35:18 UTC (rev 533) +++ trunk/src/GRHydro_Boundaries.F90 2013-05-29 16:35:22 UTC (rev 534) @@ -40,9 +40,10 @@ DECLARE_CCTK_FUNCTIONS integer, dimension(3) :: sym - integer :: ierr = 0 + integer :: ierr CCTK_INT :: GRHydro_UseGeneralCoordinates, general_coordinates + ierr = 0 sym(1) = 1 sym(2) = 1 sym(3) = 1 @@ -173,12 +174,14 @@ DECLARE_CCTK_FUNCTIONS integer, dimension(3) :: sw - integer :: ierr = 0 + integer :: ierr CCTK_INT :: GRHydro_UseGeneralCoordinates, general_coordinates CCTK_INT, parameter :: faces=CCTK_ALL_FACES CCTK_INT, parameter :: ione=1 + ierr = 0 + sw = GRHydro_stencil general_coordinates = GRHydro_UseGeneralCoordinates(cctkGH) @@ -375,12 +378,13 @@ DECLARE_CCTK_FUNCTIONS integer, dimension(3) :: sw - integer :: ierr = 0 + integer :: ierr CCTK_INT :: GRHydro_UseGeneralCoordinates, general_coordinates CCTK_INT, parameter :: faces=CCTK_ALL_FACES CCTK_INT, parameter :: ione=1 + ierr = 0 sw = GRHydro_stencil general_coordinates = GRHydro_UseGeneralCoordinates(cctkGH) @@ -513,12 +517,13 @@ DECLARE_CCTK_FUNCTIONS integer, dimension(3) :: sw - integer :: ierr = 0 + integer :: ierr CCTK_INT :: GRHydro_UseGeneralCoordinates, general_coordinates CCTK_INT, parameter :: faces=CCTK_ALL_FACES CCTK_INT, parameter :: ione=1 + ierr = 0 sw = GRHydro_stencil general_coordinates = GRHydro_UseGeneralCoordinates(cctkGH) @@ -648,11 +653,12 @@ DECLARE_CCTK_FUNCTIONS integer, dimension(3) :: sw - integer :: ierr = 0 + integer :: ierr CCTK_INT, parameter :: faces=CCTK_ALL_FACES CCTK_INT, parameter :: ione=1 + ierr = 0 sw = GRHydro_stencil

1 0

[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

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits May 2013