]>
Commit | Line | Data |
---|---|---|
71d90ac4 | 1 | `/* Implementation of the RRSPACING intrinsic |
f1717362 | 2 | Copyright (C) 2006-2016 Free Software Foundation, Inc. |
71d90ac4 | 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 | |
6bc9506f | 10 | version 3 of the License, or (at your option) any later version. |
71d90ac4 | 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 | ||
6bc9506f | 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/>. */ | |
71d90ac4 | 25 | |
71d90ac4 | 26 | #include "libgfortran.h"' |
27 | ||
28 | include(`mtype.m4')dnl | |
29 | ||
9dc59fdf | 30 | mathfunc_macro |
31 | ||
32 | `#if defined (HAVE_'real_type`) && 'hasmathfunc(fabs) && hasmathfunc(frexp)` | |
71d90ac4 | 33 | |
0a6b5f6b | 34 | extern 'real_type` rrspacing_r'kind` ('real_type` s, int p); |
35 | export_proto(rrspacing_r'kind`); | |
71d90ac4 | 36 | |
0a6b5f6b | 37 | 'real_type` |
38 | rrspacing_r'kind` ('real_type` s, int p) | |
71d90ac4 | 39 | { |
40 | int e; | |
0a6b5f6b | 41 | 'real_type` x; |
9dc59fdf | 42 | x = MATHFUNC(fabs) (s); |
71d90ac4 | 43 | if (x == 0.) |
44 | return 0.; | |
9dc59fdf | 45 | MATHFUNC(frexp) (s, &e); |
46 | #if 'hasmathfunc(ldexp)` | |
47 | return MATHFUNC(ldexp) (x, p - e); | |
bd73b228 | 48 | #else |
9dc59fdf | 49 | return MATHFUNC(scalbn) (x, p - e); |
bd73b228 | 50 | #endif |
51 | ||
71d90ac4 | 52 | } |
53 | ||
0a6b5f6b | 54 | #endif' |