LORENE
hole_bhns.C
1 /*
2  * Methods of class Hole_bhns
3  *
4  * (see file hole_bhns.h for documentation).
5  *
6  */
7 
8 /*
9  * Copyright (c) 2005-2007 Keisuke Taniguchi
10  *
11  * This file is part of LORENE.
12  *
13  * LORENE is free software; you can redistribute it and/or modify
14  * it under the terms of the GNU General Public License version 2
15  * as published by the Free Software Foundation.
16  *
17  * LORENE is distributed in the hope that it will be useful,
18  * but WITHOUT ANY WARRANTY; without even the implied warranty of
19  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
20  * GNU General Public License for more details.
21  *
22  * You should have received a copy of the GNU General Public License
23  * along with LORENE; if not, write to the Free Software
24  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
25  *
26  */
27 
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 $" ;
29 
30 /*
31  * $Id: hole_bhns.C,v 1.4 2014/10/13 08:52:59 j_novak Exp $
32  * $Log: hole_bhns.C,v $
33  * Revision 1.4 2014/10/13 08:52:59 j_novak
34  * Lorene classes and functions now belong to the namespace Lorene.
35  *
36  * Revision 1.3 2014/10/06 15:13:10 j_novak
37  * Modified #include directives to use c++ syntax.
38  *
39  * Revision 1.2 2008/05/15 19:03:00 k_taniguchi
40  * Change of some parameters.
41  *
42  * Revision 1.1 2007/06/22 01:24:16 k_taniguchi
43  * *** empty log message ***
44  *
45  *
46  * $Header: /cvsroot/Lorene/C++/Source/Hole_bhns/hole_bhns.C,v 1.4 2014/10/13 08:52:59 j_novak Exp $
47  *
48  */
49 
50 // C++ headers
51 //#include <>
52 
53 // C headers
54 #include <cmath>
55 
56 // Lorene headers
57 #include "hole_bhns.h"
58 #include "unites.h"
59 
60 
61  //----------------------//
62  // Constructors //
63  //----------------------//
64 
65 // Standard constructor
66 // --------------------
67 namespace Lorene {
68 Hole_bhns::Hole_bhns(Map& mp_i, bool kerrschild_i, bool bc_nd_i,
69  bool bc_fs_i, bool irrot_i, double massbh)
70  : Black_hole(mp_i, kerrschild_i, 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),
76  lapconf_auto(mp_i),
77  lapconf_comp(mp_i),
78  lapconf_tot(mp_i),
79  lapse_auto(mp_i),
80  lapse_tot(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()),
94  confo_auto_rs(mp_i),
95  confo_auto_bh(mp_i),
96  confo_auto(mp_i),
97  confo_comp(mp_i),
98  confo_tot(mp_i),
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),
113  taij_quad_tot(mp_i),
114  taij_quad_auto(mp_i),
115  taij_quad_comp(mp_i)
116 {
117 
118  omega_spin = 0. ;
119 
120  // The metric quantities are initialized to the flat one or zero
121  lapconf_auto_rs = 0. ;
123  lapconf_auto_bh = 1. ;
125  lapconf_auto = 1. ;
127  lapconf_comp = 0. ;
129  lapconf_tot = 1. ;
131 
132  lapse_auto = 1. ;
134  lapse_tot = 1. ;
136 
141 
147 
152 
153  confo_auto_rs = 0. ;
155  confo_auto_bh = 1. ;
157  confo_auto = 1. ;
159  confo_comp = 0. ;
161  confo_tot = 1. ;
163 
168 
176 
177  taij_quad_tot_rs = 0. ;
179  taij_quad_tot_rot = 0. ;
181  taij_quad_tot_bh = 0. ;
183  taij_quad_tot = 0. ;
185  taij_quad_auto = 0. ;
187  taij_quad_comp = 0. ;
189 
190  // Pointers of derived quantities initialized to zero :
191  set_der_0x0() ;
192 
193 }
194 
195 // Copy constructor
196 // ----------------
198  : Black_hole(hole),
202  omega_spin(hole.omega_spin),
207  lapconf_tot(hole.lapconf_tot),
208  lapse_auto(hole.lapse_auto),
209  lapse_tot(hole.lapse_tot),
216  shift_auto(hole.shift_auto),
217  shift_comp(hole.shift_comp),
218  shift_tot(hole.shift_tot),
225  confo_auto(hole.confo_auto),
226  confo_comp(hole.confo_comp),
227  confo_tot(hole.confo_tot),
232  taij_tot_rs(hole.taij_tot_rs),
234  taij_tot_bh(hole.taij_tot_bh),
235  taij_tot(hole.taij_tot),
237  taij_auto(hole.taij_auto),
238  taij_comp(hole.taij_comp),
245 {
246  set_der_0x0() ;
247 }
248 
249 // Constructor from a file
250 // -----------------------
251 Hole_bhns::Hole_bhns(Map& mp_i, FILE* fich)
252  : Black_hole(mp_i, fich),
253  lapconf_auto_rs(mp_i, *(mp_i.get_mg()), fich),
254  lapconf_auto_bh(mp_i),
255  lapconf_auto(mp_i),
256  lapconf_comp(mp_i),
257  lapconf_tot(mp_i),
258  lapse_auto(mp_i),
259  lapse_tot(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),
274  confo_auto_bh(mp_i),
275  confo_auto(mp_i),
276  confo_comp(mp_i),
277  confo_tot(mp_i),
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),
292  taij_quad_tot(mp_i),
293  taij_quad_auto(mp_i),
294  taij_quad_comp(mp_i)
295 {
296 
297  // bc_lapconf_nd, bc_lapconf_fs, irrotational, and omega_spin
298  // are read from the file
299  fread(&bc_lapconf_nd, sizeof(bool), 1, fich) ;
300  fread(&bc_lapconf_fs, sizeof(bool), 1, fich) ;
301  fread(&irrotational, sizeof(bool), 1, fich) ;
302  fread(&omega_spin, sizeof(double), 1, fich) ;
303 
304  // All other quantities are initialized to zero
305  // --------------------------------------------
306 
307  lapconf_auto_bh = 1. ;
309  lapconf_auto = 1. ;
311  lapconf_comp = 0. ;
313  lapconf_tot = 1. ;
315 
316  lapse_auto = 1. ;
318  lapse_tot = 1. ;
320 
325 
334 
335  confo_auto_bh = 1. ;
337  confo_auto = 1. ;
339  confo_comp = 0. ;
341  confo_tot = 1. ;
343 
348 
356  taij_quad_tot_rs = 0. ;
358  taij_quad_tot_rot = 0. ;
360  taij_quad_tot_bh = 0. ;
362  taij_quad_tot = 0. ;
364  taij_quad_auto = 0. ;
366  taij_quad_comp = 0. ;
368 
369  // Pointers of derived quantities initialized to zero
370  // --------------------------------------------------
371  set_der_0x0() ;
372 
373 }
374 
375 
376  //--------------------//
377  // Destructor //
378  //--------------------//
379 
381 {
382 
383  del_deriv() ;
384 
385 }
386 
387 
388  //------------------------------------------//
389  // Management of derived quantities //
390  //------------------------------------------//
391 
392 void Hole_bhns::del_deriv() const {
393 
395 
396  if (p_mass_irr_bhns != 0x0) delete p_mass_irr_bhns ;
397  if (p_spin_am_bhns != 0x0) delete p_spin_am_bhns ;
398 
399  set_der_0x0() ;
400 
401 }
402 
404 
406 
407  p_mass_irr_bhns = 0x0 ;
408  p_spin_am_bhns = 0x0 ;
409 
410 }
411 
412 
413  //--------------------//
414  // Assignment //
415  //--------------------//
416 
417 // Assignment to another Hole_bhns
418 // -------------------------------
419 void Hole_bhns::operator=(const Hole_bhns& hole) {
420 
421  // Assignment of quantities common to the derived classes of Black_hole
422  Black_hole::operator=(hole) ;
423 
424  // Assignment of proper quantities of class Hole_bhns
427  irrotational = hole.irrotational ;
428  omega_spin = hole.omega_spin ;
431  lapconf_auto = hole.lapconf_auto ;
432  lapconf_comp = hole.lapconf_comp ;
433  lapconf_tot = hole.lapconf_tot ;
434  lapse_auto = hole.lapse_auto ;
435  lapse_tot = hole.lapse_tot ;
442  shift_auto = hole.shift_auto ;
443  shift_comp = hole.shift_comp ;
444  shift_tot = hole.shift_tot ;
447  d_shift_auto = hole.d_shift_auto ;
448  d_shift_comp = hole.d_shift_comp ;
451  confo_auto = hole.confo_auto ;
452  confo_comp = hole.confo_comp ;
453  confo_tot = hole.confo_tot ;
456  d_confo_auto = hole.d_confo_auto ;
457  d_confo_comp = hole.d_confo_comp ;
458  taij_tot_rs = hole.taij_tot_rs ;
459  taij_tot_rot = hole.taij_tot_rot ;
460  taij_tot_bh = hole.taij_tot_bh ;
461  taij_tot = hole.taij_tot ;
462  taij_auto_rs = hole.taij_auto_rs ;
463  taij_auto = hole.taij_auto ;
464  taij_comp = hole.taij_comp ;
471 
472  del_deriv() ;
473 
474 }
475 
476 
478 
479  del_deriv() ;
480  return lapconf_auto_rs ;
481 
482 }
483 
485 
486  del_deriv() ;
487  return lapconf_auto_bh ;
488 
489 }
490 
492 
493  del_deriv() ;
494  return lapconf_auto ;
495 
496 }
497 
499 
500  del_deriv() ;
501  return lapconf_comp ;
502 
503 }
504 
506 
507  del_deriv() ;
508  return lapconf_tot ;
509 
510 }
511 
513 
514  del_deriv() ;
515  return lapse_auto ;
516 
517 }
518 
520 
521  del_deriv() ;
522  return lapse_tot ;
523 
524 }
525 
527 
528  del_deriv() ;
529  return shift_auto_rs ;
530 
531 }
532 
534 
535  del_deriv() ;
536  return shift_auto_bh ;
537 
538 }
539 
541 
542  del_deriv() ;
543  return shift_auto ;
544 
545 }
546 
548 
549  del_deriv() ;
550  return shift_comp ;
551 
552 }
553 
555 
556  del_deriv() ;
557  return shift_tot ;
558 
559 }
560 
562 
563  del_deriv() ;
564  return confo_auto_rs ;
565 
566 }
567 
569 
570  del_deriv() ;
571  return confo_auto_bh ;
572 
573 }
574 
576 
577  del_deriv() ;
578  return confo_auto ;
579 
580 }
581 
583 
584  del_deriv() ;
585  return confo_comp ;
586 
587 }
588 
590 
591  del_deriv() ;
592  return confo_tot ;
593 
594 }
595 
596 
597  //-----------------//
598  // Outputs //
599  //-----------------//
600 
601 // Save in a file
602 // --------------
603 void Hole_bhns::sauve(FILE* fich) const {
604 
605  Black_hole::sauve(fich) ;
606 
607  lapconf_auto_rs.sauve(fich) ;
608  shift_auto_rs.sauve(fich) ;
609  confo_auto_rs.sauve(fich) ;
610 
611  fwrite(&bc_lapconf_nd, sizeof(bool), 1, fich) ;
612  fwrite(&bc_lapconf_fs, sizeof(bool), 1, fich) ;
613  fwrite(&irrotational, sizeof(bool), 1, fich) ;
614  fwrite(&omega_spin, sizeof(double), 1, fich) ;
615 
616 }
617 
618 // Printing
619 // --------
620 ostream& Hole_bhns::operator>>(ostream& ost) const {
621 
622  using namespace Unites ;
623 
624  // Black_hole::operator>>(ost) ;
625 
626  ost << endl ;
627  ost << "Black hole in a BHNS binary" << endl ;
628  ost << "---------------------------" << endl ;
629 
630  int nt = mp.get_mg()->get_nt(1) ;
631 
632  ost << "Irreducible mass of BH : "
633  << mass_irr_bhns() / msol << " [Mo]" << endl ;
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 : "
641  << lapse_tot.val_grid_point(1,0,nt-1,0) << endl ;
642  ost << "Conformal factor on the AH : "
643  << confo_tot.val_grid_point(1,0,nt-1,0) << endl ;
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 ;
650 
651  return ost ;
652 
653 }
654 
655  //--------------------------------//
656  // Computational routines //
657  //--------------------------------//
658 
659 void Hole_bhns::relax_bhns(const Hole_bhns& hole_prev,
660  double relax_met, int mer, int fmer_met) {
661 
662  double relax_met_jm1 = 1. - relax_met ;
663 
664  if ( (mer != 0) && (mer % fmer_met == 0)) {
665 
666  lapconf_auto_rs = relax_met * lapconf_auto_rs
667  + relax_met_jm1 * hole_prev.lapconf_auto_rs ;
668 
669  shift_auto_rs = relax_met * shift_auto_rs
670  + relax_met_jm1 * hole_prev.shift_auto_rs ;
671 
672  confo_auto_rs = relax_met * confo_auto_rs
673  + relax_met_jm1 * hole_prev.confo_auto_rs ;
674 
675  }
676 
677  del_deriv() ;
678 
679 }
680 
681 }
Map & mp
Mapping associated with the black hole.
Definition: blackhole.h:80
Sym_tensor taij_auto
Part of the extrinsic curvature tensor generated by the black hole.
Definition: hole_bhns.h:216
Scalar lapconf_auto_bh
Part of the lapconf function from the analytic background.
Definition: hole_bhns.h:92
Vector & set_shift_comp()
Read/write of the shift vector generated by the companion star.
Definition: hole_bhns.C:547
Scalar taij_quad_comp
Part of the scalar from the companion star.
Definition: hole_bhns.h:241
void operator=(const Hole_bhns &)
Assignment to another Hole_bhns.
Definition: hole_bhns.C:419
Vector d_confo_auto
Derivative of the conformal factor generated by the black hole.
Definition: hole_bhns.h:182
virtual void del_deriv() const
Deletes all the derived quantities.
Definition: blackhole.C:205
Sym_tensor taij_tot_rot
Part of the extrinsic curvature tensor from the rotation shift vector.
Definition: hole_bhns.h:195
void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
Definition: hole_bhns.C:403
Scalar confo_auto
Conformal factor generated by the black hole.
Definition: hole_bhns.h:163
Lorene prototypes.
Definition: app_hor.h:64
virtual void sauve(FILE *) const
Save in a file.
Definition: blackhole.C:267
double mass_bh
Gravitational mass of BH.
Definition: blackhole.h:88
Sym_tensor taij_tot_bh
Part of the extrinsic curvature tensor from the analytic background.
Definition: hole_bhns.h:200
Standard units of space, time and mass.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Definition: map.h:765
Scalar & set_confo_tot()
Read/write of the total conformal factor.
Definition: hole_bhns.C:589
Tensor field of valence 0 (or component of a tensorial field).
Definition: scalar.h:387
Vector d_lapconf_auto_bh
Derivative of the part of the lapconf function from the analytic background.
Definition: hole_bhns.h:117
Base class for coordinate mappings.
Definition: map.h:670
Vector shift_auto
Shift vector generated by the black hole.
Definition: hole_bhns.h:132
Scalar & set_lapse_tot()
Read/write of the total lapse function.
Definition: hole_bhns.C:519
Vector d_lapconf_auto
Derivative of the lapconf function generated by the black hole.
Definition: hole_bhns.h:120
Tensor d_shift_auto
Derivative of the shift vector generated by the black hole.
Definition: hole_bhns.h:151
Sym_tensor taij_comp
Part of the extrinsic curvature tensor generated by the companion star.
Definition: hole_bhns.h:221
Vector shift_tot
Total shift vector ;.
Definition: hole_bhns.h:138
Vector & set_shift_auto_rs()
Read/write of the shift vector generated by the black hole.
Definition: hole_bhns.C:526
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field. ...
Definition: scalar.C:784
Scalar & set_confo_comp()
Read/write of the conformal factor generated by the companion star.
Definition: hole_bhns.C:582
virtual void sauve(FILE *) const
Save in a binary file.
Definition: tensor.C:906
bool irrotational
true for an irrotational black hole, false for a corotating one
Definition: hole_bhns.h:83
void operator=(const Black_hole &)
Assignment to another Black_hole.
Definition: blackhole.C:236
Tensor d_shift_comp
Derivative of the shift vector generated by the companion star.
Definition: hole_bhns.h:154
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 .
Definition: hole_bhns.C:659
Tensor field of valence 1.
Definition: vector.h:188
Scalar & set_confo_auto_bh()
Read/write of the conformal factor generated by the black hole.
Definition: hole_bhns.C:568
Sym_tensor taij_auto_rs
Part of the extrinsic curvature tensor numericalty computed for the black hole.
Definition: hole_bhns.h:211
double * p_spin_am_bhns
Irreducible mass of BH.
Definition: hole_bhns.h:248
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
Definition: scalar.h:637
Scalar lapconf_comp
Lapconf function generated by the companion star.
Definition: hole_bhns.h:98
Hole_bhns(Map &mp_i, bool kerrschild_i, bool bc_lapconf_nd, bool bc_lapconf_fs, bool irrotational, double massbh)
Standard constructor.
Definition: hole_bhns.C:68
Tensor d_shift_auto_bh
Derivative of the part of the shift vector from the analytic background.
Definition: hole_bhns.h:148
Scalar & set_lapse_auto()
Read/write of the lapse function generated by the black hole.
Definition: hole_bhns.C:512
Vector & set_shift_auto()
Read/write of the shift vector generated by the black hole.
Definition: hole_bhns.C:540
Scalar confo_auto_rs
Part of the conformal factor from the numerical computation.
Definition: hole_bhns.h:157
virtual void sauve(FILE *) const
Save in a file.
Definition: scalar.C:686
Scalar lapse_tot
Total lapse function.
Definition: hole_bhns.h:107
Vector d_lapconf_comp
Derivative of the lapconf function generated by the companion star.
Definition: hole_bhns.h:123
Scalar lapse_auto
Lapse function of the "black hole" part.
Definition: hole_bhns.h:104
Vector shift_auto_bh
Part of the shift vector from the analytic background.
Definition: hole_bhns.h:129
Scalar & set_lapconf_auto_rs()
Read/write of the lapconf function generated by the black hole.
Definition: hole_bhns.C:477
virtual void sauve(FILE *) const
Save in a file.
Definition: hole_bhns.C:603
Scalar confo_comp
Conformal factor generated by the companion star.
Definition: hole_bhns.h:166
Vector & set_shift_tot()
Read/write of the total shift vector.
Definition: hole_bhns.C:554
Scalar taij_quad_auto
Part of the scalar from the black hole.
Definition: hole_bhns.h:238
Vector shift_auto_rs
Part of the shift vector from the numerical computation.
Definition: hole_bhns.h:126
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<)
Definition: hole_bhns.C:620
Base class for black holes.
Definition: blackhole.h:74
Vector shift_comp
Shift vector generated by the companion star.
Definition: hole_bhns.h:135
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.
Definition: hole_bhns.C:575
Vector & set_shift_auto_bh()
Read/write of the shift vector generated by the black hole.
Definition: hole_bhns.C:533
Scalar & set_lapconf_comp()
Read/write of the lapconf function generated by the companion star.
Definition: hole_bhns.C:498
Vector d_lapconf_auto_rs
Derivative of the part of the lapconf function from the numerical computation.
Definition: hole_bhns.h:112
Tensor d_shift_auto_rs
Derivative of the part of the shift vector from the numerical computation.
Definition: hole_bhns.h:143
Vector d_confo_auto_bh
Derivative of the part of the conformal factor from the analytic background.
Definition: hole_bhns.h:179
Class for black holes in black hole-neutron star binaries.
Definition: hole_bhns.h:65
Sym_tensor taij_tot
Total extrinsic curvature tensor generated by shift_tot , lapse_tot , and confo_tot ...
Definition: hole_bhns.h:206
Scalar lapconf_auto
Lapconf function generated by the black hole.
Definition: hole_bhns.h:95
Scalar confo_tot
Total conformal factor.
Definition: hole_bhns.h:169
Scalar & set_lapconf_auto()
Read/write of the lapconf function generated by the black hole.
Definition: hole_bhns.C:491
Scalar taij_quad_tot
Total scalar generated by .
Definition: hole_bhns.h:235
double omega_spin
Spin angular velocity of the black hole.
Definition: hole_bhns.h:86
Vector d_confo_comp
Derivative of the conformal factor generated by the companion star.
Definition: hole_bhns.h:185
Sym_tensor taij_tot_rs
Part of the extrinsic curvature tensor from the numerical computation.
Definition: hole_bhns.h:190
Scalar taij_quad_tot_rot
Part of the scalar from the rotation shift vector.
Definition: hole_bhns.h:227
Scalar lapconf_auto_rs
Part of the lapconf function from the numerical computation.
Definition: hole_bhns.h:89
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 ...
Definition: hole_bhns.h:78
Vector d_confo_auto_rs
Derivative of the part of the conformal factor from the numerical computation.
Definition: hole_bhns.h:174
Scalar & set_lapconf_tot()
Read/write of the total lapconf function.
Definition: hole_bhns.C:505
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
Definition: tensor.C:497
bool bc_lapconf_nd
true for the Neumann type BC for the lapconf function, false for the Dirichlet type BH ...
Definition: hole_bhns.h:73
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Definition: grilles.h:457
Scalar & set_confo_auto_rs()
Read/write of the conformal factor generated by the black hole.
Definition: hole_bhns.C:561
Scalar taij_quad_tot_bh
Part of the scalar from the analytic background.
Definition: hole_bhns.h:230
void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
Definition: blackhole.C:218
virtual ~Hole_bhns()
Destructor.
Definition: hole_bhns.C:380
Scalar & set_lapconf_auto_bh()
Read/write of the lapconf function generated by the black hole.
Definition: hole_bhns.C:484
virtual void del_deriv() const
Deletes all the derived quantities.
Definition: hole_bhns.C:392
Scalar lapconf_tot
Total lapconf function.
Definition: hole_bhns.h:101
Scalar taij_quad_tot_rs
Part of the scalar from the numerical computation.
Definition: hole_bhns.h:224
Scalar confo_auto_bh
Part of the conformal factor from the analytic background.
Definition: hole_bhns.h:160