< prev index next >
src/java.base/share/native/libfdlibm/k_rem_pio2.c
Print this page
*** 196,206 ****
j = jv-jx; m = jx+jk;
for(i=0;i<=m;i++,j++) f[i] = (j<0)? zero : (double) ipio2[j];
/* compute q[0],q[1],...q[jk] */
for (i=0;i<=jk;i++) {
! for(j=0,fw=0.0;j<=jx;j++) fw += x[j]*f[jx+i-j]; q[i] = fw;
}
jz = jk;
recompute:
/* distill q[] into iq[] reversingly */
--- 196,208 ----
j = jv-jx; m = jx+jk;
for(i=0;i<=m;i++,j++) f[i] = (j<0)? zero : (double) ipio2[j];
/* compute q[0],q[1],...q[jk] */
for (i=0;i<=jk;i++) {
! for (j=0,fw=0.0;j<=jx;j++) {
! fw += x[j]*f[jx+i-j]; q[i] = fw;
! }
}
jz = jk;
recompute:
/* distill q[] into iq[] reversingly */
< prev index next >