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
28char 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// --------------------
67namespace Lorene {
68Hole_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),
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)
245{
246 set_der_0x0() ;
247}
248
249// Constructor from a file
250// -----------------------
251Hole_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
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// -------------------------------
420
421 // Assignment of quantities common to the derived classes of Black_hole
423
424 // Assignment of proper quantities of class Hole_bhns
428 omega_spin = hole.omega_spin ;
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 ;
451 confo_auto = hole.confo_auto ;
452 confo_comp = hole.confo_comp ;
453 confo_tot = hole.confo_tot ;
458 taij_tot_rs = hole.taij_tot_rs ;
460 taij_tot_bh = hole.taij_tot_bh ;
461 taij_tot = hole.taij_tot ;
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// --------------
603void 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// --------
620ostream& 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
659void 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}
Base class for black holes.
Definition blackhole.h:74
Map & mp
Mapping associated with the black hole.
Definition blackhole.h:80
virtual double rad_ah() const
Radius of the apparent horizon.
void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
Definition blackhole.C:218
virtual void sauve(FILE *) const
Save in a file.
Definition blackhole.C:267
void operator=(const Black_hole &)
Assignment to another Black_hole.
Definition blackhole.C:236
double mass_bh
Gravitational mass of BH.
Definition blackhole.h:88
virtual void del_deriv() const
Deletes all the derived quantities.
Definition blackhole.C:205
Class for black holes in black hole-neutron star binaries.
Definition hole_bhns.h:65
Scalar confo_auto
Conformal factor generated by the black hole.
Definition hole_bhns.h:163
Sym_tensor taij_tot_bh
Part of the extrinsic curvature tensor from the analytic background.
Definition hole_bhns.h:200
Scalar lapconf_auto
Lapconf function generated by the black hole.
Definition hole_bhns.h:95
Vector shift_tot
Total shift vector ;.
Definition hole_bhns.h:138
Scalar & set_confo_comp()
Read/write of the conformal factor generated by the companion star.
Definition hole_bhns.C:582
Sym_tensor taij_comp
Part of the extrinsic curvature tensor generated by the companion star.
Definition hole_bhns.h:221
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
Scalar & set_lapconf_auto_bh()
Read/write of the lapconf function generated by the black hole.
Definition hole_bhns.C:484
Vector d_lapconf_auto_rs
Derivative of the part of the lapconf function from the numerical computation.
Definition hole_bhns.h:112
Sym_tensor taij_auto
Part of the extrinsic curvature tensor generated by the black hole.
Definition hole_bhns.h:216
Scalar confo_auto_bh
Part of the conformal factor from the analytic background.
Definition hole_bhns.h:160
Tensor d_shift_auto_rs
Derivative of the part of the shift vector from the numerical computation.
Definition hole_bhns.h:143
Sym_tensor taij_tot_rs
Part of the extrinsic curvature tensor from the numerical computation.
Definition hole_bhns.h:190
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<)
Definition hole_bhns.C:620
Scalar taij_quad_auto
Part of the scalar from the black hole.
Definition hole_bhns.h:238
virtual ~Hole_bhns()
Destructor.
Definition hole_bhns.C:380
virtual double mass_irr_bhns() const
Irreducible mass of the black hole.
Vector & set_shift_auto()
Read/write of the shift vector generated by the black hole.
Definition hole_bhns.C:540
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
Vector shift_auto_rs
Part of the shift vector from the numerical computation.
Definition hole_bhns.h:126
double omega_spin
Spin angular velocity of the black hole.
Definition hole_bhns.h:86
Vector & set_shift_auto_rs()
Read/write of the shift vector generated by the black hole.
Definition hole_bhns.C:526
Scalar & set_confo_auto_rs()
Read/write of the conformal factor generated by the black hole.
Definition hole_bhns.C:561
virtual void del_deriv() const
Deletes all the derived quantities.
Definition hole_bhns.C:392
Tensor d_shift_auto_bh
Derivative of the part of the shift vector from the analytic background.
Definition hole_bhns.h:148
Scalar & set_lapconf_auto_rs()
Read/write of the lapconf function generated by the black hole.
Definition hole_bhns.C:477
Scalar & set_confo_tot()
Read/write of the total conformal factor.
Definition hole_bhns.C:589
void operator=(const Hole_bhns &)
Assignment to another Hole_bhns.
Definition hole_bhns.C:419
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
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_lapconf_auto_bh
Derivative of the part of the lapconf function from the analytic background.
Definition hole_bhns.h:117
Scalar & set_confo_auto()
Read/write of the conformal factor generated by the black hole.
Definition hole_bhns.C:575
Scalar & set_lapconf_tot()
Read/write of the total lapconf function.
Definition hole_bhns.C:505
Vector d_confo_auto_rs
Derivative of the part of the conformal factor from the numerical computation.
Definition hole_bhns.h:174
Scalar confo_auto_rs
Part of the conformal factor from the numerical computation.
Definition hole_bhns.h:157
Scalar taij_quad_tot_rot
Part of the scalar from the rotation shift vector.
Definition hole_bhns.h:227
Vector & set_shift_tot()
Read/write of the total shift vector.
Definition hole_bhns.C:554
Scalar taij_quad_tot
Total scalar generated by .
Definition hole_bhns.h:235
Scalar lapse_tot
Total lapse function.
Definition hole_bhns.h:107
Vector & set_shift_auto_bh()
Read/write of the shift vector generated by the black hole.
Definition hole_bhns.C:533
Scalar & set_lapse_auto()
Read/write of the lapse function generated by the black hole.
Definition hole_bhns.C:512
Sym_tensor taij_tot_rot
Part of the extrinsic curvature tensor from the rotation shift vector.
Definition hole_bhns.h:195
Vector & set_shift_comp()
Read/write of the shift vector generated by the companion star.
Definition hole_bhns.C:547
Vector d_lapconf_comp
Derivative of the lapconf function generated by the companion star.
Definition hole_bhns.h:123
Sym_tensor taij_tot
Total extrinsic curvature tensor generated by shift_tot , lapse_tot , and confo_tot .
Definition hole_bhns.h:206
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 confo_comp
Conformal factor generated by the companion star.
Definition hole_bhns.h:166
Scalar lapconf_comp
Lapconf function generated by the companion star.
Definition hole_bhns.h:98
Scalar lapconf_auto_bh
Part of the lapconf function from the analytic background.
Definition hole_bhns.h:92
Vector d_confo_auto
Derivative of the conformal factor generated by the black hole.
Definition hole_bhns.h:182
virtual void sauve(FILE *) const
Save in a file.
Definition hole_bhns.C:603
Vector shift_auto
Shift vector generated by the black hole.
Definition hole_bhns.h:132
void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
Definition hole_bhns.C:403
Scalar lapconf_auto_rs
Part of the lapconf function from the numerical computation.
Definition hole_bhns.h:89
Scalar & set_lapconf_auto()
Read/write of the lapconf function generated by the black hole.
Definition hole_bhns.C:491
Vector d_confo_comp
Derivative of the conformal factor generated by the companion star.
Definition hole_bhns.h:185
Scalar taij_quad_comp
Part of the scalar from the companion star.
Definition hole_bhns.h:241
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
Scalar & set_lapconf_comp()
Read/write of the lapconf function generated by the companion star.
Definition hole_bhns.C:498
Scalar lapconf_tot
Total lapconf function.
Definition hole_bhns.h:101
Tensor d_shift_comp
Derivative of the shift vector generated by the companion star.
Definition hole_bhns.h:154
Vector d_confo_auto_bh
Derivative of the part of the conformal factor from the analytic background.
Definition hole_bhns.h:179
Vector shift_comp
Shift vector generated by the companion star.
Definition hole_bhns.h:135
Scalar confo_tot
Total conformal factor.
Definition hole_bhns.h:169
Scalar & set_confo_auto_bh()
Read/write of the conformal factor generated by the black hole.
Definition hole_bhns.C:568
double * p_spin_am_bhns
Irreducible mass of BH.
Definition hole_bhns.h:248
Sym_tensor taij_auto_rs
Part of the extrinsic curvature tensor numericalty computed for the black hole.
Definition hole_bhns.h:211
Scalar taij_quad_tot_rs
Part of the scalar from the numerical computation.
Definition hole_bhns.h:224
Scalar taij_quad_tot_bh
Part of the scalar from the analytic background.
Definition hole_bhns.h:230
bool irrotational
true for an irrotational black hole, false for a corotating one
Definition hole_bhns.h:83
Scalar & set_lapse_tot()
Read/write of the total lapse function.
Definition hole_bhns.C:519
Base class for coordinate mappings.
Definition map.h:670
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Definition map.h:765
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Definition grilles.h:457
Tensor field of valence 0 (or component of a tensorial field).
Definition scalar.h:387
virtual void sauve(FILE *) const
Save in a file.
Definition scalar.C:686
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
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
Tensor field of valence 1.
Definition vector.h:188
virtual void sauve(FILE *) const
Save in a binary file.
Definition tensor.C:906
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
Definition tensor.C:497
Lorene prototypes.
Definition app_hor.h:64
Standard units of space, time and mass.