diff options
Diffstat (limited to 'libgfortran/generated/rrspacing_r16.c')
-rw-r--r-- | libgfortran/generated/rrspacing_r16.c | 58 |
1 files changed, 58 insertions, 0 deletions
diff --git a/libgfortran/generated/rrspacing_r16.c b/libgfortran/generated/rrspacing_r16.c new file mode 100644 index 000000000..470541e72 --- /dev/null +++ b/libgfortran/generated/rrspacing_r16.c @@ -0,0 +1,58 @@ +/* Implementation of the RRSPACING intrinsic + Copyright 2006, 2007, 2009, 2010 Free Software Foundation, Inc. + Contributed by Steven G. Kargl <kargl@gcc.gnu.org> + +This file is part of the GNU Fortran 95 runtime library (libgfortran). + +Libgfortran is free software; you can redistribute it and/or +modify it under the terms of the GNU General Public +License as published by the Free Software Foundation; either +version 3 of the License, or (at your option) any later version. + +Libgfortran is distributed in the hope that it will be useful, +but WITHOUT ANY WARRANTY; without even the implied warranty of +MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +GNU General Public License for more details. + +Under Section 7 of GPL version 3, you are granted additional +permissions described in the GCC Runtime Library Exception, version +3.1, as published by the Free Software Foundation. + +You should have received a copy of the GNU General Public License and +a copy of the GCC Runtime Library Exception along with this program; +see the files COPYING3 and COPYING.RUNTIME respectively. If not, see +<http://www.gnu.org/licenses/>. */ + +#include "libgfortran.h" + + + +#if defined(GFC_REAL_16_IS_FLOAT128) +#define MATHFUNC(funcname) funcname ## q +#else +#define MATHFUNC(funcname) funcname ## l +#endif + +#if defined (HAVE_GFC_REAL_16) && (defined(GFC_REAL_16_IS_FLOAT128) || defined(HAVE_FABSL)) && (defined(GFC_REAL_16_IS_FLOAT128) || defined(HAVE_FREXPL)) + +extern GFC_REAL_16 rrspacing_r16 (GFC_REAL_16 s, int p); +export_proto(rrspacing_r16); + +GFC_REAL_16 +rrspacing_r16 (GFC_REAL_16 s, int p) +{ + int e; + GFC_REAL_16 x; + x = MATHFUNC(fabs) (s); + if (x == 0.) + return 0.; + MATHFUNC(frexp) (s, &e); +#if (defined(GFC_REAL_16_IS_FLOAT128) || defined(HAVE_LDEXPL)) + return MATHFUNC(ldexp) (x, p - e); +#else + return MATHFUNC(scalbn) (x, p - e); +#endif + +} + +#endif |