Revision 1245 trunk/src/linad99/combv.cpp

combv.cpp (revision 1245)
240 240
static dvariable gammlnguts(const prevariable _z)
241 241
{
242 242
  double  z = value(_z);
243
  double zdot=1.0;
244
  const double lpi =1.1447298858494001741434272;
245
  const double pi =3.1415926535897932384626432;
243
  //double zdot=1.0;
244
  //const double lpi =1.1447298858494001741434272;
245
  //const double pi =3.1415926535897932384626432;
246 246
  const double lpp =0.9189385332046727417803297;
247 247
  int n=7;
248 248
  const double c[9]={0.99999999999980993,
......
265 265
    xdot-=c[i]*square(zinv);
266 266
  }
267 267
  double t=z+n+0.5;
268
  double tdot=zdot;
269 268
  //return lpp + (z+0.5)*log(t) -t + log(x);
270 269
  double ans= lpp + (z+0.5)*log(t) -t + log(x);
270
  //double tdot=zdot;
271 271
  //double ansdot=zdot*log(t) + (z+0.5)/t*tdot -tdot +xdot/x;
272 272
  // since tdot=1.0
273 273
  // since zdot=1.0

Also available in: Unified diff