14 #ifndef __GyotoNumericalMetricLoreneMetric_H_ 15 #define __GyotoNumericalMetricLoreneMetric_H_ 21 namespace Metric {
class NumericalMetricLorene; }
22 class FactoryMessenger;
37 #ifdef GYOTO_USE_XERCES 55 bool bosonstarcircular_;
57 int has_acceleration_vector_;
58 int specify_marginalorbits_;
65 Lorene::Scalar** lapse_tab_;
66 Lorene::Vector** shift_tab_;
67 Lorene::Sym_tensor** gamcov_tab_;
68 Lorene::Sym_tensor** gamcon_tab_;
69 Lorene::Sym_tensor** kij_tab_;
72 Lorene::Valeur** nssurf_tab_;
73 Lorene::Vector** vsurf_tab_;
74 Lorene::Vector** accel_tab_;
75 Lorene::Scalar** lorentz_tab_;
76 Lorene::Valeur** hor_tab_;
94 virtual void setMetricSource();
96 void directory(std::string const &dir) ;
97 std::string directory() const ;
98 double initialTime() const ;
99 void initialTime(double t0);
100 double horizon() const ;
101 void horizon(double t0);
102 double rico() const ;
103 void rico(double r0);
104 bool hasSurface() const;
105 void hasSurface(bool s);
106 bool hasAccelerationVector() const;
107 void hasAccelerationVector(bool aa);
108 bool bosonstarcircular() const;
109 void bosonstarcircular(bool);
110 bool specifyMarginalOrbits() const;
111 void specifyMarginalOrbits(bool s);
114 std::vector<double> refineIntegStep() const;
115 void refineIntegStep(std::vector<double> const&);
117 Lorene::Vector** getShift_tab() const;
118 Lorene::Scalar** getLapse_tab() const;
119 Lorene::Sym_tensor** getGamcon_tab() const;
120 Lorene::Sym_tensor** getGamcov_tab() const;
121 double* getTimes() const;
122 int getNbtimes() const;
123 Lorene::Valeur** getNssurf_tab() const;
124 Lorene::Vector** getVsurf_tab() const;
125 Lorene::Vector** getAccel_tab() const;
126 Lorene::Scalar** getLorentz_tab() const;
127 Lorene::Valeur** getHor_tab() const;
128 double getRms() const;
129 double getRmb() const;
130 void setLapse_tab(Lorene::Scalar* lapse, int ii);
131 void setShift_tab(Lorene::Vector* shift, int ii);
132 void setGamcov_tab(Lorene::Sym_tensor* gamcov, int ii);
133 void setGamcon_tab(Lorene::Sym_tensor* gamcon, int ii);
134 void setKij_tab(Lorene::Sym_tensor* kij, int ii);
135 void setTimes(double time,int ii);
137 virtual double getSpecificAngularMomentum(double rr) const;
138 virtual double getPotential(double const pos[4], double l_cst) const;
143 virtual int myrk4(double tt, const double coord[7], double h, double res[7]) const;
144 virtual int myrk4(Worldline* line, state_t const &coord,
145 double h, state_t &res) const;
150 int myrk4_adaptive(Gyoto::Worldline* line, state_t const &coord, double lastnorm, double normref, state_t &coordnew, double h0, double& h1, double h1max) const;
152 int myrk4_adaptive(double tt, const double coor[7], double lastnorm, double normref, double coornew[7], const double cst[2], double& tdot_used, double h0, double& h1, double& hused, double h1max) const;
158 void reverseR(double tt, double coord[7]) const;
163 void computeNBeta(const double coord[4],double &NN,double beta[3]) const;
169 double gmunu(const double x[4], int mu, int nu) const ;
171 double gmunu(const double x[3], int indice_time, int mu, int nu) const ;
176 double gmunu_up_dr(const double x[4], int mu, int nu) const ;
178 double gmunu_up_dr(const double x[3], int indice_time, int mu, int nu) const ;
180 double christoffel(const double coord[4], const int alpha, const int mu,
181 const int nu) const ;
182 double christoffel(const double coord[4],
184 const int mu, const int nu,
185 const int indice_time) const;
186 virtual int christoffel(double dst[4][4][4],
187 const double coord[4]) const;
188 int christoffel(double dst[4][4][4],
189 const double coord[4],
190 const int indice_time) const;
194 double christoffel3(const double coord[3], const int indice_time,
195 const int ii, const int jj, const int kk) const ;
197 void setParticleProperties(Worldline * line, const double coord[8]) const;
206 double Interpol3rdOrder(double tt, int indice_time, double values[4]) const;
211 double computeHorizon(const double* pos) const;
212 double computeHorizon(const double* pos, int indice) const;
218 virtual int diff(state_t const &coord, state_t &res) const;
219 int diff(double tt, const double y[7], double res[7]) const ;
220 virtual int diff(const double y[7], double res[7], int indice_time) const ;
241 void circularVelocity(double const coor[4], double vel[3],
243 void circularVelocity(double const coor[4], double vel[3],
244 double dir, int indice_time) const ;
#define GYOTO_OBJECT
Declare class::properties and class::getProperties()
Definition: GyotoObject.h:84
Base class for metrics.
Definition: GyotoMetric.h:158
Reference-counting pointers.
#define GYOTO_OBJECT_THREAD_SAFETY
Declare virtual bool isThreadSafe() const.
Definition: GyotoObject.h:99
Definition: GyotoNumericalMetricLorene.h:26
virtual double gmunu(double const x[4], int mu, int nu) const
Metric coefficients.
Timelike or null geodesics.
Base class for metric description.
Namespace for the Gyoto library.
Definition: GyotoAstrobj.h:43
Pointers performing reference counting.
Definition: GyotoProperty.h:45
Definition: GyotoNumericalMetricLorene.h:47