------- Comment #4 from pinskia at gcc dot gnu dot org 2006-01-23 18:42 ------- Further reduced: double bessel_Kn_scaled_small_x(int n) { int k; double k_term, sum1; for(k=1; k<=n-1; k++) { k_term *= -1/(k * (n-k)); sum1 += k_term; } return sum1; }
-- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=25911