none/tests/x86/x87trigOOR.c

   1
   2 /* Tests out of range handling for FSIN, FCOS, FSINCOS and FPTAN.  Be
   3    careful with the inline assembly -- this program is compiled as
   4    both a 32-bit and 64-bit test. */
   5
   6 #include <stdio.h>
   7 #include <string.h>
   8 #include <assert.h>
   9
  10 typedef  unsigned short int      UShort;
  11 typedef  unsigned int            UInt;
  12 typedef  double                  Double;
  13 typedef  unsigned long long int  ULong;
  14
  15 typedef  struct { Double arg; Double st0; Double st1; UShort fpusw; }  Res;
  16
  17 #define SHIFT_C3   14
  18 #define SHIFT_C2   10
  19 #define SHIFT_C1   9
  20 #define SHIFT_C0   8
  21
  22
  23 #define my_offsetof(type,memb) ((int)(unsigned long int)&((type*)0)->memb)
  24
  25 void do_fsin ( /*OUT*/Res* r, double d )
  26 {
  27    assert(my_offsetof(Res,arg) == 0);
  28    assert(my_offsetof(Res,st0) == 8);
  29    assert(my_offsetof(Res,st1) == 16);
  30    assert(my_offsetof(Res,fpusw) == 24);
  31    memset(r, 0, sizeof(*r));
  32    r->arg = d;
  33    __asm__ __volatile__(
  34      "finit"              "\n\t"
  35      "fldpi"              "\n\t"
  36      "fldl 0(%0)"         "\n\t" // .arg
  37      "fsin"               "\n\t"
  38      "fstsw %%ax"         "\n\t"
  39      "fstpl 8(%0)"        "\n\t" // .st0
  40      "fstpl 16(%0)"       "\n\t" // .st1
  41      "movw %%ax, 24(%0)"  "\n\t" // .fpusw
  42      "finit"              "\n"
  43      : : "r"(r) : "eax","cc","memory"
  44    );
  45 }
  46
  47 void do_fcos ( /*OUT*/Res* r, double d )
  48 {
  49    assert(my_offsetof(Res,arg) == 0);
  50    assert(my_offsetof(Res,st0) == 8);
  51    assert(my_offsetof(Res,st1) == 16);
  52    assert(my_offsetof(Res,fpusw) == 24);
  53    memset(r, 0, sizeof(*r));
  54    r->arg = d;
  55    __asm__ __volatile__(
  56      "finit"              "\n\t"
  57      "fldpi"              "\n\t"
  58      "fldl 0(%0)"         "\n\t" // .arg
  59      "fcos"               "\n\t"
  60      "fstsw %%ax"         "\n\t"
  61      "fstpl 8(%0)"        "\n\t" // .st0
  62      "fstpl 16(%0)"       "\n\t" // .st1
  63      "movw %%ax, 24(%0)"  "\n\t" // .fpusw
  64      "finit"              "\n"
  65      : : "r"(r) : "eax","cc","memory"
  66    );
  67 }
  68
  69 void do_fsincos ( /*OUT*/Res* r, double d )
  70 {
  71    assert(my_offsetof(Res,arg) == 0);
  72    assert(my_offsetof(Res,st0) == 8);
  73    assert(my_offsetof(Res,st1) == 16);
  74    assert(my_offsetof(Res,fpusw) == 24);
  75    memset(r, 0, sizeof(*r));
  76    r->arg = d;
  77    __asm__ __volatile__(
  78      "finit"              "\n\t"
  79      "fldpi"              "\n\t"
  80      "fldl 0(%0)"         "\n\t" // .arg
  81      "fsincos"            "\n\t"
  82      "fstsw %%ax"         "\n\t"
  83      "fstpl 8(%0)"        "\n\t" // .st0
  84      "fstpl 16(%0)"       "\n\t" // .st1
  85      "movw %%ax, 24(%0)"  "\n\t" // .fpusw
  86      "finit"              "\n"
  87      : : "r"(r) : "eax","cc","memory"
  88    );
  89 }
  90
  91 void do_fptan ( /*OUT*/Res* r, double d )
  92 {
  93    assert(my_offsetof(Res,arg) == 0);
  94    assert(my_offsetof(Res,st0) == 8);
  95    assert(my_offsetof(Res,st1) == 16);
  96    assert(my_offsetof(Res,fpusw) == 24);
  97    memset(r, 0, sizeof(*r));
  98    r->arg = d;
  99    __asm__ __volatile__(
 100      "finit"              "\n\t"
 101      "fldpi"              "\n\t"
 102      "fldl 0(%0)"         "\n\t" // .arg
 103      "fptan"              "\n\t"
 104      "fstsw %%ax"         "\n\t"
 105      "fstpl 8(%0)"        "\n\t" // .st0
 106      "fstpl 16(%0)"       "\n\t" // .st1
 107      "movw %%ax, 24(%0)"  "\n\t" // .fpusw
 108      "finit"              "\n"
 109      : : "r"(r) : "eax","cc","memory"
 110    );
 111 }
 112
 113
 114 void try ( char* name, void(*fn)(Res*,double), double d )
 115 {
 116    Res r;
 117    fn(&r, d);
 118    // Mask out all except C2 (range)
 119    r.fpusw &= (1 << SHIFT_C2);
 120    printf("%s  %16e --> %16e %16e %04x\n",
 121           name, r.arg, r.st0, r.st1, (UInt)r.fpusw);
 122 }
 123
 124 int main ( void )
 125 {
 126    Double limit = 9223372036854775808.0; // 2^63
 127
 128    char* names[4] = { "fsin   ", "fcos   ", "fsincos", "fptan  " };
 129    void(*fns[4])(Res*,double) = { do_fsin, do_fcos, do_fsincos, do_fptan };
 130
 131    int i;
 132    for (i = 0; i < 4; i++) {
 133       char* name = names[i];
 134       void (*fn)(Res*,double) = fns[i];
 135
 136       try( name, fn,   0.0   );
 137       try( name, fn,   0.123 );
 138       try( name, fn,  -0.456 );
 139       try( name, fn,  37.0   );
 140       try( name, fn, -53.0   );
 141       printf("\n");
 142
 143       try( name, fn, limit * 0.900000 );
 144       try( name, fn, limit * 0.999999 );
 145       try( name, fn, limit * 1.000000 );
 146       try( name, fn, limit * 1.000001 );
 147       try( name, fn, limit * 1.100000 );
 148       printf("\n");
 149
 150       try( name, fn, -limit * 0.900000 );
 151       try( name, fn, -limit * 0.999999 );
 152       try( name, fn, -limit * 1.000000 );
 153       try( name, fn, -limit * 1.000001 );
 154       try( name, fn, -limit * 1.100000 );
 155       printf("\n");
 156    }
 157
 158    return 0;
 159 }