1 /* Software floating-point emulation.
2 frexpl(x, exp)
3 Copyright (C) 1999-2017 Free Software Foundation, Inc.
4 This file is part of the GNU C Library.
5 Contributed by Jakub Jelinek (jj@ultra.linux.cz).
6
7 The GNU C Library is free software; you can redistribute it and/or
8 modify it under the terms of the GNU Lesser General Public
9 License as published by the Free Software Foundation; either
10 version 2.1 of the License, or (at your option) any later version.
11
12 The GNU C Library is distributed in the hope that it will be useful,
13 but WITHOUT ANY WARRANTY; without even the implied warranty of
14 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 Lesser General Public License for more details.
16
17 You should have received a copy of the GNU Lesser General Public
18 License along with the GNU C Library; if not, see
19 <http://www.gnu.org/licenses/>. */
20
21 /*
22 * for non-zero x
23 * x = frexpl(arg,&exp);
24 * return a long double fp quantity x such that 0.5 <= |x| <1.0
25 * and the corresponding binary exponent "exp". That is
26 * arg = x*2^exp.
27 * If arg is inf, 0.0, or NaN, then frexpl(arg,&exp) returns arg
28 * with *exp=0.
29 */
30
31 #include "soft-fp.h"
32 #include "quad.h"
33
__frexpl(long double arg,int * exp)34 long double __frexpl(long double arg, int *exp)
35 {
36 FP_DECL_EX;
37 FP_DECL_Q(A);
38 long double r;
39
40 *exp = 0;
41 FP_UNPACK_Q(A, arg);
42 if (A_c != FP_CLS_NORMAL)
43 return arg;
44 *exp = A_e + 1;
45 A_e = -1;
46 FP_PACK_Q(r, A);
47
48 return r;
49 }
50
51 weak_alias (__frexpl, frexpl)
52