lh | 9ed821d | 2023-04-07 01:36:19 -0700 | [diff] [blame^] | 1 | /* Copyright (C) 1998-2015 Free Software Foundation, Inc. |
| 2 | This file is part of the GNU C Library. |
| 3 | Contributed by Ulrich Drepper <drepper@cygnus.com>, 1998. |
| 4 | |
| 5 | The GNU C Library is free software; you can redistribute it and/or |
| 6 | modify it under the terms of the GNU Lesser General Public |
| 7 | License as published by the Free Software Foundation; either |
| 8 | version 2.1 of the License, or (at your option) any later version. |
| 9 | |
| 10 | The GNU C Library is distributed in the hope that it will be useful, |
| 11 | but WITHOUT ANY WARRANTY; without even the implied warranty of |
| 12 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
| 13 | Lesser General Public License for more details. |
| 14 | |
| 15 | You should have received a copy of the GNU Lesser General Public |
| 16 | License along with the GNU C Library; if not, see |
| 17 | <http://www.gnu.org/licenses/>. */ |
| 18 | |
| 19 | #include <math.h> |
| 20 | #include <math_private.h> |
| 21 | #include <float.h> |
| 22 | |
| 23 | long double |
| 24 | __ieee754_exp10l (long double arg) |
| 25 | { |
| 26 | if (isfinite (arg) && arg < LDBL_MIN_10_EXP - LDBL_DIG - 10) |
| 27 | return LDBL_MIN * LDBL_MIN; |
| 28 | else |
| 29 | /* This is a very stupid and inprecise implementation. It'll get |
| 30 | replaced sometime (soon?). */ |
| 31 | return __ieee754_expl (M_LN10l * arg); |
| 32 | } |
| 33 | strong_alias (__ieee754_exp10l, __exp10l_finite) |