28 char star_bhns_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star_bhns/star_bhns.C,v 1.4 2014/10/13 08:53:40 j_novak Exp $" ;
59 #include "star_bhns.h"
65 Cmp raccord_c1(
const Cmp& uu,
int l1) ;
74 :
Star(mp_i, nzet_i, eos_i),
76 irrotational(irrot_i),
78 d_psi(mp_i, COV, mp_i.get_bvect_cart()),
79 wit_w(mp_i, CON, mp_i.get_bvect_cart()),
81 bsn(mp_i, CON, mp_i.get_bvect_cart()),
90 d_lapconf_auto(mp_i, COV, mp_i.get_bvect_cart()),
91 d_lapconf_comp(mp_i, COV, mp_i.get_bvect_cart()),
92 shift_auto(mp_i, CON, mp_i.get_bvect_cart()),
93 shift_comp(mp_i, CON, mp_i.get_bvect_cart()),
94 shift_tot(mp_i, CON, mp_i.get_bvect_cart()),
95 d_shift_auto(mp_i, 2, CON, mp_i.get_bvect_cart()),
96 d_shift_comp(mp_i, 2, CON, mp_i.get_bvect_cart()),
100 d_confo_auto(mp_i, COV, mp_i.get_bvect_cart()),
101 d_confo_comp(mp_i, COV, mp_i.get_bvect_cart()),
103 taij_auto(mp_i, CON, mp_i.get_bvect_cart()),
104 taij_quad_auto(mp_i),
105 flat(mp_i, mp_i.get_bvect_cart()),
109 ssjm1_wshift(mp_i, CON, mp_i.get_bvect_cart())
181 irrotational(star.irrotational),
189 pot_centri(star.pot_centri),
190 lapconf_auto(star.lapconf_auto),
191 lapconf_comp(star.lapconf_comp),
192 lapconf_tot(star.lapconf_tot),
193 lapse_auto(star.lapse_auto),
194 lapse_tot(star.lapse_tot),
195 d_lapconf_auto(star.d_lapconf_auto),
196 d_lapconf_comp(star.d_lapconf_comp),
197 shift_auto(star.shift_auto),
198 shift_comp(star.shift_comp),
199 shift_tot(star.shift_tot),
200 d_shift_auto(star.d_shift_auto),
201 d_shift_comp(star.d_shift_comp),
202 confo_auto(star.confo_auto),
203 confo_comp(star.confo_comp),
204 confo_tot(star.confo_tot),
205 d_confo_auto(star.d_confo_auto),
206 d_confo_comp(star.d_confo_comp),
208 taij_auto(star.taij_auto),
209 taij_quad_auto(star.taij_quad_auto),
211 ssjm1_lapconf(star.ssjm1_lapconf),
212 ssjm1_confo(star.ssjm1_confo),
213 ssjm1_khi(star.ssjm1_khi),
214 ssjm1_wshift(star.ssjm1_wshift)
222 :
Star(mp_i, eos_i, fich),
225 d_psi(mp_i, COV, mp_i.get_bvect_cart()),
226 wit_w(mp_i, CON, mp_i.get_bvect_cart()),
228 bsn(mp_i, CON, mp_i.get_bvect_cart()),
232 lapconf_auto(mp_i, *(mp_i.get_mg()), fich),
237 d_lapconf_auto(mp_i, COV, mp_i.get_bvect_cart()),
238 d_lapconf_comp(mp_i, COV, mp_i.get_bvect_cart()),
239 shift_auto(mp_i, mp_i.get_bvect_cart(), fich),
240 shift_comp(mp_i, CON, mp_i.get_bvect_cart()),
241 shift_tot(mp_i, CON, mp_i.get_bvect_cart()),
242 d_shift_auto(mp_i, 2, CON, mp_i.get_bvect_cart()),
243 d_shift_comp(mp_i, 2, CON, mp_i.get_bvect_cart()),
244 confo_auto(mp_i, *(mp_i.get_mg()), fich),
247 d_confo_auto(mp_i, COV, mp_i.get_bvect_cart()),
248 d_confo_comp(mp_i, COV, mp_i.get_bvect_cart()),
250 taij_auto(mp_i, CON, mp_i.get_bvect_cart()),
251 taij_quad_auto(mp_i),
252 flat(mp_i, mp_i.get_bvect_cart()),
253 ssjm1_lapconf(mp_i, *(mp_i.get_mg()), fich),
254 ssjm1_confo(mp_i, *(mp_i.get_mg()), fich),
255 ssjm1_khi(mp_i, *(mp_i.get_mg()), fich),
256 ssjm1_wshift(mp_i, mp_i.get_bvect_cart(), fich)
502 ost <<
"Neutron star in a BHNS binary" << endl ;
503 ost <<
"-----------------------------" << endl ;
505 ost <<
"Coordinate radius R_eq_tow : "
507 ost <<
"Coordinate radius R_eq_opp : "
508 <<
ray_eq() / km <<
" [km]" << endl ;
509 ost <<
"Coordinate radius R_eq_orb : "
511 ost <<
"Coordinate radius R_pole : "
512 <<
ray_pole() / km <<
" [km]" << endl ;
513 ost <<
"Central enthalph H_ent : "
515 ost <<
"Lapse function at the center of NS : "
517 ost <<
"Conformal factor at the center of NS : "
519 ost <<
"shift(1) at the center of NS : "
520 <<
shift_tot(1).val_grid_point(0,0,0,0) << endl ;
521 ost <<
"shift(2) at the center of NS : "
522 <<
shift_tot(2).val_grid_point(0,0,0,0) << endl ;
523 ost <<
"shift(3) at the center of NS : "
524 <<
shift_tot(3).val_grid_point(0,0,0,0) << endl ;
557 for (
int i=1; i<=3; i++) {
558 v_orb.
set(i) = www(i).val_grid_point(0,0,0,0) ;
566 for (
int i=1; i<=3; i++)
571 d_psi = d_psi0 + v_orb ;
573 for (
int i=1; i<=3; i++) {
574 if (
d_psi(i).get_etat() == ETATZERO)
584 for (
int i=1; i<=3; i++) {
586 d_psi_i.
va.
set_base( d_psi0(i).get_spectral_va().base ) ;
587 d_psi_i = raccord_c1(d_psi_i,
nzet) ;
595 double relax_met,
int mer,
int fmer_met) {
597 double relax_ent_jm1 = 1. - relax_ent ;
598 double relax_met_jm1 = 1. - relax_met ;
600 ent = relax_ent *
ent + relax_ent_jm1 * star_prev.
ent ;
602 if ( (mer != 0) && (mer % fmer_met == 0)) {