floatunsisf.c revision 214152
1214152Sed//===-- lib/floatunsisf.c - uint -> single-precision conversion ---*- C -*-===// 2214152Sed// 3214152Sed// The LLVM Compiler Infrastructure 4214152Sed// 5214152Sed// This file is distributed under the University of Illinois Open Source 6214152Sed// License. See LICENSE.TXT for details. 7214152Sed// 8214152Sed//===----------------------------------------------------------------------===// 9214152Sed// 10214152Sed// This file implements unsigned integer to single-precision conversion for the 11214152Sed// compiler-rt library in the IEEE-754 default round-to-nearest, ties-to-even 12214152Sed// mode. 13214152Sed// 14214152Sed//===----------------------------------------------------------------------===// 15214152Sed 16214152Sed#define SINGLE_PRECISION 17214152Sed#include "fp_lib.h" 18214152Sed 19214152Sedfp_t __floatunsisf(unsigned int a) { 20214152Sed 21214152Sed const int aWidth = sizeof a * CHAR_BIT; 22214152Sed 23214152Sed // Handle zero as a special case to protect clz 24214152Sed if (a == 0) return fromRep(0); 25214152Sed 26214152Sed // Exponent of (fp_t)a is the width of abs(a). 27214152Sed const int exponent = (aWidth - 1) - __builtin_clz(a); 28214152Sed rep_t result; 29214152Sed 30214152Sed // Shift a into the significand field, rounding if it is a right-shift 31214152Sed if (exponent <= significandBits) { 32214152Sed const int shift = significandBits - exponent; 33214152Sed result = (rep_t)a << shift ^ implicitBit; 34214152Sed } else { 35214152Sed const int shift = exponent - significandBits; 36214152Sed result = (rep_t)a >> shift ^ implicitBit; 37214152Sed rep_t round = (rep_t)a << (typeWidth - shift); 38214152Sed if (round > signBit) result++; 39214152Sed if (round == signBit) result += result & 1; 40214152Sed } 41214152Sed 42214152Sed // Insert the exponent 43214152Sed result += (rep_t)(exponent + exponentBias) << significandBits; 44214152Sed return fromRep(result); 45214152Sed} 46