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 146 #ifndef DOXYGEN_INTERNAL 158 int y_derivs(
double r,
size_t nv,
const ubvector &vals,
164 int H_derivs(
double r,
size_t nv,
const ubvector &vals,
176 double eval_k2(
double beta,
double yR);
191 std::shared_ptr<o2scl::table_units<> >
tab;
200 void calc_y(
double &yR,
double &beta,
double &k2,
double &lambda_km5,
201 double &lambda_cgs,
bool tabulate=
false);
205 void calc_H(
double &yR,
double &beta,
double &k2,
double &lambda_km5,
210 #ifndef DOXYGEN_NO_O2NS
void calc_H(double &yR, double &beta, double &k2, double &lambda_km5, double &lambda_cgs)
Compute the love number using H.
double schwarz_km
Schwarzchild radius in km (set in constructor)
int H_derivs(double r, size_t nv, const ubvector &vals, ubvector &ders)
The derivatives and .
void set_ODE(o2scl::ode_iv_solve<> &ois_new)
Set ODE integrator.
void calc_y(double &yR, double &beta, double &k2, double &lambda_km5, double &lambda_cgs, bool tabulate=false)
Compute the love number using y.
double eval_k2(double beta, double yR)
Compute using the analytic expression.
int y_derivs(double r, size_t nv, const ubvector &vals, ubvector &ders)
The derivative .
Determination of the neutron star Love number.
o2scl::ode_iv_solve def_ois
The default ODE integrator.
double eps
The first radial point in (default 0.02)
o2scl::ode_iv_solve * oisp
The ODE integrator.
o2scl::table_units results
A table containing the solution to the differential equation(s)
std::shared_ptr< o2scl::table_units<> > tab
Pointer to the input profile.