58 #include "et_bin_nsbh.h" 59 #include "graphique.h" 62 double Et_bin_nsbh::compute_angul()
const {
65 Cmp dx_mu (
nnn().dsdx() /
nnn());
66 Cmp dx_loggam (
loggam().dsdx()) ;
67 double part_dx = dx_mu(0,0,0,0) + dx_loggam(0,0,0,0) ;
77 G_square_cmp.std_base_scal() ;
78 double dG_square = G_square_cmp.dsdx()(0,0,0,0) ;
82 G_single_cmp.std_base_scal() ;
83 double dG_single = G_single_cmp.dsdx()(0,0,0,0) ;
87 G_const_cmp.std_base_scal() ;
88 double dG_const = G_const_cmp.dsdx()(0,0,0,0) ;
91 double G_square = G_square_cmp (0,0,0,0) ;
92 double G_single = G_single_cmp (0,0,0,0) ;
93 double G_const = G_const_cmp (0,0,0,0) ;
96 double a_coef = dG_square + 2*G_square*part_dx ;
97 double b_coef = dG_single + 2*G_single*part_dx ;
98 double c_coef = dG_const + 2*G_const *part_dx;
100 double determinant = b_coef*b_coef - 4*a_coef*c_coef ;
105 if (determinant <0) {
110 double sol_un = (-b_coef -
sqrt(determinant))/2/a_coef ;
111 double sol_deux = (-b_coef +
sqrt(determinant))/2/a_coef ;
113 bool signe_un = (sol_un >0) ?
true :
false ;
114 bool signe_deux = (sol_deux >0) ?
true :
false ;
117 if (signe_un == signe_deux) {
121 res = (signe_un) ? sol_un : sol_deux ;
Coord xa
Absolute x coordinate.
Cmp sqrt(const Cmp &)
Square root.
Tenseur nnn
Total lapse function.
Tenseur shift
Total shift vector.
Tenseur confpsi
Total conformal factor $$.
Map & mp
Mapping associated with the star.
Cmp pow(const Cmp &, int)
Power .
Coord ya
Absolute y coordinate.
Tenseur loggam
Logarithm of the Lorentz factor between the fluid and the co-orbiting observer.