Sat Apr 20 17:07:17 1996 Ulrich Drepper <drepper@cygnus.com>
[kopensolaris-gnu/glibc.git] / math / test-math.c
1 #include <ansidecl.h>
2 #include <stdio.h>
3 #include <math.h>
4 #include <stdlib.h>
5 #include <errno.h>
6 #include <string.h>
7
8 void print_trig_stuff __P ((void));
9
10 int
11 DEFUN_VOID(main)
12 {
13   CONST char str[] = "123.456";
14   double x,h,li,lr,a,lrr;
15
16   x = atof (str);
17
18   printf ("%g %g\n", x, pow (10.0, 3.0));
19
20   x = sinh(2.0);
21
22   printf("sinh(2.0) = %g\n", x);
23
24   x = sinh(3.0);
25
26   printf("sinh(3.0) = %g\n", x);
27
28   h = hypot(2.0,3.0);
29
30   printf("h=%g\n", h);
31
32   a = atan2(3.0, 2.0);
33
34   printf("atan2(3,2) = %g\n", a);
35
36   lr = pow(h,4.0);
37
38   printf("pow(%g,4.0) = %g\n", h, lr);
39
40   lrr = lr;
41
42   li = 4.0 * a;
43
44   lr = lr / exp(a*5.0);
45
46   printf("%g / exp(%g * 5) = %g\n", lrr, a, lr);
47
48   lrr = li;
49
50   li += 5.0 * log(h);
51
52   printf("%g + 5*log(%g) = %g\n", lrr, h, li);
53
54   printf("cos(%g) = %g,  sin(%g) = %g\n", li, cos(li), li, sin(li));
55
56   x = drem(10.3435,6.2831852);
57
58   printf("drem(10.3435,6.2831852) = %g\n", x);
59
60   x = drem(-10.3435,6.2831852);
61
62   printf("drem(-10.3435,6.2831852) = %g\n", x);
63
64   x = drem(-10.3435,-6.2831852);
65
66   printf("drem(-10.3435,-6.2831852) = %g\n", x);
67
68   x = drem(10.3435,-6.2831852);
69
70   printf("drem(10.3435,-6.2831852) = %g\n", x);
71
72
73   printf("x%8.6gx\n", .5);
74   printf("x%-8.6gx\n", .5);
75   printf("x%6.6gx\n", .5);
76
77   {
78     double x = atof ("-1e-17-");
79     printf ("%g %c= %g %s!\n",
80             x,
81             x == -1e-17 ? '=' : '!',
82             -1e-17,
83             x == -1e-17 ? "Worked" : "Failed");
84   }
85
86   print_trig_stuff ();
87
88   return 0;
89 }
90 \f
91
92 #define PI 3.14159265358979323846264338327
93
94 const double RAD[5] = { 0, PI/2, PI, (3*PI)/2, 2*PI };
95 const int    DEG[5] = { 0, 90, 180, 360 };
96
97 #define PRINT_IT_1_ARG(_func, _arg, _value) \
98     (_value) = (_func)((_arg)); \
99     if (errno) { \
100       errno = 0; \
101       printf("%s = ERROR %s\n", #_func, strerror(errno)); \
102     } else \
103       printf("%s(%g) = %g\n", #_func, _arg, (_value)); \
104
105 #define PRINT_IT_2_ARG(_func, _arg1, _arg2, _value) \
106     (_value) = (_func)((_arg1),(_arg2)); \
107     if (errno) { \
108       errno = 0; \
109       printf("%s = ERROR %s\n", #_func, strerror(errno)); \
110     } else \
111       printf("%s(%g, %g) = %g\n", #_func, _arg1, _arg2, (_value)); \
112
113 void
114 DEFUN_VOID (print_trig_stuff)
115 {
116   double value, arg1, arg2;
117   int i;
118
119   puts ("\n\nMath Test");
120
121   errno = 0;                    /* automatically reset on error condition */
122   for (i=0; i<4; i++)
123     {
124       PRINT_IT_1_ARG (sin, RAD[i], value);
125       PRINT_IT_1_ARG (cos, RAD[i], value);
126       PRINT_IT_1_ARG (tan, RAD[i], value);
127       PRINT_IT_1_ARG (asin, RAD[i], value);
128       PRINT_IT_1_ARG (acos, RAD[i], value);
129       PRINT_IT_1_ARG (atan, RAD[i], value);
130       PRINT_IT_2_ARG (atan2, RAD[i], -RAD[i % 4], value);
131     }
132
133   arg1 = 16;
134   arg2 = 3;
135   PRINT_IT_1_ARG (exp, arg1, value);
136   PRINT_IT_1_ARG (log, arg1, value);
137   PRINT_IT_1_ARG (log10, arg1, value);
138   PRINT_IT_2_ARG (pow, arg1, arg2, value);
139   PRINT_IT_1_ARG (sqrt, arg1, value);
140   PRINT_IT_1_ARG (cbrt, arg1, value);
141   PRINT_IT_2_ARG (hypot, arg1, arg2, value);
142   PRINT_IT_1_ARG (expm1, arg1, value);
143   PRINT_IT_1_ARG (log1p, arg1, value);
144   PRINT_IT_1_ARG (sinh, arg1, value);
145   PRINT_IT_1_ARG (cosh, arg1, value);
146   PRINT_IT_1_ARG (tanh, arg1, value);
147   PRINT_IT_1_ARG (asinh, arg1, value);
148   PRINT_IT_1_ARG (acosh, arg1, value);
149   PRINT_IT_1_ARG (atanh, arg1, value);
150 }