l4/pkg/valgrind/src/valgrind-3.6.0-svn/memcheck/tests/x86/more_x86_fp.c

   1
   2 /* Derived from: */
   3
   4 /*
   5  *  x86 CPU test
   6  *
   7  *  Copyright (c) 2003 Fabrice Bellard
   8  *
   9  *  This program is free software; you can redistribute it and/or modify
  10  *  it under the terms of the GNU General Public License as published by
  11  *  the Free Software Foundation; either version 2 of the License, or
  12  *  (at your option) any later version.
  13  *
  14  *  This program is distributed in the hope that it will be useful,
  15  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
  16  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  17  *  GNU General Public License for more details.
  18  *
  19  *  You should have received a copy of the GNU General Public License
  20  *  along with this program; if not, write to the Free Software
  21  *  Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  22  */
  23
  24
  25 #include <stdlib.h>
  26 #include <stdio.h>
  27 #include <string.h>
  28 #include <inttypes.h>
  29 #include <math.h>
  30
  31 /**********************************************/
  32
  33 void test_fops(double a, double b)
  34 {
  35     printf("a=%f b=%f a+b=%f\n", a, b, a + b);
  36     printf("a=%f b=%f a-b=%f\n", a, b, a - b);
  37     printf("a=%f b=%f a*b=%f\n", a, b, a * b);
  38     printf("a=%f b=%f a/b=%f\n", a, b, a / b);
  39     printf("a=%f b=%f fmod(a, b)=%f\n", a, b, fmod(a, b));
  40     printf("a=%f sqrt(a)=%f\n", a, sqrt(a));
  41     printf("a=%f sin(a)=%f\n", a, sin(a));
  42     printf("a=%f cos(a)=%f\n", a, cos(a));
  43     printf("a=%f tan(a)=%f\n", a, tan(a));
  44     printf("a=%f log(a)=%f\n", a, log(a));
  45     printf("a=%f exp(a)=%f\n", a, exp(a));
  46     printf("a=%f b=%f atan2(a, b)=%f\n", a, b, atan2(a, b));
  47     /* just to test some op combining */
  48     printf("a=%f asin(sin(a))=%f\n", a, asin(sin(a)));
  49     printf("a=%f acos(cos(a))=%f\n", a, acos(cos(a)));
  50     printf("a=%f atan(tan(a))=%f\n", a, atan(tan(a)));
  51 }
  52 #define CC_C    0x0001
  53 #define CC_P    0x0004
  54 #define CC_A    0x0010
  55 #define CC_Z    0x0040
  56 #define CC_S    0x0080
  57 #define CC_O    0x0800
  58
  59
  60 void test_fcmp(double a, double b)
  61 {
  62     printf("(%f<%f)=%d\n",
  63            a, b, a < b);
  64     printf("(%f<=%f)=%d\n",
  65            a, b, a <= b);
  66     printf("(%f==%f)=%d\n",
  67            a, b, a == b);
  68     printf("(%f>%f)=%d\n",
  69            a, b, a > b);
  70     printf("(%f<=%f)=%d\n",
  71            a, b, a >= b);
  72     {
  73         unsigned int eflags;
  74         /* test f(u)comi instruction */
  75         asm("fcomi %2, %1\n"
  76             "pushf\n"
  77             "pop %0\n"
  78             : "=r" (eflags)
  79             : "t" (a), "u" (b));
  80         printf("fcomi(%f %f)=%08x\n", a, b, eflags & (CC_Z | CC_P | CC_C));
  81     }
  82 }
  83
  84 void test_fcvt(double a)
  85 {
  86     float fa;
  87     long double la;
  88     int16_t fpuc;
  89     int i;
  90     int64_t lla;
  91     int ia;
  92     int16_t wa;
  93     double ra;
  94
  95     fa = a;
  96     la = a;
  97     printf("(float)%f = %f\n", a, fa);
  98     printf("(long double)%f = %Lf\n", a, la);
  99     printf("a=%016llx\n", *(long long *)&a);
 100     printf("la=%016llx %04x\n", *(long long *)&la,
 101            *(unsigned short *)((char *)(&la) + 8));
 102
 103     /* test all roundings */
 104     asm volatile ("fstcw %0" : "=m" (fpuc));
 105     for(i=0;i<4;i++) {
 106       int16_t tmp = (fpuc & ~0x0c00) | (i << 10);
 107         asm volatile ("fldcw %0" : : "m" (tmp));
 108         asm volatile ("fist %0" : "=m" (wa) : "t" (a));
 109         asm volatile ("fistl %0" : "=m" (ia) : "t" (a));
 110         asm volatile ("fistpll %0" : "=m" (lla) : "t" (a) : "st");
 111         asm volatile ("frndint ; fstl %0" : "=m" (ra) : "t" (a));
 112         asm volatile ("fldcw %0" : : "m" (fpuc));
 113         printf("(short)a = %d\n", wa);
 114         printf("(int)a = %d\n", ia);
 115         printf("(int64_t)a = %lld\n", lla);
 116         printf("rint(a) = %f\n", ra);
 117     }
 118 }
 119
 120 #define TEST(N) \
 121     asm("fld" #N : "=t" (a)); \
 122     printf("fld" #N "= %f\n", a);
 123
 124 void test_fconst(void)
 125 {
 126     double a;
 127     TEST(1);
 128     TEST(l2t);
 129     TEST(l2e);
 130     TEST(pi);
 131     TEST(lg2);
 132     TEST(ln2);
 133     TEST(z);
 134 }
 135
 136 void test_fbcd(double a)
 137 {
 138     unsigned short bcd[5];
 139     double b;
 140
 141     asm("fbstp %0" : "=m" (bcd[0]) : "t" (a) : "st");
 142     asm("fbld %1" : "=t" (b) : "m" (bcd[0]));
 143     printf("a=%f bcd=%04x%04x%04x%04x%04x b=%f\n",
 144            a, bcd[4], bcd[3], bcd[2], bcd[1], bcd[0], b);
 145 }
 146
 147 #define TEST_ENV(env, save, restore)\
 148 {\
 149     memset((env), 0xaa, sizeof(*(env)));\
 150     for(i=0;i<5;i++)\
 151         asm volatile ("fldl %0" : : "m" (dtab[i]));\
 152     asm(save " %0\n" : : "m" (*(env)));\
 153     asm(restore " %0\n": : "m" (*(env)));\
 154     for(i=0;i<5;i++)\
 155         asm volatile ("fstpl %0" : "=m" (rtab[i]));\
 156     for(i=0;i<5;i++)\
 157         printf("res[%d]=%f\n", i, rtab[i]);\
 158     printf("fpuc=%04x fpus=%04x fptag=%04x\n",\
 159            (env)->fpuc,\
 160            (env)->fpus & 0xff00,\
 161            (env)->fptag);\
 162 }
 163
 164 void test_fenv(void)
 165 {
 166     struct __attribute__((packed)) {
 167         uint16_t fpuc;
 168         uint16_t dummy1;
 169         uint16_t fpus;
 170         uint16_t dummy2;
 171         uint16_t fptag;
 172         uint16_t dummy3;
 173         uint32_t ignored[4];
 174         long double fpregs[8];
 175     } float_env32;
 176     struct __attribute__((packed)) {
 177         uint16_t fpuc;
 178         uint16_t fpus;
 179         uint16_t fptag;
 180         uint16_t ignored[4];
 181         long double fpregs[8];
 182     } float_env16;
 183     double dtab[8];
 184     double rtab[8];
 185     int i;
 186
 187     for(i=0;i<8;i++)
 188         dtab[i] = i + 1;
 189
 190     TEST_ENV(&float_env16, "data16/fnstenv", "data16/fldenv");
 191     TEST_ENV(&float_env16, "data16/fnsave", "data16/frstor");
 192     TEST_ENV(&float_env32, "fnstenv", "fldenv");
 193     TEST_ENV(&float_env32, "fnsave", "frstor");
 194
 195     /* test for ffree */
 196     for(i=0;i<5;i++)
 197         asm volatile ("fldl %0" : : "m" (dtab[i]));
 198     asm volatile("ffree %st(2)");
 199     asm volatile ("fnstenv %0\n" : : "m" (float_env32));
 200     asm volatile ("fninit");
 201     printf("fptag=%04x\n", float_env32.fptag);
 202 }
 203
 204
 205 #define TEST_FCMOV(a, b, eflags, CC)\
 206 {\
 207     double res;\
 208     asm("push %3\n"\
 209         "popf\n"\
 210         "fcmov" CC " %2, %0\n"\
 211         : "=t" (res)\
 212         : "0" (a), "u" (b), "g" (eflags));\
 213     printf("fcmov%s eflags=0x%04x-> %f\n", \
 214            CC, eflags, res);\
 215 }
 216
 217 void test_fcmov(void)
 218 {
 219     double a, b;
 220     int eflags, i;
 221
 222     a = 1.0;
 223     b = 2.0;
 224     for(i = 0; i < 4; i++) {
 225         eflags = 0;
 226         if (i & 1)
 227             eflags |= CC_C;
 228         if (i & 2)
 229             eflags |= CC_Z;
 230         TEST_FCMOV(a, b, eflags, "b");
 231         TEST_FCMOV(a, b, eflags, "e");
 232         TEST_FCMOV(a, b, eflags, "be");
 233         TEST_FCMOV(a, b, eflags, "nb");
 234         TEST_FCMOV(a, b, eflags, "ne");
 235         TEST_FCMOV(a, b, eflags, "nbe");
 236     }
 237     TEST_FCMOV(a, b, 0, "u");
 238     TEST_FCMOV(a, b, CC_P, "u");
 239     TEST_FCMOV(a, b, 0, "nu");
 240     TEST_FCMOV(a, b, CC_P, "nu");
 241 }
 242
 243 void test_floats(void)
 244 {
 245     test_fops(2, 3);
 246     test_fops(1.4, -5);
 247     test_fcmp(2, -1);
 248     test_fcmp(2, 2);
 249     test_fcmp(2, 3);
 250     test_fcvt(0.5);
 251     test_fcvt(-0.5);
 252     test_fcvt(1.0/7.0);
 253     test_fcvt(-1.0/9.0);
 254     test_fcvt(32768);
 255     test_fcvt(-1e20);
 256     test_fconst();
 257 }
 258
 259 int main ( void )
 260 {
 261   test_floats();
 262   return 0;
 263 }