1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24#ifndef __MATH_EMU_SOFT_FP_H__
25#define __MATH_EMU_SOFT_FP_H__
26
27#include <asm/sfp-machine.h>
28
29
30#ifndef __BYTE_ORDER
31#include <endian.h>
32#endif
33
34#define _FP_WORKBITS 3
35#define _FP_WORK_LSB ((_FP_W_TYPE)1 << 3)
36#define _FP_WORK_ROUND ((_FP_W_TYPE)1 << 2)
37#define _FP_WORK_GUARD ((_FP_W_TYPE)1 << 1)
38#define _FP_WORK_STICKY ((_FP_W_TYPE)1 << 0)
39
40#ifndef FP_RND_NEAREST
41# define FP_RND_NEAREST 0
42# define FP_RND_ZERO 1
43# define FP_RND_PINF 2
44# define FP_RND_MINF 3
45#ifndef FP_ROUNDMODE
46# define FP_ROUNDMODE FP_RND_NEAREST
47#endif
48#endif
49
50
51#ifndef FP_EX_INVALID
52#define FP_EX_INVALID 0
53#endif
54#ifndef FP_EX_OVERFLOW
55#define FP_EX_OVERFLOW 0
56#endif
57#ifndef FP_EX_UNDERFLOW
58#define FP_EX_UNDERFLOW
59#endif
60#ifndef FP_EX_DIVZERO
61#define FP_EX_DIVZERO 0
62#endif
63#ifndef FP_EX_INEXACT
64#define FP_EX_INEXACT 0
65#endif
66#ifndef FP_EX_DENORM
67#define FP_EX_DENORM 0
68#endif
69
70#ifdef _FP_DECL_EX
71#define FP_DECL_EX \
72 int _fex = 0; \
73 _FP_DECL_EX
74#else
75#define FP_DECL_EX int _fex = 0
76#endif
77
78#ifndef FP_INIT_ROUNDMODE
79#define FP_INIT_ROUNDMODE do {} while (0)
80#endif
81
82#ifndef FP_HANDLE_EXCEPTIONS
83#define FP_HANDLE_EXCEPTIONS do {} while (0)
84#endif
85
86
87#ifndef FP_DENORM_ZERO
88#define FP_DENORM_ZERO 0
89#endif
90
91#ifndef FP_INHIBIT_RESULTS
92
93
94
95
96
97#define FP_INHIBIT_RESULTS 0
98#endif
99
100#ifndef FP_TRAPPING_EXCEPTIONS
101#define FP_TRAPPING_EXCEPTIONS 0
102#endif
103
104#define FP_SET_EXCEPTION(ex) \
105 _fex |= (ex)
106
107#define FP_UNSET_EXCEPTION(ex) \
108 _fex &= ~(ex)
109
110#define FP_CUR_EXCEPTIONS \
111 (_fex)
112
113#define FP_CLEAR_EXCEPTIONS \
114 _fex = 0
115
116#define _FP_ROUND_NEAREST(wc, X) \
117do { \
118 if ((_FP_FRAC_LOW_##wc(X) & 15) != _FP_WORK_ROUND) \
119 _FP_FRAC_ADDI_##wc(X, _FP_WORK_ROUND); \
120} while (0)
121
122#define _FP_ROUND_ZERO(wc, X) 0
123
124#define _FP_ROUND_PINF(wc, X) \
125do { \
126 if (!X##_s && (_FP_FRAC_LOW_##wc(X) & 7)) \
127 _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB); \
128} while (0)
129
130#define _FP_ROUND_MINF(wc, X) \
131do { \
132 if (X##_s && (_FP_FRAC_LOW_##wc(X) & 7)) \
133 _FP_FRAC_ADDI_##wc(X, _FP_WORK_LSB); \
134} while (0)
135
136#define _FP_ROUND(wc, X) \
137do { \
138 if (_FP_FRAC_LOW_##wc(X) & 7) \
139 FP_SET_EXCEPTION(FP_EX_INEXACT); \
140 switch (FP_ROUNDMODE) \
141 { \
142 case FP_RND_NEAREST: \
143 _FP_ROUND_NEAREST(wc,X); \
144 break; \
145 case FP_RND_ZERO: \
146 _FP_ROUND_ZERO(wc,X); \
147 break; \
148 case FP_RND_PINF: \
149 _FP_ROUND_PINF(wc,X); \
150 break; \
151 case FP_RND_MINF: \
152 _FP_ROUND_MINF(wc,X); \
153 break; \
154 } \
155} while (0)
156
157#define FP_CLS_NORMAL 0
158#define FP_CLS_ZERO 1
159#define FP_CLS_INF 2
160#define FP_CLS_NAN 3
161
162#define _FP_CLS_COMBINE(x,y) (((x) << 2) | (y))
163
164#include <math-emu/op-1.h>
165#include <math-emu/op-2.h>
166#include <math-emu/op-4.h>
167#include <math-emu/op-8.h>
168#include <math-emu/op-common.h>
169
170
171#define UWtype _FP_W_TYPE
172#define W_TYPE_SIZE _FP_W_TYPE_SIZE
173
174typedef int SItype __attribute__((mode(SI)));
175typedef int DItype __attribute__((mode(DI)));
176typedef unsigned int USItype __attribute__((mode(SI)));
177typedef unsigned int UDItype __attribute__((mode(DI)));
178#if _FP_W_TYPE_SIZE == 32
179typedef unsigned int UHWtype __attribute__((mode(HI)));
180#elif _FP_W_TYPE_SIZE == 64
181typedef USItype UHWtype;
182#endif
183
184#ifndef umul_ppmm
185#include <stdlib/longlong.h>
186#endif
187
188#endif
189
