blob: 0f473aa1aadd6b180e49979d7fdf034242ec356d [file] [log] [blame]
Stephen Canonb6d4e2e2010-07-03 00:56:03 +00001//===-- lib/floatunsidf.c - uint -> double-precision conversion ---*- C -*-===//
Stephen Canon09009c52010-07-02 23:05:46 +00002//
3// The LLVM Compiler Infrastructure
4//
Howard Hinnant9ad441f2010-11-16 22:13:33 +00005// This file is dual licensed under the MIT and the University of Illinois Open
6// Source Licenses. See LICENSE.TXT for details.
Stephen Canon09009c52010-07-02 23:05:46 +00007//
8//===----------------------------------------------------------------------===//
9//
10// This file implements unsigned integer to double-precision conversion for the
11// compiler-rt library in the IEEE-754 default round-to-nearest, ties-to-even
12// mode.
13//
14//===----------------------------------------------------------------------===//
15
16#define DOUBLE_PRECISION
17#include "fp_lib.h"
18
Stephen Canonb6d4e2e2010-07-03 00:56:03 +000019fp_t __floatunsidf(unsigned int a) {
Stephen Canon09009c52010-07-02 23:05:46 +000020
21 const int aWidth = sizeof a * CHAR_BIT;
22
23 // Handle zero as a special case to protect clz
24 if (a == 0) return fromRep(0);
25
26 // Exponent of (fp_t)a is the width of abs(a).
27 const int exponent = (aWidth - 1) - __builtin_clz(a);
28 rep_t result;
29
Stephen Canon4d055d52010-08-17 19:13:45 +000030 // Shift a into the significand field and clear the implicit bit.
31 const int shift = significandBits - exponent;
32 result = (rep_t)a << shift ^ implicitBit;
Stephen Canon09009c52010-07-02 23:05:46 +000033
34 // Insert the exponent
35 result += (rep_t)(exponent + exponentBias) << significandBits;
36 return fromRep(result);
37}