[Commits] [svn:einsteintoolkit] incoming/EOS_Omni/ (Rev. 7)

cott at tapir.caltech.edu cott at tapir.caltech.edu
Wed Aug 18 13:41:07 CDT 2010


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)
 



More information about the Commits mailing list