1  /* Copyright (C) 1998-2023 Free Software Foundation, Inc.
       2     This file is part of the GNU C Library.
       3  
       4     The GNU C Library is free software; you can redistribute it and/or
       5     modify it under the terms of the GNU Lesser General Public
       6     License as published by the Free Software Foundation; either
       7     version 2.1 of the License, or (at your option) any later version.
       8  
       9     The GNU C Library is distributed in the hope that it will be useful,
      10     but WITHOUT ANY WARRANTY; without even the implied warranty of
      11     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
      12     Lesser General Public License for more details.
      13  
      14     You should have received a copy of the GNU Lesser General Public
      15     License along with the GNU C Library; if not, see
      16     <https://www.gnu.org/licenses/>.  */
      17  
      18  #include <math.h>
      19  #include <math_private.h>
      20  #include <float.h>
      21  #include <libm-alias-finite.h>
      22  
      23  double
      24  __ieee754_exp10 (double arg)
      25  {
      26    if (isfinite (arg) && arg < DBL_MIN_10_EXP - DBL_DIG - 10)
      27      return DBL_MIN * DBL_MIN;
      28    else
      29      /* This is a very stupid and inprecise implementation.  It'll get
      30         replaced sometime (soon?).  */
      31      return __ieee754_exp (M_LN10 * arg);
      32  }
      33  libm_alias_finite (__ieee754_exp10, __exp10)