4  * SPDX-License-Identifier: BSD-3-Clause
10  * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
64 	u_int m0, m1, m2, m3, s;  in fpu_shr()  local
68 	if (rsh <= 0 || (fp->fp_class != FPC_NUM && !ISNAN(fp)))  in fpu_shr()
72 	m0 = fp->fp_mant[0];  in fpu_shr()
73 	m1 = fp->fp_mant[1];  in fpu_shr()
74 	m2 = fp->fp_mant[2];  in fpu_shr()
75 	m3 = fp->fp_mant[3];  in fpu_shr()
80 		if ((m0 | m1 | m2 | m3) == 0)  in fpu_shr()
83 		fp->fp_mant[0] = 0;  in fpu_shr()
84 		fp->fp_mant[1] = 0;  in fpu_shr()
85 		fp->fp_mant[2] = 0;  in fpu_shr()
86 		fp->fp_mant[3] = 0;  in fpu_shr()
88 		if ((m0 | m1 | m2 | m3) == 0)  in fpu_shr()
89 			fp->fp_class = FPC_ZERO;  in fpu_shr()
92 			fp->fp_sticky = 1;  in fpu_shr()
97 	s = fp->fp_sticky;  in fpu_shr()
99 		s |= m3 | m2 | m1;  in fpu_shr()
100 		m3 = m0, m2 = 0, m1 = 0, m0 = 0;  in fpu_shr()
103 		m3 = m1, m2 = m0, m1 = 0, m0 = 0;  in fpu_shr()
106 		m3 = m2, m2 = m1, m1 = m0, m0 = 0;  in fpu_shr()
111 		lsh = 32 - rsh;  in fpu_shr()
114 		m2 = (m2 >> rsh) | (m1 << lsh);  in fpu_shr()
115 		m1 = (m1 >> rsh) | (m0 << lsh);  in fpu_shr()
118 	fp->fp_mant[0] = m0;  in fpu_shr()
119 	fp->fp_mant[1] = m1;  in fpu_shr()
120 	fp->fp_mant[2] = m2;  in fpu_shr()
121 	fp->fp_mant[3] = m3;  in fpu_shr()
122 	fp->fp_sticky = s;  in fpu_shr()
131  * Internally, this may use a `supernormal' -- a number whose fp_mant
132  * is greater than or equal to 2.0 -- so as a side effect you can hand it
133  * a supernormal and it will fix it (provided fp->fp_mant[3] == 0).
138 	u_int m0, m1, m2, m3, top, sup, nrm;  in fpu_norm()  local
141 	exp = fp->fp_exp;  in fpu_norm()
142 	m0 = fp->fp_mant[0];  in fpu_norm()
143 	m1 = fp->fp_mant[1];  in fpu_norm()
144 	m2 = fp->fp_mant[2];  in fpu_norm()
145 	m3 = fp->fp_mant[3];  in fpu_norm()
147 	/* Handle severe subnormals with 32-bit moves. */  in fpu_norm()
149 		if (m1)  in fpu_norm()
150 			m0 = m1, m1 = m2, m2 = m3, m3 = 0, exp -= 32;  in fpu_norm()
152 			m0 = m2, m1 = m3, m2 = 0, m3 = 0, exp -= 2 * 32;  in fpu_norm()
154 			m0 = m3, m1 = 0, m2 = 0, m3 = 0, exp -= 3 * 32;  in fpu_norm()
156 			fp->fp_class = FPC_ZERO;  in fpu_norm()
172 		lsh = 32 - rsh;  in fpu_norm()
174 		m2 = (m2 >> rsh) | (m1 << lsh);  in fpu_norm()
175 		m1 = (m1 >> rsh) | (m0 << lsh);  in fpu_norm()
184 		exp -= lsh;  in fpu_norm()
185 		rsh = 32 - lsh;  in fpu_norm()
186 		m0 = top | (m1 >> rsh);  in fpu_norm()
187 		m1 = (m1 << lsh) | (m2 >> rsh);  in fpu_norm()
192 	fp->fp_exp = exp;  in fpu_norm()
193 	fp->fp_mant[0] = m0;  in fpu_norm()
194 	fp->fp_mant[1] = m1;  in fpu_norm()
195 	fp->fp_mant[2] = m2;  in fpu_norm()
196 	fp->fp_mant[3] = m3;  in fpu_norm()
208 	fe->fe_cx |= FPSCR_VXSNAN;  in fpu_newnan()
209 	fp = &fe->fe_f3;  in fpu_newnan()
210 	fp->fp_class = FPC_QNAN;  in fpu_newnan()
211 	fp->fp_sign = 0;  in fpu_newnan()
212 	fp->fp_mant[0] = FP_1 - 1;  in fpu_newnan()
213 	fp->fp_mant[1] = fp->fp_mant[2] = fp->fp_mant[3] = ~0;  in fpu_newnan()