26 #ifndef __GyotoKerrBL_H_
27 #define __GyotoKerrBL_H_
30 namespace Metric {
class KerrBL; }
36 #ifdef GYOTO_USE_XERCES
42 #define GYOTO_KERRBL_DEFAULT_DIFFTOL 1e-2
78 void
spin(const double spin);
84 void horizonSecurity(double drhor);
85 double horizonSecurity() const;
86 void genericIntegrator(bool);
87 bool genericIntegrator() const ;
93 void
gmunu(double g[4][4], const double * pos) const ;
94 double
gmunu(const double * const x, int mu, int nu) const ;
99 void
gmunu_up(double gup[4][4], const double * pos) const ;
100 double
gmunu_up(const double * const x, int mu, int nu) const ;
103 int christoffel(double dst[4][4][4], const double * pos) const ;
106 const double u1[4], const double u2[4]) const ;
108 void
nullifyCoord(double coord[8], double & tdot2) const;
114 double dir=1.) const ;
117 virtual void
MakeCoord(const double coordin[8], const double cst[5], double coordout[8]) const ;
121 void
MakeMomentum(const double coordin[8], const double cst[5], double coordout[8]) const;
124 virtual void
setParameter(std::string, std::string, std::string);
125 #ifdef GYOTO_USE_XERCES
136 int
myrk4(Worldline * line, const double coordin[8], double h, double res[8]) const;
139 int
myrk4(const double coor[8], const double cst[5], double h, double res[8]) const;
140 int
myrk4_adaptive(Gyoto::Worldline* line, const double coor[8], double lastnorm, double normref, double coor1[8], double h0, double& h1, double h1max) const;
146 int
CheckCons(const double coor_init[8], const double cst[5], double coor_fin[8]) const;
153 void
Normalize4v(double coord[8], const double part_mass) const;
161 virtual int
diff(const double y[8], const double cst[5],
162 double res[8]) const ;
165 virtual void
computeCst(const double coord[8], double cst[5]) const;
bool generic_integrator_
which integrator to use
Definition: GyotoKerrBL.h:67
Pointers performing reference counting.
Definition: GyotoSmartPointer.h:51
void gmunu(double g[4][4], const double *pos) const
Metric coefficients.
void nullifyCoord(double coord[8], double &tdot2) const
Set tdot (coord[4]) such that coord is light-like and return other possible tdot. ...
double drhor_
horizon security
Definition: GyotoKerrBL.h:66
Metric around a Kerr black-hole in Boyer-Lindquist coordinates.
Definition: GyotoKerrBL.h:48
virtual int isStopCondition(double const *const coord) const
Check whether integration should stop.
double getRmb() const
Returns prograde marginally bound orbit.
Base class for metrics.
Definition: GyotoMetric.h:147
double spin_
Angular momentum parameter.
Definition: GyotoKerrBL.h:54
double a2_
spin_*spin_
Definition: GyotoKerrBL.h:55
double difftol_
Numerical tuning parameter.
Definition: GyotoKerrBL.h:64
virtual void fillElement(FactoryMessenger *fmp)
called from Factory
virtual double christoffel(const double coord[8], const int alpha, const int mu, const int nu) const
Chistoffel symbol.
double a3_
a2_*spin_
Definition: GyotoKerrBL.h:56
Timelike or null geodesics.
double a4_
a2_*a2_
Definition: GyotoKerrBL.h:57
Base class for metric description.
virtual void setParameter(std::string, std::string, std::string)
Set parameter by name.
double difftol() const
Get difftol_.
void setParticleProperties(Worldline *line, const double *coord) const
Set Metric-specific constants of motion. Used e.g. in KerrBL.
int myrk4_adaptive(Gyoto::Worldline *line, const double coor[8], double lastnorm, double normref, double coor1[8], double h0, double &h1, double h1max) const
Internal-use adaptive RK4 proxy.
double rsink_
numerical horizon
Definition: GyotoKerrBL.h:65
Namespace for the Gyoto library.
Definition: GyotoAstrobj.h:42
int myrk4(Worldline *line, const double coordin[8], double h, double res[8]) const
RK4 integrator.
virtual void circularVelocity(double const pos[4], double vel[4], double dir=1.) const
Yield circular velocity at a given position.
virtual int diff(const double y[8], const double cst[5], double res[8]) const
Used in RK4 proxies.
virtual void MakeCoord(const double coordin[8], const double cst[5], double coordout[8]) const
Inverse function of MakeMomentumAndCst.
double getRms() const
Returns prograde marginally stable orbit.
void MakeMomentum(const double coordin[8], const double cst[5], double coordout[8]) const
Transforms from Boyer-Lindquist coordinates [t,r,th,phi,tdot,rdot,thdot,phidot] to [t...
int CheckCons(const double coor_init[8], const double cst[5], double coor_fin[8]) const
Ensure conservation of the constants of motion.
virtual void computeCst(const double coord[8], double cst[5]) const
virtual KerrBL * clone() const
Virtual copy constructor.
double ScalarProd(const double pos[4], const double u1[4], const double u2[4]) const
Scalar product.
void gmunu_up(double gup[4][4], const double *pos) const
gμ,ν
void Normalize4v(double coord[8], const double part_mass) const
Normalize 4-velocity.
double spin() const
Returns spin.
virtual int diff(const double y[8], double res[8]) const
F function such as dy/dtau=F(y,cst)