[Commits] [svn:einsteintoolkit] GRHydro/trunk/ (Rev. 605)

rhaas at tapir.caltech.edu rhaas at tapir.caltech.edu
Tue Apr 15 14:48:28 CDT 2014


User: rhaas
Date: 2014/04/15 02:48 PM

Modified:
 /trunk/
  interface.ccl
 /trunk/src/
  GRHydro_Prim2Con.c

Log:
 GRHydro: start changes for improved EOS and HLLE handling

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_Prim2Con.c
Delta lines: +16 -16
===================================================================
--- trunk/src/GRHydro_Prim2Con.c	2014-04-15 19:48:25 UTC (rev 604)
+++ trunk/src/GRHydro_Prim2Con.c	2014-04-15 19:48:28 UTC (rev 605)
@@ -79,11 +79,11 @@
        for(int k=0;k<cctk_ash[2];k++) 
 	 for(int j=0;j<cctk_ash[1];j++) {
 	   int i = CCTK_GFINDEX3D(cctkGH,0,j,k);
-	   EOS_Omni_press(*GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,cctk_ash[0],
-			  &(rhominus[i]),&(epsminus[i]),NULL,NULL,&(pressminus[i]),
+	   EOS_Omni_press_cs2(*GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,cctk_ash[0],
+			  &(rhominus[i]),&(epsminus[i]),NULL,NULL,&(pressminus[i]),&(cs2minus[i]),
 			  &(keyerr[i]),&anyerr);
-	   EOS_Omni_press(*GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,cctk_ash[0],
-			  &(rhoplus[i]),&(epsplus[i]),NULL,NULL,&(pressplus[i]),
+	   EOS_Omni_press_cs2(*GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,cctk_ash[0],
+			  &(rhoplus[i]),&(epsplus[i]),NULL,NULL,&(pressplus[i]),&(cs2plus[i]),
 			  &(keyerr[i]),&anyerr);
 	 }
      free(keyerr);
@@ -108,9 +108,9 @@
        for(int k=0;k<cctk_ash[2];k++) 
 	 for(int j=0;j<cctk_ash[1];j++) {
 	   int i = CCTK_GFINDEX3D(cctkGH,0,j,k);
-	   EOS_Omni_press(*GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,cctk_ash[0],
+	   EOS_Omni_press_cs2(*GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,cctk_ash[0],
 			  &rhominus[i],&epsminus[i],&tempminus[i],&Y_e_minus[i],
-			  &pressminus[i],
+			  &pressminus[i],&cs2minus[i],
 			  &keyerr[i],&anyerr);
 	 }
        
@@ -137,9 +137,9 @@
        for(int k=0;k<cctk_ash[2];k++) 
 	 for(int j=0;j<cctk_ash[1];j++) {
 	   int i = CCTK_GFINDEX3D(cctkGH,0,j,k);
-	   EOS_Omni_press(*GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,cctk_ash[0],
+	   EOS_Omni_press_cs2(*GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,cctk_ash[0],
 			  &rhoplus[i],&epsplus[i],&tempplus[i],&Y_e_plus[i],
-			  &pressplus[i],
+			  &pressplus[i],&cs2plus[i],
 			  &keyerr[i],&anyerr);
 	 }
 
@@ -185,9 +185,9 @@
        for(int k=0;k<cctk_ash[2];k++) 
 	 for(int j=0;j<cctk_ash[1];j++) {
 	   int i = CCTK_GFINDEX3D(cctkGH,0,j,k);
-	   EOS_Omni_press(*GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,cctk_ash[0],
+	   EOS_Omni_press_cs2(*GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,cctk_ash[0],
 			  &rhominus[i],&epsminus[i],&tempminus[i],&Y_e_minus[i],
-			  &pressminus[i],
+			  &pressminus[i],&cs2minus[i],
 			  &keyerr[i],&anyerr);
 	 }
 
@@ -212,9 +212,9 @@
 		 int lkeyerr[1];
 		 int lanyerr = 0;
 		 int lkeytemp = 1;
-		 EOS_Omni_press(*GRHydro_eos_handle,lkeytemp,GRHydro_eos_rf_prec,ln,
+		 EOS_Omni_press_cs2(*GRHydro_eos_handle,lkeytemp,GRHydro_eos_rf_prec,ln,
 				&rhominus[i],&epsminus[i],&tempminus[i],
-				&Y_e_minus[i],&pressminus[i],lkeyerr,&lanyerr);
+				&Y_e_minus[i],&pressminus[i],&cs2minus[i],lkeyerr,&lanyerr);
 		 if(lanyerr !=0) {
 		   CCTK_VWarn(0, __LINE__, __FILE__, CCTK_THORNSTRING,
 			      "Aborting! keyerr=%d, r=%15.6E, t=%15.6E, ye=%15.6E",
@@ -233,9 +233,9 @@
        for(int k=0;k<cctk_ash[2];k++) 
 	 for(int j=0;j<cctk_ash[1];j++) {
 	   int i = CCTK_GFINDEX3D(cctkGH,0,j,k);
-	   EOS_Omni_press(*GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,cctk_ash[0],
+	   EOS_Omni_press_cs2(*GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,cctk_ash[0],
 			  &rhoplus[i],&epsplus[i],&tempplus[i],&Y_e_plus[i],
-			  &pressplus[i],
+			  &pressplus[i],&cs2plus[i],
 			  &keyerr[i],&anyerr);
 	 }
 
@@ -259,9 +259,9 @@
 		 int lkeyerr[1];
 		 int lanyerr = 0;
 		 int lkeytemp = 1;
-		 EOS_Omni_press(*GRHydro_eos_handle,lkeytemp,GRHydro_eos_rf_prec,ln,
+		 EOS_Omni_press_cs2(*GRHydro_eos_handle,lkeytemp,GRHydro_eos_rf_prec,ln,
 				&rhoplus[i],&epsplus[i],&tempplus[i],
-				&Y_e_plus[i],&pressplus[i],lkeyerr,&lanyerr);
+				&Y_e_plus[i],&pressplus[i],&cs2plus[i],lkeyerr,&lanyerr);
 		 if(lanyerr !=0) {
 		   CCTK_VWarn(0, __LINE__, __FILE__, CCTK_THORNSTRING,
 			      "Aborting! keyerr=%d, r=%15.6E, t=%15.6E, ye=%15.6E",

Directory: /trunk/
==================

File [modified]: interface.ccl
Delta lines: +17 -2
===================================================================
--- trunk/interface.ccl	2014-04-15 19:48:25 UTC (rev 604)
+++ trunk/interface.ccl	2014-04-15 19:48:28 UTC (rev 605)
@@ -286,7 +286,22 @@
 
 USES FUNCTION EOS_Omni_press
 
+void FUNCTION EOS_Omni_press_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,  \
+			     CCTK_REAL INOUT ARRAY temp, \
+			     CCTK_REAL IN ARRAY ye,      \
+			     CCTK_REAL OUT ARRAY press,  \
+			     CCTK_REAL OUT ARRAY cs2,  \
+			     CCTK_INT OUT ARRAY keyerr,  \
+			     CCTK_INT OUT anyerr)
 
+USES FUNCTION EOS_Omni_press_cs2
+
+
 void FUNCTION EOS_Omni_DPressByDEps(CCTK_INT IN eoskey,      \
 			     CCTK_INT IN havetemp,           \
 			     CCTK_REAL IN rf_precision,       \
@@ -455,8 +470,8 @@
 real GRHydro_prim_bext type = GF Timelevels = 1 tags='Prolongation="None" checkpoint="no"'
 {
   rhoplus, velxplus, velyplus, velzplus, pressplus, epsplus, w_lorentzplus, 
-  rhominus, velxminus, velyminus, velzminus, pressminus, epsminus, 
-  w_lorentzminus
+  cs2plus, rhominus, velxminus, velyminus, velzminus, pressminus, epsminus, 
+  w_lorentzminus, cs2minus
 } "Primitive variables extended to the cell boundaries"
 
 int GRHydro_scalars type = SCALAR tags='checkpoint="no"'



More information about the Commits mailing list