15  * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
21   These vectors were originally generated on a sun3 with a 68881 using
24   Each function is called with a variety of interesting arguments.
27   useful to test - eg sin(1e30) is pointless - the arg has to be
29   left to evaluate with - any number would be just as precise as any
74 				(f) = -__builtin_nansf("");	\
79 				(f) = -__builtin_nanf("");	\
89 				(d) = -__builtin_nans("");	\
94 				(d) = -__builtin_nan("");	\
117   thedouble(p->qs[id].msw, p->qs[id].lsw, &correct.value);  in ffcheck_id()
119   int error_bit = p->error_bit;  in ffcheck_id()
127    * On i386, just returning a value converts signaling nans to quiet  in ffcheck_id()
144 	   name,  p->line, mag,  in ffcheck_id()
153   if (p->qs[id].merror != merror)  in ffcheck_id()
157 	   name, p->line, merror, p->qs[id].merror);  in ffcheck_id()
160   if (p->qs[id].errno_val != errno)  in ffcheck_id()
163 	   name, p->line, errno, p->qs[id].errno_val);  in ffcheck_id()
200   thedouble(p->qs[id].msw,p->qs[id].lsw, &correct_double.value);  in fffcheck_id()
205   int error_bit = p->error_bit;  in fffcheck_id()
228 	   name,  p->line, mag,  in fffcheck_id()
238   if (p->qs[id].merror != merror)  in fffcheck_id()
242 	   name, p->line, merror, p->qs[id].merror);  in fffcheck_id()
245   if (p->qs[id].errno_val != errno)  in fffcheck_id()
248 	   name, p->line, errno, p->qs[id].errno_val);  in fffcheck_id()
270     while (len--)  in volatile_memcpy()
311   if (reduce && p->error_bit < mag)  in frontline()
313     fprintf(f, "{%2d,", p->error_bit);  in frontline()
330     translate_to(f, thedouble(p->qs[0].msw, p->qs[0].lsw, NULL));  in frontline()
335   fprintf(f,"0x%08lx, 0x%08lx", (unsigned long) p->qs[1].msw, (unsigned long) p->qs[1].lsw);  in frontline()
341     fprintf(f,"0x%08lx, 0x%08lx", (unsigned long) p->qs[2].msw, (unsigned long) p->qs[2].lsw);  in frontline()
345   	  thedouble(p->qs[1].msw, p->qs[1].lsw, NULL));  in frontline()
349     fprintf(f,", %g", thedouble(p->qs[2].msw,p->qs[2].lsw, NULL));  in frontline()
402     p->qs[2].msw = result2_double.parts.msw;  in finish2()
403     p->qs[2].lsw = result2_double.parts.lsw;  in finish2()
456     p->qs[2].msw = result2_double.parts.msw;  in ffinish2()
457     p->qs[2].lsw = result2_double.parts.lsw;  in ffinish2()
470 	return -(double) FLT_MAX <= arg && arg <= (double) FLT_MAX;  in in_float_range()
494 	struct { uint32_t a, b; };  in run_vector_1()  member
498       for (k = -.2; k < .2; k+= 0.00132)  in run_vector_1()
503 		(unsigned long) d.a, (unsigned long) d.b, (unsigned long) d4.a, (unsigned long) d4.b);  in run_vector_1()
507       for (k = -1.2; k < 1.2; k+= 0.01)  in run_vector_1()
512 		(unsigned long) d.a, (unsigned long) d.b, (unsigned long) d4.a, (unsigned long) d4.b);  in run_vector_1()
515       for (k = -M_PI *2; k < M_PI *2; k+= M_PI/2)  in run_vector_1()
520 		(unsigned long) d.a, (unsigned long) d.b, (unsigned long) d4.a, (unsigned long) d4.b);  in run_vector_1()
524       for (k = -30; k < 30; k+= 1.7)  in run_vector_1()
529 		(unsigned long) d.a, (unsigned long) d.b, (unsigned long) d4.a, (unsigned long) d4.b);  in run_vector_1()
539   while (p->line)  in run_vector_1()
542     thedouble(p->qs[1].msw, p->qs[1].lsw, &arg1);  in run_vector_1()
544     thedouble(p->qs[2].msw, p->qs[2].lsw, &arg2);  in run_vector_1()
551     line(p->line);  in run_vector_1()
560       /* Double function returning a double */  in run_vector_1()
572       /* Double function returning a double */  in run_vector_1()
675 #if TEST_PART == 0 || TEST_PART == -1  in test_math()
686 #if TEST_PART > 0 || TEST_PART == -1  in test_math()
690 #if TEST_PART == 0 || TEST_PART == -1  in test_math()
709 #if TEST_PART == 1 || TEST_PART == -1  in test_math()
736 #if TEST_PART == 2 || TEST_PART == -1  in test_math()
771 float cosf (float a) { return cos((double)a); }
772 float sinf (float  a) { return sin((double)a); }
773 float log1pf (float a) { return log1p((double)a); }
774 float tanf (float a) { return tan((double)a); }
775 float ceilf (float a) { return ceil(a); }
776 float floorf (float a) { return floor(a); }
782 float fmodf(a,b) float a,b; { return fmod(a,b); }
783 float hypotf(a,b) float a,b; { return hypot(a,b); }
785 float acosf(a) float a; { return acos(a); }
786 float acoshf(a) float a; { return acosh(a); }
787 float asinf(a) float a; { return asin(a); }
788 float asinhf(a) float a; { return asinh(a); }
789 float atanf(a) float a; { return atan(a); }
790 float atanhf(a) float a; { return atanh(a); }
792 float coshf(a) float a; { return cosh(a); }
793 float erff(a) float a; { return erf(a); }
794 float erfcf(a) float a; { return erfc(a); }
795 float expf(a) float a; { return exp(a); }
796 float fabsf(a) float a; { return fabs(a); }
798 float gammaf(a) float a; { return gamma(a); }
799 float j0f(a) float a; { return j0(a); }
800 float j1f(a) float a; { return j1(a); }
801 float log10f(a) float a; { return log10(a); }
803 float logf(a) float a; { return log(a); }
805 float sinhf(a) float a; { return sinh(a); }
806 float sqrtf(a) float a; { return sqrt(a); }
808 float tanhf(a) float a; { return tanh(a); }
809 float y0f(a) float a; { return y0(a); }
810 float y1f(a) float a; { return y1(a); }