Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 145)
by knarf＠cct.lsu.edu 16 Aug '10

16 Aug '10

User: knarf Date: 2010/08/16 01:27 PM Modified: /trunk/src/ GRHydro_Eigenproblem.F90 Log: more #define File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Eigenproblem.F90 Delta lines: +2 -0 =================================================================== --- trunk/src/GRHydro_Eigenproblem.F90 2010-08-16 18:23:38 UTC (rev 144) +++ trunk/src/GRHydro_Eigenproblem.F90 2010-08-16 18:27:37 UTC (rev 145) @@ -58,10 +58,12 @@ CCTK_INT handle CCTK_REAL dpdrho,dpdeps,press +#if !USE_EOS_OMNI #ifdef _EOS_BASE_INC_ #undef _EOS_BASE_INC_ #endif #include "EOS_Base.inc" +#endif #if USE_EOS_OMNI ! begin EOS Omni vars

1 0

[svn:einsteintoolkit] GRHydro_InitData/trunk/src/ (Rev. 110)
by knarf＠cct.lsu.edu 16 Aug '10

16 Aug '10

User: knarf Date: 2010/08/16 01:24 PM Modified: /trunk/src/ GRHydro_C2P2C.F90, GRHydro_Con2Prim.F90, GRHydro_SimpleWave.F90 Log: use eos handles File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_C2P2C.F90 Delta lines: +2 -4 =================================================================== --- trunk/src/GRHydro_C2P2C.F90 2010-08-16 13:45:31 UTC (rev 109) +++ trunk/src/GRHydro_C2P2C.F90 2010-08-16 18:24:10 UTC (rev 110) @@ -55,7 +55,6 @@ #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: poly_eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -63,7 +62,6 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - poly_eoskey = GRHydro_poly_eoskey ! end EOS Omni vars #endif @@ -108,10 +106,10 @@ GRHydro_rho_min = 1.e-10 #if USE_EOS_OMNI - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(GRHydro_eos_handle,keytemp,n,& GRHydro_rho_min,1.0d0,xtemp,xye,pmin,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(GRHydro_eos_handle,keytemp,n,& GRHydro_rho_min,xeps,xtemp,xye,pmin,epsmin,keyerr,anyerr) #else pmin = EOS_Pressure(GRHydro_eos_handle, GRHydro_rho_min, 1.0d0) File [modified]: GRHydro_Con2Prim.F90 Delta lines: +2 -4 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2010-08-16 13:45:31 UTC (rev 109) +++ trunk/src/GRHydro_Con2Prim.F90 2010-08-16 18:24:10 UTC (rev 110) @@ -62,7 +62,6 @@ #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: poly_eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -70,7 +69,6 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - poly_eoskey = GRHydro_poly_eoskey ! end EOS Omni vars #endif @@ -110,10 +108,10 @@ epsnegative = .false. #if USE_EOS_OMNI - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(GRHydro_eos_handle,keytemp,n,& GRHydro_rho_min,1.0d0,xtemp,xye,pmin,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(GRHydro_eos_handle,keytemp,n,& GRHydro_rho_min,xeps,xtemp,xye,pmin,epsmin,keyerr,anyerr) #else pmin = EOS_Pressure(GRHydro_eos_handle, GRHydro_rho_min, 1.0d0) File [modified]: GRHydro_SimpleWave.F90 Delta lines: +1 -3 =================================================================== --- trunk/src/GRHydro_SimpleWave.F90 2010-08-16 13:45:31 UTC (rev 109) +++ trunk/src/GRHydro_SimpleWave.F90 2010-08-16 18:24:10 UTC (rev 110) @@ -59,7 +59,6 @@ #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: poly_eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -67,7 +66,6 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - poly_eoskey = GRHydro_poly_eoskey ! end EOS Omni vars #endif @@ -151,7 +149,7 @@ velx(i,1,1) = 0.d0 w_lorentz(i,1,1) = 1.d0 #if USE_EOS_OMNI - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,n,& rho(i,1,1),1.0d0,xtemp,xye,press(i,1,1),keyerr,anyerr) #else press(i,1,1) = EOS_Pressure(GRHydro_polytrope_handle, rho(i,1,1), 1.0d0)

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 144)
by knarf＠cct.lsu.edu 16 Aug '10

16 Aug '10

