Displaying 1 result from an estimated 1 matches for "gnumabs".
Did you mean:
gnueabi
2002 Oct 25
0
qgamma precision (PR#2214)
...gnum_float e0 = pgamma (x0, alpha, scale, lower_tail, log_p) - p;
if (e0 != 0 && lower_tail && !log_p) {
gnum_float d0 = dgamma (x0, alpha, scale, log_p);
if (d0) {
gnum_float x1 = x0 - e0 / d0;
gnum_float e1 = pgamma (x1, alpha, scale, lower_tail, log_p) - p;
if (gnumabs (e1) < gnumabs (e0))
x0 = x1;
}
}
return x0;
}
(Just read "gnum_float" as "double" and "gnumabs" as "fabs".)
The test for lower_tail is because (a) that is where I use it, and (b) the d0
calculation is probably wrong otherwise.
With th...