61 #include "blackhole.h" 62 #include "utilitaires.h" 66 double ff(
double,
const double) ;
82 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
94 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
100 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
106 double hh = 1./double(nn) ;
111 double x1, x2, x3, x4, x5 ;
119 for (
int i=0; i<mm; i++) {
121 x1 = hh * double(4*i) ;
122 x2 = hh * double(4*i+1) ;
123 x3 = hh * double(4*i+2) ;
124 x4 = hh * double(4*i+3) ;
125 x5 = hh * double(4*i+4) ;
127 integ += (hh/45.) * (14.*ff(x1,cc) + 64.*ff(x2,cc)
128 + 24.*ff(x3,cc) + 64.*ff(x4,cc)
133 rah = 2. *
exp(integ) ;
141 double ff(
double xx,
const double cc) {
143 double tcc2 = cc*cc/16. ;
144 double tmp =
sqrt(1. - xx + tcc2*
pow(xx, 4.)) ;
146 double resu = (-1. + tcc2 *
pow(xx, 3.)) / tmp / (1. + tmp) ;
Cmp exp(const Cmp &)
Exponential.
Cmp sqrt(const Cmp &)
Square root.
double rah_iso(bool neumann, bool first) const
Computes a radius of apparent horizon (excised surface) in isotropic coordinates. ...
Cmp pow(const Cmp &, int)
Power .