26#ifndef __BIN_NS_NCP_H_
27#define __BIN_NS_NCP_H_
153 Map& mp2,
int nzet2,
const Eos& eos2,
int irrot2,
154 int relat,
int conf_flat,
const Metrique& flat1,
const Metrique& flat2,
170 const Metrique& flat1,
const Metrique& flat2, FILE* fich) ;
195 { assert( (i==1) || (i==2) );
211 { assert( (i==1) || (i==2) );
229 void sauve(FILE *)
const ;
341 void orbit(
double fact_omeg_min,
double fact_omeg_max,
double& xgg1,
Extended description of the class for Doc++ documentation.
double * p_mass_kom
Total Komar mass of the system.
double * p_virial_gb
Virial theorem error by E.Gourgoulhon and S.Bonazzola.
double virial_gb() const
Estimates the relative error on the virial theorem calculated by E.Gourgoulhon and S....
Et_bin_ncp * et[2]
Array of the two stars (to perform loops on the stars): {\tt et[0]} contains the address of {\tt star...
double virial() const
Estimates the relative error on the virial theorem (for a relativistic one, it returns $|1 - M_{\rm K...
const Tbl & mom_constr() const
Estimates the relative error on the momentum constraint equation by comparing ${\overline\nabla}_j K^...
double separation() const
Returns the coordinate separation of the two stellar centers [{\tt r_unit}].
void operator=(const Bin_ns_ncp &)
Assignment to another {\tt Bin_ns_ncp}.
Et_bin_ncp & set(int i)
Read/write of the star no. i.
double ham_constr() const
Estimates the relative error on the Hamiltonian constraint equation by comparing $\underline\Delta\ln...
const Tbl & angu_mom() const
Total angular momentum.
double get_x_axe() const
Returns the absolute coordinate X of the rotation axis [{\tt r_unit}].
double virial_fus() const
Estimates the relative error on the virial theorem calculated by J.L.Friedman, K.Uryu,...
double x_axe
Absolute X coordinate of the rotation axis.
double omega
Angular velocity with respect to an asymptotically inertial observer.
ostream & operator>>(ostream &) const
Operator >> (function called by the operator <<).
Bin_ns_ncp(Map &mp1, int nzet1, const Eos &eos1, int irrot1, Map &mp2, int nzet2, const Eos &eos2, int irrot2, int relat, int conf_flat, const Metrique &flat1, const Metrique &flat2, const Tenseur_sym &source1, const Tenseur_sym &source2)
Standard constructor.
void write_global(ostream &) const
Write global quantities in a formatted file.
void analytical_shift()
Sets some analytical template for the shift vector (via the members {\tt w_shift} and {\tt khi_shift}...
Tbl * p_angu_mom
Total angular momentum of the system.
void del_deriv() const
Destructor.
double * p_mass_adm
Total ADM mass of the system.
void fait_decouple()
Calculates {tt decouple} which is used to obtain {\tt a_car_auto} by the formula : {\tt a_car_auto} =...
Bin_ns_ncp(Map &mp1, const Eos &eos1, Map &mp2, const Eos &eos2, const Metrique &flat1, const Metrique &flat2, FILE *fich)
Copy constructor.
double mass_kom() const
Total Komar mass.
double * p_virial_fus
Virial theorem error by J.L.Friedman, K.Uryu, and M.Shibata.
void analytical_omega()
Sets the orbital angular velocity to some 2-PN analytical value (Keplerian value in the Newtonian cas...
friend ostream & operator<<(ostream &, const Bin_ns_ncp &)
Save in a file.
double * p_total_ener
Total energy of the system.
Tbl * p_mom_constr
Relative error on the momentum constraint.
const Et_bin_ncp & operator()(int i) const
Returns a reference to the star no. i.
double * p_ham_constr
Relative error on the Hamiltonian constraint.
void set_der_0x0() const
Sets to {\tt 0x0} all the pointers on derived quantities.
double total_ener() const
Total energy (excluding the rest mass energy).
Et_bin_ncp star2
Second star ncp of the system.
double get_omega() const
Returns the orbital angular velocity [{\tt f_unit}].
double mass_adm() const
Total ADM mass.
double * p_virial
Virial theorem error.
void orbit(double fact_omeg_min, double fact_omeg_max, double &xgg1, double &xgg2)
Computes the orbital angular velocity {\tt omega} and the position of the rotation axis {\tt x_axe}.
void display_poly(ostream &) const
Display in polytropic units.
Et_bin_ncp star1
First star ncp of the system.
double & set_omega()
Sets the orbital angular velocity [{\tt f_unit}].
double & set_x_axe()
Sets the absolute coordinate X of the rotation axis [{\tt r_unit}].
Equation of state base class.
Base class for coordinate mappings.
Class intended to describe tensors with a symmetry on the two last indices *** DEPRECATED : use class...