LORENE
et_bin_global.C
1/*
2 * Methods for computing global quantities within the class Etoile_bin
3 *
4 * (see file etoile.h for documentation)
5 */
6
7/*
8 * Copyright (c) 2000-2001 Eric Gourgoulhon
9 *
10 * This file is part of LORENE.
11 *
12 * LORENE is free software; you can redistribute it and/or modify
13 * it under the terms of the GNU General Public License as published by
14 * the Free Software Foundation; either version 2 of the License, or
15 * (at your option) any later version.
16 *
17 * LORENE is distributed in the hope that it will be useful,
18 * but WITHOUT ANY WARRANTY; without even the implied warranty of
19 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
20 * GNU General Public License for more details.
21 *
22 * You should have received a copy of the GNU General Public License
23 * along with LORENE; if not, write to the Free Software
24 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
25 *
26 */
27
28
29char 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 $" ;
30
31/*
32 * $Id: et_bin_global.C,v 1.12 2014/10/13 08:52:55 j_novak Exp $
33 * $Log: et_bin_global.C,v $
34 * Revision 1.12 2014/10/13 08:52:55 j_novak
35 * Lorene classes and functions now belong to the namespace Lorene.
36 *
37 * Revision 1.11 2004/12/30 17:24:42 f_limousin
38 * Remove #include "et_bin_ncp.h"
39 *
40 * Revision 1.10 2004/12/30 17:22:25 f_limousin
41 * Remove all functions belonging to the class Et_bin_ncp
42 *
43 * Revision 1.9 2003/06/20 14:02:44 f_limousin
44 * pow(gamma(), 1./2.) is changed to pow(a_car(), 3./2.)
45 *
46 * Revision 1.8 2003/02/06 17:25:56 f_limousin
47 * Add global quantities for the class et_bin_ncp
48 *
49 * Revision 1.7 2003/01/17 13:32:27 f_limousin
50 * Add comments
51 *
52 * Revision 1.6 2002/12/17 09:48:21 k_taniguchi
53 * Suppress every desaliasing multiplications using "%".
54 *
55 * Revision 1.5 2002/12/16 17:32:47 k_taniguchi
56 * Suppress the things I did in the previous version.
57 *
58 * Revision 1.4 2002/12/16 16:59:39 k_taniguchi
59 * Set some Cmp to the state of "std_base_scal()".
60 *
61 * Revision 1.3 2002/12/16 14:36:39 k_taniguchi
62 * Introduce a new Cmp for the calculation of gravitational mass.
63 *
64 * Revision 1.2 2002/12/10 15:45:25 k_taniguchi
65 * Change the multiplication "*" to "%".
66 *
67 * Revision 1.1.1.1 2001/11/20 15:19:28 e_gourgoulhon
68 * LORENE
69 *
70 * Revision 2.4 2000/07/06 10:02:22 eric
71 * *** empty log message ***
72 *
73 * Revision 2.3 2000/07/06 09:40:37 eric
74 * Ajout de la fonction xa_barycenter().
75 *
76 * Revision 2.2 2000/02/02 09:23:12 eric
77 * 1ere version operationnelle dans le cas relativiste.
78 *
79 * Revision 2.1 2000/02/01 16:00:13 eric
80 * Le calcul de mass_b est implemente dans le cas relativiste.
81 *
82 * Revision 2.0 2000/01/31 15:57:21 eric
83 * *** empty log message ***
84 *
85 *
86 * $Header: /cvsroot/Lorene/C++/Source/Etoile/et_bin_global.C,v 1.12 2014/10/13 08:52:55 j_novak Exp $
87 *
88 */
89
90// Headers C
91
92// Headers Lorene
93#include "etoile.h"
94
95 //--------------------------//
96 // Baryon mass //
97 //--------------------------//
98
99namespace Lorene {
100double Etoile_bin::mass_b() const {
101
102 if (p_mass_b == 0x0) { // a new computation is required
103
104 if (relativistic) {
105
106 Cmp sqrt_acar = sqrt(a_car()) ;
107 sqrt_acar.std_base_scal() ;
108
109 // See Eq (69) from Gourgoulhon et al. (2001)
110 // Cmp dens = a_car() % sqrt_acar % gam_euler() % nbar() ;
111
112 Cmp dens = a_car() * sqrt( a_car() ) * gam_euler() * nbar() ;
113
114 dens.std_base_scal() ;
115
116 p_mass_b = new double( dens.integrale() ) ;
117
118 }
119 else{
120 assert(nbar.get_etat() == ETATQCQ) ;
121
122 p_mass_b = new double( nbar().integrale() ) ;
123
124 }
125
126 }
127
128 return *p_mass_b ;
129
130}
131
132
133 //----------------------------//
134 // Gravitational mass //
135 //----------------------------//
136
137double Etoile_bin::mass_g() const {
138
139 if (p_mass_g == 0x0) { // a new computation is required
140
141 if (relativistic) {
142
143
144 Cmp sqrt_acar = sqrt(a_car()) ;
145 sqrt_acar.std_base_scal() ;
146 Cmp dens = a_car() * sqrt( a_car() ) * nnn()
147 * ( ener_euler() + s_euler() ) ;
148
149 dens.std_base_scal() ;
150
151 p_mass_g = new double( dens.integrale() ) ;
152
153 }
154 else{
155 p_mass_g = new double( mass_b() ) ; // in the Newtonian case
156 // M_g = M_b
157 }
158 }
159
160 return *p_mass_g ;
161
162}
163
164
165
166 //----------------------------------//
167 // X coordinate of the barycenter //
168 //----------------------------------//
169
171
172 if (p_xa_barycenter == 0x0) { // a new computation is required
173
174 Cmp xxa(mp) ;
175 xxa = mp.xa ; // Absolute X coordinate
176 xxa.std_base_scal() ;
177
178 Cmp dens = a_car() * sqrt( a_car() ) * gam_euler() * nbar() * xxa ;
179
180 dens.std_base_scal() ;
181
182 p_xa_barycenter = new double( dens.integrale() / mass_b() ) ;
183
184 }
185
186 return *p_xa_barycenter ;
187
188}
189
190
191
192}
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
Definition cmp.h:446
void std_base_scal()
Sets the spectral bases of the Valeur va to the standard ones for a scalar.
Definition cmp.C:644
double integrale() const
Computes the integral over all space of *this .
Definition cmp_integ.C:55
virtual double mass_b() const
Baryon mass.
double * p_xa_barycenter
Absolute coordinate X of the barycenter of the baryon density.
Definition etoile.h:1006
virtual double xa_barycenter() const
Absolute coordinate X of the barycenter of the baryon density, defined according to the formula.
virtual double mass_g() const
Gravitational mass.
double * p_mass_b
Baryon mass.
Definition etoile.h:547
double * p_mass_g
Gravitational mass.
Definition etoile.h:548
Tenseur nnn
Total lapse function.
Definition etoile.h:509
Tenseur nbar
Baryon density in the fluid frame.
Definition etoile.h:459
Tenseur gam_euler
Lorentz factor between the fluid and Eulerian observers.
Definition etoile.h:471
Map & mp
Mapping associated with the star.
Definition etoile.h:429
bool relativistic
Indicator of relativity: true for a relativistic star, false for a Newtonian one.
Definition etoile.h:437
Tenseur ener_euler
Total energy density in the Eulerian frame.
Definition etoile.h:465
Tenseur s_euler
Trace of the stress tensor in the Eulerian frame.
Definition etoile.h:468
Tenseur a_car
Total conformal factor .
Definition etoile.h:515
Coord xa
Absolute x coordinate.
Definition map.h:730
int get_etat() const
Returns the logical state.
Definition tenseur.h:707
Cmp sqrt(const Cmp &)
Square root.
Definition cmp_math.C:220
Lorene prototypes.
Definition app_hor.h:64