Commits February 2013

commits@lists.einsteintoolkit.org

6 participants
52 discussions

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

11 Feb '13

User: rhaas Date: 2013/02/11 10:35 AM Modified: /trunk/ param.ccl /trunk/src/ EOS_Omni_Module.F90, EOS_Omni_MultiVarCalls.F90, EOS_Omni_SingleVarCalls.F90, EOS_Omni_Startup.F90 Log: remove poly_gamma_ini EOS_Omni used this to compute (initially) conversion factors from code (assumed to e G=c=M_sun=1) to cgs units. It then goes back and forth between these units even for polytropic and gamma law eos. This is very confusing since poly_gamma_ini since poly_gamma_ini must mast poly_gamma and gl_gamma for a "usual" run with a Gamma law eos. File Changes: Directory: /trunk/src/ ====================== File [modified]: EOS_Omni_Module.F90 Delta lines: +1 -5 =================================================================== --- trunk/src/EOS_Omni_Module.F90 2012-12-29 21:37:30 UTC (rev 72) +++ trunk/src/EOS_Omni_Module.F90 2013-02-11 16:35:39 UTC (rev 73) @@ -26,10 +26,6 @@ ! These values are set by EOS_Omni_Startup - real*8 :: poly_k_cgs = 0.0d0 - real*8 :: gl_k_cgs = 0.0d0 + real*8 :: hybrid_k2 = 0.0d0 - real*8 :: hybrid_k1_cgs = 0.0d0 - real*8 :: hybrid_k2_cgs = 0.0d0 - end module EOS_Omni_Module File [modified]: EOS_Omni_MultiVarCalls.F90 Delta lines: +32 -32 =================================================================== --- trunk/src/EOS_Omni_MultiVarCalls.F90 2012-12-29 21:37:30 UTC (rev 72) +++ trunk/src/EOS_Omni_MultiVarCalls.F90 2013-02-11 16:35:39 UTC (rev 73) @@ -104,7 +104,7 @@ integer :: i character(256) :: warnstring real*8 :: hybrid_local_gamma - real*8 :: hybrid_local_k_cgs + real*8 :: hybrid_local_k real*8 :: hybrid_dp_poly,hybrid_dp_th1,hybrid_dp_th2 ! temporary vars for nuc_eos real*8 :: xrho,xye,xtemp,xenr,xent @@ -120,23 +120,23 @@ ! polytropic EOS if(keytemp.eq.1) then do i=1,npoints - eps(i) = press_gf * poly_k_cgs * & - (rho(i)*inv_rho_gf)**(poly_gamma) / & + eps(i) = poly_k * & + rho(i)**(poly_gamma) / & (poly_gamma - 1.0d0) / rho(i) enddo endif do i=1,npoints - dpdrhoe(i) = press_gf * poly_k_cgs * & - poly_gamma * inv_rho_gf * & - (rho(i)*inv_rho_gf) ** (poly_gamma - 1.d0) + dpdrhoe(i) = poly_k * & + poly_gamma * & + rho(i) ** (poly_gamma - 1.d0) dpderho(i) = 0.0d0 enddo case (2) ! gamma-law EOS if(keytemp.eq.1) then do i=1,npoints - eps(i) = press_gf * gl_k_cgs * & - (rho(i)*inv_rho_gf)**(gl_gamma) / & + eps(i) = gl_k * & + rho(i)**(gl_gamma) / & (gl_gamma - 1.0d0) / rho(i) enddo endif @@ -151,23 +151,23 @@ do i=1,npoints if(rho(i).gt.hybrid_rho_nuc) then hybrid_local_gamma = hybrid_gamma2 - hybrid_local_k_cgs = hybrid_k2_cgs + hybrid_local_k = hybrid_k2 else hybrid_local_gamma = hybrid_gamma1 - hybrid_local_k_cgs = hybrid_k1_cgs + hybrid_local_k = hybrid_k1 endif - hybrid_dp_poly = hybrid_local_gamma * press_gf * & - hybrid_local_k_cgs * rho(i)**(hybrid_local_gamma - 1.0d0) * & - inv_rho_gf**hybrid_local_gamma + hybrid_dp_poly = hybrid_local_gamma * & + hybrid_local_k * rho(i)**(hybrid_local_gamma - 1.0d0) + - hybrid_dp_th1 = - hybrid_local_gamma * press_gf * hybrid_local_k_cgs * & + hybrid_dp_th1 = - hybrid_local_gamma * hybrid_local_k * & (hybrid_gamma_th - 1.d0) / (hybrid_local_gamma - 1.d0) * & - rho(i)**(hybrid_local_gamma - 1.d0) * inv_rho_gf**hybrid_local_gamma + rho(i)**(hybrid_local_gamma - 1.d0) hybrid_dp_th2 = (hybrid_gamma_th - 1.d0) * eps(i) & - (hybrid_gamma_th - 1.d0) * (hybrid_local_gamma - hybrid_gamma1) / & (hybrid_gamma1 - 1.d0) / (hybrid_gamma2 - 1.d0) * & - press_gf * hybrid_k1_cgs * inv_rho_gf**hybrid_gamma1 * & + hybrid_k1 * & hybrid_rho_nuc**(hybrid_gamma1 - 1.d0) dpdrhoe(i) = hybrid_dp_poly + hybrid_dp_th1 + hybrid_dp_th2 @@ -224,7 +224,7 @@ integer :: i character(256) :: warnstring real*8 :: hybrid_local_gamma - real*8 :: hybrid_local_k_cgs + real*8 :: hybrid_local_k real*8 :: hybrid_dp_poly,hybrid_dp_th1,hybrid_dp_th2 ! temporary vars for nuc_eos real*8 :: xrho,xye,xtemp,xenr,xent @@ -240,23 +240,23 @@ ! polytropic EOS if(keytemp.eq.1) then do i=1,npoints - eps(i) = press_gf * poly_k_cgs * & - (rho(i)*inv_rho_gf)**(poly_gamma) / & + eps(i) = poly_k * & + rho(i)**(poly_gamma) / & (poly_gamma - 1.0d0) / rho(i) enddo endif do i=1,npoints depsdpress(i) = 1.0d0/(poly_gamma - 1.0d0)/rho(i) - depsdrho(i) = depsdpress(i) * press_gf * poly_k_cgs * & - poly_gamma * inv_rho_gf * & - (rho(i)*inv_rho_gf) ** (poly_gamma - 1.d0) + depsdrho(i) = depsdpress(i) * poly_k * & + poly_gamma * & + rho(i) ** (poly_gamma - 1.d0) enddo case (2) ! gamma-law EOS if(keytemp.eq.1) then do i=1,npoints - eps(i) = press_gf * gl_k_cgs * & - (rho(i)*inv_rho_gf)**(gl_gamma) / & + eps(i) = gl_k * & + rho(i)**(gl_gamma) / & (gl_gamma - 1.0d0) / rho(i) enddo endif @@ -270,23 +270,23 @@ do i=1,npoints if(rho(i).gt.hybrid_rho_nuc) then hybrid_local_gamma = hybrid_gamma2 - hybrid_local_k_cgs = hybrid_k2_cgs + hybrid_local_k = hybrid_k2 else hybrid_local_gamma = hybrid_gamma1 - hybrid_local_k_cgs = hybrid_k1_cgs + hybrid_local_k = hybrid_k1 endif - hybrid_dp_poly = hybrid_local_gamma * press_gf * & - hybrid_local_k_cgs * rho(i)**(hybrid_local_gamma - 1.0d0) * & - inv_rho_gf**hybrid_local_gamma + hybrid_dp_poly = hybrid_local_gamma * & + hybrid_local_k * rho(i)**(hybrid_local_gamma - 1.0d0) + - hybrid_dp_th1 = - hybrid_local_gamma * press_gf * hybrid_local_k_cgs * & + hybrid_dp_th1 = - hybrid_local_gamma * hybrid_local_k * & (hybrid_gamma_th - 1.d0) / (hybrid_local_gamma - 1.d0) * & - rho(i)**(hybrid_local_gamma - 1.d0) * inv_rho_gf**hybrid_local_gamma + rho(i)**(hybrid_local_gamma - 1.d0) hybrid_dp_th2 = (hybrid_gamma_th - 1.d0) * eps(i) & - (hybrid_gamma_th - 1.d0) * (hybrid_local_gamma - hybrid_gamma1) / & (hybrid_gamma1 - 1.d0) / (hybrid_gamma2 - 1.d0) * & - press_gf * hybrid_k1_cgs * inv_rho_gf**hybrid_gamma1 * & + hybrid_k1 * & hybrid_rho_nuc**(hybrid_gamma1 - 1.d0) xdpdrhoe = hybrid_dp_poly + hybrid_dp_th1 + hybrid_dp_th2 File [modified]: EOS_Omni_SingleVarCalls.F90 Delta lines: +48 -48 =================================================================== --- trunk/src/EOS_Omni_SingleVarCalls.F90 2012-12-29 21:37:30 UTC (rev 72) +++ trunk/src/EOS_Omni_SingleVarCalls.F90 2013-02-11 16:35:39 UTC (rev 73) @@ -27,7 +27,7 @@ ! local vars integer :: i character(256) :: warnstring - real*8 :: hybrid_local_gamma, hybrid_local_k_cgs, & + real*8 :: hybrid_local_gamma, hybrid_local_k, & hybrid_p_poly, hybrid_p_th real*8,parameter :: zero = 0.0d0 ! temporary vars for nuc_eos @@ -43,21 +43,21 @@ ! polytropic EOS if(keytemp.eq.1) then do i=1,npoints - eps(i) = press_gf * poly_k_cgs * & - (rho(i)*inv_rho_gf)**(poly_gamma) / & + eps(i) = poly_k * & + rho(i)**(poly_gamma) / & (poly_gamma - 1.0d0) / rho(i) enddo endif do i=1,npoints - press(i) = press_gf * poly_k_cgs * & - (rho(i)*inv_rho_gf)**poly_gamma + press(i) = poly_k * & + rho(i)**poly_gamma enddo case (2) ! gamma-law EOS if(keytemp.eq.1) then do i=1,npoints - eps(i) = press_gf * gl_k_cgs * & - (rho(i)*inv_rho_gf)**(gl_gamma) / & + eps(i) = gl_k * & + rho(i)**(gl_gamma) / & (gl_gamma - 1.0d0) / rho(i) enddo endif @@ -70,19 +70,19 @@ do i=1,npoints if(rho(i).gt.hybrid_rho_nuc) then hybrid_local_gamma = hybrid_gamma2 - hybrid_local_k_cgs = hybrid_k2_cgs + hybrid_local_k = hybrid_k2 else hybrid_local_gamma = hybrid_gamma1 - hybrid_local_k_cgs = hybrid_k1_cgs + hybrid_local_k = hybrid_k1 endif - hybrid_p_poly = press_gf * hybrid_local_k_cgs * & - (rho(i) * inv_rho_gf)**hybrid_local_gamma - hybrid_p_th = - press_gf * hybrid_local_k_cgs * (hybrid_gamma_th - 1.d0) / & - (hybrid_local_gamma - 1.0d0) * (rho(i) * inv_rho_gf)**hybrid_local_gamma + & + hybrid_p_poly = hybrid_local_k * & + rho(i)**hybrid_local_gamma + hybrid_p_th = - hybrid_local_k * (hybrid_gamma_th - 1.d0) / & + (hybrid_local_gamma - 1.0d0) * rho(i)**hybrid_local_gamma + & (hybrid_gamma_th - 1.0d0) * rho(i) * eps(i) - & (hybrid_gamma_th - 1.d0) * (hybrid_local_gamma - hybrid_gamma1) / & (hybrid_gamma1 - 1.d0) / (hybrid_gamma2 - 1.d0) * & - press_gf * hybrid_k1_cgs * inv_rho_gf**hybrid_gamma1 * & + hybrid_k1 * & hybrid_rho_nuc**(hybrid_gamma1 - 1.d0) * rho(i) hybrid_p_th = max(zero, hybrid_p_th) press(i) = hybrid_p_poly + hybrid_p_th @@ -140,7 +140,7 @@ ! local vars integer :: i character(256) :: warnstring - real*8 :: hybrid_local_gamma, hybrid_local_k_cgs, & + real*8 :: hybrid_local_gamma, hybrid_local_k, & hybrid_dp_poly, hybrid_dp_th1, hybrid_dp_th2 real*8,parameter :: zero = 0.0d0 ! temporary vars for nuc_eos @@ -156,22 +156,22 @@ ! polytropic EOS if(keytemp.eq.1) then do i=1,npoints - eps(i) = press_gf * poly_k_cgs * & - (rho(i)*inv_rho_gf)**(poly_gamma) / & + eps(i) = poly_k * & + rho(i)**(poly_gamma) / & (poly_gamma - 1.0d0) / rho(i) enddo endif do i=1,npoints - dpdrhoe(i) = press_gf * poly_k_cgs * & - poly_gamma * inv_rho_gf * & - (rho(i)*inv_rho_gf) ** (poly_gamma - 1.d0) + dpdrhoe(i) = poly_k * & + poly_gamma * & + rho(i) ** (poly_gamma - 1.d0) enddo case (2) ! gamma-law EOS if(keytemp.eq.1) then do i=1,npoints - eps(i) = press_gf * gl_k_cgs * & - (rho(i)*inv_rho_gf)**(gl_gamma) / & + eps(i) = gl_k * & + rho(i)**(gl_gamma) / & (gl_gamma - 1.0d0) / rho(i) enddo endif @@ -184,23 +184,23 @@ do i=1,npoints if(rho(i).gt.hybrid_rho_nuc) then hybrid_local_gamma = hybrid_gamma2 - hybrid_local_k_cgs = hybrid_k2_cgs + hybrid_local_k = hybrid_k2 else hybrid_local_gamma = hybrid_gamma1 - hybrid_local_k_cgs = hybrid_k1_cgs + hybrid_local_k = hybrid_k1 endif - hybrid_dp_poly = hybrid_local_gamma * press_gf * & - hybrid_local_k_cgs * rho(i)**(hybrid_local_gamma - 1.0d0) * & - inv_rho_gf**hybrid_local_gamma + hybrid_dp_poly = hybrid_local_gamma * & + hybrid_local_k * rho(i)**(hybrid_local_gamma - 1.0d0) + - hybrid_dp_th1 = - hybrid_local_gamma * press_gf * hybrid_local_k_cgs * & + hybrid_dp_th1 = - hybrid_local_gamma * hybrid_local_k * & (hybrid_gamma_th - 1.d0) / (hybrid_local_gamma - 1.d0) * & - rho(i)**(hybrid_local_gamma - 1.d0) * inv_rho_gf**hybrid_local_gamma + rho(i)**(hybrid_local_gamma - 1.d0) hybrid_dp_th2 = (hybrid_gamma_th - 1.d0) * eps(i) & - (hybrid_gamma_th - 1.d0) * (hybrid_local_gamma - hybrid_gamma1) / & (hybrid_gamma1 - 1.d0) / (hybrid_gamma2 - 1.d0) * & - press_gf * hybrid_k1_cgs * inv_rho_gf**hybrid_gamma1 * & + hybrid_k1 * & hybrid_rho_nuc**(hybrid_gamma1 - 1.d0) dpdrhoe(i) = hybrid_dp_poly + max(0.0d0,hybrid_dp_th1 + hybrid_dp_th2) @@ -269,8 +269,8 @@ ! polytropic EOS if(keytemp.eq.1) then do i=1,npoints - eps(i) = press_gf * poly_k_cgs * & - (rho(i)*inv_rho_gf)**(poly_gamma) / & + eps(i) = poly_k * & + rho(i)**(poly_gamma) / & (poly_gamma - 1.0d0) / rho(i) enddo endif @@ -281,8 +281,8 @@ ! gamma-law EOS if(keytemp.eq.1) then do i=1,npoints - eps(i) = press_gf * gl_k_cgs * & - (rho(i)*inv_rho_gf)**(gl_gamma) / & + eps(i) = gl_k * & + rho(i)**(gl_gamma) / & (gl_gamma - 1.0d0) / rho(i) enddo endif @@ -348,7 +348,7 @@ integer :: i character(256) :: warnstring real*8 :: xpress,xdpdrhoe,xdpderho - real*8 :: hybrid_local_gamma, hybrid_local_k_cgs, & + real*8 :: hybrid_local_gamma, hybrid_local_k, & hybrid_p_poly, hybrid_p_th real*8,parameter :: zero = 0.0d0 ! temporary vars for nuc_eos @@ -364,14 +364,14 @@ ! polytropic EOS if(keytemp.eq.1) then do i=1,npoints - eps(i) = press_gf * poly_k_cgs * & - (rho(i)*inv_rho_gf)**(poly_gamma) / & + eps(i) = poly_k * & + rho(i)**(poly_gamma) / & (poly_gamma - 1.0d0) / rho(i) enddo endif do i=1,npoints - xpress = press_gf*poly_k_cgs * & - (rho(i)*inv_rho_gf)**(poly_gamma) + xpress = poly_k * & + rho(i)**(poly_gamma) cs2(i) = poly_gamma * xpress / rho(i) / & (1 + eps(i) + xpress/rho(i)) enddo @@ -379,8 +379,8 @@ ! gamma-law EOS if(keytemp.eq.1) then do i=1,npoints - eps(i) = press_gf * gl_k_cgs * & - (rho(i)*inv_rho_gf)**(gl_gamma) / & + eps(i) = gl_k * & + rho(i)**(gl_gamma) / & (gl_gamma - 1.0d0) / rho(i) enddo endif @@ -396,20 +396,20 @@ do i=1,npoints if(rho(i).gt.hybrid_rho_nuc) then hybrid_local_gamma = hybrid_gamma2 - hybrid_local_k_cgs = hybrid_k2_cgs + hybrid_local_k = hybrid_k2 else hybrid_local_gamma = hybrid_gamma1 - hybrid_local_k_cgs = hybrid_k1_cgs + hybrid_local_k = hybrid_k1 endif ! first calculate the pressure - hybrid_p_poly = press_gf * hybrid_local_k_cgs * & - (rho(i) * inv_rho_gf)**hybrid_local_gamma - hybrid_p_th = - press_gf * hybrid_local_k_cgs * (hybrid_gamma_th - 1.d0) / & - (hybrid_local_gamma - 1.0d0) * (rho(i) * inv_rho_gf)**hybrid_local_gamma + & + hybrid_p_poly = hybrid_local_k * & + rho(i)**hybrid_local_gamma + hybrid_p_th = - hybrid_local_k * (hybrid_gamma_th - 1.d0) / & + (hybrid_local_gamma - 1.0d0) * rho(i)**hybrid_local_gamma + & (hybrid_gamma_th - 1.0d0) * rho(i) * eps(i) - & (hybrid_gamma_th - 1.d0) * (hybrid_local_gamma - hybrid_gamma1) / & (hybrid_gamma1 - 1.d0) / (hybrid_gamma2 - 1.d0) * & - press_gf * hybrid_k1_cgs * inv_rho_gf**hybrid_gamma1 * & + hybrid_k1 * & hybrid_rho_nuc**(hybrid_gamma1 - 1.d0) * rho(i) hybrid_p_th = max(zero, hybrid_p_th) xpress = hybrid_p_poly + hybrid_p_th File [modified]: EOS_Omni_Startup.F90 Delta lines: +2 -8 =================================================================== --- trunk/src/EOS_Omni_Startup.F90 2012-12-29 21:37:30 UTC (rev 72) +++ trunk/src/EOS_Omni_Startup.F90 2013-02-11 16:35:39 UTC (rev 73) @@ -11,13 +11,7 @@ DECLARE_CCTK_PARAMETERS DECLARE_CCTK_ARGUMENTS - poly_k_cgs = poly_k * rho_gf**poly_gamma_ini / press_gf + hybrid_k2 = hybrid_k1 * & + hybrid_rho_nuc**(hybrid_gamma1-hybrid_gamma2) - gl_k_cgs = gl_k * rho_gf**poly_gamma_ini / press_gf - - hybrid_k1_cgs = hybrid_k1 * rho_gf**poly_gamma_ini / press_gf - - hybrid_k2_cgs = hybrid_k1_cgs * & - (hybrid_rho_nuc * inv_rho_gf)**(hybrid_gamma1-hybrid_gamma2) - end subroutine EOS_Omni_Startup Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +0 -5 =================================================================== --- trunk/param.ccl 2012-12-29 21:37:30 UTC (rev 72) +++ trunk/param.ccl 2013-02-11 16:35:39 UTC (rev 73) @@ -9,11 +9,6 @@ : :: "" } 2.0 -REAL poly_gamma_ini "Initial Adiabatic Index for poly EOS" STEERABLE=RECOVER -{ - : :: "" -} 2.0 - REAL poly_k "Polytropic constant in c=G=Msun=1" STEERABLE=RECOVER { : :: ""

1 0

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

11 Feb '13

User: rhaas Date: 2013/02/11 10:30 AM Added: /trunk/test/ test_one_hybrid.par /trunk/test/test_one_hybrid/ dens_3D_diagonal.xg, dens_maximum.xg, dens_minimum.xg, dens_norm1.xg, dens_norm2.xg, dens_x_[16][16].xg, dens_y_[16][16].xg, dens_z_[16][16].xg, eps_3D_diagonal.xg, eps_maximum.xg, eps_minimum.xg, eps_norm1.xg, eps_norm2.xg, eps_x_[16][16].xg, eps_y_[16][16].xg, eps_z_[16][16].xg, press_3D_diagonal.xg, press_maximum.xg, press_minimum.xg, press_norm1.xg, press_norm2.xg, press_x_[16][16].xg, press_y_[16][16].xg, press_z_[16][16].xg, rho_3D_diagonal.xg, rho_maximum.xg, rho_minimum.xg, rho_norm1.xg, rho_norm2.xg, rho_x_[16][16].xg, rho_y_[16][16].xg, rho_z_[16][16].xg, scon[0]_3D_diagonal.xg, scon[0]_maximum.xg, scon[0]_minimum.xg, scon[0]_norm1.xg, scon[0]_norm2.xg, scon[0]_x_[16][16].xg, scon[0]_y_[16][16].xg, scon[0]_z_[16][16].xg, scon[1]_3D_diagonal.xg, scon[1]_maximum.xg, scon[1]_minimum.xg, scon[1]_norm1.xg, scon[1]_norm2.xg, scon[1]_x_[16][16].xg, scon[1]_y_[16][16].xg, scon[1]_z_[16][16].xg, scon[2]_3D_diagonal.xg, scon[2]_maximum.xg, scon[2]_minimum.xg, scon[2]_norm1.xg, scon[2]_norm2.xg, scon[2]_x_[16][16].xg, scon[2]_y_[16][16].xg, scon[2]_z_[16][16].xg, tau_3D_diagonal.xg, tau_maximum.xg, tau_minimum.xg, tau_norm1.xg, tau_norm2.xg, tau_x_[16][16].xg, tau_y_[16][16].xg, tau_z_[16][16].xg, test_one_hybrid.par, vel[0]_3D_diagonal.xg, vel[0]_maximum.xg, vel[0]_minimum.xg, vel[0]_norm1.xg, vel[0]_norm2.xg, vel[0]_x_[16][16].xg, vel[0]_y_[16][16].xg, vel[0]_z_[16][16].xg, vel[1]_3D_diagonal.xg, vel[1]_maximum.xg, vel[1]_minimum.xg, vel[1]_norm1.xg, vel[1]_norm2.xg, vel[1]_x_[16][16].xg, vel[1]_y_[16][16].xg, vel[1]_z_[16][16].xg, vel[2]_3D_diagonal.xg, vel[2]_maximum.xg, vel[2]_minimum.xg, vel[2]_norm1.xg, vel[2]_norm2.xg, vel[2]_x_[16][16].xg, vel[2]_y_[16][16].xg, vel[2]_z_[16][16].xg, w_lorentz_3D_diagonal.xg, w_lorentz_maximum.xg, w_lorentz_minimum.xg, w_lorentz_norm1.xg, w_lorentz_norm2.xg, w_lorentz_x_[16][16].xg, w_lorentz_y_[16][16].xg, w_lorentz_z_[16][16].xg Log: GRHydro: add test case for hybrid EOS From: Roland Haas <rhaas(a)tapir.caltech.edu> File Changes: Directory: /trunk/test/ ======================= File [added]: test_one_hybrid.par Delta lines: <Unavailable> Directory: /trunk/test/test_one_hybrid/ ======================================= File [added]: dens_3D_diagonal.xg Delta lines: <Unavailable> File [added]: dens_maximum.xg Delta lines: <Unavailable> File [added]: dens_minimum.xg Delta lines: <Unavailable> File [added]: dens_norm1.xg Delta lines: <Unavailable> File [added]: dens_norm2.xg Delta lines: <Unavailable> File [added]: dens_x_[16][16].xg Delta lines: <Unavailable> File [added]: dens_y_[16][16].xg Delta lines: <Unavailable> File [added]: dens_z_[16][16].xg Delta lines: <Unavailable> File [added]: eps_3D_diagonal.xg Delta lines: <Unavailable> File [added]: eps_maximum.xg Delta lines: <Unavailable> File [added]: eps_minimum.xg Delta lines: <Unavailable> File [added]: eps_norm1.xg Delta lines: <Unavailable> File [added]: eps_norm2.xg Delta lines: <Unavailable> File [added]: eps_x_[16][16].xg Delta lines: <Unavailable> File [added]: eps_y_[16][16].xg Delta lines: <Unavailable> File [added]: eps_z_[16][16].xg Delta lines: <Unavailable> File [added]: press_3D_diagonal.xg Delta lines: <Unavailable> File [added]: press_maximum.xg Delta lines: <Unavailable> File [added]: press_minimum.xg Delta lines: <Unavailable> File [added]: press_norm1.xg Delta lines: <Unavailable> File [added]: press_norm2.xg Delta lines: <Unavailable> File [added]: press_x_[16][16].xg Delta lines: <Unavailable> File [added]: press_y_[16][16].xg Delta lines: <Unavailable> File [added]: press_z_[16][16].xg Delta lines: <Unavailable> File [added]: rho_3D_diagonal.xg Delta lines: <Unavailable> File [added]: rho_maximum.xg Delta lines: <Unavailable> File [added]: rho_minimum.xg Delta lines: <Unavailable> File [added]: rho_norm1.xg Delta lines: <Unavailable> File [added]: rho_norm2.xg Delta lines: <Unavailable> File [added]: rho_x_[16][16].xg Delta lines: <Unavailable> File [added]: rho_y_[16][16].xg Delta lines: <Unavailable> File [added]: rho_z_[16][16].xg Delta lines: <Unavailable> File [added]: scon[0]_3D_diagonal.xg Delta lines: <Unavailable> File [added]: scon[0]_maximum.xg Delta lines: <Unavailable> File [added]: scon[0]_minimum.xg Delta lines: <Unavailable> File [added]: scon[0]_norm1.xg Delta lines: <Unavailable> File [added]: scon[0]_norm2.xg Delta lines: <Unavailable> File [added]: scon[0]_x_[16][16].xg Delta lines: <Unavailable> File [added]: scon[0]_y_[16][16].xg Delta lines: <Unavailable> File [added]: scon[0]_z_[16][16].xg Delta lines: <Unavailable> File [added]: scon[1]_3D_diagonal.xg Delta lines: <Unavailable> File [added]: scon[1]_maximum.xg Delta lines: <Unavailable> File [added]: scon[1]_minimum.xg Delta lines: <Unavailable> File [added]: scon[1]_norm1.xg Delta lines: <Unavailable> File [added]: scon[1]_norm2.xg Delta lines: <Unavailable> File [added]: scon[1]_x_[16][16].xg Delta lines: <Unavailable> File [added]: scon[1]_y_[16][16].xg Delta lines: <Unavailable> File [added]: scon[1]_z_[16][16].xg Delta lines: <Unavailable> File [added]: scon[2]_3D_diagonal.xg Delta lines: <Unavailable> File [added]: scon[2]_maximum.xg Delta lines: <Unavailable> File [added]: scon[2]_minimum.xg Delta lines: <Unavailable> File [added]: scon[2]_norm1.xg Delta lines: <Unavailable> File [added]: scon[2]_norm2.xg Delta lines: <Unavailable> File [added]: scon[2]_x_[16][16].xg Delta lines: <Unavailable> File [added]: scon[2]_y_[16][16].xg Delta lines: <Unavailable> File [added]: scon[2]_z_[16][16].xg Delta lines: <Unavailable> File [added]: tau_3D_diagonal.xg Delta lines: <Unavailable> File [added]: tau_maximum.xg Delta lines: <Unavailable> File [added]: tau_minimum.xg Delta lines: <Unavailable> File [added]: tau_norm1.xg Delta lines: <Unavailable> File [added]: tau_norm2.xg Delta lines: <Unavailable> File [added]: tau_x_[16][16].xg Delta lines: <Unavailable> File [added]: tau_y_[16][16].xg Delta lines: <Unavailable> File [added]: tau_z_[16][16].xg Delta lines: <Unavailable> File [added]: test_one_hybrid.par Delta lines: <Unavailable> File [added]: vel[0]_3D_diagonal.xg Delta lines: <Unavailable> File [added]: vel[0]_maximum.xg Delta lines: <Unavailable> File [added]: vel[0]_minimum.xg Delta lines: <Unavailable> File [added]: vel[0]_norm1.xg Delta lines: <Unavailable> File [added]: vel[0]_norm2.xg Delta lines: <Unavailable> File [added]: vel[0]_x_[16][16].xg Delta lines: <Unavailable> File [added]: vel[0]_y_[16][16].xg Delta lines: <Unavailable> File [added]: vel[0]_z_[16][16].xg Delta lines: <Unavailable> File [added]: vel[1]_3D_diagonal.xg Delta lines: <Unavailable> File [added]: vel[1]_maximum.xg Delta lines: <Unavailable> File [added]: vel[1]_minimum.xg Delta lines: <Unavailable> File [added]: vel[1]_norm1.xg Delta lines: <Unavailable> File [added]: vel[1]_norm2.xg Delta lines: <Unavailable> File [added]: vel[1]_x_[16][16].xg Delta lines: <Unavailable> File [added]: vel[1]_y_[16][16].xg Delta lines: <Unavailable> File [added]: vel[1]_z_[16][16].xg Delta lines: <Unavailable> File [added]: vel[2]_3D_diagonal.xg Delta lines: <Unavailable> File [added]: vel[2]_maximum.xg Delta lines: <Unavailable> File [added]: vel[2]_minimum.xg Delta lines: <Unavailable> File [added]: vel[2]_norm1.xg Delta lines: <Unavailable> File [added]: vel[2]_norm2.xg Delta lines: <Unavailable> File [added]: vel[2]_x_[16][16].xg Delta lines: <Unavailable> File [added]: vel[2]_y_[16][16].xg Delta lines: <Unavailable> File [added]: vel[2]_z_[16][16].xg Delta lines: <Unavailable> File [added]: w_lorentz_3D_diagonal.xg Delta lines: <Unavailable> File [added]: w_lorentz_maximum.xg Delta lines: <Unavailable> File [added]: w_lorentz_minimum.xg Delta lines: <Unavailable> File [added]: w_lorentz_norm1.xg Delta lines: <Unavailable> File [added]: w_lorentz_norm2.xg Delta lines: <Unavailable> File [added]: w_lorentz_x_[16][16].xg Delta lines: <Unavailable> File [added]: w_lorentz_y_[16][16].xg Delta lines: <Unavailable> File [added]: w_lorentz_z_[16][16].xg Delta lines: <Unavailable>

1 0

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

11 Feb '13

User: rhaas Date: 2013/02/11 10:30 AM Modified: /trunk/test/ GRHydro_test_shock.par, GRHydro_test_shock_ppm.par, GRHydro_test_tov_ppm_ML.par, GRHydro_test_tov_ppm_ML_disable_internal_excision.par, GRHydro_test_tov_ppm_no_trp_ML.par, test_tov_carpet_refined_nosync.par, tov_carpetevolutionmask.par, tov_carpetevolutionmask2.par, tov_slowsector.par Log: GRHydro: remove poly_gamma_ini from test cases From: Roland Haas <roland.haas(a)physics.gatech.edu> File Changes: Directory: /trunk/test/ ======================= File [modified]: GRHydro_test_shock.par Delta lines: +0 -1 =================================================================== --- trunk/test/GRHydro_test_shock.par 2013-02-11 16:30:23 UTC (rev 481) +++ trunk/test/GRHydro_test_shock.par 2013-02-11 16:30:26 UTC (rev 482) @@ -64,7 +64,6 @@ eos_omni::poly_gamma = 2.0 -eos_omni::poly_gamma_ini = 2.0 eos_omni::poly_k = 80.0 eos_omni::gl_gamma = 1.66666666666666 eos_omni::gl_k = 80.0 File [modified]: GRHydro_test_shock_ppm.par Delta lines: +0 -1 =================================================================== --- trunk/test/GRHydro_test_shock_ppm.par 2013-02-11 16:30:23 UTC (rev 481) +++ trunk/test/GRHydro_test_shock_ppm.par 2013-02-11 16:30:26 UTC (rev 482) @@ -60,7 +60,6 @@ GRHydro::bound = "flat" eos_omni::poly_gamma = 2.0 -eos_omni::poly_gamma_ini = 2.0 eos_omni::poly_k = 80.0 eos_omni::gl_gamma = 1.66666666666666 eos_omni::gl_k = 80.0 File [modified]: GRHydro_test_tov_ppm_ML.par Delta lines: +0 -1 =================================================================== --- trunk/test/GRHydro_test_tov_ppm_ML.par 2013-02-11 16:30:23 UTC (rev 481) +++ trunk/test/GRHydro_test_tov_ppm_ML.par 2013-02-11 16:30:26 UTC (rev 482) @@ -88,7 +88,6 @@ tovsolver::TOV_Combine_Method = "maximum" eos_omni::poly_gamma = 2.0 -eos_omni::poly_gamma_ini = 2.0 eos_omni::poly_k = 100.0 IOBasic::outInfo_every = 1 File [modified]: GRHydro_test_tov_ppm_ML_disable_internal_excision.par Delta lines: +0 -1 =================================================================== --- trunk/test/GRHydro_test_tov_ppm_ML_disable_internal_excision.par 2013-02-11 16:30:23 UTC (rev 481) +++ trunk/test/GRHydro_test_tov_ppm_ML_disable_internal_excision.par 2013-02-11 16:30:26 UTC (rev 482) @@ -90,7 +90,6 @@ tovsolver::TOV_Combine_Method = "maximum" eos_omni::poly_gamma = 2.0 -eos_omni::poly_gamma_ini = 2.0 eos_omni::poly_k = 100.0 IOBasic::outInfo_every = 1 File [modified]: GRHydro_test_tov_ppm_no_trp_ML.par Delta lines: +0 -1 =================================================================== --- trunk/test/GRHydro_test_tov_ppm_no_trp_ML.par 2013-02-11 16:30:23 UTC (rev 481) +++ trunk/test/GRHydro_test_tov_ppm_no_trp_ML.par 2013-02-11 16:30:26 UTC (rev 482) @@ -89,7 +89,6 @@ tovsolver::TOV_Combine_Method = "maximum" eos_omni::poly_gamma = 2.0 -eos_omni::poly_gamma_ini = 2.0 eos_omni::poly_k = 100.0 IOBasic::outInfo_every = 1 File [modified]: test_tov_carpet_refined_nosync.par Delta lines: +0 -1 =================================================================== --- trunk/test/test_tov_carpet_refined_nosync.par 2013-02-11 16:30:23 UTC (rev 481) +++ trunk/test/test_tov_carpet_refined_nosync.par 2013-02-11 16:30:26 UTC (rev 482) @@ -84,7 +84,6 @@ GRHydro::rho_abs_min = 1.e-7 GRHydro::initial_atmosphere_factor = 0.999 EOS_Omni::poly_gamma = 2.0 -EOS_Omni::poly_gamma_ini = 2.0 EOS_Omni::gl_gamma = 2.0 EOS_Omni::poly_k = 100 EOS_Omni::gl_k = 100 File [modified]: tov_carpetevolutionmask.par Delta lines: +0 -1 =================================================================== --- trunk/test/tov_carpetevolutionmask.par 2013-02-11 16:30:23 UTC (rev 481) +++ trunk/test/tov_carpetevolutionmask.par 2013-02-11 16:30:26 UTC (rev 482) @@ -91,7 +91,6 @@ 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 File [modified]: tov_carpetevolutionmask2.par Delta lines: +0 -1 =================================================================== --- trunk/test/tov_carpetevolutionmask2.par 2013-02-11 16:30:23 UTC (rev 481) +++ trunk/test/tov_carpetevolutionmask2.par 2013-02-11 16:30:26 UTC (rev 482) @@ -92,7 +92,6 @@ tovsolver::TOV_Combine_Method = "maximum" eos_omni::gl_gamma=2.0 -eos_omni::poly_gamma_ini=2.0 eos_omni::gl_k=100. Carpet::prolongation_order_space = 3 File [modified]: tov_slowsector.par Delta lines: +0 -1 =================================================================== --- trunk/test/tov_slowsector.par 2013-02-11 16:30:23 UTC (rev 481) +++ trunk/test/tov_slowsector.par 2013-02-11 16:30:26 UTC (rev 482) @@ -78,7 +78,6 @@ GRHydro::initial_atmosphere_factor = 0.999 # these do not match TOVSolver to make the star actually move EOS_Omni::poly_gamma = 2.0 -EOS_Omni::poly_gamma_ini = 2.0 EOS_Omni::gl_gamma = 2.0 EOS_Omni::poly_k = 100 EOS_Omni::gl_k = 100

1 0

[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

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits February 2013