star_bhns_global.C
/*
* Methods of class Star_bhns to compute global quantities
*
* (see file star_bhns.h for documentation).
*
*/
/*
* Copyright (c) 2005,2007 Keisuke Taniguchi
*
* 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 version 2
* as published by the Free Software Foundation.
*
* 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 star_bhns_global_C[] = "$Header: /cvsroot/Lorene/C++/Source/Star_bhns/star_bhns_global.C,v 1.3 2014/10/13 08:53:40 j_novak Exp $" ;
/*
* $Id: star_bhns_global.C,v 1.3 2014/10/13 08:53:40 j_novak Exp $
* $Log: star_bhns_global.C,v $
* Revision 1.3 2014/10/13 08:53:40 j_novak
* Lorene classes and functions now belong to the namespace Lorene.
*
* Revision 1.2 2008/05/15 19:15:20 k_taniguchi
* Change of a parameter.
*
* Revision 1.1 2007/06/22 01:31:24 k_taniguchi
* *** empty log message ***
*
*
* $Header: /cvsroot/Lorene/C++/Source/Star_bhns/star_bhns_global.C,v 1.3 2014/10/13 08:53:40 j_novak Exp $
*
*/
// C++ headers
//#include <>
// C headers
//#include <>
// Lorene headers
#include "star_bhns.h"
#include "unites.h"
#include "utilitaires.h"
//-------------------------------//
// Baryon mass //
//-------------------------------//
namespace Lorene {
double Star_bhns::mass_b() const {
using namespace Unites ;
if (p_mass_b == 0x0) { // a new computation is required
Scalar dens = pow(confo_tot, 6.) * nbar * gam_euler ;
dens.std_spectral_base() ;
p_mass_b = new double( dens.integrale() ) ;
}
return *p_mass_b ;
/*
cout << "Star_bhns::mass_b() is not available !!!" << endl
<< " --> Use Star_bhns::mass_b_bhns(kerrschild, mass_bh, sepa)"
<< endl ;
abort() ;
*/
}
double Star_bhns::mass_b_bhns(bool kerrschild, const double& mass_bh,
const double& sepa) const {
using namespace Unites ;
if (p_mass_b_bhns == 0x0) { // a new computation is required
Scalar tmp(mp) ;
if (kerrschild) {
Scalar xx(mp) ;
xx = mp.x ;
xx.std_spectral_base() ;
Scalar yy(mp) ;
yy = mp.y ;
yy.std_spectral_base() ;
Scalar zz(mp) ;
zz = mp.z ;
zz.std_spectral_base() ;
double yns = mp.get_ori_y() ;
Scalar rr(mp) ;
rr = sqrt( (xx+sepa)*(xx+sepa) + (yy+yns)*(yy+yns) + zz*zz ) ;
rr.std_spectral_base() ;
tmp = sqrt(1. + 2.*ggrav*mass_bh/rr) ;
}
else { // Isotropic coordinates with the maximal slicing
tmp = 1. ;
}
tmp.std_spectral_base() ;
Scalar dens = pow(confo_tot, 6.) * nbar * gam_euler * tmp ;
dens.std_spectral_base() ;
p_mass_b_bhns = new double( dens.integrale() ) ;
}
return *p_mass_b_bhns ;
}
//--------------------------------------//
// Gravitational mass //
//--------------------------------------//
double Star_bhns::mass_g() const {
cout << "Star_bhns::mass_g() is not available !!!" << endl
<< " --> Use Star_bhns::mass_g_bhns()"
<< endl ;
abort() ;
}
double Star_bhns::mass_g_bhns() const {
// This mass is valid only for an isolated spherical star
if (p_mass_g_bhns == 0x0) { // a new computation is required
Scalar dens = lapconf_tot * pow(confo_tot, 5.)
* (ener_euler + s_euler) ;
dens.std_spectral_base() ;
p_mass_g_bhns = new double( dens.integrale() ) ;
}
return *p_mass_g_bhns ;
}
}
Lorene::Scalar
Tensor field of valence 0 (or component of a tensorial field).
Definition scalar.h:387
Lorene::Scalar::integrale
double integrale() const
Computes the integral over all space of *this .
Definition scalar_integ.C:61
Lorene::Scalar::std_spectral_base
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
Definition scalar.C:784
Lorene::Star_bhns::confo_tot
Scalar confo_tot
Total conformal factor.
Definition star_bhns.h:163
Lorene::Star_bhns::lapconf_tot
Scalar lapconf_tot
Total lapconf function.
Definition star_bhns.h:119
Lorene::Star_bhns::p_mass_b_bhns
double * p_mass_b_bhns
Baryon mass.
Definition star_bhns.h:225
Lorene::Star_bhns::mass_g
virtual double mass_g() const
Gravitational mass.
Definition star_bhns_global.C:138
Lorene::Star_bhns::mass_b
virtual double mass_b() const
Baryon mass.
Definition star_bhns_global.C:63
Lorene::Star_bhns::p_mass_g_bhns
double * p_mass_g_bhns
Gravitational mass.
Definition star_bhns.h:226
Lorene::Star::p_mass_b
double * p_mass_b
Baryon mass.
Definition star.h:268
Lorene::Star::nbar
Scalar nbar
Baryon density in the fluid frame.
Definition star.h:192
Lorene::Star::ener_euler
Scalar ener_euler
Total energy density in the Eulerian frame.
Definition star.h:198
Lorene::Star::gam_euler
Scalar gam_euler
Lorentz factor between the fluid and Eulerian observers.
Definition star.h:204
Lorene::Star::s_euler
Scalar s_euler
Trace of the stress scalar in the Eulerian frame.
Definition star.h:201
Lorene::Star::mp
Map & mp
Mapping associated with the star.
Definition star.h:180
Lorene::sqrt
Cmp sqrt(const Cmp &)
Square root.
Definition cmp_math.C:220
Lorene::pow
Cmp pow(const Cmp &, int)
Power .
Definition cmp_math.C:348
Lorene
Lorene prototypes.
Definition app_hor.h:64
Unites
Standard units of space, time and mass.