31char star_rot_isco_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star/star_rot_isco.C,v 1.5 2014/10/13 08:53:39 j_novak Exp $" ;
62#include "utilitaires.h"
65double funct_star_rot_isco(
double,
const Param& ) ;
112 2 / (1 + r *
bdlog ) ;
128 orbit.std_spectral_base() ;
173 par_ms.add_scalar(orbit) ;
186 " number of iterations used in zerosec to locate the ISCO : "
216 (
double(2) *
M_PI) ) ;
239 (
double(2) *
M_PI) ) ;
326double funct_star_rot_isco(
double xi,
const Param&
par){
Time evolution with partial storage (*** under development ***).
Coord r
r coordinate centered on the grid
virtual double val_r(int l, double xi, double theta, double pphi) const =0
Returns the value of the radial coordinate r for a given in a given domain.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tensor field of valence 0 (or component of a tensorial field).
const Scalar & dsdr() const
Returns of *this .
const Valeur & get_spectral_va() const
Returns va (read only version)
virtual double lspec_isco() const
Angular momentum of a particle on the ISCO.
double * p_r_isco
Circumferential radius of the ISCO.
virtual double espec_isco() const
Energy of a particle on the ISCO.
virtual double r_isco(ostream *ost=0x0) const
Circumferential radius of the innermost stable circular orbit (ISCO).
double * p_lspec_isco
Specific angular momentum of a particle on the ISCO.
Scalar bbb
Metric factor B.
Scalar nphi
Metric coefficient .
double * p_f_isco
Orbital frequency of the ISCO.
virtual double f_isco() const
Orbital frequency at the innermost stable circular orbit (ISCO).
double * p_espec_isco
Specific energy of a particle on the ISCO.
double * p_f_eq
Orbital frequency at the equator.
virtual double f_eq() const
Orbital frequency at the equator.
Scalar nn
Lapse function N .
double ray_eq() const
Coordinate radius at , [r_unit].
Map & mp
Mapping associated with the star.
int nzet
Number of domains of *mp occupied by the star.
Values and coefficients of a (real-value) function.
Cmp sqrt(const Cmp &)
Square root.
double zerosec(double(*f)(double, const Param &), const Param &par, double a, double b, double precis, int nitermax, int &niter, bool abort=true)
Finding the zero a function.
void annule_domain(int l)
Sets the Tensor to zero in a given domain.