Revision 1258 trunk/src/df1b2separable/df1b2f20.cpp
df1b2f20.cpp (revision 1258)  

19  19 
ADUNCONST(df1b2variable,zz) 
20  20 
df1b2variable u; 
21  21 
double z = value(_zz); 
22 
double zdot=1.0; 

23 
double z2dot=0.0; 

24 
double z3dot=0.0; 

25 
const double lpi =1.1447298858494001741434272; 

26 
const double pi =3.1415926535897932384626432; 

22 
//double zdot=1.0; 

27  23 
const double lpp =0.9189385332046727417803297; 
28  24 
int n=7; 
29  25 
const double c[9]={0.99999999999980993, 
...  ...  
51  47 
x3dot=6.0*c[i]*fourth(zinv); 
52  48 
} 
53  49 
double t=z+n+0.5; 
54 
double tdot=zdot; 

55  50 
//return lpp + (z+0.5)*log(t) t + log(x); 
56  51 
double ans= lpp + (z+0.5)*log(t) t + log(x); 
52 
//double tdot=zdot; 

57  53 
//double dfx=zdot*log(t) + (z+0.5)/t*tdot tdot +xdot/x; 
58  54 
// since tdot=1.0 
59  55 
// since zdot=1.0 
Also available in: Unified diff