GremlinEq
specialradial.h
1 //---------------------------------------------------------------------------
2 //
3 // $Id: specialradial.h,v 1.1 2013/05/17 15:30:12 sahughes Exp $
4 //
5 //---------------------------------------------------------------------------
6 
7 
8 /* specialradial.h
9  * Functions to calculate properties of the radial continued fraction
10  * at special points
11  */
12 
13 #ifndef SPECIALRADIAL_H_SEEN
14 #define SPECIALRADIAL_H_SEEN
15 
16 #include "teukolskydefs.h"
17 
18 /* estimate of the precision of the results of the special functions, as
19  * a multiple of their term sizes */
20 #define SPECIAL_TERM_NOISE 3e-15
21 
22 /* calculate the value of the function at nu = -1/2
23  * If term is not NULL, the address it points to will be set to the magnitude
24  * of one of the terms in the final addition */
25 Real radial_half_value(Real epsilon, Real q, int m, Real lambda, Real *term);
26 
27 /* calculate the slope of the function at nu = 1
28  * If term is not NULL, the address it points to will be set to the magnitude
29  * of one of the terms in the final addition */
30 Real radial_int_slope(Real epsilon, Real q, int m, Real lambda, Real *term);
31 
32 #endif /* SPECIALRADIAL_H_SEEN */