29 char blackhole_rah_iso_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Black_hole/blackhole_rah_iso.C,v 1.4 2014/10/13 08:52:46 j_novak Exp $" ;
58 #include "blackhole.h" 59 #include "utilitaires.h" 63 double ff(
double,
const double) ;
79 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
91 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
97 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
103 double hh = 1./double(nn) ;
108 double x1, x2, x3, x4, x5 ;
116 for (
int i=0; i<mm; i++) {
118 x1 = hh * double(4*i) ;
119 x2 = hh * double(4*i+1) ;
120 x3 = hh * double(4*i+2) ;
121 x4 = hh * double(4*i+3) ;
122 x5 = hh * double(4*i+4) ;
124 integ += (hh/45.) * (14.*ff(x1,cc) + 64.*ff(x2,cc)
125 + 24.*ff(x3,cc) + 64.*ff(x4,cc)
130 rah = 2. *
exp(integ) ;
138 double ff(
double xx,
const double cc) {
140 double tcc2 = cc*cc/16. ;
141 double tmp =
sqrt(1. - xx + tcc2*
pow(xx, 4.)) ;
143 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 .