Revision 1258 trunk/src/df1b2-separable/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