2 * IEEE754 floating point arithmetic
3 * double precision: MIN{,A}.f
4 * MIN : Scalar Floating-Point Minimum
5 * MINA: Scalar Floating-Point argument with Minimum Absolute Value
7 * MIN.D : FPR[fd] = minNum(FPR[fs],FPR[ft])
8 * MINA.D: 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 "ieee754dp.h"
21 union ieee754dp ieee754dp_fmin(union ieee754dp x, union ieee754dp 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 ieee754dp_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 ieee754dp_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 ieee754dp_zero(xs | ys);
97 case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_DNORM):
100 case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_DNORM):
104 case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_NORM):
108 /* Finally get to do some computation */
110 assert(xm & DP_HIDDEN_BIT);
111 assert(ym & DP_HIDDEN_BIT);
119 /* Signs of inputs are the same, let's compare exponents */
121 /* Inputs are both positive */
127 /* Inputs are both negative */
134 /* Signs and exponents of inputs are equal, let's compare mantissas */
136 /* Inputs are both positive, with equal signs and exponents */
141 /* Inputs are both negative, with equal signs and exponents */
147 union ieee754dp ieee754dp_fmina(union ieee754dp x, union ieee754dp y)
160 switch (CLPAIR(xc, yc)) {
161 case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_SNAN):
162 case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_SNAN):
163 case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_SNAN):
164 case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_SNAN):
165 case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_SNAN):
166 return ieee754dp_nanxcpt(y);
168 case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_SNAN):
169 case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_QNAN):
170 case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_ZERO):
171 case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_NORM):
172 case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_DNORM):
173 case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_INF):
174 return ieee754dp_nanxcpt(x);
181 * The case of both inputs quiet NaNs
183 case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_QNAN):
187 * The cases of exactly one input quiet NaN (numbers
188 * are here preferred as returned values to NaNs)
190 case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_QNAN):
191 case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_QNAN):
192 case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_QNAN):
193 case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_QNAN):
196 case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_ZERO):
197 case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_NORM):
198 case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_DNORM):
199 case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_INF):
203 * Infinity and zero handling
205 case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_INF):
206 return ieee754dp_inf(xs | ys);
208 case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_ZERO):
209 case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_NORM):
210 case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_DNORM):
211 case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_ZERO):
212 case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_ZERO):
215 case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_INF):
216 case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_INF):
217 case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_INF):
218 case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_NORM):
219 case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_DNORM):
222 case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_ZERO):
223 return ieee754dp_zero(xs | ys);
225 case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_DNORM):
228 case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_DNORM):
232 case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_NORM):
236 /* Finally get to do some computation */
238 assert(xm & DP_HIDDEN_BIT);
239 assert(ym & DP_HIDDEN_BIT);
241 /* Compare exponent */
247 /* Compare mantissa */