28 char et_bin_nsbh_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Etoile/et_bin_nsbh.C,v 1.12 2014/10/13 08:52:56 j_novak Exp $" ;
81 #include "et_bin_nsbh.h"
84 #include "utilitaires.h"
97 :
Etoile_bin(mp_i, nzet_i, relat, eos_i, irrot, ref_triad_i),
100 d_n_auto(mp_i, 1, COV, ref_triad_i),
101 d_n_comp(mp_i, 1, COV, ref_triad_i),
105 d_confpsi_auto(mp_i, 1, COV, ref_triad_i),
106 d_confpsi_comp(mp_i, 1, COV, ref_triad_i),
107 taij_auto(mp_i, 2, CON, ref_triad_i) ,
108 taij_comp(mp_i, 2, CON, ref_triad_i),
109 taij_tot(mp_i, 2, CON, ref_triad_i) ,
110 tkij_auto(mp_i, 2, CON, ref_triad_i),
111 tkij_tot(mp_i, 2, CON, ref_triad_i),
113 ssjm1_confpsi(mp_i) {
153 d_n_auto(et.d_n_auto),
154 d_n_comp(et.d_n_comp),
156 confpsi_auto(et.confpsi_auto),
157 confpsi_comp(et.confpsi_comp),
158 d_confpsi_auto(et.d_confpsi_auto),
159 d_confpsi_comp(et.d_confpsi_comp),
160 taij_auto(et.taij_auto),
161 taij_comp(et.taij_comp),
162 taij_tot(et.taij_tot),
163 tkij_auto(et.tkij_auto),
164 tkij_tot(et.tkij_tot),
165 ssjm1_lapse(et.ssjm1_lapse),
166 ssjm1_confpsi(et.ssjm1_confpsi) {
174 const Base_vect& ref_triad_i, FILE* fich,
bool old)
178 d_n_auto(mp_i, 1, COV, ref_triad_i),
179 d_n_comp(mp_i, 1, COV, ref_triad_i ),
183 d_confpsi_auto(mp_i, 1, COV, ref_triad_i),
184 d_confpsi_comp(mp_i, 1, COV, ref_triad_i),
185 taij_auto(mp_i, 2, CON, ref_triad_i),
186 taij_comp(mp_i, 2, CON, ref_triad_i),
187 taij_tot(mp_i, 2, CON, ref_triad_i),
188 tkij_auto(mp_i, 2, CON, ref_triad_i),
189 tkij_tot(mp_i, 2, CON, ref_triad_i),
191 ssjm1_confpsi(mp_i) {
194 Cmp n_from_file (mp_i, *(mp_i.
get_mg()), fich) ;
198 Cmp psi_from_file (mp_i, *(mp_i.
get_mg()), fich) ;
203 Tenseur shift_auto_file (
mp, ref_triad_i, fich) ;
226 Cmp ssjm1_lapse_file(mp_i, *(mp_i.
get_mg()), fich) ;
229 Cmp ssjm1_confpsi_file(mp_i, *(mp_i.
get_mg()), fich) ;
241 Et_bin_nsbh::~Et_bin_nsbh(){
335 ost <<
"Neutron star in a binary system" << endl ;
336 ost <<
"-------------------------------" << endl ;
339 ost <<
"irrotational configuration" << endl ;
342 ost <<
"corotating configuration" << endl ;
345 ost <<
"Absolute abscidia of the stellar center: " <<
348 ost <<
"Absolute abscidia of the barycenter of the baryon density : " <<
353 double d_tilde = 2 * d_ns / r_0 ;
355 ost <<
"d_tilde : " << d_tilde << endl ;
357 ost <<
"Orientation with respect to the absolute frame : " <<
360 ost <<
"Central value of gam_euler : "
363 ost <<
"Central u_euler (U^X, U^Y, U^Z) [c] : "
364 <<
u_euler(0)(0, 0, 0, 0) <<
" "
365 <<
u_euler(1)(0, 0, 0, 0) <<
" "
366 <<
u_euler(2)(0, 0, 0, 0) << endl ;
369 ost <<
"Central d_psi (X, Y, Z) [c] : "
370 <<
d_psi(0)(0, 0, 0, 0) <<
" "
371 <<
d_psi(1)(0, 0, 0, 0) <<
" "
372 <<
d_psi(2)(0, 0, 0, 0) << endl ;
374 ost <<
"Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : "
375 <<
wit_w(0)(0, 0, 0, 0) <<
" "
376 <<
wit_w(1)(0, 0, 0, 0) <<
" "
377 <<
wit_w(2)(0, 0, 0, 0) << endl ;
379 ost <<
"Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : "
384 ost <<
"Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : "
389 double r_surf =
mp.
val_r(0,1.,M_PI/4,M_PI/4) ;
391 ost <<
"Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : "
392 <<
wit_w(0).val_point(r_surf,M_PI/4,M_PI/4) <<
" "
393 <<
wit_w(1).val_point(r_surf,M_PI/4,M_PI/4) <<
" "
394 <<
wit_w(2).val_point(r_surf,M_PI/4,M_PI/4) << endl ;
396 ost <<
"Central value of loggam : "
397 <<
loggam()(0, 0, 0, 0) << endl ;
400 ost <<
"Central value of lapse(N) auto : "
401 <<
n_auto()(0, 0, 0, 0) << endl ;
403 ost <<
"Central value of confpsi auto : "
406 ost <<
"Central value of shift (N^X, N^Y, N^Z) [c] : "
407 <<
shift(0)(0, 0, 0, 0) <<
" "
408 <<
shift(1)(0, 0, 0, 0) <<
" "
409 <<
shift(2)(0, 0, 0, 0) << endl ;
411 ost <<
" ... shift_auto part of it [c] : "
416 ost <<
" ... shift_comp part of it [c] : "
421 ost <<
" ... w_shift (NB: components in the star Cartesian frame) [c] : "
423 <<
w_shift(0)(0, 0, 0, 0) <<
" "
424 <<
w_shift(1)(0, 0, 0, 0) <<
" "
425 <<
w_shift(2)(0, 0, 0, 0) << endl ;
427 ost <<
"Central value of khi_shift [km c] : "
428 <<
khi_shift()(0, 0, 0, 0) / km << endl ;
430 ost << endl <<
"Central value of (B^X, B^Y, B^Z)/N [c] : "
431 <<
bsn(0)(0, 0, 0, 0) <<
" "
432 <<
bsn(1)(0, 0, 0, 0) <<
" "
433 <<
bsn(2)(0, 0, 0, 0) << endl ;
436 "Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : "
437 <<
d_n_auto(0)(0, 0, 0, 0) * km <<
" "
438 <<
d_n_auto(1)(0, 0, 0, 0) * km <<
" "
439 <<
d_n_auto(2)(0, 0, 0, 0) * km << endl ;
442 "Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : "