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, see <http://www.gnu.org/licenses/>.
  21  */
  22
  23
  24 #include <stdlib.h>
  25 #include <stdio.h>
  26 #include <string.h>
  27 #include <inttypes.h>
  28 #include <math.h>
  29
  30 /**********************************************/
  31
  32 void test_fops(double a, double b)
  33 {
  34     printf("a=%f b=%f a+b=%f\n", a, b, a + b);
  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 fmod(a, b)=%f\n", a, b, fmod(a, b));
  39     printf("a=%f sqrt(a)=%f\n", a, sqrt(a));
  40     printf("a=%f sin(a)=%f\n", a, sin(a));
  41     printf("a=%f cos(a)=%f\n", a, cos(a));
  42     printf("a=%f tan(a)=%f\n", a, tan(a));
  43     printf("a=%f log(a)=%f\n", a, log(a));
  44     printf("a=%f exp(a)=%f\n", a, exp(a));
  45     printf("a=%f b=%f atan2(a, b)=%f\n", a, b, atan2(a, b));
  46     /* just to test some op combining */
  47     printf("a=%f asin(sin(a))=%f\n", a, asin(sin(a)));
  48     printf("a=%f acos(cos(a))=%f\n", a, acos(cos(a)));
  49     printf("a=%f atan(tan(a))=%f\n", a, atan(tan(a)));
  50 }
  51 #define CC_C    0x0001
  52 #define CC_P    0x0004
  53 #define CC_A    0x0010
  54 #define CC_Z    0x0040
  55 #define CC_S    0x0080
  56 #define CC_O    0x0800
  57
  58
  59 void test_fcmp(double a, double b)
  60 {
  61     printf("(%f<%f)=%d\n",
  62            a, b, a < b);
  63     printf("(%f<=%f)=%d\n",
  64            a, b, a <= b);
  65     printf("(%f==%f)=%d\n",
  66            a, b, a == b);
  67     printf("(%f>%f)=%d\n",
  68            a, b, a > b);
  69     printf("(%f<=%f)=%d\n",
  70            a, b, a >= b);
  71     {
  72         unsigned int eflags;
  73         /* test f(u)comi instruction */
  74         asm("fcomi %2, %1\n"
  75             "pushf\n"
  76             "pop %0\n"
  77             : "=r" (eflags)
  78             : "t" (a), "u" (b));
  79         printf("fcomi(%f %f)=%08x\n", a, b, eflags & (CC_Z | CC_P | CC_C));
  80     }
  81 }
  82
  83 void test_fcvt(double a)
  84 {
  85     float fa;
  86     long double la;
  87     int16_t fpuc;
  88     int i;
  89     int64_t lla;
  90     int ia;
  91     int16_t wa;
  92     double ra;
  93
  94     fa = a;
  95     la = a;
  96     printf("(float)%f = %f\n", a, fa);
  97     printf("(long double)%f = %Lf\n", a, la);
  98     printf("a=%016llx\n", *(long long *)&a);
  99     printf("la=%016llx %04x\n", *(long long *)&la,
 100            *(unsigned short *)((char *)(&la) + 8));
 101
 102     /* test all roundings */
 103     asm volatile ("fstcw %0" : "=m" (fpuc));
 104     for(i=0;i<4;i++) {
 105       int16_t tmp = (fpuc & ~0x0c00) | (i << 10);
 106         asm volatile ("fldcw %0" : : "m" (tmp));
 107         asm volatile ("fists %0" : "=m" (wa) : "t" (a));
 108         asm volatile ("fistl %0" : "=m" (ia) : "t" (a));
 109         asm volatile ("fistpll %0" : "=m" (lla) : "t" (a) : "st");
 110         asm volatile ("frndint ; fstl %0" : "=m" (ra) : "t" (a));
 111         asm volatile ("fldcw %0" : : "m" (fpuc));
 112         printf("(short)a = %d\n", wa);
 113         printf("(int)a = %d\n", ia);
 114         printf("(int64_t)a = %lld\n", lla);
 115         printf("rint(a) = %f\n", ra);
 116     }
 117 }
 118
 119 #define TEST(N) \
 120     asm("fld" #N : "=t" (a)); \
 121     printf("fld" #N "= %f\n", a);
 122
 123 void test_fconst(void)
 124 {
 125     double a;
 126     TEST(1);
 127     TEST(l2t);
 128     TEST(l2e);
 129     TEST(pi);
 130     TEST(lg2);
 131     TEST(ln2);
 132     TEST(z);
 133 }
 134
 135 void test_fbcd(double a)
 136 {
 137     unsigned short bcd[5];
 138     double b;
 139
 140     asm("fbstp %0" : "=m" (bcd[0]) : "t" (a) : "st");
 141     asm("fbld %1" : "=t" (b) : "m" (bcd[0]));
 142     printf("a=%f bcd=%04x%04x%04x%04x%04x b=%f\n",
 143            a, bcd[4], bcd[3], bcd[2], bcd[1], bcd[0], b);
 144 }
 145
 146 #define TEST_ENV(env, save, restore)\
 147 {\
 148     memset((env), 0xaa, sizeof(*(env)));\
 149     for(i=0;i<5;i++)\
 150         asm volatile ("fldl %0" : : "m" (dtab[i]));\
 151     asm(save " %0\n" : : "m" (*(env)));\
 152     asm(restore " %0\n": : "m" (*(env)));\
 153     for(i=0;i<5;i++)\
 154         asm volatile ("fstpl %0" : "=m" (rtab[i]));\
 155     for(i=0;i<5;i++)\
 156         printf("res[%d]=%f\n", i, rtab[i]);\
 157     printf("fpuc=%04x fpus=%04x fptag=%04x\n",\
 158            (env)->fpuc,\
 159            (env)->fpus & 0xff00,\
 160            (env)->fptag);\
 161 }
 162
 163 void test_fenv(void)
 164 {
 165     struct __attribute__((packed)) {
 166         uint16_t fpuc;
 167         uint16_t dummy1;
 168         uint16_t fpus;
 169         uint16_t dummy2;
 170         uint16_t fptag;
 171         uint16_t dummy3;
 172         uint32_t ignored[4];
 173         long double fpregs[8];
 174     } float_env32;
 175     struct __attribute__((packed)) {
 176         uint16_t fpuc;
 177         uint16_t fpus;
 178         uint16_t fptag;
 179         uint16_t ignored[4];
 180         long double fpregs[8];
 181     } float_env16;
 182     double dtab[8];
 183     double rtab[8];
 184     int i;
 185
 186     for(i=0;i<8;i++)
 187         dtab[i] = i + 1;
 188
 189     TEST_ENV(&float_env16, "data16/fnstenv", "data16/fldenv");
 190     TEST_ENV(&float_env16, "data16/fnsave", "data16/frstor");
 191     TEST_ENV(&float_env32, "fnstenv", "fldenv");
 192     TEST_ENV(&float_env32, "fnsave", "frstor");
 193
 194     /* test for ffree */
 195     for(i=0;i<5;i++)
 196         asm volatile ("fldl %0" : : "m" (dtab[i]));
 197     asm volatile("ffree %st(2)");
 198     asm volatile ("fnstenv %0\n" : : "m" (float_env32));
 199     asm volatile ("fninit");
 200     printf("fptag=%04x\n", float_env32.fptag);
 201 }
 202
 203
 204 #define TEST_FCMOV(a, b, eflags, CC)\
 205 {\
 206     double res;\
 207     asm("push %3\n"\
 208         "popf\n"\
 209         "fcmov" CC " %2, %0\n"\
 210         : "=t" (res)\
 211         : "0" (a), "u" (b), "g" (eflags));\
 212     printf("fcmov%s eflags=0x%04x-> %f\n", \
 213            CC, eflags, res);\
 214 }
 215
 216 void test_fcmov(void)
 217 {
 218     double a, b;
 219     int eflags, i;
 220
 221     a = 1.0;
 222     b = 2.0;
 223     for(i = 0; i < 4; i++) {
 224         eflags = 0;
 225         if (i & 1)
 226             eflags |= CC_C;
 227         if (i & 2)
 228             eflags |= CC_Z;
 229         TEST_FCMOV(a, b, eflags, "b");
 230         TEST_FCMOV(a, b, eflags, "e");
 231         TEST_FCMOV(a, b, eflags, "be");
 232         TEST_FCMOV(a, b, eflags, "nb");
 233         TEST_FCMOV(a, b, eflags, "ne");
 234         TEST_FCMOV(a, b, eflags, "nbe");
 235     }
 236     TEST_FCMOV(a, b, 0, "u");
 237     TEST_FCMOV(a, b, CC_P, "u");
 238     TEST_FCMOV(a, b, 0, "nu");
 239     TEST_FCMOV(a, b, CC_P, "nu");
 240 }
 241
 242 void test_floats(void)
 243 {
 244     test_fops(2, 3);
 245     test_fops(1.4, -5);
 246     test_fcmp(2, -1);
 247     test_fcmp(2, 2);
 248     test_fcmp(2, 3);
 249     test_fcvt(0.5);
 250     test_fcvt(-0.5);
 251     test_fcvt(1.0/7.0);
 252     test_fcvt(-1.0/9.0);
 253     test_fcvt(32768);
 254     test_fcvt(-1e20);
 255     test_fconst();
 256 }
 257
 258 int main ( void )
 259 {
 260   test_floats();
 261   return 0;
 262 }