Lines Matching +full:0 +full:x0
82 static const unsigned thresh[2] = { 0x3fc90000, 0x3fc40000 };
89 * y[i*stridey] := cos( x[i*stridex] ), for i = 0..n.
102 double x0, x1, x2, *py0 = 0, *py1 = 0, *py2, *xsave, *ysave; in __vcos() local
103 unsigned hx0, hx1, hx2, xsb0, xsb1 = 0, xsb2; in __vcos()
111 biguns = 0; in __vcos()
113 x0 = *x; /* 'x0' may be used uninitialized */ in __vcos()
119 hx0 = xsb0 & ~0x80000000; /* mask off sign bit */ in __vcos()
120 if (hx0 > 0x3fe921fb) { in __vcos()
125 if (hx0 < 0x3e400000) { in __vcos()
131 i = 0; in __vcos()
132 if (--n <= 0) in __vcos()
136 x0 = *x; in __vcos()
141 if (--n <= 0) in __vcos()
146 hx1 = xsb1 & ~0x80000000; in __vcos()
147 if (hx1 > 0x3fe921fb) in __vcos()
152 if (hx1 < 0x3e400000) in __vcos()
159 if (--n <= 0) in __vcos()
168 if (--n <= 0) in __vcos()
173 hx2 = xsb2 & ~0x80000000; in __vcos()
174 if (hx2 > 0x3fe921fb) in __vcos()
179 if (hx2 < 0x3e400000) in __vcos()
186 if (--n <= 0) in __vcos()
194 * 0x3fc40000 = 5/32 ~ 0.15625 in __vcos()
198 i = (hx0 - 0x3fc40000) >> 31; in __vcos()
199 i |= ((hx1 - 0x3fc40000) >> 30) & 2; in __vcos()
200 i |= ((hx2 - 0x3fc40000) >> 29) & 4; in __vcos()
207 case 0: /* All are > 5/32 */ in __vcos()
208 j0 = (xsb0 + 0x4000) & 0xffff8000; in __vcos()
209 j1 = (xsb1 + 0x4000) & 0xffff8000; in __vcos()
210 j2 = (xsb2 + 0x4000) & 0xffff8000; in __vcos()
214 LO(&t0) = 0; in __vcos()
215 LO(&t1) = 0; in __vcos()
216 LO(&t2) = 0; in __vcos()
217 x0 -= t0; in __vcos()
220 z0 = x0 * x0; in __vcos()
226 w0 = x0 * (one + z0 * (pp1 + z0 * pp2)); in __vcos()
229 j0 = (((j0 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
230 j1 = (((j1 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
231 j2 = (((j2 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
249 j1 = (xsb1 + 0x4000) & 0xffff8000; in __vcos()
250 j2 = (xsb2 + 0x4000) & 0xffff8000; in __vcos()
253 LO(&t1) = 0; in __vcos()
254 LO(&t2) = 0; in __vcos()
257 z0 = x0 * x0; in __vcos()
266 j1 = (((j1 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
267 j2 = (((j2 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
280 j0 = (xsb0 + 0x4000) & 0xffff8000; in __vcos()
281 j2 = (xsb2 + 0x4000) & 0xffff8000; in __vcos()
284 LO(&t0) = 0; in __vcos()
285 LO(&t2) = 0; in __vcos()
286 x0 -= t0; in __vcos()
288 z0 = x0 * x0; in __vcos()
294 w0 = x0 * (one + z0 * (pp1 + z0 * pp2)); in __vcos()
297 j0 = (((j0 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
298 j2 = (((j2 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
311 j2 = (xsb2 + 0x4000) & 0xffff8000; in __vcos()
313 LO(&t2) = 0; in __vcos()
315 z0 = x0 * x0; in __vcos()
324 j2 = (((j2 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
334 j0 = (xsb0 + 0x4000) & 0xffff8000; in __vcos()
335 j1 = (xsb1 + 0x4000) & 0xffff8000; in __vcos()
338 LO(&t0) = 0; in __vcos()
339 LO(&t1) = 0; in __vcos()
340 x0 -= t0; in __vcos()
342 z0 = x0 * x0; in __vcos()
348 w0 = x0 * (one + z0 * (pp1 + z0 * pp2)); in __vcos()
351 j0 = (((j0 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
352 j1 = (((j1 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
365 j1 = (xsb1 + 0x4000) & 0xffff8000; in __vcos()
367 LO(&t1) = 0; in __vcos()
369 z0 = x0 * x0; in __vcos()
378 j1 = (((j1 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
388 j0 = (xsb0 + 0x4000) & 0xffff8000; in __vcos()
390 LO(&t0) = 0; in __vcos()
391 x0 -= t0; in __vcos()
392 z0 = x0 * x0; in __vcos()
398 w0 = x0 * (one + z0 * (pp1 + z0 * pp2)); in __vcos()
401 j0 = (((j0 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
411 z0 = x0 * x0; in __vcos()
428 i = 0; in __vcos()
429 } while (--n > 0); /* END MAIN LOOP */ in __vcos()
432 * CLEAN UP last 0, 1, or 2 elts. in __vcos()
434 if (i > 0) /* Clean up elts at tail. i < 3. */ in __vcos()
442 if (hx1 < 0x3fc40000) in __vcos()
452 j1 = (xsb1 + 0x4000) & 0xffff8000; in __vcos()
454 LO(&t1) = 0; in __vcos()
459 j1 = (((j1 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
467 if (hx0 < 0x3fc40000) in __vcos()
469 z0 = x0 * x0; in __vcos()
477 j0 = (xsb0 + 0x4000) & 0xffff8000; in __vcos()
479 LO(&t0) = 0; in __vcos()
480 x0 -= t0; in __vcos()
481 z0 = x0 * x0; in __vcos()
483 w0 = x0 * (one + z0 * (pp1 + z0 * pp2)); in __vcos()
484 j0 = (((j0 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
519 biguns = 0; in __vcos()
527 biguns = 0; in __vcos()
530 biguns = 0; in __vcos()
545 hx &= ~0x80000000; in __vcos()
546 if (hx > 0x413921fb) /* (1.6471e+06) Too big: leave it. */ in __vcos()
548 if (hx >= 0x7ff00000) /* Inf or NaN */ in __vcos()
550 x0 = *x; in __vcos()
551 *y = x0 - x0; in __vcos()
557 i = 0; in __vcos()
558 if (--n <= 0) in __vcos()
562 x0 = *x; in __vcos()
567 if (--n <= 0) in __vcos()
573 hx &= ~0x80000000; in __vcos()
574 if (hx > 0x413921fb) in __vcos()
576 if (hx >= 0x7ff00000) in __vcos()
586 if (--n <= 0) in __vcos()
595 if (--n <= 0) in __vcos()
601 hx &= ~0x80000000; in __vcos()
602 if (hx > 0x413921fb) in __vcos()
604 if (hx >= 0x7ff00000) in __vcos()
614 if (--n <= 0) in __vcos()
621 n0 = (int) (x0 * invpio2 + half[xsb0]); in __vcos()
630 a0 = x0 - fn0 * pio2_1; in __vcos()
636 x0 = a0 - w0; in __vcos()
639 y0 = (a0 - x0) - w0; in __vcos()
642 a0 = x0; in __vcos()
648 x0 = a0 - w0; in __vcos()
651 y0 = (a0 - x0) - w0; in __vcos()
654 a0 = x0; in __vcos()
660 x0 = a0 - w0; in __vcos()
663 y0 = (a0 - x0) - w0; in __vcos()
666 xsb0 = HI(&x0); in __vcos()
667 i = ((xsb0 & ~0x80000000) - thresh[n0&1]) >> 31; in __vcos()
669 i |= (((xsb1 & ~0x80000000) - thresh[n1&1]) >> 30) & 2; in __vcos()
671 i |= (((xsb2 & ~0x80000000) - thresh[n2&1]) >> 29) & 4; in __vcos()
677 case 0: in __vcos()
678 j0 = (xsb0 + 0x4000) & 0xffff8000; in __vcos()
679 j1 = (xsb1 + 0x4000) & 0xffff8000; in __vcos()
680 j2 = (xsb2 + 0x4000) & 0xffff8000; in __vcos()
684 LO(&t0) = 0; in __vcos()
685 LO(&t1) = 0; in __vcos()
686 LO(&t2) = 0; in __vcos()
687 x0 = (x0 - t0) + y0; in __vcos()
690 z0 = x0 * x0; in __vcos()
696 w0 = x0 * (one + z0 * (pp1 + z0 * pp2)); in __vcos()
699 j0 = (((j0 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
700 j1 = (((j1 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
701 j2 = (((j2 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
724 j1 = (xsb1 + 0x4000) & 0xffff8000; in __vcos()
725 j2 = (xsb2 + 0x4000) & 0xffff8000; in __vcos()
728 LO(&t1) = 0; in __vcos()
729 LO(&t2) = 0; in __vcos()
730 x0_or_one[0] = x0; in __vcos()
731 x0_or_one[2] = -x0; in __vcos()
732 y0_or_zero[0] = y0; in __vcos()
736 z0 = x0 * x0; in __vcos()
745 j1 = (((j1 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
746 j2 = (((j2 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
764 j0 = (xsb0 + 0x4000) & 0xffff8000; in __vcos()
766 j2 = (xsb2 + 0x4000) & 0xffff8000; in __vcos()
769 LO(&t0) = 0; in __vcos()
770 LO(&t2) = 0; in __vcos()
771 x1_or_one[0] = x1; in __vcos()
773 x0 = (x0 - t0) + y0; in __vcos()
774 y1_or_zero[0] = y1; in __vcos()
777 z0 = x0 * x0; in __vcos()
783 w0 = x0 * (one + z0 * (pp1 + z0 * pp2)); in __vcos()
786 j0 = (((j0 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
787 j2 = (((j2 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
807 j2 = (xsb2 + 0x4000) & 0xffff8000; in __vcos()
809 LO(&t2) = 0; in __vcos()
810 x0_or_one[0] = x0; in __vcos()
811 x0_or_one[2] = -x0; in __vcos()
812 x1_or_one[0] = x1; in __vcos()
814 y0_or_zero[0] = y0; in __vcos()
816 y1_or_zero[0] = y1; in __vcos()
819 z0 = x0 * x0; in __vcos()
828 j2 = (((j2 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
842 j0 = (xsb0 + 0x4000) & 0xffff8000; in __vcos()
843 j1 = (xsb1 + 0x4000) & 0xffff8000; in __vcos()
847 LO(&t0) = 0; in __vcos()
848 LO(&t1) = 0; in __vcos()
849 x2_or_one[0] = x2; in __vcos()
851 x0 = (x0 - t0) + y0; in __vcos()
853 y2_or_zero[0] = y2; in __vcos()
855 z0 = x0 * x0; in __vcos()
861 w0 = x0 * (one + z0 * (pp1 + z0 * pp2)); in __vcos()
864 j0 = (((j0 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
865 j1 = (((j1 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
884 j1 = (xsb1 + 0x4000) & 0xffff8000; in __vcos()
887 LO(&t1) = 0; in __vcos()
888 x0_or_one[0] = x0; in __vcos()
889 x0_or_one[2] = -x0; in __vcos()
890 x2_or_one[0] = x2; in __vcos()
892 y0_or_zero[0] = y0; in __vcos()
895 y2_or_zero[0] = y2; in __vcos()
897 z0 = x0 * x0; in __vcos()
906 j1 = (((j1 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
920 j0 = (xsb0 + 0x4000) & 0xffff8000; in __vcos()
924 LO(&t0) = 0; in __vcos()
925 x1_or_one[0] = x1; in __vcos()
927 x2_or_one[0] = x2; in __vcos()
929 x0 = (x0 - t0) + y0; in __vcos()
930 y1_or_zero[0] = y1; in __vcos()
932 y2_or_zero[0] = y2; in __vcos()
934 z0 = x0 * x0; in __vcos()
940 w0 = x0 * (one + z0 * (pp1 + z0 * pp2)); in __vcos()
943 j0 = (((j0 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
960 x0_or_one[0] = x0; in __vcos()
961 x0_or_one[2] = -x0; in __vcos()
962 x1_or_one[0] = x1; in __vcos()
964 x2_or_one[0] = x2; in __vcos()
966 y0_or_zero[0] = y0; in __vcos()
968 y1_or_zero[0] = y1; in __vcos()
970 y2_or_zero[0] = y2; in __vcos()
972 z0 = x0 * x0; in __vcos()
992 i = 0; in __vcos()
993 } while (--n > 0); in __vcos()
995 if (i > 0) in __vcos()
1020 if ((xsb1 & ~0x80000000) < thresh[n1&1]) in __vcos()
1023 x1_or_one[0] = x1; in __vcos()
1025 y1_or_zero[0] = y1; in __vcos()
1035 j1 = (xsb1 + 0x4000) & 0xffff8000; in __vcos()
1037 LO(&t1) = 0; in __vcos()
1042 j1 = (((j1 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
1051 n0 = (int) (x0 * invpio2 + half[xsb0]); in __vcos()
1054 a0 = x0 - fn0 * pio2_1; in __vcos()
1056 x0 = a0 - w0; in __vcos()
1057 y0 = (a0 - x0) - w0; in __vcos()
1058 a0 = x0; in __vcos()
1060 x0 = a0 - w0; in __vcos()
1061 y0 = (a0 - x0) - w0; in __vcos()
1062 a0 = x0; in __vcos()
1064 x0 = a0 - w0; in __vcos()
1065 y0 = (a0 - x0) - w0; in __vcos()
1066 xsb0 = HI(&x0); in __vcos()
1067 if ((xsb0 & ~0x80000000) < thresh[n0&1]) in __vcos()
1070 x0_or_one[0] = x0; in __vcos()
1071 x0_or_one[2] = -x0; in __vcos()
1072 y0_or_zero[0] = y0; in __vcos()
1074 z0 = x0 * x0; in __vcos()
1082 j0 = (xsb0 + 0x4000) & 0xffff8000; in __vcos()
1084 LO(&t0) = 0; in __vcos()
1085 x0 = (x0 - t0) + y0; in __vcos()
1086 z0 = x0 * x0; in __vcos()
1088 w0 = x0 * (one + z0 * (pp1 + z0 * pp2)); in __vcos()
1089 j0 = (((j0 & ~0x80000000) - 0x3fc40000) >> 13) & ~0x3; in __vcos()
1100 __vlibm_vcos_big(nsave, xsave, sxsave, ysave, sysave, 0x413921fb); in __vcos()