blob: 996b15ba07a073d314b05275d3785203ce96d3f8 [file] [log] [blame]
Thomas Gleixner9d5a6342019-05-31 01:09:55 -07001// SPDX-License-Identifier: GPL-2.0-only
Linus Torvalds1da177e2005-04-16 15:20:36 -07002/* IEEE754 floating point arithmetic
3 * double precision: common utilities
4 */
5/*
6 * MIPS floating point support
7 * Copyright (C) 1994-2000 Algorithmics Ltd.
Linus Torvalds1da177e2005-04-16 15:20:36 -07008 */
9
Linus Torvalds1da177e2005-04-16 15:20:36 -070010#include "ieee754dp.h"
11
Ralf Baechle2209bcb2014-04-16 01:31:11 +020012union ieee754dp ieee754dp_fint(int x)
Linus Torvalds1da177e2005-04-16 15:20:36 -070013{
Atsushi Nemotoecf52d32006-06-01 01:00:03 +090014 u64 xm;
15 int xe;
16 int xs;
Linus Torvalds1da177e2005-04-16 15:20:36 -070017
Ralf Baechle9e8bad12014-04-19 00:36:32 +020018 ieee754_clearcx();
Linus Torvalds1da177e2005-04-16 15:20:36 -070019
Linus Torvalds1da177e2005-04-16 15:20:36 -070020 if (x == 0)
21 return ieee754dp_zero(0);
22 if (x == 1 || x == -1)
23 return ieee754dp_one(x < 0);
24 if (x == 10 || x == -10)
25 return ieee754dp_ten(x < 0);
26
27 xs = (x < 0);
28 if (xs) {
29 if (x == (1 << 31))
30 xm = ((unsigned) 1 << 31); /* max neg can't be safely negated */
31 else
32 xm = -x;
33 } else {
34 xm = x;
35 }
36
Linus Torvalds1da177e2005-04-16 15:20:36 -070037 /* normalize - result can never be inexact or overflow */
Ralf Baechlead8fb5532014-04-22 15:51:55 +020038 xe = DP_FBITS;
39 while ((xm >> DP_FBITS) == 0) {
Linus Torvalds1da177e2005-04-16 15:20:36 -070040 xm <<= 1;
41 xe--;
42 }
43 return builddp(xs, xe + DP_EBIAS, xm & ~DP_HIDDEN_BIT);
Linus Torvalds1da177e2005-04-16 15:20:36 -070044}