]>
Commit | Line | Data |
---|---|---|
cc6d3bde | 1 | /* Implementation of the SPACING intrinsic |
748086b7 | 2 | Copyright 2006, 2007, 2009 Free Software Foundation, Inc. |
cc6d3bde SK |
3 | Contributed by Steven G. Kargl <kargl@gcc.gnu.org> |
4 | ||
5 | This file is part of the GNU Fortran 95 runtime library (libgfortran). | |
6 | ||
7 | Libgfortran is free software; you can redistribute it and/or | |
8 | modify it under the terms of the GNU General Public | |
9 | License as published by the Free Software Foundation; either | |
748086b7 | 10 | version 3 of the License, or (at your option) any later version. |
cc6d3bde SK |
11 | |
12 | Libgfortran 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 | |
15 | GNU General Public License for more details. | |
16 | ||
748086b7 JJ |
17 | Under Section 7 of GPL version 3, you are granted additional |
18 | permissions described in the GCC Runtime Library Exception, version | |
19 | 3.1, as published by the Free Software Foundation. | |
20 | ||
21 | You should have received a copy of the GNU General Public License and | |
22 | a copy of the GCC Runtime Library Exception along with this program; | |
23 | see the files COPYING3 and COPYING.RUNTIME respectively. If not, see | |
24 | <http://www.gnu.org/licenses/>. */ | |
cc6d3bde | 25 | |
cc6d3bde SK |
26 | #include "libgfortran.h" |
27 | ||
28 | ||
08fd13d4 FXC |
29 | |
30 | #if defined(GFC_REAL_16_IS_FLOAT128) | |
31 | #define MATHFUNC(funcname) funcname ## q | |
32 | #else | |
33 | #define MATHFUNC(funcname) funcname ## l | |
34 | #endif | |
35 | ||
36 | #if defined (HAVE_GFC_REAL_16) && (defined(GFC_WITH_QUAD_LIB) || defined(HAVE_FREXPL)) | |
cc6d3bde SK |
37 | |
38 | extern GFC_REAL_16 spacing_r16 (GFC_REAL_16 s, int p, int emin, GFC_REAL_16 tiny); | |
39 | export_proto(spacing_r16); | |
40 | ||
41 | GFC_REAL_16 | |
42 | spacing_r16 (GFC_REAL_16 s, int p, int emin, GFC_REAL_16 tiny) | |
43 | { | |
44 | int e; | |
45 | if (s == 0.) | |
46 | return tiny; | |
08fd13d4 | 47 | MATHFUNC(frexp) (s, &e); |
cc6d3bde SK |
48 | e = e - p; |
49 | e = e > emin ? e : emin; | |
08fd13d4 FXC |
50 | #if (defined(GFC_WITH_QUAD_LIB) || defined(HAVE_LDEXPL)) |
51 | return MATHFUNC(ldexp) (1., e); | |
4010308c | 52 | #else |
08fd13d4 | 53 | return MATHFUNC(scalbn) (1., e); |
4010308c | 54 | #endif |
cc6d3bde SK |
55 | } |
56 | ||
57 | #endif |