* remove option to have different gamma and K for different

  stars. That's physically impossible.
* adjust test par files (will do so with GRHydro next)
* modify the way solutions are added so that atmosphere cuts
  are applied after solutions are added. Otherwise one of the
  stars will have the atmosphere added to it. This breaks symmetry.



git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TOVSolver/trunk@137 1bdb13ef-5d69-4035-bb54-08abeb3aa7f1

1 files changed, 5 insertions, 0 deletions

diff --git a/src/tov.c b/src/tov.c
index 7fbbcd8..fa65e58 100644
--- a/src/tov.c
+++ b/src/tov.c
@@ -829,6 +829,10 @@ void TOV_C_Exact(CCTK_ARGUMENTS)
                          (r_to_star[star_i] * r_to_star[star_i] + 1.0e-30)) /
                         my_psi4;
           rho[i3D] += rho_point[star_i];
+	  if( fabs(x[i3D]-15.0) <= 1.0e-10 ||
+	      fabs(x[i3D]+15.0) <= 1.0e-10 ) {
+	    fprintf(stderr,"%22.14E %22.14E\n",x[i3D],rho[i3D]);
+	  }
           eps[i3D] += eps_point[star_i];
           press[i3D] += press_point[star_i];
           /* we still have to know if we are inside one star - and which */
@@ -904,6 +908,7 @@ void TOV_C_Exact(CCTK_ARGUMENTS)
                                          (1.0 + eps[i3D]) +
                               press[i3D]*(w_lorentz[i3D]*w_lorentz[i3D]-1.0) ) -
                  dens[i3D];
+      abort();
     }
   /* if used, recalculate the derivatives of the conformal factor */
   if (*conformal_state > 1)