diff options
-rw-r--r-- | src/tov.c | 5 | ||||
-rw-r--r-- | test/test_one_boost_max.par | 4 | ||||
-rw-r--r-- | test/test_one_static_max.par | 4 | ||||
-rw-r--r-- | test/test_tov_carpet.par | 4 | ||||
-rw-r--r-- | test/test_two_av.par | 8 | ||||
-rw-r--r-- | test/test_two_max.par | 8 |
6 files changed, 19 insertions, 14 deletions
@@ -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) diff --git a/test/test_one_boost_max.par b/test/test_one_boost_max.par index d949733..6dd2846 100644 --- a/test/test_one_boost_max.par +++ b/test/test_one_boost_max.par @@ -41,8 +41,8 @@ ADMBase::initial_lapse = "tov" ADMBase::initial_shift = "tov" tovsolver::TOV_Rho_Central[0] = 1.28e-3 -tovsolver::TOV_Gamma[0] = 2.0 -tovsolver::TOV_K[0] = 100.0 +tovsolver::TOV_Gamma = 2.0 +tovsolver::TOV_K = 100.0 tovsolver::TOV_Velocity_x[0] = 0.1 tovsolver::TOV_Velocity_y[0] = 0.2 tovsolver::TOV_Velocity_z[0] = 0.3 diff --git a/test/test_one_static_max.par b/test/test_one_static_max.par index c65b2e7..6cb4b5b 100644 --- a/test/test_one_static_max.par +++ b/test/test_one_static_max.par @@ -41,8 +41,8 @@ ADMBase::initial_lapse = "tov" ADMBase::initial_shift = "tov" tovsolver::TOV_Rho_Central[0] = 1.28e-3 -tovsolver::TOV_Gamma[0] = 2.0 -tovsolver::TOV_K[0] = 100.0 +tovsolver::TOV_Gamma = 2.0 +tovsolver::TOV_K = 100.0 tovsolver::TOV_Combine_Method = "maximum" IO::out_dir = $parfile diff --git a/test/test_tov_carpet.par b/test/test_tov_carpet.par index ab4c26e..4073448 100644 --- a/test/test_tov_carpet.par +++ b/test/test_tov_carpet.par @@ -52,8 +52,8 @@ ADMBase::initial_lapse = "tov" ADMBase::initial_shift = "tov" tovsolver::TOV_Rho_Central[0] = 1.28e-3 -tovsolver::TOV_Gamma[0] = 2.0 -tovsolver::TOV_K[0] = 100.0 +tovsolver::TOV_Gamma = 2.0 +tovsolver::TOV_K = 100.0 tovsolver::TOV_Combine_Method = "maximum" Carpet::prolongation_order_space = 3 diff --git a/test/test_two_av.par b/test/test_two_av.par index 0b5725e..8f074b0 100644 --- a/test/test_two_av.par +++ b/test/test_two_av.par @@ -45,13 +45,13 @@ tovsolver::TOV_Num_Radial = 200000 tovsolver::TOV_Combine_Method = "average" tovsolver::TOV_Rho_Central[0] = 0.16e-3 -tovsolver::TOV_Gamma[0] = 2.0 -tovsolver::TOV_K[0] = 100.0 +tovsolver::TOV_Gamma = 2.0 +tovsolver::TOV_K = 100.0 tovsolver::TOV_Position_x[0] = -15.0 tovsolver::TOV_Rho_Central[1] = 0.32e-3 -tovsolver::TOV_Gamma[1] = 2.0 -tovsolver::TOV_K[1] = 100.0 +tovsolver::TOV_Gamma = 2.0 +tovsolver::TOV_K = 100.0 tovsolver::TOV_Position_x[1] = 15.0 IO::out_fileinfo="none" diff --git a/test/test_two_max.par b/test/test_two_max.par index 39761a3..4a2819a 100644 --- a/test/test_two_max.par +++ b/test/test_two_max.par @@ -45,13 +45,13 @@ tovsolver::TOV_Num_Radial = 200000 tovsolver::TOV_Combine_Method = "maximum" tovsolver::TOV_Rho_Central[0] = 0.16e-3 -tovsolver::TOV_Gamma[0] = 2.0 -tovsolver::TOV_K[0] = 100.0 +tovsolver::TOV_Gamma = 2.0 +tovsolver::TOV_K = 100.0 tovsolver::TOV_Position_x[0] = -15.0 tovsolver::TOV_Rho_Central[1] = 0.32e-3 -tovsolver::TOV_Gamma[1] = 2.0 -tovsolver::TOV_K[1] = 100.0 +tovsolver::TOV_Gamma = 2.0 +tovsolver::TOV_K = 100.0 tovsolver::TOV_Position_x[1] = 15.0 IO::out_fileinfo="none" |