28 char blackhole_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Black_hole/blackhole.C,v 1.5 2014/10/13 08:52:45 j_novak Exp $" ;
60 #include "blackhole.h"
61 #include "utilitaires.h"
73 kerrschild(kerrschild_i),
79 shift(mp_i, CON, mp_i.get_bvect_cart()),
80 shift_rs(mp_i, CON, mp_i.get_bvect_cart()),
81 shift_bh(mp_i, CON, mp_i.get_bvect_cart()),
83 taij(mp_i, CON, mp_i.get_bvect_cart()),
84 taij_rs(mp_i, CON, mp_i.get_bvect_cart()),
87 flat(mp_i, mp_i.get_bvect_cart()) {
120 kerrschild(bh.kerrschild),
123 lapconf_rs(bh.lapconf_rs),
124 lapconf_bh(bh.lapconf_bh),
127 shift_rs(bh.shift_rs),
128 shift_bh(bh.shift_bh),
132 taij_quad(bh.taij_quad),
133 taij_quad_rs(bh.taij_quad_rs),
145 lapconf_rs(mp_i, *(mp_i.get_mg()), fich),
148 shift(mp_i, CON, mp_i.get_bvect_cart()),
149 shift_rs(mp_i, mp_i.get_bvect_cart(), fich),
150 shift_bh(mp_i, CON, mp_i.get_bvect_cart()),
151 confo(mp_i, *(mp_i.get_mg()), fich),
152 taij(mp_i, CON, mp_i.get_bvect_cart()),
153 taij_rs(mp_i, CON, mp_i.get_bvect_cart()),
156 flat(mp_i, mp_i.get_bvect_cart()) {
161 fread(&
mass_bh,
sizeof(
double), 1, fich) ;
238 assert( &(bh.
mp) == &
mp ) ;
274 fwrite(&
mass_bh,
sizeof(
double), 1, fich) ;
296 ost <<
"Kerr-Schild background" << endl ;
297 ost <<
"----------------------" << endl ;
300 ost <<
"Conformally flat background" << endl ;
301 ost <<
"---------------------------" << endl ;
304 ost <<
"lapconf on the AH : "
306 ost <<
"lapse on the AH : "
308 ost <<
"shift(1) on the AH : "
309 <<
shift(1).val_grid_point(1,0,nt-1,0) << endl ;
310 ost <<
"shift(2) on the AH : "
311 <<
shift(2).val_grid_point(1,0,nt-1,0) << endl ;
312 ost <<
"shift(3) on the AH : "
313 <<
shift(3).val_grid_point(1,0,nt-1,0) << endl ;
314 ost <<
"confo on the AH : "
316 ost <<
"Gravitational mass : "
317 <<
mass_bh / msol <<
" M_sol" << endl ;
318 ost <<
"Irreducible mass : "
319 <<
mass_irr() / msol <<
" M_sol" << endl ;
321 <<
mass_adm() / msol <<
" M_sol" << endl ;
322 ost <<
"Komar mass : "
323 <<
mass_kom() / msol <<
" M_sol" << endl ;
325 double irr_gm, adm_gm, kom_gm ;
329 ost <<
"Diff. (Mirr-Mg)/Mg : " << irr_gm << endl ;
330 ost <<
"Diff. (Madm-Mg)/Mg : " << adm_gm << endl ;
331 ost <<
"Diff. (Mkom-Mg)/Mg : " << kom_gm << endl ;