et_bin_global.C
/*
* Methods for computing global quantities within the class Etoile_bin
*
* (see file etoile.h for documentation)
*/
/*
* Copyright (c) 2000-2001 Eric Gourgoulhon
*
* This file is part of LORENE.
*
* LORENE is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* LORENE is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with LORENE; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*
*/
char et_bin_global_C[] = "$Header: /cvsroot/Lorene/C++/Source/Etoile/et_bin_global.C,v 1.12 2014/10/13 08:52:55 j_novak Exp $" ;
/*
* $Id: et_bin_global.C,v 1.12 2014/10/13 08:52:55 j_novak Exp $
* $Log: et_bin_global.C,v $
* Revision 1.12 2014/10/13 08:52:55 j_novak
* Lorene classes and functions now belong to the namespace Lorene.
*
* Revision 1.11 2004/12/30 17:24:42 f_limousin
* Remove #include "et_bin_ncp.h"
*
* Revision 1.10 2004/12/30 17:22:25 f_limousin
* Remove all functions belonging to the class Et_bin_ncp
*
* Revision 1.9 2003/06/20 14:02:44 f_limousin
* pow(gamma(), 1./2.) is changed to pow(a_car(), 3./2.)
*
* Revision 1.8 2003/02/06 17:25:56 f_limousin
* Add global quantities for the class et_bin_ncp
*
* Revision 1.7 2003/01/17 13:32:27 f_limousin
* Add comments
*
* Revision 1.6 2002/12/17 09:48:21 k_taniguchi
* Suppress every desaliasing multiplications using "%".
*
* Revision 1.5 2002/12/16 17:32:47 k_taniguchi
* Suppress the things I did in the previous version.
*
* Revision 1.4 2002/12/16 16:59:39 k_taniguchi
* Set some Cmp to the state of "std_base_scal()".
*
* Revision 1.3 2002/12/16 14:36:39 k_taniguchi
* Introduce a new Cmp for the calculation of gravitational mass.
*
* Revision 1.2 2002/12/10 15:45:25 k_taniguchi
* Change the multiplication "*" to "%".
*
* Revision 1.1.1.1 2001/11/20 15:19:28 e_gourgoulhon
* LORENE
*
* Revision 2.4 2000/07/06 10:02:22 eric
* *** empty log message ***
*
* Revision 2.3 2000/07/06 09:40:37 eric
* Ajout de la fonction xa_barycenter().
*
* Revision 2.2 2000/02/02 09:23:12 eric
* 1ere version operationnelle dans le cas relativiste.
*
* Revision 2.1 2000/02/01 16:00:13 eric
* Le calcul de mass_b est implemente dans le cas relativiste.
*
* Revision 2.0 2000/01/31 15:57:21 eric
* *** empty log message ***
*
*
* $Header: /cvsroot/Lorene/C++/Source/Etoile/et_bin_global.C,v 1.12 2014/10/13 08:52:55 j_novak Exp $
*
*/
// Headers C
// Headers Lorene
#include "etoile.h"
//--------------------------//
// Baryon mass //
//--------------------------//
namespace Lorene {
double Etoile_bin::mass_b() const {
if (p_mass_b == 0x0) { // a new computation is required
if (relativistic) {
Cmp sqrt_acar = sqrt(a_car()) ;
sqrt_acar.std_base_scal() ;
// See Eq (69) from Gourgoulhon et al. (2001)
// Cmp dens = a_car() % sqrt_acar % gam_euler() % nbar() ;
Cmp dens = a_car() * sqrt( a_car() ) * gam_euler() * nbar() ;
dens.std_base_scal() ;
p_mass_b = new double( dens.integrale() ) ;
}
else{
assert(nbar.get_etat() == ETATQCQ) ;
p_mass_b = new double( nbar().integrale() ) ;
}
}
return *p_mass_b ;
}
//----------------------------//
// Gravitational mass //
//----------------------------//
double Etoile_bin::mass_g() const {
if (p_mass_g == 0x0) { // a new computation is required
if (relativistic) {
Cmp sqrt_acar = sqrt(a_car()) ;
sqrt_acar.std_base_scal() ;
Cmp dens = a_car() * sqrt( a_car() ) * nnn()
* ( ener_euler() + s_euler() ) ;
dens.std_base_scal() ;
p_mass_g = new double( dens.integrale() ) ;
}
else{
p_mass_g = new double( mass_b() ) ; // in the Newtonian case
// M_g = M_b
}
}
return *p_mass_g ;
}
//----------------------------------//
// X coordinate of the barycenter //
//----------------------------------//
double Etoile_bin::xa_barycenter() const {
if (p_xa_barycenter == 0x0) { // a new computation is required
Cmp xxa(mp) ;
xxa = mp.xa ; // Absolute X coordinate
xxa.std_base_scal() ;
Cmp dens = a_car() * sqrt( a_car() ) * gam_euler() * nbar() * xxa ;
dens.std_base_scal() ;
p_xa_barycenter = new double( dens.integrale() / mass_b() ) ;
}
return *p_xa_barycenter ;
}
}
Lorene::Cmp
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
Definition cmp.h:446
Lorene::Cmp::std_base_scal
void std_base_scal()
Sets the spectral bases of the Valeur va to the standard ones for a scalar.
Definition cmp.C:644
Lorene::Cmp::integrale
double integrale() const
Computes the integral over all space of *this .
Definition cmp_integ.C:55
Lorene::Etoile_bin::mass_b
virtual double mass_b() const
Baryon mass.
Definition et_bin_global.C:100
Lorene::Etoile_bin::p_xa_barycenter
double * p_xa_barycenter
Absolute coordinate X of the barycenter of the baryon density.
Definition etoile.h:1006
Lorene::Etoile_bin::xa_barycenter
virtual double xa_barycenter() const
Absolute coordinate X of the barycenter of the baryon density, defined according to the formula.
Definition et_bin_global.C:170
Lorene::Etoile_bin::mass_g
virtual double mass_g() const
Gravitational mass.
Definition et_bin_global.C:137
Lorene::Etoile::p_mass_b
double * p_mass_b
Baryon mass.
Definition etoile.h:547
Lorene::Etoile::p_mass_g
double * p_mass_g
Gravitational mass.
Definition etoile.h:548
Lorene::Etoile::nnn
Tenseur nnn
Total lapse function.
Definition etoile.h:509
Lorene::Etoile::nbar
Tenseur nbar
Baryon density in the fluid frame.
Definition etoile.h:459
Lorene::Etoile::gam_euler
Tenseur gam_euler
Lorentz factor between the fluid and Eulerian observers.
Definition etoile.h:471
Lorene::Etoile::mp
Map & mp
Mapping associated with the star.
Definition etoile.h:429
Lorene::Etoile::relativistic
bool relativistic
Indicator of relativity: true for a relativistic star, false for a Newtonian one.
Definition etoile.h:437
Lorene::Etoile::ener_euler
Tenseur ener_euler
Total energy density in the Eulerian frame.
Definition etoile.h:465
Lorene::Etoile::s_euler
Tenseur s_euler
Trace of the stress tensor in the Eulerian frame.
Definition etoile.h:468
Lorene::Etoile::a_car
Tenseur a_car
Total conformal factor .
Definition etoile.h:515
Lorene::sqrt
Cmp sqrt(const Cmp &)
Square root.
Definition cmp_math.C:220
Lorene
Lorene prototypes.
Definition app_hor.h:64
Lorene::integrale
virtual Tbl * integrale(const Cmp &) const =0
Computes the integral over all space of a Cmp .