diff --git a/src/k_rem_pio2.c b/src/k_rem_pio2.c index b0b5f63..e1224c3 100644 --- a/src/k_rem_pio2.c +++ b/src/k_rem_pio2.c @@ -311,7 +311,8 @@ __kernel_rem_pio2(double *x, double *y, int e0, int nx, int prec) /* 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; + for(j=0,fw=0.0;j<=jx;j++) fw += x[j]*f[jx+i-j]; + q[i] = fw; } jz = jk;