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

[cott at tapir.caltech.edu]

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)