blob: ba6c2cfd2a1742b44030677859a73df547e522db [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
Anton Korobeynikov37b97d12011-04-19 17:51:24 +000019#include "int_lib.h"
20
Chandler Carruth0193b742012-06-22 21:09:15 +000021ARM_EABI_FNALIAS(ui2d, floatunsidf)
Anton Korobeynikov37b97d12011-04-19 17:51:24 +000022
Stephen Canonb6d4e2e2010-07-03 00:56:03 +000023fp_t __floatunsidf(unsigned int a) {
Stephen Canon09009c52010-07-02 23:05:46 +000024
25 const int aWidth = sizeof a * CHAR_BIT;
26
27 // Handle zero as a special case to protect clz
28 if (a == 0) return fromRep(0);
29
30 // Exponent of (fp_t)a is the width of abs(a).
31 const int exponent = (aWidth - 1) - __builtin_clz(a);
32 rep_t result;
33
Stephen Canon4d055d52010-08-17 19:13:45 +000034 // Shift a into the significand field and clear the implicit bit.
35 const int shift = significandBits - exponent;
36 result = (rep_t)a << shift ^ implicitBit;
Stephen Canon09009c52010-07-02 23:05:46 +000037
38 // Insert the exponent
39 result += (rep_t)(exponent + exponentBias) << significandBits;
40 return fromRep(result);
41}