28 #ifndef __GyotoRotStar3_1_H_
29 #define __GyotoRotStar3_1_H_
37 namespace Metric {
class RotStar3_1; }
44 #ifdef GYOTO_USE_XERCES
82 int
myrk4(const double coord[6], double h, double res[6]) const;
90 int
myrk4_adaptive(Gyoto::Worldline* line, const double coord[8], double lastnorm, double normref, double coordnew[8], double h0, double& h1, double h1max) const;
97 int
myrk4_adaptive(const double coor[6], double lastnorm, double normref, double coornew[6], double cst[2], double& tdot_used, double h0, double& h1, double h1max, double& hused) const;
102 int
diff(const double coord[8], double res[8]) const ;
111 int
diff(const double y[6], double res[6], int) const ;
117 void
Normalize4v(const double coordin[6], double coordout[6], const double cst[2], double& tdot_used) const;
119 double
gmunu(const double * x, int mu, int nu) const ;
121 double
christoffel(const double coord[8], const int alpha, const int mu,
122 const int nu) const ;
125 const double u1[4], const double u2[4]) const ;
127 virtual void
setParameter(std::string, std::string, std::string);
128 #ifdef GYOTO_USE_XERCES
int myrk4(const double coord[6], double h, double res[6]) const
RK4 integrator.
Pointers performing reference counting.
Definition: GyotoSmartPointer.h:51
int integ_kind_
1 if RotStar3_1::myrk4(), 0 if Metric::myrk4()
Definition: GyotoRotStar3_1.h:58
Numerical metric around a rotating star in 3+1 formalism.
Definition: GyotoRotStar3_1.h:52
Base class for metrics.
Definition: GyotoMetric.h:147
Reference-counting pointers.
int myrk4_adaptive(Gyoto::Worldline *line, const double coord[8], double lastnorm, double normref, double coordnew[8], double h0, double &h1, double h1max) const
Adaptive RK4 integrator.
char * filename_
Lorene output file name.
Definition: GyotoRotStar3_1.h:56
virtual void setParameter(std::string, std::string, std::string)
Set parameter by name.
void Normalize4v(const double coordin[6], double coordout[6], const double cst[2], double &tdot_used) const
Tweak coordinates to insure conservation of cst.
double ScalarProd(const double pos[4], const double u1[4], const double u2[4]) const
Scalar product.
virtual void fillElement(FactoryMessenger *fmp)
called from Factory
Star_rot * star_
Pointer to underlying Lorene Star_rot instance.
Definition: GyotoRotStar3_1.h:57
Timelike or null geodesics.
Base class for metric description.
Namespace for the Gyoto library.
Definition: GyotoAstrobj.h:42
virtual void setParameters(Gyoto::FactoryMessenger *fmp)
Main loop in Subcontractor_t function.
double christoffel(const double coord[8], const int alpha, const int mu, const int nu) const
Chistoffel symbol.
int integKind() const
Get integ_kind_.
char const * fileName() const
Get filename_.
double gmunu(const double *x, int mu, int nu) const
Metric coefficients.
virtual int myrk4(Worldline *line, const double coord[8], double h, double res[8]) const
RK4 integrator.
int diff(const double coord[8], double res[8]) const
F function such as dy/dtau=F(y,cst)
virtual RotStar3_1 * clone() const
Cloner (uses RotStar3_1(file, integ_kind))