162306a36Sopenharmony_ci// SPDX-License-Identifier: GPL-2.0-only 262306a36Sopenharmony_ci/* IEEE754 floating point arithmetic 362306a36Sopenharmony_ci * single precision 462306a36Sopenharmony_ci */ 562306a36Sopenharmony_ci/* 662306a36Sopenharmony_ci * MIPS floating point support 762306a36Sopenharmony_ci * Copyright (C) 1994-2000 Algorithmics Ltd. 862306a36Sopenharmony_ci */ 962306a36Sopenharmony_ci 1062306a36Sopenharmony_ci#include "ieee754sp.h" 1162306a36Sopenharmony_ci 1262306a36Sopenharmony_ciint ieee754sp_tint(union ieee754sp x) 1362306a36Sopenharmony_ci{ 1462306a36Sopenharmony_ci u32 residue; 1562306a36Sopenharmony_ci int round; 1662306a36Sopenharmony_ci int sticky; 1762306a36Sopenharmony_ci int odd; 1862306a36Sopenharmony_ci 1962306a36Sopenharmony_ci COMPXSP; 2062306a36Sopenharmony_ci 2162306a36Sopenharmony_ci ieee754_clearcx(); 2262306a36Sopenharmony_ci 2362306a36Sopenharmony_ci EXPLODEXSP; 2462306a36Sopenharmony_ci FLUSHXSP; 2562306a36Sopenharmony_ci 2662306a36Sopenharmony_ci switch (xc) { 2762306a36Sopenharmony_ci case IEEE754_CLASS_SNAN: 2862306a36Sopenharmony_ci case IEEE754_CLASS_QNAN: 2962306a36Sopenharmony_ci ieee754_setcx(IEEE754_INVALID_OPERATION); 3062306a36Sopenharmony_ci return ieee754si_indef(); 3162306a36Sopenharmony_ci 3262306a36Sopenharmony_ci case IEEE754_CLASS_INF: 3362306a36Sopenharmony_ci ieee754_setcx(IEEE754_INVALID_OPERATION); 3462306a36Sopenharmony_ci return ieee754si_overflow(xs); 3562306a36Sopenharmony_ci 3662306a36Sopenharmony_ci case IEEE754_CLASS_ZERO: 3762306a36Sopenharmony_ci return 0; 3862306a36Sopenharmony_ci 3962306a36Sopenharmony_ci case IEEE754_CLASS_DNORM: 4062306a36Sopenharmony_ci case IEEE754_CLASS_NORM: 4162306a36Sopenharmony_ci break; 4262306a36Sopenharmony_ci } 4362306a36Sopenharmony_ci if (xe >= 31) { 4462306a36Sopenharmony_ci /* look for valid corner case */ 4562306a36Sopenharmony_ci if (xe == 31 && xs && xm == SP_HIDDEN_BIT) 4662306a36Sopenharmony_ci return -0x80000000; 4762306a36Sopenharmony_ci /* Set invalid. We will only use overflow for floating 4862306a36Sopenharmony_ci point overflow */ 4962306a36Sopenharmony_ci ieee754_setcx(IEEE754_INVALID_OPERATION); 5062306a36Sopenharmony_ci return ieee754si_overflow(xs); 5162306a36Sopenharmony_ci } 5262306a36Sopenharmony_ci /* oh gawd */ 5362306a36Sopenharmony_ci if (xe > SP_FBITS) { 5462306a36Sopenharmony_ci xm <<= xe - SP_FBITS; 5562306a36Sopenharmony_ci } else { 5662306a36Sopenharmony_ci if (xe < -1) { 5762306a36Sopenharmony_ci residue = xm; 5862306a36Sopenharmony_ci round = 0; 5962306a36Sopenharmony_ci sticky = residue != 0; 6062306a36Sopenharmony_ci xm = 0; 6162306a36Sopenharmony_ci } else { 6262306a36Sopenharmony_ci /* Shifting a u32 32 times does not work, 6362306a36Sopenharmony_ci * so we do it in two steps. Be aware that xe 6462306a36Sopenharmony_ci * may be -1 */ 6562306a36Sopenharmony_ci residue = xm << (xe + 1); 6662306a36Sopenharmony_ci residue <<= 31 - SP_FBITS; 6762306a36Sopenharmony_ci round = (residue >> 31) != 0; 6862306a36Sopenharmony_ci sticky = (residue << 1) != 0; 6962306a36Sopenharmony_ci xm >>= SP_FBITS - xe; 7062306a36Sopenharmony_ci } 7162306a36Sopenharmony_ci odd = (xm & 0x1) != 0x0; 7262306a36Sopenharmony_ci switch (ieee754_csr.rm) { 7362306a36Sopenharmony_ci case FPU_CSR_RN: 7462306a36Sopenharmony_ci if (round && (sticky || odd)) 7562306a36Sopenharmony_ci xm++; 7662306a36Sopenharmony_ci break; 7762306a36Sopenharmony_ci case FPU_CSR_RZ: 7862306a36Sopenharmony_ci break; 7962306a36Sopenharmony_ci case FPU_CSR_RU: /* toward +Infinity */ 8062306a36Sopenharmony_ci if ((round || sticky) && !xs) 8162306a36Sopenharmony_ci xm++; 8262306a36Sopenharmony_ci break; 8362306a36Sopenharmony_ci case FPU_CSR_RD: /* toward -Infinity */ 8462306a36Sopenharmony_ci if ((round || sticky) && xs) 8562306a36Sopenharmony_ci xm++; 8662306a36Sopenharmony_ci break; 8762306a36Sopenharmony_ci } 8862306a36Sopenharmony_ci if ((xm >> 31) != 0) { 8962306a36Sopenharmony_ci /* This can happen after rounding */ 9062306a36Sopenharmony_ci ieee754_setcx(IEEE754_INVALID_OPERATION); 9162306a36Sopenharmony_ci return ieee754si_overflow(xs); 9262306a36Sopenharmony_ci } 9362306a36Sopenharmony_ci if (round || sticky) 9462306a36Sopenharmony_ci ieee754_setcx(IEEE754_INEXACT); 9562306a36Sopenharmony_ci } 9662306a36Sopenharmony_ci if (xs) 9762306a36Sopenharmony_ci return -xm; 9862306a36Sopenharmony_ci else 9962306a36Sopenharmony_ci return xm; 10062306a36Sopenharmony_ci} 101