User: cott Date: 2010/12/04 07:19 PM
Modified: /TOVSolverHot/src/ tov.c
Log: * small mapping fix
File Changes:
Directory: /TOVSolverHot/src/ =============================
File [modified]: tov.c Delta lines: +5 -0 =================================================================== --- TOVSolverHot/src/tov.c 2010-12-05 00:44:50 UTC (rev 28) +++ TOVSolverHot/src/tov.c 2010-12-05 01:19:49 UTC (rev 29) @@ -848,6 +848,8 @@ } else { keytemp = 0; } + // initial guess + rho_point[star] = 0.001e0; EOS_Omni_RhoFromPressEpsTempEnt(eoskey, keytemp, TOVSolverHot_eos_rf_prec, @@ -1058,6 +1060,7 @@ (r_to_star[star_i] * r_to_star[star_i] + 1.0e-30)) / my_psi4; rho[i3D] += rho_point[star_i]; + eps[i3D] += eps_point[star_i]; press[i3D] += press_point[star_i]; if(eoskey==4) { @@ -1077,6 +1080,7 @@ /* Reset atmosphere according to chosen star */ if(rho[i3D] <= TOV_Atmosphere[star_i]) rho[i3D] = TOV_Atmosphere[star_i]; + }
if (TOV_Conformal_Flat_Three_Metric) @@ -1141,6 +1145,7 @@ Y_e_con[i3D] = Y_e[i3D] * dens[i3D]; }
+ } /* if used, recalculate the derivatives of the conformal factor */ if (*conformal_state > 1)