memcheck/tests/vcpu_fnfns.c

   1
   2 /* Program to check that the FP stuff underlying these common FP
   3    functions isn't too badly broken.  Carefully kludged to print the
   4    same answers on different platforms (even when run natively). */
   5
   6 #include <stdio.h>
   7 #include <math.h>
   8
   9 int main ( void )
  10 {
  11    double d;
  12    float f;
  13    int i;
  14
  15    const double tinyD = 0.0000000001;
  16    const double tinyF = 0.0001;
  17
  18    /* -------------------- any arg -------------------- */
  19
  20    d = -2.0;
  21    for (i = 0; i < 41; i++) {
  22       printf("floorD(%+20.13e) = %+20.13e\n", d, floor(d));
  23       d += 0.1-tinyD;
  24    }
  25    f = -2.0;
  26    for (i = 0; i < 41; i++) {
  27       printf("floorF(%+20.4e) = %+20.4e\n", (double)f, (double)floorf(f));
  28       f += 0.1-tinyF;
  29    }
  30
  31
  32    d = -2.0;
  33    for (i = 0; i < 41; i++) {
  34       printf(" ceilD(%+20.13e) = %+20.13e\n", d, ceil(d));
  35       d += 0.1-tinyD;
  36    }
  37    f = -2.0;
  38    for (i = 0; i < 41; i++) {
  39       printf(" ceilF(%+20.4e) = %+20.4e\n", (double)f, (double)ceilf(f));
  40       f += 0.1-tinyF;
  41    }
  42
  43
  44    d = -2.0;
  45    for (i = 0; i < 41; i++) {
  46       printf("  sinD(%+20.13e) = %+20.13e\n", d, sin(d));
  47       d += 0.1-tinyD;
  48    }
  49    f = -2.0;
  50    for (i = 0; i < 41; i++) {
  51       printf("  sinF(%+20.4e) = %+20.4e\n", (double)f, (double)sinf(f));
  52       f += 0.1-tinyF;
  53    }
  54
  55
  56    d = -2.0;
  57    for (i = 0; i < 41; i++) {
  58       printf("  cosD(%+20.13e) = %+20.13e\n", d, cos(d));
  59       d += 0.1-tinyD;
  60    }
  61    f = -2.0;
  62    for (i = 0; i < 41; i++) {
  63       printf("  cosF(%+20.4e) = %+20.4e\n", (double)f, (double)cosf(f));
  64       f += 0.1-tinyF;
  65    }
  66
  67
  68    d = -2.0;
  69    for (i = 0; i < 41; i++) {
  70       printf("  tanD(%+20.13e) = %+20.13e\n", d, tan(d));
  71       d += 0.1-tinyD;
  72    }
  73    f = -2.0;
  74    for (i = 0; i < 41; i++) {
  75       printf("  tanF(%+20.4e) = %+20.4e\n", (double)f, (double)tanf(f));
  76       f += 0.1-tinyF;
  77    }
  78
  79
  80    d = -2.0;
  81    for (i = 0; i < 41; i++) {
  82       printf("  expD(%+20.13e) = %+20.13e\n", d, exp(d));
  83       d += 0.1-tinyD;
  84    }
  85    f = -2.0;
  86    for (i = 0; i < 41; i++) {
  87       printf("  expF(%+20.4e) = %+20.4e\n", (double)f, (double)expf(f));
  88       f += 0.1-tinyF;
  89    }
  90
  91    /* -------------------- >= 0 arg -------------------- */
  92
  93    d = 0.0;
  94    for (i = 0; i < 21; i++) {
  95       printf(" sqrtD(%+20.13e) = %+20.13e\n", d, sqrt(d));
  96       d += 0.1-tinyD;
  97    }
  98    f = 0.0;
  99    for (i = 0; i < 21; i++) {
 100       printf(" sqrtF(%+20.4e) = %+20.4e\n", (double)f, (double)sqrtf(f));
 101       f += 0.1-tinyF;
 102    }
 103
 104
 105    d = 0.0;
 106    for (i = 0; i < 21; i++) {
 107       printf("  logD(%+20.13e) = %+20.13e\n", d, log(d));
 108       d += 0.1-tinyD;
 109    }
 110    f = 0.0;
 111    for (i = 0; i < 21; i++) {
 112       printf("  logF(%+20.4e) = %+20.4e\n", (double)f, (double)logf(f));
 113       f += 0.1-tinyF;
 114    }
 115
 116
 117    d = 0.0;
 118    for (i = 0; i < 21; i++) {
 119       printf("log10D(%+20.13e) = %+20.13e\n", d, log10(d));
 120       d += 0.1-tinyD;
 121    }
 122    f = 0.0;
 123    for (i = 0; i < 21; i++) {
 124       printf("log10F(%+20.4e) = %+20.4e\n", (double)f, (double)log10f(f));
 125       f += 0.1-tinyF;
 126    }
 127
 128    /* -------------------- -1 .. +1 arg -------------------- */
 129
 130    d = -1.0;
 131    for (i = 0; i < 21; i++) {
 132       printf(" asinD(%+20.13e) = %+20.13e\n", d, asin(d));
 133       d += 0.1-tinyD;
 134    }
 135    f = -1.0;
 136    for (i = 0; i < 21; i++) {
 137       printf(" asinF(%+20.4e) = %+20.4e\n", (double)f, (double)asinf(f));
 138       f += 0.1-tinyF;
 139    }
 140
 141    /* acos(double) seems very prone to accuracy loss near the end of
 142       the range (arg --> +1.0).  Hence is different from the rest to
 143       stop it getting so close to 1.0. */
 144    d = -1.0;
 145    for (i = 0; i < 21; i++) {
 146       printf(" acosD(%+20.13e) = %+20.10e\n", d, acos(d));
 147       d += 0.1 - 1000.0*tinyD;
 148    }
 149    f = -1.0;
 150    for (i = 0; i < 21; i++) {
 151       printf(" acosF(%+20.4e) = %+20.4e\n", (double)f, (double)acosf(f));
 152       f += 0.1-tinyF;
 153    }
 154
 155
 156    d = -1.0;
 157    for (i = 0; i < 21; i++) {
 158       printf(" atanD(%+20.13e) = %+20.13e\n", d, atan(d));
 159       d += 0.1-tinyD;
 160    }
 161    f = -1.0;
 162    for (i = 0; i < 21; i++) {
 163       printf(" atanF(%+20.4e) = %+20.4e\n", (double)f, (double)atanf(f));
 164       f += 0.1-tinyF;
 165    }
 166
 167
 168    d = -1.0;
 169    for (i = 0; i < 21; i++) {
 170       printf("atan2D(%+20.13e) = %+20.13e\n", d, atan2(d, 1.0));
 171       d += 0.1-tinyD;
 172    }
 173    f = -1.0;
 174    for (i = 0; i < 21; i++) {
 175       printf("atan2F(%+20.4e) = %+20.4e\n", (double)f, (double)atan2f(f,1.0));
 176       f += 0.1-tinyF;
 177    }
 178
 179    return 0;
 180 }