htswanalysis/MACS/lib/gsl/gsl-1.11/specfunc/cheb_eval.c

   1
   2 static inline int
   3 cheb_eval_e(const cheb_series * cs,
   4             const double x,
   5             gsl_sf_result * result)
   6 {
   7   int j;
   8   double d  = 0.0;
   9   double dd = 0.0;
  10
  11   double y  = (2.0*x - cs->a - cs->b) / (cs->b - cs->a);
  12   double y2 = 2.0 * y;
  13
  14   double e = 0.0;
  15
  16   for(j = cs->order; j>=1; j--) {
  17     double temp = d;
  18     d = y2*d - dd + cs->c[j];
  19     e += fabs(y2*temp) + fabs(dd) + fabs(cs->c[j]);
  20     dd = temp;
  21   }
  22
  23   {
  24     double temp = d;
  25     d = y*d - dd + 0.5 * cs->c[0];
  26     e += fabs(y*temp) + fabs(dd) + 0.5 * fabs(cs->c[0]);
  27   }
  28
  29   result->val = d;
  30   result->err = GSL_DBL_EPSILON * e + fabs(cs->c[cs->order]);
  31
  32   return GSL_SUCCESS;
  33 }
  34