diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2014-04-15 19:48:28 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2014-04-15 19:48:28 +0000 |
| commit | 42214f960ac8c1d2d62bf49d587638c7e4021664 (patch) | |
| tree | 26dafb7f630d2ed4ec27771bcdedd134e45613b1 /src/GRHydro_Prim2Con.c | |
| parent | 09bde56fc60492d2f51385680efb66c97be59782 (diff) | |
GRHydro: start changes for improved EOS and HLLE handling
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@605 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Prim2Con.c')
| -rw-r--r-- | src/GRHydro_Prim2Con.c | 32 |
1 files changed, 16 insertions, 16 deletions
diff --git a/src/GRHydro_Prim2Con.c b/src/GRHydro_Prim2Con.c index 4f206d1..16a8476 100644 --- a/src/GRHydro_Prim2Con.c +++ b/src/GRHydro_Prim2Con.c @@ -79,11 +79,11 @@ void Primitive2ConservativeC(CCTK_ARGUMENTS) 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 @@ void Primitive2ConservativeC(CCTK_ARGUMENTS) 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 @@ void Primitive2ConservativeC(CCTK_ARGUMENTS) 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 @@ void Primitive2ConservativeC(CCTK_ARGUMENTS) 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 @@ void Primitive2ConservativeC(CCTK_ARGUMENTS) 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 @@ void Primitive2ConservativeC(CCTK_ARGUMENTS) 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 @@ void Primitive2ConservativeC(CCTK_ARGUMENTS) 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", |