diff options
Diffstat (limited to 'math/test-math.c')
-rw-r--r-- | math/test-math.c | 150 |
1 files changed, 0 insertions, 150 deletions
diff --git a/math/test-math.c b/math/test-math.c deleted file mode 100644 index edd82a3..0000000 --- a/math/test-math.c +++ /dev/null @@ -1,150 +0,0 @@ -#include <stdio.h> -#include <math.h> -#include <stdlib.h> -#include <errno.h> -#include <string.h> - -void print_trig_stuff __P ((void)); - -int -main (int argc, char *argv[]) -{ - const char str[] = "123.456"; - double x,h,li,lr,a,lrr; - - if (signbit (argc < 1 ? -0.0 : 0.0)) - /* I don't make this an error for now. --drepper */ - fputs ("\n*** Your compiler has a bug. Consider upgrading\n\n", stderr); - - x = atof (str); - - printf ("%g %g\n", x, pow (10.0, 3.0)); - - x = sinh(2.0); - - printf("sinh(2.0) = %g\n", x); - - x = sinh(3.0); - - printf("sinh(3.0) = %g\n", x); - - h = hypot(2.0,3.0); - - printf("h=%g\n", h); - - a = atan2(3.0, 2.0); - - printf("atan2(3,2) = %g\n", a); - - lr = pow(h,4.0); - - printf("pow(%g,4.0) = %g\n", h, lr); - - lrr = lr; - - li = 4.0 * a; - - lr = lr / exp(a*5.0); - - printf("%g / exp(%g * 5) = %g\n", lrr, a, lr); - - lrr = li; - - li += 5.0 * log(h); - - printf("%g + 5*log(%g) = %g\n", lrr, h, li); - - printf("cos(%g) = %g, sin(%g) = %g\n", li, cos(li), li, sin(li)); - - x = drem(10.3435,6.2831852); - - printf("drem(10.3435,6.2831852) = %g\n", x); - - x = drem(-10.3435,6.2831852); - - printf("drem(-10.3435,6.2831852) = %g\n", x); - - x = drem(-10.3435,-6.2831852); - - printf("drem(-10.3435,-6.2831852) = %g\n", x); - - x = drem(10.3435,-6.2831852); - - printf("drem(10.3435,-6.2831852) = %g\n", x); - - - printf("x%8.6gx\n", .5); - printf("x%-8.6gx\n", .5); - printf("x%6.6gx\n", .5); - - { - double x = atof ("-1e-17-"); - printf ("%g %c= %g %s!\n", - x, - x == -1e-17 ? '=' : '!', - -1e-17, - x == -1e-17 ? "Worked" : "Failed"); - } - - print_trig_stuff (); - - return 0; -} - - -const double RAD[5] = { 0, M_PI/2, M_PI, (3*M_PI)/2, 2*M_PI }; - -#define PRINT_IT_1_ARG(_func, _arg, _value) \ - (_value) = (_func) ((_arg)); \ - if (errno) { \ - errno = 0; \ - printf ("%s = ERROR %s\n", #_func, strerror (errno)); \ - } else \ - printf ("%s(%g) = %g\n", #_func, _arg, (_value)); \ - -#define PRINT_IT_2_ARG(_func, _arg1, _arg2, _value) \ - (_value) = (_func) ((_arg1),(_arg2)); \ - if (errno) { \ - errno = 0; \ - printf ("%s = ERROR %s\n", #_func, strerror (errno)); \ - } else \ - printf ("%s(%g, %g) = %g\n", #_func, _arg1, _arg2, (_value)); \ - -void -print_trig_stuff(void) -{ - double value, arg1, arg2; - int i; - - puts ("\n\nMath Test"); - - errno = 0; /* automatically reset on error condition */ - for (i=0; i<4; i++) - { - PRINT_IT_1_ARG (sin, RAD[i], value); - PRINT_IT_1_ARG (cos, RAD[i], value); - PRINT_IT_1_ARG (tan, RAD[i], value); - PRINT_IT_1_ARG (asin, RAD[i], value); - PRINT_IT_1_ARG (acos, RAD[i], value); - PRINT_IT_1_ARG (atan, RAD[i], value); - PRINT_IT_2_ARG (atan2, RAD[i], -RAD[i % 4], value); - } - - arg1 = 16; - arg2 = 3; - PRINT_IT_1_ARG (exp, arg1, value); - PRINT_IT_1_ARG (log, arg1, value); - PRINT_IT_1_ARG (log10, arg1, value); - PRINT_IT_2_ARG (pow, arg1, arg2, value); - PRINT_IT_1_ARG (sqrt, arg1, value); - PRINT_IT_1_ARG (cbrt, arg1, value); - PRINT_IT_2_ARG (hypot, arg1, arg2, value); - PRINT_IT_1_ARG (expm1, arg1, value); - PRINT_IT_1_ARG (log1p, arg1, value); - PRINT_IT_1_ARG (sinh, arg1, value); - PRINT_IT_1_ARG (cosh, arg1, value); - PRINT_IT_1_ARG (tanh, arg1, value); - PRINT_IT_1_ARG (asinh, arg1, value); - PRINT_IT_1_ARG (acosh, arg1, value); - PRINT_IT_1_ARG (atanh, arg1, value); -} |