User: knarf Date: 2010/08/16 01:23 PM Modified: /trunk/src/ GRHydro_Con2Prim.F90, GRHydro_Eigenproblem.F90, GRHydro_Eigenproblem_Marquina.F90, GRHydro_EoSChangeGamma.F90, GRHydro_Prim2Con.F90 Log: more #defines, use integer handles if given, remove eos_key usage File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: +32 -42 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2010-08-16 17:52:56 UTC (rev 143) +++ trunk/src/GRHydro_Con2Prim.F90 2010-08-16 18:23:38 UTC (rev 144) @@ -68,7 +68,6 @@ #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: poly_eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -76,7 +75,6 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - poly_eoskey = GRHydro_poly_eoskey ! end EOS Omni vars #endif @@ -96,9 +94,9 @@ #if USE_EOS_OMNI ! this is a poly call - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,n,& GRHydro_rho_min,xeps,xtemp,xye,pmin,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,n,& GRHydro_rho_min,xeps,xtemp,xye,pmin,epsmin,keyerr,anyerr) #else pmin = EOS_Pressure(GRHydro_polytrope_handle, GRHydro_rho_min, 1.0d0) @@ -160,10 +158,10 @@ vel(i,j,k,:) = 0.d0 w_lorentz(i,j,k) = 1.d0 #if USE_EOS_OMNI - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,n,& rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,n,& rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),eps(i,j,k),keyerr,anyerr) #else press(i,j,k) = EOS_Pressure(GRHydro_polytrope_handle, rho(i,j,k), eps(i,j,k)) @@ -214,10 +212,10 @@ vel(i,j,k,:) = 0.d0 w_lorentz(i,j,k) = 1.d0 #if USE_EOS_OMNI - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,n,& rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,n,& rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),eps(i,j,k),keyerr,anyerr) #else press(i,j,k) = EOS_Pressure(GRHydro_polytrope_handle, rho(i,j,k), eps(i,j,k)) @@ -297,14 +295,12 @@ #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 real*8 :: xpress = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - eoskey = GRHydro_eoskey ! end EOS Omni vars #endif @@ -321,7 +317,7 @@ !!$ Set initial guess for pressure: #if USE_EOS_OMNI - call EOS_Omni_press(eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rho,epsilon,xtemp,xye,xpress,keyerr,anyerr) pold = max(1.d-10,xpress) #else @@ -356,7 +352,7 @@ !!$ Calculate the function #if USE_EOS_OMNI - call EOS_Omni_press(eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rho,epsilon,xtemp,xye,xpress,keyerr,anyerr) f = pold - xpress @@ -408,10 +404,10 @@ end if #if USE_EOS_OMNI - call EOS_Omni_DPressByDRho(eoskey,keytemp,n,& + call EOS_Omni_DPressByDRho(handle,keytemp,n,& rho,epsilon,xtemp,xye,dpressbydrho,keyerr,anyerr) - call EOS_Omni_DPressByDEps(eoskey,keytemp,n,& + call EOS_Omni_DPressByDEps(handle,keytemp,n,& rho,epsilon,xtemp,xye,dpressbydeps,keyerr,anyerr) #else dpressbydrho = EOS_DPressByDRho(handle,rho,epsilon) @@ -437,7 +433,7 @@ udens)/udens #if USE_EOS_OMNI - call EOS_Omni_press(eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rho,epsilon,xtemp,xye,xpress,keyerr,anyerr) f = pnew - xpress @@ -452,10 +448,10 @@ do i=1,GRHydro_polish #if USE_EOS_OMNI - call EOS_Omni_DPressByDRho(eoskey,keytemp,n,& + call EOS_Omni_DPressByDRho(handle,keytemp,n,& rho,epsilon,xtemp,xye,dpressbydrho,keyerr,anyerr) - call EOS_Omni_DPressByDEps(eoskey,keytemp,n,& + call EOS_Omni_DPressByDEps(handle,keytemp,n,& rho,epsilon,xtemp,xye,dpressbydeps,keyerr,anyerr) #else dpressbydrho = EOS_DPressByDRho(handle,rho,epsilon) @@ -479,7 +475,7 @@ udens)/udens #if USE_EOS_OMNI - call EOS_Omni_press(eoskey,keytemp,n, & + call EOS_Omni_press(handle,keytemp,n, & rho,epsilon,xtemp,xye,xpress,keyerr,anyerr) f = pold - xpress @@ -579,8 +575,6 @@ #if USE_EOS_OMNI ! begin EOS omni - CCTK_INT :: eoskey = 0 - CCTK_INT :: poly_eoskey = 0 CCTK_INT :: keyerr(1) = 0 CCTK_INT :: anyerr = 0 CCTK_INT :: keytemp = 0 @@ -588,18 +582,16 @@ CCTK_REAL :: xye = 0.0d0 CCTK_REAL :: xeps = 0.0d0 CCTK_REAL :: xtemp = 0.0d0 - eoskey = GRHydro_eoskey - poly_eoskey = GRHydro_poly_eoskey ! end EOS omni #endif #if USE_EOS_OMNI ! this is a poly call - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,n,& GRHydro_rho_min,1.0d0,xtemp,xye,pmin,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,n,& GRHydro_rho_min,epsmin,xtemp,xye,pmin,epsmin,keyerr,anyerr) #else pmin=EOS_Pressure(GRHydro_polytrope_handle, GRHydro_rho_min, 1.0d0) @@ -937,7 +929,6 @@ #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: poly_eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -945,7 +936,6 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - poly_eoskey = GRHydro_poly_eoskey ! end EOS Omni vars #endif @@ -964,10 +954,10 @@ rhoold = max(GRHydro_rho_min,rho) #if USE_EOS_OMNI -! call EOS_Omni_press(poly_eoskey,keytemp,n,& +! call EOS_Omni_press(handle,keytemp,n,& ! rhoold,xeps,xtemp,xye,xpress,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(handle,keytemp,n,& rhoold,xeps,xtemp,xye,press,xeps,keyerr,anyerr) enthalpy = 1.0d0 + xeps + press / rhoold @@ -979,7 +969,7 @@ w_lorentz = sqrt(1.d0 + s2 / ( (udens*enthalpy)**2 )) #if USE_EOS_OMNI - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rhoold,epsilon,xtemp,xye,press,keyerr,anyerr) #else press = EOS_Pressure(handle, rhoold, epsilon) @@ -1017,10 +1007,10 @@ udens = rhonew dens = sqrtdet * rhonew #if USE_EOS_OMNI - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rhonew,1.0d0,xtemp,xye,press,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(handle,keytemp,n,& rhonew,xeps,xtemp,xye,press,epsilon,keyerr,anyerr) #else press = EOS_Pressure(handle, rhonew, 1.d0) @@ -1042,7 +1032,7 @@ !!$ 2D_Polytrope it gives the right answers. #if USE_EOS_OMNI - call EOS_Omni_DPressByDrho(poly_eoskey,keytemp,n,& + call EOS_Omni_DPressByDrho(handle,keytemp,n,& rhonew,epsilon,xtemp,xye,denthalpy,keyerr,anyerr) denthalpy = denthalpy / rhonew #else @@ -1069,10 +1059,10 @@ !!$ Recalculate primitive variables and function #if USE_EOS_OMNI - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rhonew,epsilon,xtemp,xye,press,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(handle,keytemp,n,& rhonew,xeps,xtemp,xye,press,epsilon,keyerr,anyerr) enthalpy = 1.0d0 + epsilon + press / rhonew @@ -1100,10 +1090,10 @@ dens = sqrtdet * rhonew #if USE_EOS_OMNI - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rhonew,1.0d0,xtemp,xye,press,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(handle,keytemp,n,& rhonew,xeps,xtemp,xye,press,epsilon,keyerr,anyerr) #else press = EOS_Pressure(handle, rhonew, 1.d0) @@ -1124,10 +1114,10 @@ 50 rho = rhonew #if USE_EOS_OMNI - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rhonew,epsilon,xtemp,xye,press,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(handle,keytemp,n,& rhonew,xeps,xtemp,xye,press,epsilon,keyerr,anyerr) enthalpy = 1.0d0 + epsilon + press / rhonew @@ -1138,10 +1128,10 @@ #endif #if USE_EOS_OMNI - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rhonew,epsilon,xtemp,xye,press,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(handle,keytemp,n,& rhonew,xeps,xtemp,xye,press,epsilon,keyerr,anyerr) #else press = EOS_Pressure(handle, rhonew, epsilon) @@ -1166,9 +1156,9 @@ rho = GRHydro_rho_min dens = sqrtdet * rho #if USE_EOS_OMNI - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rhonew,1.d0,xtemp,xye,press,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(handle,keytemp,n,& rhonew,xeps,xtemp,xye,press,epsilon,keyerr,anyerr) #else press = EOS_Pressure(handle, rhonew, 1.d0) File [modified]: GRHydro_Eigenproblem.F90 Delta lines: +10 -16 =================================================================== --- trunk/src/GRHydro_Eigenproblem.F90 2010-08-16 17:52:56 UTC (rev 143) +++ trunk/src/GRHydro_Eigenproblem.F90 2010-08-16 18:23:38 UTC (rev 144) @@ -66,7 +66,6 @@ #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -74,7 +73,6 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - eoskey = GRHydro_eoskey ! end EOS Omni vars #endif @@ -86,15 +84,15 @@ #if USE_EOS_OMNI -! call EOS_Omni_cs2(eoskey,keytemp,n,& +! call EOS_Omni_cs2(handle,keytemp,n,& ! rho,eps,xtemp,xye,cs2,keyerr,anyerr) - call EOS_Omni_press(eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rho,eps,xtemp,xye,press,keyerr,anyerr) - call EOS_Omni_DPressByDEps(eoskey,keytemp,n,& + call EOS_Omni_DPressByDEps(handle,keytemp,n,& rho,eps,xtemp,xye,dpdeps,keyerr,anyerr) - call EOS_Omni_DPressByDRho(eoskey,keytemp,n,& + call EOS_Omni_DPressByDRho(handle,keytemp,n,& rho,eps,xtemp,xye,dpdrho,keyerr,anyerr) cs2 = (dpdrho + press * dpdeps / (rho**2))/ & @@ -192,7 +190,6 @@ #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -200,7 +197,6 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - eoskey = GRHydro_eoskey ! end EOS Omni vars #endif @@ -210,18 +206,18 @@ !!$ Set required fluid quantities #if USE_EOS_OMNI - call EOS_Omni_press(eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rho,eps,xtemp,xye,press,keyerr,anyerr) - call EOS_Omni_DPressByDEps(eoskey,keytemp,n,& + call EOS_Omni_DPressByDEps(handle,keytemp,n,& rho,eps,xtemp,xye,dpdeps,keyerr,anyerr) - call EOS_Omni_DPressByDRho(eoskey,keytemp,n,& + call EOS_Omni_DPressByDRho(handle,keytemp,n,& rho,eps,xtemp,xye,dpdrho,keyerr,anyerr) cs2 = (dpdrho + press * dpdeps / (rho**2))/ & (1.0d0 + eps + press/rho) -! call EOS_Omni_cs2(eoskey,keytemp,n,& +! call EOS_Omni_cs2(handle,keytemp,n,& ! rho,eps,xtemp,xye,cs2,keyerr,anyerr) #else press = EOS_Pressure(handle,rho,eps) @@ -577,7 +573,6 @@ #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -585,7 +580,6 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - eoskey = GRHydro_eoskey ! end EOS Omni vars #endif @@ -594,10 +588,10 @@ !!$ Set required fluid quantities #if USE_EOS_OMNI - call EOS_Omni_press(eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rho,eps,xtemp,xye,press,keyerr,anyerr) - call EOS_Omni_cs2(eoskey,keytemp,n,& + call EOS_Omni_cs2(handle,keytemp,n,& rho,eps,xtemp,xye,cs2,keyerr,anyerr) #else press = EOS_Pressure(handle,rho,eps) File [modified]: GRHydro_Eigenproblem_Marquina.F90 Delta lines: +10 -10 =================================================================== --- trunk/src/GRHydro_Eigenproblem_Marquina.F90 2010-08-16 17:52:56 UTC (rev 143) +++ trunk/src/GRHydro_Eigenproblem_Marquina.F90 2010-08-16 18:23:38 UTC (rev 144) @@ -92,15 +92,16 @@ !!$ Warning, warning. Nasty hack follows +#if !USE_EOS_OMNI #ifdef _EOS_BASE_INC_ #undef _EOS_BASE_INC_ #endif #include "EOS_Base.inc" +#endif #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -108,7 +109,6 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - eoskey = GRHydro_eoskey ! end EOS Omni vars #endif @@ -120,13 +120,13 @@ !!$ Set required fluid quantities invrhol = one / rhol #if USE_EOS_OMNI - call EOS_Omni_press(eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rhol,epsl,xtemp,xye,pressl,keyerr,anyerr) -! call EOS_Omni_cs2(eoskey,keytemp,n,& +! call EOS_Omni_cs2(handle,keytemp,n,& ! rhol,epsl,xtemp,xye,cs2l,keyerr,anyerr) - call EOS_Omni_DPressByDEps(eoskey,keytemp,n,& + call EOS_Omni_DPressByDEps(handle,keytemp,n,& rhol,epsl,xtemp,xye,dpdepsl,keyerr,anyerr) - call EOS_Omni_DPressByDRho(eoskey,keytemp,n,& + call EOS_Omni_DPressByDRho(handle,keytemp,n,& rhol,epsl,xtemp,xye,dpdrhol,keyerr,anyerr) cs2l = (dpdrhol + pressl * dpdepsl * invrhol**2)/ & (1.0d0 + epsl + pressl*invrhol) @@ -174,13 +174,13 @@ !!$ Set required fluid quantities invrhor = one / rhor #if USE_EOS_OMNI - call EOS_Omni_press(eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& rhor,epsr,xtemp,xye,pressr,keyerr,anyerr) -! call EOS_Omni_cs2(eoskey,keytemp,n,& +! call EOS_Omni_cs2(handle,keytemp,n,& ! rhor,epsr,xtemp,xye,cs2r,keyerr,anyerr) - call EOS_Omni_DPressByDEps(eoskey,keytemp,n,& + call EOS_Omni_DPressByDEps(handle,keytemp,n,& rhor,epsr,xtemp,xye,dpdepsr,keyerr,anyerr) - call EOS_Omni_DPressByDRho(eoskey,keytemp,n,& + call EOS_Omni_DPressByDRho(handle,keytemp,n,& rhor,epsr,xtemp,xye,dpdrhor,keyerr,anyerr) cs2r = (dpdrhor + pressr * dpdepsr * invrhor**2)/ & (1.0d0 + epsr + pressr*invrhor) File [modified]: GRHydro_EoSChangeGamma.F90 Delta lines: +6 -12 =================================================================== --- trunk/src/GRHydro_EoSChangeGamma.F90 2010-08-16 17:52:56 UTC (rev 143) +++ trunk/src/GRHydro_EoSChangeGamma.F90 2010-08-16 18:23:38 UTC (rev 144) @@ -62,7 +62,6 @@ #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: poly_eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -70,12 +69,11 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - poly_eoskey = GRHydro_poly_eoskey ! end EOS Omni vars - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,n,& 1.0d0,1.0d0,xtemp,xye,xpress,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,n,& 1.0d0,1.0d0,xtemp,xye,xpress,xeps,keyerr,anyerr) local_Gamma = 1.0d0 + xpress/xeps @@ -161,7 +159,6 @@ #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: poly_eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -169,12 +166,11 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - poly_eoskey = GRHydro_poly_eoskey ! end EOS Omni vars - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,n,& 1.0d0,1.0d0,xtemp,xye,xpress,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,n,& 1.0d0,1.0d0,xtemp,xye,xpress,xeps,keyerr,anyerr) local_Gamma = 1.0d0 + xpress/xeps @@ -267,7 +263,6 @@ #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: poly_eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -275,12 +270,11 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - poly_eoskey = GRHydro_poly_eoskey ! end EOS Omni vars - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,n,& 1.0d0,1.0d0,xtemp,xye,xpress,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,n,& 1.0d0,1.0d0,xtemp,xye,xpress,xeps,keyerr,anyerr) local_Gamma = 1.0d0 + xpress/xeps File [modified]: GRHydro_Prim2Con.F90 Delta lines: +7 -7 =================================================================== --- trunk/src/GRHydro_Prim2Con.F90 2010-08-16 17:52:56 UTC (rev 143) +++ trunk/src/GRHydro_Prim2Con.F90 2010-08-16 18:23:38 UTC (rev 144) @@ -129,12 +129,13 @@ deps, dpress, w, vlowx, vlowy, vlowz CCTK_INT :: handle +#if !USE_EOS_OMNI #include "EOS_Base.inc" +#endif #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -142,7 +143,6 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - eoskey = GRHydro_eoskey ! end EOS Omni vars #endif @@ -151,7 +151,7 @@ *dvely*dvelz)) #if USE_EOS_OMNI - call EOS_Omni_press(eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& drho,deps,xtemp,xye,dpress,keyerr,anyerr) #else dpress = EOS_Pressure(handle, drho, deps) @@ -341,15 +341,16 @@ CCTK_INT :: handle character(len=256) NaN_WarnLine +#if !USE_EOS_OMNI #ifdef _EOS_BASE_INC_ #undef _EOS_BASE_INC_ #endif #include "EOS_Base.inc" +#endif #if USE_EOS_OMNI ! begin EOS Omni vars integer :: n = 1 - integer :: poly_eoskey = 0 integer :: keytemp = 0 integer :: anyerr = 0 integer :: keyerr(1) = 0 @@ -357,7 +358,6 @@ real*8 :: xeps = 0.0d0 real*8 :: xtemp = 0.0d0 real*8 :: xye = 0.0d0 - poly_eoskey = GRHydro_poly_eoskey ! end EOS Omni vars #endif @@ -380,10 +380,10 @@ endif #if USE_EOS_OMNI - call EOS_Omni_press(poly_eoskey,keytemp,n,& + call EOS_Omni_press(handle,keytemp,n,& drho,xeps,xtemp,xye,dpress,keyerr,anyerr) - call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + call EOS_Omni_EpsFromPress(handle,keytemp,n,& drho,xeps,xtemp,xye,dpress,deps,keyerr,anyerr) #else if (handle .ge. 0) then

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 143)
by knarf＠cct.lsu.edu 16 Aug '10

