< 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 >