29#ifndef __ET_ROT_BIFLUID_H_
30#define __ET_ROT_BIFLUID_H_
118#include "eos_bifluid.h"
124 Cmp prolonge_c1(
const Cmp& uu,
const int nzet) ;
315 double precis = 1.e-14) ;
328 double precis = 1.e-14) ;
365 virtual void sauve(FILE *)
const ;
368 virtual ostream&
operator>>(ostream& )
const ;
422 virtual double mass_b()
const ;
423 virtual double mass_g()
const ;
430 virtual double grv2()
const ;
443 virtual double grv3(ostream* ost = 0x0)
const ;
445 virtual double r_circ2()
const ;
446 virtual double area2()
const ;
448 virtual double aplat2()
const ;
574 double omega20,
const Tbl& ent_limit,
575 const Tbl& ent2_limit,
const Itbl& icontrol,
576 const Tbl& control,
Tbl& diff,
577 int mer_mass,
double mbar1_wanted,
double mbar2_wanted,
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
2-fluids equation of state base class.
Class for two-fluid rotating relativistic stars.
const Tenseur & get_K_pp() const
Returns the coefficient Kpp.
virtual void del_hydro_euler()
Sets to ETATNONDEF (undefined state) the hydrodynamical quantities relative to the Eulerian observer.
double * p_angu_mom_1_part2_2
To compute Xn (2nd version)
void equilibrium_bi(double ent_c, double ent_c2, double omega0, double omega20, const Tbl &ent_limit, const Tbl &ent2_limit, const Itbl &icontrol, const Tbl &control, Tbl &diff, int mer_mass, double mbar1_wanted, double mbar2_wanted, double aexp_mass)
Computes an equilibrium configuration.
Tenseur gam_euler2
Lorentz factor between the fluid 2 and Eulerian observers
Tenseur enerps_euler
the combination : useful because in the Newtonian limit .
Tenseur K_np
Coefficient Knp.
virtual ~Et_rot_bifluid()
Destructor.
double * p_mass_b1
Baryon mass of fluid 1.
virtual double angu_mom_1_part1_1() const
To compute In (1st version)
double * p_mass_b2
Baryon mass of fluid 2.
virtual double grv3(ostream *ost=0x0) const
Error on the virial identity GRV3.
virtual const Itbl & l_surf() const
Description of the surface of fluid 1: returns a 2-D Itbl containing the values of the domain index...
Tenseur j_euler12_1
To compute Ip (1st version)
const Tenseur & get_ent2() const
Returns the enthalpy field for fluid 2.
virtual double mom_quad_Bo() const
Part of the quadrupole moment.
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
void set_enthalpies(const Cmp &, const Cmp &)
Sets both enthalpy profiles.
double * p_angu_mom_1_part1_2
To compute In (2nd version)
Tenseur j_euler2
To compute Jp.
Tenseur sphph_euler
The component of the stress tensor .
const Tenseur & get_nbar2() const
Returns the proper baryon density for fluid 2.
double * p_angu_mom_2_part2_1
To compute Xp (1st version)
virtual double area2() const
Surface area for fluid 2.
Tenseur K_pp
Coefficient Kpp.
void operator=(const Et_rot_bifluid &)
Assignment to another Et_rot_bifluid.
virtual void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
const Eos_bifluid & get_eos() const
Returns the equation of state.
virtual double angu_mom_2_part1_2() const
To compute Ip (2nd version)
Tbl * p_xi_surf2
Description of the surface of fluid 2: 2-D Tbl containing the values of the radial coordinate on the...
virtual double angu_mom_2_part1_1() const
To compute Ip (1st version)
double * p_angu_mom_1
Angular momentum of fluid 1.
virtual double angu_mom_2() const
Angular momentum of fluid 2.
Tenseur j_euler11_1
To compute In (1st version)
virtual void del_deriv() const
Deletes all the derived quantities.
virtual double angu_mom_2_part2_2() const
To compute Xp (2nd version)
double * p_ray_eq2
Coordinate radius at , .
virtual void sauve(FILE *) const
Save in a file.
const Tenseur & get_K_np() const
Returns the coefficient Knp.
double * p_angu_mom_1_part2_1
To compute Xn (1st version)
double ray_pole2() const
Coordinate radius for fluid 2 at [r_unit].
void equilibrium_spher_bi(double ent_c, double ent_c2, double precis=1.e-14)
Computes a spherical static configuration.
const Tenseur & get_K_nn() const
Returns the coefficient Knn.
virtual double mass_b() const
Total Baryon mass.
Tenseur j_euler1
To compute Jn.
Tenseur j_euler
Total angular momentum (flat-space!) 3-vector , which is related to of the "3+1" decomposition,...
virtual void equation_of_state()
Computes the proper baryon and energy densities, as well as pressure and the coefficients Knn,...
const Tenseur & get_gam_euler2() const
Returns the Lorentz factor between the fluid 2 and Eulerian observers.
double get_omega2() const
Returns the rotation angular velocity of fluid 2([f_unit] )
const Tbl & xi_surf2() const
Description of the surface of fluid 2: returns a 2-D Tbl containing the values of the radial coordi...
double * p_ray_pole2
Coordinate radius at .
Tenseur K_nn
Coefficient Knn.
virtual double mass_g() const
Gravitational mass.
virtual double angu_mom_2_part2_1() const
To compute Xp (1st version)
double * p_aplat2
Flatening r_pole/r_eq of fluid no.2.
double * p_angu_mom_2_part2_2
To compute Xp (2nd version)
void equil_spher_regular(double ent_c, double ent_c2, double precis=1.e-14)
Computes a spherical static configuration.
double * p_area2
Surface area of fluid no.2.
double * p_angu_mom_2
Angular momentum of fluid 2.
Tenseur delta_car
The "relative velocity" (squared) of the two fluids.
const Itbl & l_surf2() const
Description of the surface of fluid 2: returns a 2-D Itbl containing the values of the domain index...
double omega2
Rotation angular velocity for fluid 2 ([f_unit] )
const Eos_bifluid & eos
Equation of state for two-fluids model.
double * p_angu_mom_2_part1_2
To compute Ip (2nd version)
double mass_b1() const
Baryon mass of fluid 1.
virtual double angu_mom_1_part1_2() const
To compute In (2nd version)
virtual double angu_mom_1_part2_1() const
To compute Xn (1st version)
Tenseur uuu2
Norm of the (fluid no.2) 3-velocity with respect to the eulerian observer.
double ray_eq2() const
Coordinate radius for fluid 2 at , [r_unit].
virtual double r_circ2() const
Circumferential radius for fluid 2.
double mass_b2() const
Baryon mass of fluid 2.
virtual double mom_quad() const
Quadrupole moment.
virtual double aplat2() const
Flatening r_pole/r_eq for fluid 2.
virtual double grv2() const
Error on the virial identity GRV2.
double * p_angu_mom_1_part1_1
To compute In (1st version)
double * p_r_circ2
Circumferential radius of fluid no.2.
Tenseur ent2
Log-enthalpy for the second fluid.
double * p_ray_eq2_pi
Coordinate radius at , .
virtual double angu_mom_1_part2_2() const
To compute Xn (2nd version)
const Tenseur & get_uuu2() const
Returns the norm of the fluid 2 3-velocity with respect to the eulerian frame.
virtual void partial_display(ostream &) const
Printing of some informations, excluding all global quantities.
double * p_angu_mom_2_part1_1
To compute Ip (1st version)
Itbl * p_l_surf2
Description of the surface of fluid 2: 2-D Itbl containing the values of the domain index l on the su...
Tenseur nbar2
Baryon density in the fluid frame, for fluid 2.
virtual double angu_mom_1() const
Angular momentum of fluid 1.
virtual double angu_mom() const
Angular momentum.
double * p_ray_eq2_pis2
Coordinate radius at , .
double ray_eq2_pis2() const
Coordinate radius for fluid 2 at , [r_unit].
double ray_eq2_pi() const
Coordinate radius for fluid 2 at , [r_unit].
const Tenseur & get_delta_car() const
Returns the "relative velocity" (squared) of the two fluids.
virtual void hydro_euler()
Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid fram...
virtual double mean_radius2() const
Mean radius for fluid 2.
virtual double mom_quad_old() const
Part of the quadrupole moment.
Class for isolated rotating stars *** DEPRECATED : use class Star_rot instead ***.
Basic integer array class.
Base class for coordinate mappings.
Tensor handling *** DEPRECATED : use class Tensor instead ***.