16 Aug '10

User: knarf Date: 2010/08/16 12:52 PM Modified: /trunk/src/ GRHydro_Con2Prim.F90, GRHydro_EOS.c, GRHydro_EoSChangeGamma.F90 Log: some more #defines for EOS_Omni, use handle-function to get EOS_Handle from EOS_Omni File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: +15 -0 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2010-07-30 03:20:39 UTC (rev 142) +++ trunk/src/GRHydro_Con2Prim.F90 2010-08-16 17:52:56 UTC (rev 143) @@ -48,10 +48,12 @@ DECLARE_CCTK_ARGUMENTS DECLARE_CCTK_PARAMETERS +#if !USE_EOS_OMNI #ifdef _EOS_BASE_INC_ #undef _EOS_BASE_INC_ #endif #include "EOS_Base.inc" +#endif integer :: i, j, k, itracer, nx, ny, nz CCTK_REAL :: uxx, uxy, uxz, uyy, uyz, uzz, det, psi4pt, pmin, epsmin @@ -284,10 +286,12 @@ character(len=200) warnline logical epsnegative +#if !USE_EOS_OMNI #ifdef _EOS_BASE_INC_ #undef _EOS_BASE_INC_ #endif #include "EOS_Base.inc" +#endif #if USE_EOS_OMNI @@ -565,10 +569,12 @@ CCTK_REAL :: local_min_tracer +#if !USE_EOS_OMNI #ifdef _EOS_BASE_INC_ #undef _EOS_BASE_INC_ #endif #include "EOS_Base.inc" +#endif #if USE_EOS_OMNI @@ -920,10 +926,12 @@ enthalpy, denthalpy, sqrtdet, invsqrtdet, invfac, GRHydro_C2P_failed character(len=200) warnline +#if !USE_EOS_OMNI #ifdef _EOS_BASE_INC_ #undef _EOS_BASE_INC_ #endif #include "EOS_Base.inc" +#endif #if USE_EOS_OMNI @@ -1157,8 +1165,15 @@ rho = GRHydro_rho_min dens = sqrtdet * rho +#if USE_EOS_OMNI + call EOS_Omni_press(poly_eoskey,keytemp,n,& + rhonew,1.d0,xtemp,xye,press,keyerr,anyerr) + call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + rhonew,xeps,xtemp,xye,press,epsilon,keyerr,anyerr) +#else press = EOS_Pressure(handle, rhonew, 1.d0) epsilon = EOS_SpecificIntEnergy(handle, rhonew, press) +#endif ! w_lorentz=1, so the expression for tau reduces to: tau = sqrtdet * (rho+rho*epsilon) - dens usx = 0.d0 File [modified]: GRHydro_EOS.c Delta lines: +19 -29 =================================================================== --- trunk/src/GRHydro_EOS.c 2010-07-30 03:20:39 UTC (rev 142) +++ trunk/src/GRHydro_EOS.c 2010-08-16 17:52:56 UTC (rev 143) @@ -15,7 +15,9 @@ #include <stdio.h> #include <stdlib.h> +#if !USE_EOS_OMNI #include "EOS_Base.h" +#endif /*@@ @routine GRHydro_EOSHandle @@ -38,34 +40,22 @@ DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; -#ifndef USE_EOS_OMNI - int ierr = 0; - char *infoline; - - *GRHydro_eos_handle = EOS_Handle(GRHydro_eos_table); - - if (*GRHydro_eos_handle < 0) { - ierr = *GRHydro_eos_handle; - CCTK_WARN(0, "Table handle is incorrect"); - } - - *GRHydro_polytrope_handle = EOS_Handle("2D_Polytrope"); - - if (*GRHydro_polytrope_handle < 0) { - ierr = *GRHydro_polytrope_handle; - CCTK_WARN(0, "GRHydro requires the 2D_Polytrope or the 2D_HybridPolytrope EoSs to be registered, for the atmosphere"); - } - - infoline = (char *)malloc(100*sizeof(char)); - if (!infoline) { - CCTK_WARN(0, "Failed to malloc an array of 100 chars in GRHydro_EOS."); - } - else { - sprintf(infoline,"GRHydro will use the %s equation of state.", - GRHydro_eos_table); - CCTK_INFO(infoline); - } - free(infoline); - infoline = NULL; +#ifdef USE_EOS_OMNI + CCTK_INFO("Trying to get EOS handles"); + EOS_Omni_GetHandle("2D_Polytrope"); + CCTK_INFO("Trying to get EOS handles"); + if (!(*GRHydro_eos_handle = EOS_Omni_GetHandle(GRHydro_eos_table))) + CCTK_WARN(0, "Cannot get EOS handle, aborting."); + if (!(*GRHydro_polytrope_handle = EOS_Omni_GetHandle("2D_Polytrope"))) + CCTK_WARN(0, "Cannot get polytrope EOS handle, aborting."); + CCTK_VInfo(CCTK_THORNSTRING, "GRHydro will use the %s equation of state.", + GRHydro_eos_table); +#else + if ((*GRHydro_eos_handle = EOS_Handle(GRHydro_eos_table)) < 0) + CCTK_WARN(0, "Cannot get EOS handle, aborting."); + if ((*GRHydro_polytrope_handle = EOS_Handle("2D_Polytrope")) < 0); + CCTK_WARN(0, "Cannot get polytrope EOS handle, aborting."); + CCTK_VInfo(CCTK_THORNSTRING, "GRHydro will use the %s equation of state.", + GRHydro_eos_table); #endif } File [modified]: GRHydro_EoSChangeGamma.F90 Delta lines: +8 -0 =================================================================== --- trunk/src/GRHydro_EoSChangeGamma.F90 2010-07-30 03:20:39 UTC (rev 142) +++ trunk/src/GRHydro_EoSChangeGamma.F90 2010-08-16 17:52:56 UTC (rev 143) @@ -54,7 +54,9 @@ CCTK_REAL :: local_gamma +#if !USE_EOS_OMNI #include "EOS_Base.inc" +#endif !!$ Set up the fluid constants #if USE_EOS_OMNI @@ -135,7 +137,9 @@ subroutine GRHydro_EoSChangeK(CCTK_ARGUMENTS) +#if !USE_EOS_OMNI USE EOS_Polytrope_Scalars +#endif implicit none @@ -149,7 +153,9 @@ CCTK_REAL, dimension(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)) :: Q +#if !USE_EOS_OMNI #include "EOS_Base.inc" +#endif !!$ Set up the fluid constants #if USE_EOS_OMNI @@ -253,7 +259,9 @@ character(len=100) infoline +#if !USE_EOS_OMNI #include "EOS_Base.inc" +#endif !!$ Set up the fluid constants #if USE_EOS_OMNI

