1+
12/* @(#)k_rem_pio2.c 1.3 95/01/18 */
23/*
34 * ====================================================
@@ -310,7 +311,7 @@ __kernel_rem_pio2(double *x, double *y, int e0, int nx, int prec)
310311
311312 /* compute q[0],q[1],...q[jk] */
312313 for (i = 0 ;i <=jk ;i ++ ) {
313- for (j = 0 ,fw = zero ;j <=jx ;j ++ ) fw += x [j ]* f [jx + i - j ]; q [i ] = fw ;
314+ for (j = 0 ,fw = 0.0 ;j <=jx ;j ++ ) fw += x [j ]* f [jx + i - j ]; q [i ] = fw ;
314315 }
315316
316317 jz = jk ;
@@ -369,7 +370,7 @@ __kernel_rem_pio2(double *x, double *y, int e0, int nx, int prec)
369370
370371 for (i = jz + 1 ;i <=jz + k ;i ++ ) { /* add q[jz+1] to q[jz+k] */
371372 f [jx + i ] = (double ) ipio2 [jv + i ];
372- for (j = 0 ,fw = zero ;j <=jx ;j ++ ) fw += x [j ]* f [jx + i - j ];
373+ for (j = 0 ,fw = 0.0 ;j <=jx ;j ++ ) fw += x [j ]* f [jx + i - j ];
373374 q [i ] = fw ;
374375 }
375376 jz += k ;
@@ -378,7 +379,7 @@ __kernel_rem_pio2(double *x, double *y, int e0, int nx, int prec)
378379 }
379380
380381 /* chop off zero terms */
381- if (z == zero ) {
382+ if (z == 0.0 ) {
382383 jz -= 1 ; q0 -= 24 ;
383384 while (iq [jz ]== 0 ) { jz -- ; q0 -= 24 ;}
384385 } else { /* break z into 24-bit if necessary */
@@ -399,20 +400,20 @@ __kernel_rem_pio2(double *x, double *y, int e0, int nx, int prec)
399400
400401 /* compute PIo2[0,...,jp]*q[jz,...,0] */
401402 for (i = jz ;i >=0 ;i -- ) {
402- for (fw = zero ,k = 0 ;k <=jp && k <=jz - i ;k ++ ) fw += PIo2 [k ]* q [i + k ];
403+ for (fw = 0.0 ,k = 0 ;k <=jp && k <=jz - i ;k ++ ) fw += PIo2 [k ]* q [i + k ];
403404 fq [jz - i ] = fw ;
404405 }
405406
406407 /* compress fq[] into y[] */
407408 switch (prec ) {
408409 case 0 :
409- fw = zero ;
410+ fw = 0.0 ;
410411 for (i = jz ;i >=0 ;i -- ) fw += fq [i ];
411412 y [0 ] = (ih == 0 )? fw : - fw ;
412413 break ;
413414 case 1 :
414415 case 2 :
415- fw = zero ;
416+ fw = 0.0 ;
416417 for (i = jz ;i >=0 ;i -- ) fw += fq [i ];
417418 STRICT_ASSIGN (double ,fw ,fw );
418419 y [0 ] = (ih == 0 )? fw : - fw ;
@@ -431,7 +432,7 @@ __kernel_rem_pio2(double *x, double *y, int e0, int nx, int prec)
431432 fq [i ] += fq [i - 1 ]- fw ;
432433 fq [i - 1 ] = fw ;
433434 }
434- for (fw = zero ,i = jz ;i >=2 ;i -- ) fw += fq [i ];
435+ for (fw = 0.0 ,i = jz ;i >=2 ;i -- ) fw += fq [i ];
435436 if (ih == 0 ) {
436437 y [0 ] = fq [0 ]; y [1 ] = fq [1 ]; y [2 ] = fw ;
437438 } else {
0 commit comments