60 #include "star_bhns.h" 65 const double& sepa,
double omega,
66 double x_rot,
double y_rot) {
88 bsn.
set(1) = - omega * (ya - y_rot) ;
89 bsn.
set(2) = omega * (xa - x_rot) ;
95 bsn.
set(1) = omega * (ya - y_rot) ;
96 bsn.
set(2) = - omega * (xa - x_rot) ;
106 for (
int i=1; i<=3; i++) {
124 double mass = ggrav * mass_bh ;
139 rbh =
sqrt( (xx+sepa)*(xx+sepa) + (yy+yns)*(yy+yns) + zz*zz ) ;
144 ll.set(1) = (xx+sepa) / rbh ;
145 ll.
set(2) = (yy+yns) / rbh ;
146 ll.
set(3) = zz / rbh ;
154 llbsn = ll(1)%
bsn(1) + ll(2)%
bsn(2) + ll(3)%
bsn(3) ;
158 tmp1 = 2. * msr % llbsn % llbsn ;
Scalar psi4
Fourth power of the total conformal factor.
Coord xa
Absolute x coordinate.
void annule_domain(int l)
Sets the Tensor to zero in a given domain.
Vector shift_tot
Total shift vector.
Scalar lapconf_tot
Total lapconf function.
Cmp log(const Cmp &)
Neperian logarithm.
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
Map & mp
Mapping associated with the star.
virtual void del_deriv() const
Deletes all the derived quantities.
Cmp sqrt(const Cmp &)
Square root.
Vector bsn
3-vector shift, divided by N , of the rotating coordinates, .
double get_ori_y() const
Returns the y coordinate of the origin.
Standard units of space, time and mass.
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).
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field. ...
double get_rot_phi() const
Returns the angle between the x –axis and X –axis.
Tensor field of valence 1.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
Scalar confo_tot
Total conformal factor.
void kinema_bhns(bool kerrschild, const double &mass_bh, const double &sepa, double omega, double x_rot, double y_rot)
Computes the quantities bsn and pot_centri .
int get_nzone() const
Returns the number of domains.
Scalar gam0
Lorentz factor between the co-orbiting observer and the Eulerian one.
Active physical coordinates and mapping derivatives.
Scalar pot_centri
Centrifugal potential.
Coord ya
Absolute y coordinate.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping, i.e.
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
Coord y
y coordinate centered on the grid
Coord x
x coordinate centered on the grid
Scalar & set(int)
Read/write access to a component.
Coord z
z coordinate centered on the grid