User: cott Date: 2010/08/18 01:41 PM
Modified: /EOS_Omni/ interface.ccl /EOS_Omni/src/ EOS_Omni_SingleVarCalls.F90
Log: * change interface, add scalar double precision variable
rf_precision
that allows the application thorn to specify the relative precision between iterations at which EOS root finding should be stopped.
File Changes:
Directory: /EOS_Omni/ =====================
File [modified]: interface.ccl Delta lines: +6 -0 =================================================================== --- EOS_Omni/interface.ccl 2010-08-16 17:52:05 UTC (rev 6) +++ EOS_Omni/interface.ccl 2010-08-18 18:41:07 UTC (rev 7) @@ -10,6 +10,7 @@
void FUNCTION EOS_Omni_press(CCTK_INT IN eoskey, \ CCTK_INT IN havetemp, \ + CCTK_REAL IN rf_precision, \ CCTK_INT IN npoints, \ CCTK_REAL IN ARRAY rho, \ CCTK_REAL INOUT ARRAY eps, \ @@ -24,6 +25,7 @@
void FUNCTION EOS_Omni_DPressByDEps(CCTK_INT IN eoskey, \ CCTK_INT IN havetemp, \ + CCTK_REAL IN rf_precision, \ CCTK_INT IN npoints, \ CCTK_REAL IN ARRAY rho, \ CCTK_REAL INOUT ARRAY eps, \ @@ -38,6 +40,7 @@
void FUNCTION EOS_Omni_DPressByDRho(CCTK_INT IN eoskey, \ CCTK_INT IN havetemp, \ + CCTK_REAL IN rf_precision, \ CCTK_INT IN npoints, \ CCTK_REAL IN ARRAY rho, \ CCTK_REAL INOUT ARRAY eps, \ @@ -52,6 +55,7 @@
void FUNCTION EOS_Omni_cs2(CCTK_INT IN eoskey, \ CCTK_INT IN havetemp, \ + CCTK_REAL IN rf_precision, \ CCTK_INT IN npoints, \ CCTK_REAL IN ARRAY rho, \ CCTK_REAL INOUT ARRAY eps, \ @@ -66,6 +70,7 @@
void FUNCTION EOS_Omni_EpsFromPress(CCTK_INT IN eoskey, \ CCTK_INT IN havetemp, \ + CCTK_REAL IN rf_precision, \ CCTK_INT IN npoints, \ CCTK_REAL IN ARRAY rho, \ CCTK_REAL INOUT ARRAY eps, \ @@ -81,6 +86,7 @@
void FUNCTION EOS_Omni_RestMassDensityFromEpsPress(CCTK_INT IN eoskey, \ CCTK_INT IN havetemp, \ + CCTK_REAL IN rf_precision, \ CCTK_INT IN npoints, \ CCTK_REAL OUT ARRAY rho, \ CCTK_REAL IN ARRAY eps, \
Directory: /EOS_Omni/src/ =========================
File [modified]: EOS_Omni_SingleVarCalls.F90 Delta lines: +13 -7 =================================================================== --- EOS_Omni/src/EOS_Omni_SingleVarCalls.F90 2010-08-16 17:52:05 UTC (rev 6) +++ EOS_Omni/src/EOS_Omni_SingleVarCalls.F90 2010-08-18 18:41:07 UTC (rev 7) @@ -9,7 +9,7 @@ ! 3 --- hybrid EOS ! 4 --- finite-T microphysical NSE EOS
-subroutine EOS_Omni_EOS_Press(eoskey,keytemp,npoints,& +subroutine EOS_Omni_EOS_Press(eoskey,keytemp,rf_precision,npoints,& rho,eps,temp,ye,press,keyerr,anyerr)
use EOS_Omni_Module @@ -19,6 +19,7 @@ CCTK_INT, intent(in) :: eoskey,keytemp,npoints CCTK_INT, intent(out) :: keyerr(npoints) CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rf_precision CCTK_REAL, intent(in) :: rho(npoints),ye(npoints) CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) CCTK_REAL, intent(out) :: press(npoints) @@ -92,7 +93,7 @@
-subroutine EOS_Omni_EOS_DPressByDRho(eoskey,keytemp,npoints,& +subroutine EOS_Omni_EOS_DPressByDRho(eoskey,keytemp,rf_precision,npoints,& rho,eps,temp,ye,dpdrhoe,keyerr,anyerr)
use EOS_Omni_Module @@ -102,6 +103,7 @@ CCTK_INT, intent(in) :: eoskey,keytemp,npoints CCTK_INT, intent(out) :: keyerr(npoints) CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rf_precision CCTK_REAL, intent(in) :: rho(npoints),ye(npoints) CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) CCTK_REAL, intent(out) :: dpdrhoe(npoints) @@ -179,7 +181,7 @@
end subroutine EOS_Omni_EOS_DPressByDRho
-subroutine EOS_Omni_EOS_DPressByDEps(eoskey,keytemp,npoints,& +subroutine EOS_Omni_EOS_DPressByDEps(eoskey,keytemp,rf_precision,npoints,& rho,eps,temp,ye,dpdepsrho,keyerr,anyerr)
use EOS_Omni_Module @@ -189,6 +191,7 @@ CCTK_INT, intent(in) :: eoskey,keytemp,npoints CCTK_INT, intent(out) :: keyerr(npoints) CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rf_precision CCTK_REAL, intent(in) :: rho(npoints),ye(npoints) CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) CCTK_REAL, intent(out) :: dpdepsrho(npoints) @@ -238,7 +241,7 @@
end subroutine EOS_Omni_EOS_DPressByDEps
-subroutine EOS_Omni_EOS_cs2(eoskey,keytemp,npoints,& +subroutine EOS_Omni_EOS_cs2(eoskey,keytemp,rf_precision,npoints,& rho,eps,temp,ye,cs2,keyerr,anyerr)
use EOS_Omni_Module @@ -248,6 +251,7 @@ CCTK_INT, intent(in) :: eoskey,keytemp,npoints CCTK_INT, intent(out) :: keyerr(npoints) CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rf_precision CCTK_REAL, intent(in) :: rho(npoints),ye(npoints) CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) CCTK_REAL, intent(out) :: cs2(npoints) @@ -328,7 +332,7 @@ end subroutine EOS_Omni_EOS_cs2
-subroutine EOS_Omni_EOS_eps_from_press(eoskey,keytemp,npoints,& +subroutine EOS_Omni_EOS_eps_from_press(eoskey,keytemp,rf_precision,npoints,& rho,eps,temp,ye,press,xeps,keyerr,anyerr)
use EOS_Omni_Module @@ -338,6 +342,7 @@ CCTK_INT, intent(in) :: eoskey,keytemp,npoints CCTK_INT, intent(out) :: keyerr(npoints) CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rf_precision CCTK_REAL, intent(in) :: rho(npoints),ye(npoints),press(npoints) CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) CCTK_REAL, intent(out) :: xeps(npoints) @@ -377,8 +382,8 @@ end subroutine EOS_Omni_EOS_eps_from_press
-subroutine EOS_Omni_EOS_RestMassDensityFromEpsPress(eoskey,keytemp,npoints,& - rho,eps,temp,ye,press,keyerr,anyerr) +subroutine EOS_Omni_EOS_RestMassDensityFromEpsPress(eoskey,keytemp,rf_precision,& + npoints,rho,eps,temp,ye,press,keyerr,anyerr)
use EOS_Omni_Module implicit none @@ -387,6 +392,7 @@ CCTK_INT, intent(in) :: eoskey,keytemp,npoints CCTK_INT, intent(out) :: keyerr(npoints) CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rf_precision CCTK_REAL, intent(in) :: ye(npoints),press(npoints),eps(npoints) CCTK_REAL, intent(inout) :: rho(npoints),temp(npoints)