math/test-math.c

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