------- Comment #3 from pinskia at gcc dot gnu dot org 2006-01-23 18:37 ------- Reduced testcase: double bessel_Kn_scaled_small_x(int n, double x) { int k; double y = 0.25 * x * x; double ex = exp(x); double k_term; double term1, sum1, ln_pre1; double term2, sum2, pre2; for(k=1; k<=n-1; k++) { k_term *= -y/(k * (n-k)); sum1 += k_term; } term1 = 0.5 * exp(ln_pre1) * sum1; return ex * (term1 + term2); }
-- pinskia at gcc dot gnu dot org changed: What |Removed |Added ---------------------------------------------------------------------------- Status|UNCONFIRMED |NEW Ever Confirmed|0 |1 Last reconfirmed|0000-00-00 00:00:00 |2006-01-23 18:37:52 date| | http://gcc.gnu.org/bugzilla/show_bug.cgi?id=25911