1 0

[svn:einsteintoolkit] incoming/EOS_Omni/ (Rev. 6)
by knarf＠cct.lsu.edu 16 Aug '10

16 Aug '10

User: knarf Date: 2010/08/16 12:52 PM Added: /EOS_Omni/src/ EOS_Omni_Handles.c Removed: /EOS_Omni/src/ EOS_Omni_Names.cpp Modified: /EOS_Omni/ interface.ccl /EOS_Omni/src/ make.code.defn Log: actually provide the handle function, change to c for simplicity (from c++) File Changes: Directory: /EOS_Omni/ ===================== File [modified]: interface.ccl Delta lines: +2 -0 =================================================================== --- EOS_Omni/interface.ccl 2010-08-16 14:12:22 UTC (rev 5) +++ EOS_Omni/interface.ccl 2010-08-16 17:52:05 UTC (rev 6) @@ -6,6 +6,8 @@ CCTK_INT FUNCTION EOS_Omni_GetHandle(CCTK_STRING IN name) +PROVIDES FUNCTION EOS_Omni_GetHandle with EOS_Omni_GetHandle_ LANGUAGE C + void FUNCTION EOS_Omni_press(CCTK_INT IN eoskey, \ CCTK_INT IN havetemp, \ CCTK_INT IN npoints, \ Directory: /EOS_Omni/src/ ========================= File [added]: EOS_Omni_Handles.c Delta lines: +15 -0 =================================================================== --- EOS_Omni/src/EOS_Omni_Handles.c (rev 0) +++ EOS_Omni/src/EOS_Omni_Handles.c 2010-08-16 17:52:05 UTC (rev 6) @@ -0,0 +1,15 @@ +#include <cctk.h> +#include <cctk_Arguments.h> + +CCTK_INT EOS_Omni_GetHandle_(CCTK_STRING name) +{ + CCTK_INFO("GetHandle"); + if (CCTK_EQUALS(name, "2D_Polytrope")) + return 1; + if (CCTK_EQUALS(name, "Ideal_Fluid")) + return 2; + if (CCTK_EQUALS(name, "Hybrid")) + return 3; + return 0; +} + File [removed]: EOS_Omni_Names.cpp Delta lines: +0 -15 =================================================================== --- EOS_Omni/src/EOS_Omni_Names.cpp 2010-08-16 14:12:22 UTC (rev 5) +++ EOS_Omni/src/EOS_Omni_Names.cpp 2010-08-16 17:52:05 UTC (rev 6) @@ -1,15 +0,0 @@ -#include <cctk.h> -#include <cctk_Arguments.h> -#include <cctk_Parameters.h> - -CCTK_INT EOS_Omni_GetHandle(CCTK_STRING name) -{ - if (CCTK_EQUALS(name, "2D_Polytrope")) - return 1; - if (CCTK_EQUALS(name, "Ideal_Fluid")) - return 2; - if (CCTK_EQUALS(name, "Hybrid")) - return 3; - return 0; -} - File [modified]: make.code.defn Delta lines: +1 -1 =================================================================== --- EOS_Omni/src/make.code.defn 2010-08-16 14:12:22 UTC (rev 5) +++ EOS_Omni/src/make.code.defn 2010-08-16 17:52:05 UTC (rev 6) @@ -1,7 +1,7 @@ # Main make.code.defn file for thorn EOS_Omni # Source files in this directory -SRCS = EOS_Omni_Module.F90 EOS_Omni_Startup.F90 EOS_Omni_SingleVarCalls.F90 EOS_Omni_Names.cpp +SRCS = EOS_Omni_Module.F90 EOS_Omni_Startup.F90 EOS_Omni_SingleVarCalls.F90 EOS_Omni_Handles.c # Subdirectories containing source files SUBDIRS =

1 0

[svn:einsteintoolkit] incoming/EOS_Omni/ (Rev. 5)
by knarf＠cct.lsu.edu 16 Aug '10

16 Aug '10

User: knarf Date: 2010/08/16 09:12 AM Added: /EOS_Omni/src/ EOS_Omni_Names.cpp Modified: /EOS_Omni/ interface.ccl /EOS_Omni/src/ make.code.defn Log: provide function to get EOS key/handle by string/name: EOS_Omni_GetHandle File Changes: Directory: /EOS_Omni/ ===================== File [modified]: interface.ccl Delta lines: +2 -0 =================================================================== --- EOS_Omni/interface.ccl 2010-08-09 14:20:16 UTC (rev 4) +++ EOS_Omni/interface.ccl 2010-08-16 14:12:22 UTC (rev 5) @@ -4,6 +4,8 @@ inherits: +CCTK_INT FUNCTION EOS_Omni_GetHandle(CCTK_STRING IN name) + void FUNCTION EOS_Omni_press(CCTK_INT IN eoskey, \ CCTK_INT IN havetemp, \ CCTK_INT IN npoints, \ Directory: /EOS_Omni/src/ ========================= File [added]: EOS_Omni_Names.cpp Delta lines: +15 -0 =================================================================== --- EOS_Omni/src/EOS_Omni_Names.cpp (rev 0) +++ EOS_Omni/src/EOS_Omni_Names.cpp 2010-08-16 14:12:22 UTC (rev 5) @@ -0,0 +1,15 @@ +#include <cctk.h> +#include <cctk_Arguments.h> +#include <cctk_Parameters.h> + +CCTK_INT EOS_Omni_GetHandle(CCTK_STRING name) +{ + if (CCTK_EQUALS(name, "2D_Polytrope")) + return 1; + if (CCTK_EQUALS(name, "Ideal_Fluid")) + return 2; + if (CCTK_EQUALS(name, "Hybrid")) + return 3; + return 0; +} + File [modified]: make.code.defn Delta lines: +1 -1 =================================================================== --- EOS_Omni/src/make.code.defn 2010-08-09 14:20:16 UTC (rev 4) +++ EOS_Omni/src/make.code.defn 2010-08-16 14:12:22 UTC (rev 5) @@ -1,7 +1,7 @@ # Main make.code.defn file for thorn EOS_Omni # Source files in this directory -SRCS = EOS_Omni_Module.F90 EOS_Omni_Startup.F90 EOS_Omni_SingleVarCalls.F90 +SRCS = EOS_Omni_Module.F90 EOS_Omni_Startup.F90 EOS_Omni_SingleVarCalls.F90 EOS_Omni_Names.cpp # Subdirectories containing source files SUBDIRS =

1 0

[svn:einsteintoolkit] GRHydro_InitData/trunk/src/ (Rev. 109)
by knarf＠cct.lsu.edu 16 Aug '10

16 Aug '10

User: knarf Date: 2010/08/16 08:45 AM Modified: /trunk/src/ GRHydro_C2P2C.F90, GRHydro_Con2Prim.F90, GRHydro_ReadConformalData.F90, GRHydro_SimpleWave.F90 Log: add USE_EOS_OMNI checks File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_C2P2C.F90 Delta lines: +25 -0 =================================================================== --- trunk/src/GRHydro_C2P2C.F90 2010-07-12 14:50:33 UTC (rev 108) +++ trunk/src/GRHydro_C2P2C.F90 2010-08-16 13:45:31 UTC (rev 109) @@ -34,10 +34,12 @@ DECLARE_CCTK_ARGUMENTS DECLARE_CCTK_PARAMETERS +#if !USE_EOS_OMNI #ifdef _EOS_BASE_INC_ #undef _EOS_BASE_INC_ #endif #include "EOS_Base.inc" +#endif integer didit,i,j,k,nx,ny,nz CCTK_REAL det @@ -50,6 +52,21 @@ CCTK_INT C2P_failed logical epsnegative +#if USE_EOS_OMNI +! begin EOS Omni vars + integer :: n = 1 + integer :: poly_eoskey = 0 + integer :: keytemp = 0 + integer :: anyerr = 0 + integer :: keyerr(1) = 0 + real*8 :: xpress = 0.0d0 + real*8 :: xeps = 0.0d0 + real*8 :: xtemp = 0.0d0 + real*8 :: xye = 0.0d0 + poly_eoskey = GRHydro_poly_eoskey +! end EOS Omni vars +#endif + call CCTK_WARN(1,"This test works only with Ideal_Fluid EoS") nx = cctk_lsh(1) @@ -90,8 +107,16 @@ epsnegative = .false. GRHydro_rho_min = 1.e-10 +#if USE_EOS_OMNI + call EOS_Omni_press(poly_eoskey,keytemp,n,& + GRHydro_rho_min,1.0d0,xtemp,xye,pmin,keyerr,anyerr) + + call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + GRHydro_rho_min,xeps,xtemp,xye,pmin,epsmin,keyerr,anyerr) +#else pmin = EOS_Pressure(GRHydro_eos_handle, GRHydro_rho_min, 1.0d0) epsmin = EOS_SpecificIntEnergy(GRHydro_eos_handle, GRHydro_rho_min, pmin) +#endif C2P_failed = 0 write(*,*) 'C2P2C test: initial values.' File [modified]: GRHydro_Con2Prim.F90 Delta lines: +25 -0 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2010-07-12 14:50:33 UTC (rev 108) +++ trunk/src/GRHydro_Con2Prim.F90 2010-08-16 13:45:31 UTC (rev 109) @@ -52,11 +52,28 @@ CCTK_INT C2P_failed logical epsnegative +#ifndef USE_EOS_OMNI #ifdef _EOS_BASE_INC_ #undef _EOS_BASE_INC_ #endif #include "EOS_Base.inc" +#endif +#if USE_EOS_OMNI +! begin EOS Omni vars + integer :: n = 1 + integer :: poly_eoskey = 0 + integer :: keytemp = 0 + integer :: anyerr = 0 + integer :: keyerr(1) = 0 + real*8 :: xpress = 0.0d0 + real*8 :: xeps = 0.0d0 + real*8 :: xtemp = 0.0d0 + real*8 :: xye = 0.0d0 + poly_eoskey = GRHydro_poly_eoskey +! end EOS Omni vars +#endif + call CCTK_WARN(1,"For this test, remember to use a polytropic EoS and to set eos_gamma = 2.0 and eos_k = 100.0") nx = cctk_lsh(1) @@ -92,8 +109,16 @@ epsnegative = .false. +#if USE_EOS_OMNI + call EOS_Omni_press(poly_eoskey,keytemp,n,& + GRHydro_rho_min,1.0d0,xtemp,xye,pmin,keyerr,anyerr) + + call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + GRHydro_rho_min,xeps,xtemp,xye,pmin,epsmin,keyerr,anyerr) +#else pmin = EOS_Pressure(GRHydro_eos_handle, GRHydro_rho_min, 1.0d0) epsmin = EOS_SpecificIntEnergy(GRHydro_eos_handle, GRHydro_rho_min, pmin) +#endif C2P_failed = 0 write(*,*) 'Con2Prim test: converting to primitive variables' File [modified]: GRHydro_ReadConformalData.F90 Delta lines: +24 -4 =================================================================== --- trunk/src/GRHydro_ReadConformalData.F90 2010-07-12 14:50:33 UTC (rev 108) +++ trunk/src/GRHydro_ReadConformalData.F90 2010-08-16 13:45:31 UTC (rev 109) @@ -51,9 +51,23 @@ CCTK_INT :: i,j,k,handle,ierr - CCTK_REAL :: eos_k, eos_gamma, rho_min, det +! CCTK_REAL :: eos_k, eos_gamma + CCTK_REAL :: rho_min, det - +#if USE_EOS_OMNI +! begin EOS Omni vars + integer :: n = 1 + integer :: poly_eoskey = 0 + integer :: keytemp = 0 + integer :: anyerr = 0 + integer :: keyerr(1) = 0 + real*8 :: xpress = 0.0d0 + real*8 :: xeps = 0.0d0 + real*8 :: xtemp = 0.0d0 + real*8 :: xye = 0.0d0 + poly_eoskey = GRHydro_poly_eoskey +! end EOS Omni vars +#endif ! only gxx has been read in; copy it into gyy and gzz as well gyy = gxx @@ -131,8 +145,8 @@ ! if (handle < 0) call CCTK_WARN(0,"For this hack you need to compile with EOS_Polytrope") - eos_k = EOS_Pressure(GRHydro_eos_handle, 1.0, 1.0) - eos_gamma = 1.0 + eos_k / EOS_SpecificIntEnergy(GRHydro_eos_handle,1.0,1.0) +! eos_k = EOS_Pressure(GRHydro_eos_handle, 1.0, 1.0) +! eos_gamma = 1.0 + eos_k / EOS_SpecificIntEnergy(GRHydro_eos_handle,1.0,1.0) ! press = eos_k * rho**eos_gamma @@ -140,7 +154,13 @@ do i=1,cctk_lsh(1) do j=1,cctk_lsh(2) do k=1,cctk_lsh(3) +#ifdef USE_EOS_OMNI + call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + rho(i,j,k),xeps,xtemp,xye,press(i,j,k),eps(i,j,k),keyerr,anyerr) + +#else eps(i,j,k) = EOS_SpecificIntEnergy(GRHydro_eos_handle,rho(i,j,k),press(i,j,k)) +#endif end do end do end do File [modified]: GRHydro_SimpleWave.F90 Delta lines: +23 -1 =================================================================== --- trunk/src/GRHydro_SimpleWave.F90 2010-07-12 14:50:33 UTC (rev 108) +++ trunk/src/GRHydro_SimpleWave.F90 2010-08-16 13:45:31 UTC (rev 109) @@ -46,14 +46,31 @@ DECLARE_CCTK_PARAMETERS DECLARE_CCTK_FUNCTIONS +#ifndef USE_EOS_OMNI #ifdef _EOS_BASE_INC_ #undef _EOS_BASE_INC_ #endif #include "EOS_Base.inc" - +#endif + CCTK_INT :: i, j, k, nx, ny, nz CCTK_REAL :: dr, k1, k2, k3, k4, in_data, old_data, source_data, new_data, c_0, det, pi +#if USE_EOS_OMNI +! begin EOS Omni vars + integer :: n = 1 + integer :: poly_eoskey = 0 + integer :: keytemp = 0 + integer :: anyerr = 0 + integer :: keyerr(1) = 0 + real*8 :: xpress = 0.0d0 + real*8 :: xeps = 0.0d0 + real*8 :: xtemp = 0.0d0 + real*8 :: xye = 0.0d0 + poly_eoskey = GRHydro_poly_eoskey +! end EOS Omni vars +#endif + call CCTK_INFO("Setting initial data for a simple wave as Anile Miller Motta") call CCTK_WARN(1, "The simple-wave initial-data routine works only for unigrid and on one node.") @@ -133,7 +150,12 @@ eps(i,1,1) = rho_abs_min**(1.d0/3.d0) velx(i,1,1) = 0.d0 w_lorentz(i,1,1) = 1.d0 +#if USE_EOS_OMNI + call EOS_Omni_press(poly_eoskey,keytemp,n,& + rho(i,1,1),1.0d0,xtemp,xye,press(i,1,1),keyerr,anyerr) +#else press(i,1,1) = EOS_Pressure(GRHydro_polytrope_handle, rho(i,1,1), 1.0d0) +#endif ! polytrope only (initial data) end if

1 0

[svn:einsteintoolkit] manifest/trunk/ (Rev. 36)
by schnetter＠cct.lsu.edu 10 Aug '10

10 Aug '10

User: eschnett Date: 2010/08/10 10:15 PM Modified: /trunk/ einsteintoolkit.th Log: Undo adding CactusTest and CactusWave thorns; not all test cases pass. File Changes: Directory: /trunk/ ================== File [modified]: einsteintoolkit.th Delta lines: +0 -41 =================================================================== --- trunk/einsteintoolkit.th 2010-08-11 01:37:59 UTC (rev 35) +++ trunk/einsteintoolkit.th 2010-08-11 03:15:13 UTC (rev 36) @@ -98,51 +98,17 @@ CactusPUGH/PUGHInterp CactusPUGH/PUGHReduce CactusPUGH/PUGHSlab - CactusPUGHIO/IOHDF5Util CactusPUGHIO/IOHDF5 -#CactusTest/TestAllTypes -CactusTest/TestArrays -CactusTest/TestComplex -CactusTest/TestCoordinates -CactusTest/TestFortranDependencies1 -CactusTest/TestFortranDependencies2 -CactusTest/TestFpointerNULL -CactusTest/TestFreeF90 -CactusTest/TestGlobalReduce -CactusTest/TestInclude1 -CactusTest/TestInclude2 -CactusTest/TestLocalReduce -CactusTest/TestLoop -CactusTest/TestPar -CactusTest/TestReduce -CactusTest/TestSchedule -CactusTest/TestStrings -CactusTest/TestTable -CactusTest/TestTimers -CactusTest/TestTypes - CactusUtils/Formaline CactusUtils/NaNChecker CactusUtils/Nice CactusUtils/TerminationTrigger CactusUtils/TimerReport -CactusWave/IDScalarWave -CactusWave/IDScalarWaveC -CactusWave/IDScalarWaveCXX -CactusWave/IDScalarWaveElliptic -CactusWave/WaveBinarySource -CactusWave/WaveToyC -CactusWave/WaveToyCXX -CactusWave/WaveToyExtra -CactusWave/WaveToyF77 -CactusWave/WaveToyF90 -CactusWave/WaveToyFreeF90 - # The Einstein Toolkit !TARGET = $ARR !TYPE = svn @@ -293,10 +259,3 @@ Carpet/CarpetSlab Carpet/CarpetTracker Carpet/LoopControl - - - -# Private thorns -!TARGET = $ARR -!TYPE = ignore -!CHECKOUT =

1 0

[svn:einsteintoolkit] manifest/trunk/ (Rev. 35)
by schnetter＠cct.lsu.edu 10 Aug '10

10 Aug '10

User: eschnett Date: 2010/08/10 08:37 PM Modified: /trunk/ einsteintoolkit.th Log: Add CactusTest arrangement, which contains many good test cases. Add CactusWave barrangement, which contains our example application. Add (empty) example for adding local thorns. File Changes: Directory: /trunk/ ================== File [modified]: einsteintoolkit.th Delta lines: +41 -0 =================================================================== --- trunk/einsteintoolkit.th 2010-07-12 14:48:59 UTC (rev 34) +++ trunk/einsteintoolkit.th 2010-08-11 01:37:59 UTC (rev 35) @@ -98,17 +98,51 @@ CactusPUGH/PUGHInterp CactusPUGH/PUGHReduce CactusPUGH/PUGHSlab + CactusPUGHIO/IOHDF5Util CactusPUGHIO/IOHDF5 +#CactusTest/TestAllTypes +CactusTest/TestArrays +CactusTest/TestComplex +CactusTest/TestCoordinates +CactusTest/TestFortranDependencies1 +CactusTest/TestFortranDependencies2 +CactusTest/TestFpointerNULL +CactusTest/TestFreeF90 +CactusTest/TestGlobalReduce +CactusTest/TestInclude1 +CactusTest/TestInclude2 +CactusTest/TestLocalReduce +CactusTest/TestLoop +CactusTest/TestPar +CactusTest/TestReduce +CactusTest/TestSchedule +CactusTest/TestStrings +CactusTest/TestTable +CactusTest/TestTimers +CactusTest/TestTypes + CactusUtils/Formaline CactusUtils/NaNChecker CactusUtils/Nice CactusUtils/TerminationTrigger CactusUtils/TimerReport +CactusWave/IDScalarWave +CactusWave/IDScalarWaveC +CactusWave/IDScalarWaveCXX +CactusWave/IDScalarWaveElliptic +CactusWave/WaveBinarySource +CactusWave/WaveToyC +CactusWave/WaveToyCXX +CactusWave/WaveToyExtra +CactusWave/WaveToyF77 +CactusWave/WaveToyF90 +CactusWave/WaveToyFreeF90 + # The Einstein Toolkit !TARGET = $ARR !TYPE = svn @@ -259,3 +293,10 @@ Carpet/CarpetSlab Carpet/CarpetTracker Carpet/LoopControl + + + +# Private thorns +!TARGET = $ARR +!TYPE = ignore +!CHECKOUT =

1 0

[svn:einsteintoolkit] AHFinderDirect/trunk/ (Rev. 1551)
by knarf＠cct.lsu.edu 10 Aug '10

10 Aug '10

User: knarf Date: 2010/08/10 09:12 AM Modified: /trunk/ schedule.ccl Log: only whitespace change: replace tabs by spaces for consistent intending File Changes: Directory: /trunk/ ================== File [modified]: schedule.ccl Delta lines: +241 -241 =================================================================== --- trunk/schedule.ccl 2010-08-10 14:10:23 UTC (rev 1550) +++ trunk/schedule.ccl 2010-08-10 14:12:17 UTC (rev 1551) @@ -7,66 +7,66 @@ # # setup # -schedule AHFinderDirect_setup at CCTK_BASEGRID \ +schedule AHFinderDirect_setup at CCTK_BASEGRID \ after SpatialCoordinates - { - lang: C - options: global - } "setup data structures" + { + lang: C + options: global + } "setup data structures" # # checkpoint/recover support # schedule AHFinderDirect_recover at CCTK_POST_RECOVER_VARIABLES - { - lang: C - options: global - } "import horizon data from Cactus variables" + { + lang: C + options: global + } "import horizon data from Cactus variables" # # find horizons # if (run_at_CCTK_ANALYSIS != 0) - { - schedule AHFinderDirect_find_horizons at CCTK_ANALYSIS - { - lang: C - options: global - } "find apparent horizon(s) after this time step" - } + { + schedule AHFinderDirect_find_horizons at CCTK_ANALYSIS + { + lang: C + options: global + } "find apparent horizon(s) after this time step" + } if (run_at_CCTK_POSTSTEP != 0) - { - schedule AHFinderDirect_find_horizons at CCTK_POSTSTEP - { - lang: C - options: global - } "find apparent horizon(s) after this time step" - } + { + schedule AHFinderDirect_find_horizons at CCTK_POSTSTEP + { + lang: C + options: global + } "find apparent horizon(s) after this time step" + } if (run_at_CCTK_POSTINITIAL != 0) - { - schedule AHFinderDirect_find_horizons at CCTK_POSTINITIAL - { - lang: C - options: global - } "find apparent horizon(s) after this time step" - } + { + schedule AHFinderDirect_find_horizons at CCTK_POSTINITIAL + { + lang: C + options: global + } "find apparent horizon(s) after this time step" + } if (run_at_CCTK_POSTPOSTINITIAL != 0) - { - schedule AHFinderDirect_find_horizons at CCTK_POSTPOSTINITIAL - { - lang: C - options: global - } "find apparent horizon(s) after this time step" - } + { + schedule AHFinderDirect_find_horizons at CCTK_POSTPOSTINITIAL + { + lang: C + options: global + } "find apparent horizon(s) after this time step" + } if (run_at_CCTK_POST_RECOVER_VARIABLES != 0) - { - schedule AHFinderDirect_find_horizons at CCTK_POST_RECOVER_VARIABLES \ + { + schedule AHFinderDirect_find_horizons at CCTK_POST_RECOVER_VARIABLES \ after MoL_PostStep - { - lang: C - options: global - } "find apparent horizon(s) after this time step" - } + { + lang: C + options: global + } "find apparent horizon(s) after this time step" + } # # find horizons, set excision mask, maybe announce centroid to other thorns @@ -75,245 +75,245 @@ # AHFinderDirect_do_masks would be more apt # # We use two Cactus Groups here: -# group_for_mask_stuff -# group_where_mask_is_set -# AHFinderDirect_maybe_do_masks +# group_for_mask_stuff +# group_where_mask_is_set +# AHFinderDirect_maybe_do_masks # this way other thorns can schedule routines which need to run immediately # before/after we set the mask in group_for_mask_stuff before/after # group_where_mask_is_set, without having to know in what higher-level # Cactus schedule bins/groups this thorn runs. # if (run_at_CCTK_ANALYSIS != 0) - { - schedule AHFinderDirect_import_mask at CCTK_ANALYSIS \ - before AHFinderDirect_find_horizons - { - lang: C - options: global loop-local - } "import the excision mask" + { + schedule AHFinderDirect_import_mask at CCTK_ANALYSIS \ + before AHFinderDirect_find_horizons + { + lang: C + options: global loop-local + } "import the excision mask" - schedule AHFinderDirect_store at CCTK_ANALYSIS \ - after AHFinderDirect_find_horizons \ - before SphericalSurface_HasBeenSet - { - lang: C - options: global - } "store apparent horizon(s) into spherical surface(s)" + schedule AHFinderDirect_store at CCTK_ANALYSIS \ + after AHFinderDirect_find_horizons \ + before SphericalSurface_HasBeenSet + { + lang: C + options: global + } "store apparent horizon(s) into spherical surface(s)" - schedule AHFinderDirect_save at CCTK_ANALYSIS \ - after AHFinderDirect_find_horizons \ - { - lang: C - options: global - } "save apparent horizon(s) into Cactus variables" + schedule AHFinderDirect_save at CCTK_ANALYSIS \ + after AHFinderDirect_find_horizons \ + { + lang: C + options: global + } "save apparent horizon(s) into Cactus variables" - if (which_horizon_to_announce_centroid != 0) - { - schedule AHFinderDirect_announce at CCTK_ANALYSIS \ - before DriftCorrect \ - after AHFinderDirect_find_horizons - { - lang: C - options: global - } "announce horizon position(s) to other thorns" - } + if (which_horizon_to_announce_centroid != 0) + { + schedule AHFinderDirect_announce at CCTK_ANALYSIS \ + before DriftCorrect \ + after AHFinderDirect_find_horizons + { + lang: C + options: global + } "announce horizon position(s) to other thorns" + } - schedule AHFinderDirect_maybe_do_masks at CCTK_ANALYSIS \ - after AHFinderDirect_find_horizons - { - lang: C - options: global loop-local - } "set mask(s) based on apparent horizon position(s)" - } + schedule AHFinderDirect_maybe_do_masks at CCTK_ANALYSIS \ + after AHFinderDirect_find_horizons + { + lang: C + options: global loop-local + } "set mask(s) based on apparent horizon position(s)" + } if (run_at_CCTK_POSTSTEP != 0) - { - schedule AHFinderDirect_import_mask at CCTK_POSTSTEP \ - before AHFinderDirect_find_horizons - { - lang: C - options: global loop-local - } "import the excision mask" + { + schedule AHFinderDirect_import_mask at CCTK_POSTSTEP \ + before AHFinderDirect_find_horizons + { + lang: C + options: global loop-local + } "import the excision mask" - schedule AHFinderDirect_store at CCTK_POSTSTEP \ - after AHFinderDirect_find_horizons \ - before SphericalSurface_HasBeenSet - { - lang: C - options: global - } "store apparent horizon(s) into spherical surface(s)" + schedule AHFinderDirect_store at CCTK_POSTSTEP \ + after AHFinderDirect_find_horizons \ + before SphericalSurface_HasBeenSet + { + lang: C + options: global + } "store apparent horizon(s) into spherical surface(s)" - schedule AHFinderDirect_save at CCTK_POSTSTEP \ - after AHFinderDirect_find_horizons \ - { - lang: C - options: global - } "save apparent horizon(s) into Cactus variables" + schedule AHFinderDirect_save at CCTK_POSTSTEP \ + after AHFinderDirect_find_horizons \ + { + lang: C + options: global + } "save apparent horizon(s) into Cactus variables" - if (which_horizon_to_announce_centroid != 0) - { - schedule AHFinderDirect_announce at CCTK_POSTSTEP \ - before DriftCorrect \ - after AHFinderDirect_find_horizons - { - lang: C - options: global - } "announce horizon position(s) to other thorns" - } + if (which_horizon_to_announce_centroid != 0) + { + schedule AHFinderDirect_announce at CCTK_POSTSTEP \ + before DriftCorrect \ + after AHFinderDirect_find_horizons + { + lang: C + options: global + } "announce horizon position(s) to other thorns" + } - schedule AHFinderDirect_maybe_do_masks at CCTK_POSTSTEP \ - after AHFinderDirect_find_horizons - { - lang: C - options: global loop-local - } "set mask(s) based on apparent horizon position(s)" - } + schedule AHFinderDirect_maybe_do_masks at CCTK_POSTSTEP \ + after AHFinderDirect_find_horizons + { + lang: C + options: global loop-local + } "set mask(s) based on apparent horizon position(s)" + } if (run_at_CCTK_POSTINITIAL != 0) - { - schedule AHFinderDirect_import_mask at CCTK_POSTINITIAL \ - before AHFinderDirect_find_horizons - { - lang: C - options: global loop-local - } "import the excision mask" + { + schedule AHFinderDirect_import_mask at CCTK_POSTINITIAL \ + before AHFinderDirect_find_horizons + { + lang: C + options: global loop-local + } "import the excision mask" - schedule AHFinderDirect_store at CCTK_POSTINITIAL \ - after AHFinderDirect_find_horizons \ - before SphericalSurface_HasBeenSet - { - lang: C - options: global - } "store apparent horizon(s) into spherical surface(s)" + schedule AHFinderDirect_store at CCTK_POSTINITIAL \ + after AHFinderDirect_find_horizons \ + before SphericalSurface_HasBeenSet + { + lang: C + options: global + } "store apparent horizon(s) into spherical surface(s)" - schedule AHFinderDirect_save at CCTK_POSTINITIAL \ - after AHFinderDirect_find_horizons \ - { - lang: C - options: global - } "save apparent horizon(s) into Cactus variables" + schedule AHFinderDirect_save at CCTK_POSTINITIAL \ + after AHFinderDirect_find_horizons \ + { + lang: C + options: global + } "save apparent horizon(s) into Cactus variables" - if (which_horizon_to_announce_centroid != 0) - { - schedule AHFinderDirect_announce at CCTK_POSTINITIAL \ - before DriftCorrect \ - after AHFinderDirect_find_horizons - { - lang: C - options: global - } "announce horizon position(s) to other thorns" - } + if (which_horizon_to_announce_centroid != 0) + { + schedule AHFinderDirect_announce at CCTK_POSTINITIAL \ + before DriftCorrect \ + after AHFinderDirect_find_horizons + { + lang: C + options: global + } "announce horizon position(s) to other thorns" + } - schedule AHFinderDirect_maybe_do_masks at CCTK_POSTINITIAL \ - after AHFinderDirect_find_horizons - { - lang: C - options: global loop-local - } "set mask(s) based on apparent horizon position(s)" - } + schedule AHFinderDirect_maybe_do_masks at CCTK_POSTINITIAL \ + after AHFinderDirect_find_horizons + { + lang: C + options: global loop-local + } "set mask(s) based on apparent horizon position(s)" + } if (run_at_CCTK_POSTPOSTINITIAL != 0) - { - schedule AHFinderDirect_import_mask at CCTK_POSTPOSTINITIAL \ - before AHFinderDirect_find_horizons - { - lang: C - options: global loop-local - } "import the excision mask" + { + schedule AHFinderDirect_import_mask at CCTK_POSTPOSTINITIAL \ + before AHFinderDirect_find_horizons + { + lang: C + options: global loop-local + } "import the excision mask" - schedule AHFinderDirect_store at CCTK_POSTPOSTINITIAL \ - after AHFinderDirect_find_horizons \ - before SphericalSurface_HasBeenSet - { - lang: C - options: global - } "store apparent horizon(s) into spherical surface(s)" + schedule AHFinderDirect_store at CCTK_POSTPOSTINITIAL \ + after AHFinderDirect_find_horizons \ + before SphericalSurface_HasBeenSet + { + lang: C + options: global + } "store apparent horizon(s) into spherical surface(s)" - schedule AHFinderDirect_save at CCTK_POSTPOSTINITIAL \ - after AHFinderDirect_find_horizons \ - { - lang: C - options: global - } "save apparent horizon(s) into Cactus variables" + schedule AHFinderDirect_save at CCTK_POSTPOSTINITIAL \ + after AHFinderDirect_find_horizons \ + { + lang: C + options: global + } "save apparent horizon(s) into Cactus variables" - if (which_horizon_to_announce_centroid != 0) - { - schedule AHFinderDirect_announce at CCTK_POSTPOSTINITIAL \ - before DriftCorrect \ - after AHFinderDirect_find_horizons - { - lang: C - options: global - } "announce horizon position(s) to other thorns" - } + if (which_horizon_to_announce_centroid != 0) + { + schedule AHFinderDirect_announce at CCTK_POSTPOSTINITIAL \ + before DriftCorrect \ + after AHFinderDirect_find_horizons + { + lang: C + options: global + } "announce horizon position(s) to other thorns" + } - schedule AHFinderDirect_maybe_do_masks at CCTK_POSTPOSTINITIAL \ - after AHFinderDirect_find_horizons - { - lang: C - options: global loop-local - } "set mask(s) based on apparent horizon position(s)" - } + schedule AHFinderDirect_maybe_do_masks at CCTK_POSTPOSTINITIAL \ + after AHFinderDirect_find_horizons + { + lang: C + options: global loop-local + } "set mask(s) based on apparent horizon position(s)" + } if (run_at_CCTK_POST_RECOVER_VARIABLES != 0) - { - schedule AHFinderDirect_import_mask at CCTK_POST_RECOVER_VARIABLES \ - before AHFinderDirect_find_horizons - { - lang: C - options: global loop-local - } "import the excision mask" + { + schedule AHFinderDirect_import_mask at CCTK_POST_RECOVER_VARIABLES \ + before AHFinderDirect_find_horizons + { + lang: C + options: global loop-local + } "import the excision mask" - schedule AHFinderDirect_store at CCTK_POST_RECOVER_VARIABLES \ - after AHFinderDirect_find_horizons \ - before SphericalSurface_HasBeenSet - { - lang: C - options: global - } "store apparent horizon(s) into spherical surface(s)" + schedule AHFinderDirect_store at CCTK_POST_RECOVER_VARIABLES \ + after AHFinderDirect_find_horizons \ + before SphericalSurface_HasBeenSet + { + lang: C + options: global + } "store apparent horizon(s) into spherical surface(s)" - schedule AHFinderDirect_save at CCTK_POST_RECOVER_VARIABLES \ - after AHFinderDirect_find_horizons \ - { - lang: C - options: global - } "save apparent horizon(s) into Cactus variables" + schedule AHFinderDirect_save at CCTK_POST_RECOVER_VARIABLES \ + after AHFinderDirect_find_horizons \ + { + lang: C + options: global + } "save apparent horizon(s) into Cactus variables" - if (which_horizon_to_announce_centroid != 0) - { - schedule AHFinderDirect_announce at CCTK_POST_RECOVER_VARIABLES \ - before DriftCorrect \ - after AHFinderDirect_find_horizons - { - lang: C - options: global - } "announce horizon position(s) to other thorns" - } + if (which_horizon_to_announce_centroid != 0) + { + schedule AHFinderDirect_announce at CCTK_POST_RECOVER_VARIABLES \ + before DriftCorrect \ + after AHFinderDirect_find_horizons + { + lang: C + options: global + } "announce horizon position(s) to other thorns" + } - schedule AHFinderDirect_maybe_do_masks at CCTK_POST_RECOVER_VARIABLES \ - after AHFinderDirect_find_horizons - { - lang: C - options: global loop-local - } "set mask(s) based on apparent horizon position(s)" - } + schedule AHFinderDirect_maybe_do_masks at CCTK_POST_RECOVER_VARIABLES \ + after AHFinderDirect_find_horizons + { + lang: C + options: global loop-local + } "set mask(s) based on apparent horizon position(s)" + } # TODO: Don't call this, since no horizon has been found yet schedule AHFinderDirect_maybe_do_masks at CCTK_POSTREGRIDINITIAL \ after (MaskOne MaskZero) - { - lang: C - options: global loop-local - } "set mask(s) based on apparent horizon position(s)" + { + lang: C + options: global loop-local + } "set mask(s) based on apparent horizon position(s)" schedule AHFinderDirect_maybe_do_masks at CCTK_POSTREGRID \ after (MaskOne MaskZero) - { - lang: C - options: global loop-local - } "set mask(s) based on apparent horizon position(s)" + { + lang: C + options: global loop-local + } "set mask(s) based on apparent horizon position(s)" ######################################## ## ## uncomment this if you want to use the test driver -## src/util/test_patch_system.cc +## src/util/test_patch_system.cc ## ##schedule test_patch_system at CCTK_POSTINITIAL ##{

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits