60 #include "hole_bhns.h" 72 bool bc_fs_i,
bool irrot_i,
double massbh)
74 bc_lapconf_nd(bc_nd_i),
75 bc_lapconf_fs(bc_fs_i),
76 irrotational(irrot_i),
77 lapconf_auto_rs(mp_i),
78 lapconf_auto_bh(mp_i),
84 d_lapconf_auto_rs(mp_i, COV, mp_i.get_bvect_cart()),
85 d_lapconf_auto_bh(mp_i, COV, mp_i.get_bvect_cart()),
86 d_lapconf_auto(mp_i, COV, mp_i.get_bvect_cart()),
87 d_lapconf_comp(mp_i, COV, mp_i.get_bvect_cart()),
88 shift_auto_rs(mp_i, CON, mp_i.get_bvect_cart()),
89 shift_auto_bh(mp_i, CON, mp_i.get_bvect_cart()),
90 shift_auto(mp_i, CON, mp_i.get_bvect_cart()),
91 shift_comp(mp_i, CON, mp_i.get_bvect_cart()),
92 shift_tot(mp_i, CON, mp_i.get_bvect_cart()),
93 d_shift_auto_rs(mp_i, 2, CON, mp_i.get_bvect_cart()),
94 d_shift_auto_bh(mp_i, 2, 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()),
102 d_confo_auto_rs(mp_i, COV, mp_i.get_bvect_cart()),
103 d_confo_auto_bh(mp_i, COV, mp_i.get_bvect_cart()),
104 d_confo_auto(mp_i, COV, mp_i.get_bvect_cart()),
105 d_confo_comp(mp_i, COV, mp_i.get_bvect_cart()),
106 taij_tot_rs(mp_i, CON, mp_i.get_bvect_cart()),
107 taij_tot_rot(mp_i, CON, mp_i.get_bvect_cart()),
108 taij_tot_bh(mp_i, CON, mp_i.get_bvect_cart()),
109 taij_tot(mp_i, CON, mp_i.get_bvect_cart()),
110 taij_auto_rs(mp_i, CON, mp_i.get_bvect_cart()),
111 taij_auto(mp_i, CON, mp_i.get_bvect_cart()),
112 taij_comp(mp_i, CON, mp_i.get_bvect_cart()),
113 taij_quad_tot_rs(mp_i),
114 taij_quad_tot_rot(mp_i),
115 taij_quad_tot_bh(mp_i),
117 taij_quad_auto(mp_i),
202 bc_lapconf_nd(hole.bc_lapconf_nd),
203 bc_lapconf_fs(hole.bc_lapconf_fs),
204 irrotational(hole.irrotational),
205 omega_spin(hole.omega_spin),
206 lapconf_auto_rs(hole.lapconf_auto_rs),
207 lapconf_auto_bh(hole.lapconf_auto_bh),
208 lapconf_auto(hole.lapconf_auto),
209 lapconf_comp(hole.lapconf_comp),
210 lapconf_tot(hole.lapconf_tot),
211 lapse_auto(hole.lapse_auto),
212 lapse_tot(hole.lapse_tot),
213 d_lapconf_auto_rs(hole.d_lapconf_auto_rs),
214 d_lapconf_auto_bh(hole.d_lapconf_auto_bh),
215 d_lapconf_auto(hole.d_lapconf_auto),
216 d_lapconf_comp(hole.d_lapconf_comp),
217 shift_auto_rs(hole.shift_auto_rs),
218 shift_auto_bh(hole.shift_auto_bh),
219 shift_auto(hole.shift_auto),
220 shift_comp(hole.shift_comp),
221 shift_tot(hole.shift_tot),
222 d_shift_auto_rs(hole.d_shift_auto_rs),
223 d_shift_auto_bh(hole.d_shift_auto_bh),
224 d_shift_auto(hole.d_shift_auto),
225 d_shift_comp(hole.d_shift_comp),
226 confo_auto_rs(hole.confo_auto_rs),
227 confo_auto_bh(hole.confo_auto_bh),
228 confo_auto(hole.confo_auto),
229 confo_comp(hole.confo_comp),
230 confo_tot(hole.confo_tot),
231 d_confo_auto_rs(hole.d_confo_auto_rs),
232 d_confo_auto_bh(hole.d_confo_auto_bh),
233 d_confo_auto(hole.d_confo_auto),
234 d_confo_comp(hole.d_confo_comp),
235 taij_tot_rs(hole.taij_tot_rs),
236 taij_tot_rot(hole.taij_tot_rot),
237 taij_tot_bh(hole.taij_tot_bh),
238 taij_tot(hole.taij_tot),
239 taij_auto_rs(hole.taij_auto_rs),
240 taij_auto(hole.taij_auto),
241 taij_comp(hole.taij_comp),
242 taij_quad_tot_rs(hole.taij_quad_tot_rs),
243 taij_quad_tot_rot(hole.taij_quad_tot_rot),
244 taij_quad_tot_bh(hole.taij_quad_tot_bh),
245 taij_quad_tot(hole.taij_quad_tot),
246 taij_quad_auto(hole.taij_quad_auto),
247 taij_quad_comp(hole.taij_quad_comp)
256 lapconf_auto_rs(mp_i, *(mp_i.get_mg()), fich),
257 lapconf_auto_bh(mp_i),
263 d_lapconf_auto_rs(mp_i, COV, mp_i.get_bvect_cart()),
264 d_lapconf_auto_bh(mp_i, COV, mp_i.get_bvect_cart()),
265 d_lapconf_auto(mp_i, COV, mp_i.get_bvect_cart()),
266 d_lapconf_comp(mp_i, COV, mp_i.get_bvect_cart()),
267 shift_auto_rs(mp_i, mp_i.get_bvect_cart(), fich),
268 shift_auto_bh(mp_i, CON, mp_i.get_bvect_cart()),
269 shift_auto(mp_i, CON, mp_i.get_bvect_cart()),
270 shift_comp(mp_i, CON, mp_i.get_bvect_cart()),
271 shift_tot(mp_i, CON, mp_i.get_bvect_cart()),
272 d_shift_auto_rs(mp_i, 2, CON, mp_i.get_bvect_cart()),
273 d_shift_auto_bh(mp_i, 2, CON, mp_i.get_bvect_cart()),
274 d_shift_auto(mp_i, 2, CON, mp_i.get_bvect_cart()),
275 d_shift_comp(mp_i, 2, CON, mp_i.get_bvect_cart()),
276 confo_auto_rs(mp_i, *(mp_i.get_mg()), fich),
281 d_confo_auto_rs(mp_i, COV, mp_i.get_bvect_cart()),
282 d_confo_auto_bh(mp_i, COV, mp_i.get_bvect_cart()),
283 d_confo_auto(mp_i, COV, mp_i.get_bvect_cart()),
284 d_confo_comp(mp_i, COV, mp_i.get_bvect_cart()),
285 taij_tot_rs(mp_i, CON, mp_i.get_bvect_cart()),
286 taij_tot_rot(mp_i, CON, mp_i.get_bvect_cart()),
287 taij_tot_bh(mp_i, CON, mp_i.get_bvect_cart()),
288 taij_tot(mp_i, CON, mp_i.get_bvect_cart()),
289 taij_auto_rs(mp_i, CON, mp_i.get_bvect_cart()),
290 taij_auto(mp_i, CON, mp_i.get_bvect_cart()),
291 taij_comp(mp_i, CON, mp_i.get_bvect_cart()),
292 taij_quad_tot_rs(mp_i),
293 taij_quad_tot_rot(mp_i),
294 taij_quad_tot_bh(mp_i),
296 taij_quad_auto(mp_i),
399 if (p_mass_irr_bhns != 0x0)
delete p_mass_irr_bhns ;
410 p_mass_irr_bhns = 0x0 ;
617 fwrite(&
omega_spin,
sizeof(
double), 1, fich) ;
630 ost <<
"Black hole in a BHNS binary" << endl ;
631 ost <<
"---------------------------" << endl ;
635 ost <<
"Irreducible mass of BH : " 637 ost <<
"Mass in the background : " 638 <<
mass_bh / msol <<
" [Mo]" << endl ;
639 ost <<
"Radius of the apparent horison : " 640 <<
rad_ah() / km <<
" [km]" << endl ;
641 ost <<
"Spin angular velocity : " 642 <<
omega_spin * f_unit <<
" [rad/s]" << endl ;
643 ost <<
"Lapse function on the AH : " 645 ost <<
"Conformal factor on the AH : " 647 ost <<
"shift(1) on the AH : " 648 <<
shift_tot(1).val_grid_point(1,0,nt-1,0) << endl ;
649 ost <<
"shift(2) on the AH : " 650 <<
shift_tot(2).val_grid_point(1,0,nt-1,0) << endl ;
651 ost <<
"shift(3) on the AH : " 652 <<
shift_tot(3).val_grid_point(1,0,nt-1,0) << endl ;
663 double relax_met,
int mer,
int fmer_met) {
665 double relax_met_jm1 = 1. - relax_met ;
667 if ( (mer != 0) && (mer % fmer_met == 0)) {
Map & mp
Mapping associated with the black hole.
Sym_tensor taij_auto
Part of the extrinsic curvature tensor generated by the black hole.
Scalar lapconf_auto_bh
Part of the lapconf function from the analytic background.
Vector & set_shift_comp()
Read/write of the shift vector generated by the companion star.
Scalar taij_quad_comp
Part of the scalar from the companion star.
void operator=(const Hole_bhns &)
Assignment to another Hole_bhns.
Vector d_confo_auto
Derivative of the conformal factor generated by the black hole.
virtual void del_deriv() const
Deletes all the derived quantities.
Sym_tensor taij_tot_rot
Part of the extrinsic curvature tensor from the rotation shift vector.
void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
Scalar confo_auto
Conformal factor generated by the black hole.
virtual void sauve(FILE *) const
Save in a file.
double mass_bh
Gravitational mass of BH.
Sym_tensor taij_tot_bh
Part of the extrinsic curvature tensor from the analytic background.
Standard units of space, time and mass.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Scalar & set_confo_tot()
Read/write of the total conformal factor.
Tensor field of valence 0 (or component of a tensorial field).
Vector d_lapconf_auto_bh
Derivative of the part of the lapconf function from the analytic background.
Base class for coordinate mappings.
Vector shift_auto
Shift vector generated by the black hole.
Scalar & set_lapse_tot()
Read/write of the total lapse function.
Vector d_lapconf_auto
Derivative of the lapconf function generated by the black hole.
Tensor d_shift_auto
Derivative of the shift vector generated by the black hole.
Sym_tensor taij_comp
Part of the extrinsic curvature tensor generated by the companion star.
Vector shift_tot
Total shift vector ;.
Vector & set_shift_auto_rs()
Read/write of the shift vector generated by the black hole.
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field. ...
Scalar & set_confo_comp()
Read/write of the conformal factor generated by the companion star.
virtual void sauve(FILE *) const
Save in a binary file.
bool irrotational
true for an irrotational black hole, false for a corotating one
void operator=(const Black_hole &)
Assignment to another Black_hole.
Tensor d_shift_comp
Derivative of the shift vector generated by the companion star.
void relax_bhns(const Hole_bhns &hole_prev, double relax_met, int mer, int fmer_met)
Performs a relaxation on lapse_auto , shift_auto , confo_auto .
Tensor field of valence 1.
Scalar & set_confo_auto_bh()
Read/write of the conformal factor generated by the black hole.
Sym_tensor taij_auto_rs
Part of the extrinsic curvature tensor numericalty computed for the black hole.
double * p_spin_am_bhns
Irreducible mass of BH.
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
Scalar lapconf_comp
Lapconf function generated by the companion star.
Hole_bhns(Map &mp_i, bool kerrschild_i, bool bc_lapconf_nd, bool bc_lapconf_fs, bool irrotational, double massbh)
Standard constructor.
Tensor d_shift_auto_bh
Derivative of the part of the shift vector from the analytic background.
Scalar & set_lapse_auto()
Read/write of the lapse function generated by the black hole.
Vector & set_shift_auto()
Read/write of the shift vector generated by the black hole.
Scalar confo_auto_rs
Part of the conformal factor from the numerical computation.
virtual void sauve(FILE *) const
Save in a file.
Scalar lapse_tot
Total lapse function.
Vector d_lapconf_comp
Derivative of the lapconf function generated by the companion star.
Scalar lapse_auto
Lapse function of the "black hole" part.
Vector shift_auto_bh
Part of the shift vector from the analytic background.
Scalar & set_lapconf_auto_rs()
Read/write of the lapconf function generated by the black hole.
virtual void sauve(FILE *) const
Save in a file.
Scalar confo_comp
Conformal factor generated by the companion star.
Vector & set_shift_tot()
Read/write of the total shift vector.
Scalar taij_quad_auto
Part of the scalar from the black hole.
Vector shift_auto_rs
Part of the shift vector from the numerical computation.
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<)
Base class for black holes.
Vector shift_comp
Shift vector generated by the companion star.
virtual double rad_ah() const
Radius of the apparent horizon.
Scalar & set_confo_auto()
Read/write of the conformal factor generated by the black hole.
Vector & set_shift_auto_bh()
Read/write of the shift vector generated by the black hole.
Scalar & set_lapconf_comp()
Read/write of the lapconf function generated by the companion star.
Vector d_lapconf_auto_rs
Derivative of the part of the lapconf function from the numerical computation.
Tensor d_shift_auto_rs
Derivative of the part of the shift vector from the numerical computation.
Vector d_confo_auto_bh
Derivative of the part of the conformal factor from the analytic background.
Class for black holes in black hole-neutron star binaries.
Sym_tensor taij_tot
Total extrinsic curvature tensor generated by shift_tot , lapse_tot , and confo_tot ...
Scalar lapconf_auto
Lapconf function generated by the black hole.
Scalar confo_tot
Total conformal factor.
Scalar & set_lapconf_auto()
Read/write of the lapconf function generated by the black hole.
Scalar taij_quad_tot
Total scalar generated by .
double omega_spin
Spin angular velocity of the black hole.
Vector d_confo_comp
Derivative of the conformal factor generated by the companion star.
Sym_tensor taij_tot_rs
Part of the extrinsic curvature tensor from the numerical computation.
Scalar taij_quad_tot_rot
Part of the scalar from the rotation shift vector.
Scalar lapconf_auto_rs
Part of the lapconf function from the numerical computation.
virtual double mass_irr_bhns() const
Irreducible mass of the black hole.
bool bc_lapconf_fs
true for the first type BC for the lapconf function, false for the second type BH ...
Vector d_confo_auto_rs
Derivative of the part of the conformal factor from the numerical computation.
Scalar & set_lapconf_tot()
Read/write of the total lapconf function.
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
bool bc_lapconf_nd
true for the Neumann type BC for the lapconf function, false for the Dirichlet type BH ...
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Scalar & set_confo_auto_rs()
Read/write of the conformal factor generated by the black hole.
Scalar taij_quad_tot_bh
Part of the scalar from the analytic background.
void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
virtual ~Hole_bhns()
Destructor.
Scalar & set_lapconf_auto_bh()
Read/write of the lapconf function generated by the black hole.
virtual void del_deriv() const
Deletes all the derived quantities.
Scalar lapconf_tot
Total lapconf function.
Scalar taij_quad_tot_rs
Part of the scalar from the numerical computation.
Scalar confo_auto_bh
Part of the conformal factor from the analytic background.