[svn:einsteintoolkit] incoming/TOVSolverHot/src/ (Rev. 29) - Commits

4 Dec 2010

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)

    