26 #include <gsl/gsl_math.h>
28 #include <boost/numeric/ublas/vector.hpp>
30 #include <o2scl/astep_gsl.h>
31 #include <o2scl/table_units.h>
32 #include <o2scl/astep_nonadapt.h>
33 #include <o2scl/ode_rk8pd_gsl.h>
34 #include <o2scl/ode_iv_solve.h>
36 #ifndef DOXYGEN_NO_O2NS
170 typedef std::function<int(
double,
size_t,
171 const std::vector<double> &,
172 std::vector<double> &)> ode_funct2;
174 #ifndef DOXYGEN_INTERNAL
183 int y_derivs(
double r,
size_t nv,
const std::vector<double> &vals,
184 std::vector<double> &ders);
189 int H_derivs(
double r,
size_t nv,
const std::vector<double> &vals,
190 std::vector<double> &ders);
201 double eval_k2(
double beta,
double yR);
240 std::shared_ptr<o2scl::table_units<> >
tab;
250 int calc_y(
double &yR,
double &beta,
double &k2,
double &lambda_km5,
251 double &lambda_cgs,
bool tabulate=
false);
267 int calc_H(
double &yR,
double &beta,
double &k2,
double &lambda_km5,
272 #ifndef DOXYGEN_NO_O2NS
o2scl::ode_iv_solve< ode_funct2, std::vector< double > > def_ois
The default ODE integrator.
int y_derivs(double r, size_t nv, const std::vector< double > &vals, std::vector< double > &ders)
The derivative .
int calc_y(double &yR, double &beta, double &k2, double &lambda_km5, double &lambda_cgs, bool tabulate=false)
Compute the love number using y.
double delta
The radial step for resolving discontinuities in km (default )
void add_disc(double rd)
Add a discontinuity at radius rd (in km)
double eval_k2(double beta, double yR)
Compute using the analytic expression.
int show_ode
If greater than zero, show the ODE output (default 0)
bool addl_testing
Additional testing if the ODE solver fails.
std::vector< double > disc
List of discontinuities.
void set_ODE(o2scl::ode_iv_solve< ode_funct2, std::vector< double > > &ois_new)
Set ODE integrator.
int H_derivs(double r, size_t nv, const std::vector< double > &vals, std::vector< double > &ders)
The derivatives and .
o2scl::ode_iv_solve< ode_funct2, std::vector< double > > * oisp
The ODE integrator.
o2scl::table_units results
A table containing the solution to the differential equation(s)
double eps
The first radial point in (default 0.02)
bool err_nonconv
If true, call the error handler if the solution does not converge (default true)
int calc_H(double &yR, double &beta, double &k2, double &lambda_km5, double &lambda_cgs)
Compute the love number using H.
std::shared_ptr< o2scl::table_units<> > tab
Pointer to the input profile.
void clear_discs()
Remove all discontinuities.
Determination of the neutron star Love number.
double schwarz_km
Schwarzchild radius in km (set in constructor)
Documentation generated with Doxygen. Provided under the
GNU Free Documentation License (see License Information).