28 char hole_bhns_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Hole_bhns/hole_bhns.C,v 1.4 2014/10/13 08:52:59 j_novak Exp $" ;
57 #include "hole_bhns.h"
69 bool bc_fs_i,
bool irrot_i,
double massbh)
71 bc_lapconf_nd(bc_nd_i),
72 bc_lapconf_fs(bc_fs_i),
73 irrotational(irrot_i),
74 lapconf_auto_rs(mp_i),
75 lapconf_auto_bh(mp_i),
81 d_lapconf_auto_rs(mp_i, COV, mp_i.get_bvect_cart()),
82 d_lapconf_auto_bh(mp_i, COV, mp_i.get_bvect_cart()),
83 d_lapconf_auto(mp_i, COV, mp_i.get_bvect_cart()),
84 d_lapconf_comp(mp_i, COV, mp_i.get_bvect_cart()),
85 shift_auto_rs(mp_i, CON, mp_i.get_bvect_cart()),
86 shift_auto_bh(mp_i, CON, mp_i.get_bvect_cart()),
87 shift_auto(mp_i, CON, mp_i.get_bvect_cart()),
88 shift_comp(mp_i, CON, mp_i.get_bvect_cart()),
89 shift_tot(mp_i, CON, mp_i.get_bvect_cart()),
90 d_shift_auto_rs(mp_i, 2, CON, mp_i.get_bvect_cart()),
91 d_shift_auto_bh(mp_i, 2, CON, mp_i.get_bvect_cart()),
92 d_shift_auto(mp_i, 2, CON, mp_i.get_bvect_cart()),
93 d_shift_comp(mp_i, 2, CON, mp_i.get_bvect_cart()),
99 d_confo_auto_rs(mp_i, COV, mp_i.get_bvect_cart()),
100 d_confo_auto_bh(mp_i, COV, mp_i.get_bvect_cart()),
101 d_confo_auto(mp_i, COV, mp_i.get_bvect_cart()),
102 d_confo_comp(mp_i, COV, mp_i.get_bvect_cart()),
103 taij_tot_rs(mp_i, CON, mp_i.get_bvect_cart()),
104 taij_tot_rot(mp_i, CON, mp_i.get_bvect_cart()),
105 taij_tot_bh(mp_i, CON, mp_i.get_bvect_cart()),
106 taij_tot(mp_i, CON, mp_i.get_bvect_cart()),
107 taij_auto_rs(mp_i, CON, mp_i.get_bvect_cart()),
108 taij_auto(mp_i, CON, mp_i.get_bvect_cart()),
109 taij_comp(mp_i, CON, mp_i.get_bvect_cart()),
110 taij_quad_tot_rs(mp_i),
111 taij_quad_tot_rot(mp_i),
112 taij_quad_tot_bh(mp_i),
114 taij_quad_auto(mp_i),
199 bc_lapconf_nd(hole.bc_lapconf_nd),
200 bc_lapconf_fs(hole.bc_lapconf_fs),
201 irrotational(hole.irrotational),
202 omega_spin(hole.omega_spin),
203 lapconf_auto_rs(hole.lapconf_auto_rs),
204 lapconf_auto_bh(hole.lapconf_auto_bh),
205 lapconf_auto(hole.lapconf_auto),
206 lapconf_comp(hole.lapconf_comp),
207 lapconf_tot(hole.lapconf_tot),
208 lapse_auto(hole.lapse_auto),
209 lapse_tot(hole.lapse_tot),
210 d_lapconf_auto_rs(hole.d_lapconf_auto_rs),
211 d_lapconf_auto_bh(hole.d_lapconf_auto_bh),
212 d_lapconf_auto(hole.d_lapconf_auto),
213 d_lapconf_comp(hole.d_lapconf_comp),
214 shift_auto_rs(hole.shift_auto_rs),
215 shift_auto_bh(hole.shift_auto_bh),
216 shift_auto(hole.shift_auto),
217 shift_comp(hole.shift_comp),
218 shift_tot(hole.shift_tot),
219 d_shift_auto_rs(hole.d_shift_auto_rs),
220 d_shift_auto_bh(hole.d_shift_auto_bh),
221 d_shift_auto(hole.d_shift_auto),
222 d_shift_comp(hole.d_shift_comp),
223 confo_auto_rs(hole.confo_auto_rs),
224 confo_auto_bh(hole.confo_auto_bh),
225 confo_auto(hole.confo_auto),
226 confo_comp(hole.confo_comp),
227 confo_tot(hole.confo_tot),
228 d_confo_auto_rs(hole.d_confo_auto_rs),
229 d_confo_auto_bh(hole.d_confo_auto_bh),
230 d_confo_auto(hole.d_confo_auto),
231 d_confo_comp(hole.d_confo_comp),
232 taij_tot_rs(hole.taij_tot_rs),
233 taij_tot_rot(hole.taij_tot_rot),
234 taij_tot_bh(hole.taij_tot_bh),
235 taij_tot(hole.taij_tot),
236 taij_auto_rs(hole.taij_auto_rs),
237 taij_auto(hole.taij_auto),
238 taij_comp(hole.taij_comp),
239 taij_quad_tot_rs(hole.taij_quad_tot_rs),
240 taij_quad_tot_rot(hole.taij_quad_tot_rot),
241 taij_quad_tot_bh(hole.taij_quad_tot_bh),
242 taij_quad_tot(hole.taij_quad_tot),
243 taij_quad_auto(hole.taij_quad_auto),
244 taij_quad_comp(hole.taij_quad_comp)
253 lapconf_auto_rs(mp_i, *(mp_i.get_mg()), fich),
254 lapconf_auto_bh(mp_i),
260 d_lapconf_auto_rs(mp_i, COV, mp_i.get_bvect_cart()),
261 d_lapconf_auto_bh(mp_i, COV, mp_i.get_bvect_cart()),
262 d_lapconf_auto(mp_i, COV, mp_i.get_bvect_cart()),
263 d_lapconf_comp(mp_i, COV, mp_i.get_bvect_cart()),
264 shift_auto_rs(mp_i, mp_i.get_bvect_cart(), fich),
265 shift_auto_bh(mp_i, CON, mp_i.get_bvect_cart()),
266 shift_auto(mp_i, CON, mp_i.get_bvect_cart()),
267 shift_comp(mp_i, CON, mp_i.get_bvect_cart()),
268 shift_tot(mp_i, CON, mp_i.get_bvect_cart()),
269 d_shift_auto_rs(mp_i, 2, CON, mp_i.get_bvect_cart()),
270 d_shift_auto_bh(mp_i, 2, CON, mp_i.get_bvect_cart()),
271 d_shift_auto(mp_i, 2, CON, mp_i.get_bvect_cart()),
272 d_shift_comp(mp_i, 2, CON, mp_i.get_bvect_cart()),
273 confo_auto_rs(mp_i, *(mp_i.get_mg()), fich),
278 d_confo_auto_rs(mp_i, COV, mp_i.get_bvect_cart()),
279 d_confo_auto_bh(mp_i, COV, mp_i.get_bvect_cart()),
280 d_confo_auto(mp_i, COV, mp_i.get_bvect_cart()),
281 d_confo_comp(mp_i, COV, mp_i.get_bvect_cart()),
282 taij_tot_rs(mp_i, CON, mp_i.get_bvect_cart()),
283 taij_tot_rot(mp_i, CON, mp_i.get_bvect_cart()),
284 taij_tot_bh(mp_i, CON, mp_i.get_bvect_cart()),
285 taij_tot(mp_i, CON, mp_i.get_bvect_cart()),
286 taij_auto_rs(mp_i, CON, mp_i.get_bvect_cart()),
287 taij_auto(mp_i, CON, mp_i.get_bvect_cart()),
288 taij_comp(mp_i, CON, mp_i.get_bvect_cart()),
289 taij_quad_tot_rs(mp_i),
290 taij_quad_tot_rot(mp_i),
291 taij_quad_tot_bh(mp_i),
293 taij_quad_auto(mp_i),
396 if (p_mass_irr_bhns != 0x0)
delete p_mass_irr_bhns ;
407 p_mass_irr_bhns = 0x0 ;
614 fwrite(&
omega_spin,
sizeof(
double), 1, fich) ;
627 ost <<
"Black hole in a BHNS binary" << endl ;
628 ost <<
"---------------------------" << endl ;
632 ost <<
"Irreducible mass of BH : "
634 ost <<
"Mass in the background : "
635 <<
mass_bh / msol <<
" [Mo]" << endl ;
636 ost <<
"Radius of the apparent horison : "
637 <<
rad_ah() / km <<
" [km]" << endl ;
638 ost <<
"Spin angular velocity : "
639 <<
omega_spin * f_unit <<
" [rad/s]" << endl ;
640 ost <<
"Lapse function on the AH : "
642 ost <<
"Conformal factor on the AH : "
644 ost <<
"shift(1) on the AH : "
645 <<
shift_tot(1).val_grid_point(1,0,nt-1,0) << endl ;
646 ost <<
"shift(2) on the AH : "
647 <<
shift_tot(2).val_grid_point(1,0,nt-1,0) << endl ;
648 ost <<
"shift(3) on the AH : "
649 <<
shift_tot(3).val_grid_point(1,0,nt-1,0) << endl ;
660 double relax_met,
int mer,
int fmer_met) {
662 double relax_met_jm1 = 1. - relax_met ;
664 if ( (mer != 0) && (mer % fmer_met == 0)) {