summary refs log tree commit diff
path: root/soft-fp/sysdeps/i386/sfp-machine.h
blob: 047452b905d9214d837c74ebf8b73b08fb978904 (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
#define _FP_W_TYPE_SIZE		32
#define _FP_W_TYPE		unsigned long
#define _FP_WS_TYPE		signed long
#define _FP_I_TYPE		long

#define __FP_FRAC_ADD_2(rh, rl, xh, xl, yh, yl)				\
  __asm__("addl %5,%1; adcl %3,%0"					\
	  : "=r"(rh), "=r"(rl)						\
	  : "%0"(xh), "g"(yh), "%1"(xl), "g"(yl)			\
	  : "cc")

#define __FP_FRAC_ADD_4(r3,r2,r1,r0,x3,x2,x1,x0,y3,y2,y1,y0)		\
  do { 									\
    __asm__ volatile("addl %5,%1; adcl %3,%0"				\
		     : "=r"(r1), "=r"(r0)				\
		     : "%0"(x1), "g"(y1), "%1"(x0), "g"(y0)		\
		     : "cc");						\
    __asm__ volatile("adcl %5,%1; adcl %3,%0"				\
		     : "=r"(r3), "=r"(r2)				\
		     : "%0"(x3), "g"(y3), "%1"(x2), "g"(y2)		\
		     : "cc");						\
  } while (0)

#define __FP_FRAC_SUB_2(rh, rl, xh, xl, yh, yl)				\
  __asm__("subl %5,%1; sbbl %4,%0"					\
	  : "=r"(rh), "=r"(rl)						\
	  : "0"(xh), "1"(xl), "g"(yh), "g"(yl)				\
	  : "cc")

#define __FP_CLZ(r, x)							\
  do {									\
    __asm__("bsrl %1,%0" : "=r"(r) : "g"(x) : "cc");			\
    r ^= 31;								\
  } while (0)

#define _i386_mul_32_64(rh, rl, x, y)					\
  __asm__("mull %2" : "=d"(rh), "=a"(rl) : "%g"(x), "1"(y) : "cc")

#define _i386_div_64_32(q, r, nh, nl, d)				\
  __asm__ ("divl %4" : "=a"(q), "=d"(r) : "0"(nl), "1"(nh), "g"(d) : "cc")


#define _FP_MUL_MEAT_S(R,X,Y)					\
  _FP_MUL_MEAT_1_wide(_FP_WFRACBITS_S,R,X,Y,_i386_mul_32_64)
#define _FP_MUL_MEAT_D(R,X,Y)					\
  _FP_MUL_MEAT_2_wide(_FP_WFRACBITS_D,R,X,Y,_i386_mul_32_64)

#define _FP_DIV_MEAT_S(R,X,Y)	_FP_DIV_MEAT_1_udiv(S,R,X,Y)
#define _FP_DIV_MEAT_D(R,X,Y)	_FP_DIV_MEAT_2_udiv(D,R,X,Y)

#define _FP_NANFRAC_S		_FP_QNANBIT_S
#define _FP_NANFRAC_D		_FP_QNANBIT_D, 0
#define _FP_NANSIGN_S		1
#define _FP_NANSIGN_D		1

#define _FP_KEEPNANFRACP 1
/* Here is something Intel misdesigned: the specs don't define
   the case where we have two NaNs with same mantissas, but
   different sign. Different operations pick up different NaNs.
 */
#define _FP_CHOOSENAN(fs, wc, R, X, Y, OP)			\
  do {								\
    if (_FP_FRAC_GT_##wc(X, Y)					\
	|| (_FP_FRAC_EQ_##wc(X,Y) && (OP == '+' || OP == '*')))	\
      {								\
	R##_s = X##_s;						\
        _FP_FRAC_COPY_##wc(R,X);				\
      }								\
    else							\
      {								\
	R##_s = Y##_s;						\
        _FP_FRAC_COPY_##wc(R,Y);				\
      }								\
    R##_c = FP_CLS_NAN;						\
  } while (0)

#define FP_EX_INVALID           (1 << 0)
#define FP_EX_DENORM		(1 << 1)
#define FP_EX_DIVZERO           (1 << 2)
#define FP_EX_OVERFLOW          (1 << 3)
#define FP_EX_UNDERFLOW         (1 << 4)
#define FP_EX_INEXACT           (1 << 5)

#define FP_RND_NEAREST		0
#define FP_RND_ZERO		3
#define FP_RND_PINF		2
#define FP_RND_MINF		1