2 * IEEE754 floating point arithmetic
3 * single precision: MIN{,A}.f
4 * MIN : Scalar Floating-Point Minimum
5 * MINA: Scalar Floating-Point argument with Minimum Absolute Value
7 * MIN.S : FPR[fd] = minNum(FPR[fs],FPR[ft])
8 * MINA.S: FPR[fd] = maxNumMag(FPR[fs],FPR[ft])
10 * MIPS floating point support
11 * Copyright (C) 2015 Imagination Technologies, Ltd.
12 * Author: Markos Chandras <markos.chandras@imgtec.com>
14 * This program is free software; you can distribute it and/or modify it
15 * under the terms of the GNU General Public License as published by the
16 * Free Software Foundation; version 2 of the License.
19 #include "ieee754sp.h"
21 union ieee754sp
ieee754sp_fmin(union ieee754sp x
, union ieee754sp y
)
34 switch (CLPAIR(xc
, yc
)) {
35 case CLPAIR(IEEE754_CLASS_QNAN
, IEEE754_CLASS_SNAN
):
36 case CLPAIR(IEEE754_CLASS_ZERO
, IEEE754_CLASS_SNAN
):
37 case CLPAIR(IEEE754_CLASS_NORM
, IEEE754_CLASS_SNAN
):
38 case CLPAIR(IEEE754_CLASS_DNORM
, IEEE754_CLASS_SNAN
):
39 case CLPAIR(IEEE754_CLASS_INF
, IEEE754_CLASS_SNAN
):
40 return ieee754sp_nanxcpt(y
);
42 case CLPAIR(IEEE754_CLASS_SNAN
, IEEE754_CLASS_SNAN
):
43 case CLPAIR(IEEE754_CLASS_SNAN
, IEEE754_CLASS_QNAN
):
44 case CLPAIR(IEEE754_CLASS_SNAN
, IEEE754_CLASS_ZERO
):
45 case CLPAIR(IEEE754_CLASS_SNAN
, IEEE754_CLASS_NORM
):
46 case CLPAIR(IEEE754_CLASS_SNAN
, IEEE754_CLASS_DNORM
):
47 case CLPAIR(IEEE754_CLASS_SNAN
, IEEE754_CLASS_INF
):
48 return ieee754sp_nanxcpt(x
);
55 * The case of both inputs quiet NaNs
57 case CLPAIR(IEEE754_CLASS_QNAN
, IEEE754_CLASS_QNAN
):
61 * The cases of exactly one input quiet NaN (numbers
62 * are here preferred as returned values to NaNs)
64 case CLPAIR(IEEE754_CLASS_ZERO
, IEEE754_CLASS_QNAN
):
65 case CLPAIR(IEEE754_CLASS_NORM
, IEEE754_CLASS_QNAN
):
66 case CLPAIR(IEEE754_CLASS_DNORM
, IEEE754_CLASS_QNAN
):
67 case CLPAIR(IEEE754_CLASS_INF
, IEEE754_CLASS_QNAN
):
70 case CLPAIR(IEEE754_CLASS_QNAN
, IEEE754_CLASS_ZERO
):
71 case CLPAIR(IEEE754_CLASS_QNAN
, IEEE754_CLASS_NORM
):
72 case CLPAIR(IEEE754_CLASS_QNAN
, IEEE754_CLASS_DNORM
):
73 case CLPAIR(IEEE754_CLASS_QNAN
, IEEE754_CLASS_INF
):
77 * Infinity and zero handling
79 case CLPAIR(IEEE754_CLASS_INF
, IEEE754_CLASS_ZERO
):
80 case CLPAIR(IEEE754_CLASS_INF
, IEEE754_CLASS_NORM
):
81 case CLPAIR(IEEE754_CLASS_INF
, IEEE754_CLASS_DNORM
):
82 case CLPAIR(IEEE754_CLASS_NORM
, IEEE754_CLASS_ZERO
):
83 case CLPAIR(IEEE754_CLASS_DNORM
, IEEE754_CLASS_ZERO
):
86 case CLPAIR(IEEE754_CLASS_INF
, IEEE754_CLASS_INF
):
87 case CLPAIR(IEEE754_CLASS_NORM
, IEEE754_CLASS_INF
):
88 case CLPAIR(IEEE754_CLASS_DNORM
, IEEE754_CLASS_INF
):
89 case CLPAIR(IEEE754_CLASS_ZERO
, IEEE754_CLASS_INF
):
90 case CLPAIR(IEEE754_CLASS_ZERO
, IEEE754_CLASS_NORM
):
91 case CLPAIR(IEEE754_CLASS_ZERO
, IEEE754_CLASS_DNORM
):
94 case CLPAIR(IEEE754_CLASS_ZERO
, IEEE754_CLASS_ZERO
):
95 return ieee754sp_zero(xs
| ys
);
97 case CLPAIR(IEEE754_CLASS_DNORM
, IEEE754_CLASS_DNORM
):
101 case CLPAIR(IEEE754_CLASS_NORM
, IEEE754_CLASS_DNORM
):
105 case CLPAIR(IEEE754_CLASS_DNORM
, IEEE754_CLASS_NORM
):
109 /* Finally get to do some computation */
111 assert(xm
& SP_HIDDEN_BIT
);
112 assert(ym
& SP_HIDDEN_BIT
);
120 /* Signs of inputs are the same, let's compare exponents */
122 /* Inputs are both positive */
128 /* Inputs are both negative */
135 /* Signs and exponents of inputs are equal, let's compare mantissas */
137 /* Inputs are both positive, with equal signs and exponents */
142 /* Inputs are both negative, with equal signs and exponents */
148 union ieee754sp
ieee754sp_fmina(union ieee754sp x
, union ieee754sp y
)
161 switch (CLPAIR(xc
, yc
)) {
162 case CLPAIR(IEEE754_CLASS_QNAN
, IEEE754_CLASS_SNAN
):
163 case CLPAIR(IEEE754_CLASS_ZERO
, IEEE754_CLASS_SNAN
):
164 case CLPAIR(IEEE754_CLASS_NORM
, IEEE754_CLASS_SNAN
):
165 case CLPAIR(IEEE754_CLASS_DNORM
, IEEE754_CLASS_SNAN
):
166 case CLPAIR(IEEE754_CLASS_INF
, IEEE754_CLASS_SNAN
):
167 return ieee754sp_nanxcpt(y
);
169 case CLPAIR(IEEE754_CLASS_SNAN
, IEEE754_CLASS_SNAN
):
170 case CLPAIR(IEEE754_CLASS_SNAN
, IEEE754_CLASS_QNAN
):
171 case CLPAIR(IEEE754_CLASS_SNAN
, IEEE754_CLASS_ZERO
):
172 case CLPAIR(IEEE754_CLASS_SNAN
, IEEE754_CLASS_NORM
):
173 case CLPAIR(IEEE754_CLASS_SNAN
, IEEE754_CLASS_DNORM
):
174 case CLPAIR(IEEE754_CLASS_SNAN
, IEEE754_CLASS_INF
):
175 return ieee754sp_nanxcpt(x
);
182 * The case of both inputs quiet NaNs
184 case CLPAIR(IEEE754_CLASS_QNAN
, IEEE754_CLASS_QNAN
):
188 * The cases of exactly one input quiet NaN (numbers
189 * are here preferred as returned values to NaNs)
191 case CLPAIR(IEEE754_CLASS_ZERO
, IEEE754_CLASS_QNAN
):
192 case CLPAIR(IEEE754_CLASS_NORM
, IEEE754_CLASS_QNAN
):
193 case CLPAIR(IEEE754_CLASS_DNORM
, IEEE754_CLASS_QNAN
):
194 case CLPAIR(IEEE754_CLASS_INF
, IEEE754_CLASS_QNAN
):
197 case CLPAIR(IEEE754_CLASS_QNAN
, IEEE754_CLASS_ZERO
):
198 case CLPAIR(IEEE754_CLASS_QNAN
, IEEE754_CLASS_NORM
):
199 case CLPAIR(IEEE754_CLASS_QNAN
, IEEE754_CLASS_DNORM
):
200 case CLPAIR(IEEE754_CLASS_QNAN
, IEEE754_CLASS_INF
):
204 * Infinity and zero handling
206 case CLPAIR(IEEE754_CLASS_INF
, IEEE754_CLASS_INF
):
207 return ieee754sp_inf(xs
| ys
);
209 case CLPAIR(IEEE754_CLASS_INF
, IEEE754_CLASS_ZERO
):
210 case CLPAIR(IEEE754_CLASS_INF
, IEEE754_CLASS_NORM
):
211 case CLPAIR(IEEE754_CLASS_INF
, IEEE754_CLASS_DNORM
):
212 case CLPAIR(IEEE754_CLASS_NORM
, IEEE754_CLASS_ZERO
):
213 case CLPAIR(IEEE754_CLASS_DNORM
, IEEE754_CLASS_ZERO
):
216 case CLPAIR(IEEE754_CLASS_NORM
, IEEE754_CLASS_INF
):
217 case CLPAIR(IEEE754_CLASS_DNORM
, IEEE754_CLASS_INF
):
218 case CLPAIR(IEEE754_CLASS_ZERO
, IEEE754_CLASS_INF
):
219 case CLPAIR(IEEE754_CLASS_ZERO
, IEEE754_CLASS_NORM
):
220 case CLPAIR(IEEE754_CLASS_ZERO
, IEEE754_CLASS_DNORM
):
223 case CLPAIR(IEEE754_CLASS_ZERO
, IEEE754_CLASS_ZERO
):
224 return ieee754sp_zero(xs
| ys
);
226 case CLPAIR(IEEE754_CLASS_DNORM
, IEEE754_CLASS_DNORM
):
230 case CLPAIR(IEEE754_CLASS_NORM
, IEEE754_CLASS_DNORM
):
234 case CLPAIR(IEEE754_CLASS_DNORM
, IEEE754_CLASS_NORM
):
238 /* Finally get to do some computation */
240 assert(xm
& SP_HIDDEN_BIT
);
241 assert(ym
& SP_HIDDEN_BIT
);
243 /* Compare exponent */
249 /* Compare mantissa */