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);
227 std::shared_ptr<o2scl::table_units<> >
tab;
237 void calc_y(
double &yR,
double &beta,
double &k2,
double &lambda_km5,
238 double &lambda_cgs,
bool tabulate=
false);
254 void calc_H(
double &yR,
double &beta,
double &k2,
double &lambda_km5,
259 #ifndef DOXYGEN_NO_O2NS
int y_derivs(double r, size_t nv, const std::vector< double > &vals, std::vector< double > &ders)
The derivative .
double delta
The radial step for resolving discontinuities in km (default )
void calc_H(double &yR, double &beta, double &k2, double &lambda_km5, double &lambda_cgs)
Compute the love number using H.
std::vector< double > disc
List of discontinuities.
double schwarz_km
Schwarzchild radius in km (set in constructor)
void calc_y(double &yR, double &beta, double &k2, double &lambda_km5, double &lambda_cgs, bool tabulate=false)
Compute the love number using y.
o2scl::ode_iv_solve< ode_funct2, std::vector< double > > * oisp
The ODE integrator.
void set_ODE(o2scl::ode_iv_solve< ode_funct2, std::vector< double > > &ois_new)
Set ODE integrator.
double eval_k2(double beta, double yR)
Compute using the analytic expression.
Determination of the neutron star Love number.
void add_disc(double rd)
Add a discontinuity at radius rd (in km)
int H_derivs(double r, size_t nv, const std::vector< double > &vals, std::vector< double > &ders)
The derivatives and .
double eps
The first radial point in (default 0.02)
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.
void clear_discs()
Remove all discontinuities.
o2scl::ode_iv_solve< ode_funct2, std::vector< double > > def_ois
The default ODE integrator.