28 char compobj_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Compobj/compobj.C,v 1.9 2014/10/13 08:52:49 j_novak Exp $" ;
73 #include "utilitaires.h"
85 beta(map_i, CON, map_i.get_bvect_spher()) ,
86 gamma(map_i.flat_met_spher()) ,
88 mom_euler(map_i, CON, map_i.get_bvect_spher()) ,
89 stress_euler(map_i, COV, map_i.get_bvect_spher()) ,
90 kk(map_i, COV, map_i.get_bvect_spher())
114 ener_euler(co.ener_euler) ,
115 mom_euler(co.mom_euler) ,
116 stress_euler(co.stress_euler) ,
128 nn(map_i, *(map_i.get_mg()), fich) ,
129 beta(map_i, map_i.get_bvect_spher(), fich) ,
131 ener_euler(map_i, *(map_i.get_mg()), fich) ,
132 mom_euler(map_i, map_i.get_bvect_spher(), fich) ,
133 stress_euler(map_i, map_i.get_bvect_spher(), fich) ,
134 kk(map_i, COV, map_i.get_bvect_spher())
177 assert( &(co.
mp) == &
mp ) ;
212 FILE* file_out = fopen(file_name,
"w") ;
213 double total_time = 0. ;
215 fwrite_be(&total_time,
sizeof(
double), 1, file_out) ;
227 cout <<
"WRITING TO GYOTO FILE - OK: " << endl ;
233 ostream& operator<<(ostream& ost,
const Compobj& co) {
241 ost << endl <<
"Compact object (class Compobj) " << endl ;
242 ost <<
"Mapping : " <<
mp << endl ;
243 ost <<
"Central values of various fields : " << endl ;
244 ost <<
"-------------------------------- " << endl ;
246 ost <<
" metric components gamma_{ij} : " << endl
247 <<
" ( " <<
gamma.
cov()(1,1).val_grid_point(0,0,0,0) <<
" "
248 <<
gamma.
cov()(1,2).val_grid_point(0,0,0,0) <<
" "
249 <<
gamma.
cov()(1,3).val_grid_point(0,0,0,0) <<
" )" << endl
250 <<
" ( " <<
gamma.
cov()(2,1).val_grid_point(0,0,0,0) <<
" "
251 <<
gamma.
cov()(2,2).val_grid_point(0,0,0,0) <<
" "
252 <<
gamma.
cov()(2,3).val_grid_point(0,0,0,0) <<
" )" << endl
253 <<
" ( " <<
gamma.
cov()(3,1).val_grid_point(0,0,0,0) <<
" "
254 <<
gamma.
cov()(3,2).val_grid_point(0,0,0,0) <<
" "
255 <<
gamma.
cov()(3,3).val_grid_point(0,0,0,0) <<
" )" << endl ;
256 ost <<
" components of the extrinsic curvature K_{ij} : " << endl
257 <<
" ( " <<
kk(1,1).val_grid_point(0,0,0,0) <<
" "
258 <<
kk(1,2).val_grid_point(0,0,0,0) <<
" "
259 <<
kk(1,3).val_grid_point(0,0,0,0) <<
" )" << endl
260 <<
" ( " <<
kk(2,1).val_grid_point(0,0,0,0) <<
" "
261 <<
kk(2,2).val_grid_point(0,0,0,0) <<
" "
262 <<
kk(2,3).val_grid_point(0,0,0,0) <<
" )" << endl
263 <<
" ( " <<
kk(3,1).val_grid_point(0,0,0,0) <<
" "
264 <<
kk(3,2).val_grid_point(0,0,0,0) <<
" "
265 <<
kk(3,3).val_grid_point(0,0,0,0) <<
" )" << endl ;
267 ost <<
" components of the stress tensor S_{ij} / Eulerian observer : " << endl
268 <<
" ( " <<
stress_euler(1,1).val_grid_point(0,0,0,0) <<
" "
270 <<
stress_euler(1,3).val_grid_point(0,0,0,0) <<
" )" << endl
271 <<
" ( " <<
stress_euler(2,1).val_grid_point(0,0,0,0) <<
" "
273 <<
stress_euler(2,3).val_grid_point(0,0,0,0) <<
" )" << endl
274 <<
" ( " <<
stress_euler(3,1).val_grid_point(0,0,0,0) <<
" "
276 <<
stress_euler(3,3).val_grid_point(0,0,0,0) <<
" )" << endl ;
292 cout <<
"WARNING: Compobj::extrinsic_curvature() NOT TESTED !" << endl ;
300 for (
int i=1; i<=3; i++) {
301 for (
int j=i; j<=3; j++) {
302 kk.
set(i, j) = (dn(i, j) + dn(j, i))/(2*
nn) ;
320 p_adm_mass =
new double( ww.
flux(__infinity, ff) / (16.* M_PI) ) ;