]>
Commit | Line | Data |
---|---|---|
4ee9c684 | 1 | /* Implementation of the NEAREST intrinsic |
f1717362 | 2 | Copyright (C) 2003-2016 Free Software Foundation, Inc. |
4ee9c684 | 3 | Contributed by Richard Henderson <rth@redhat.com>. |
4 | ||
b417ea8c | 5 | This file is part of the GNU Fortran 95 runtime library (libgfortran). |
4ee9c684 | 6 | |
b417ea8c | 7 | Libgfortran is free software; you can redistribute it and/or |
8 | modify it under the terms of the GNU General Public | |
4ee9c684 | 9 | License as published by the Free Software Foundation; either |
6bc9506f | 10 | version 3 of the License, or (at your option) any later version. |
b417ea8c | 11 | |
12 | Libgfortran is distributed in the hope that it will be useful, | |
4ee9c684 | 13 | but WITHOUT ANY WARRANTY; without even the implied warranty of |
14 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | |
b417ea8c | 15 | GNU General Public License for more details. |
4ee9c684 | 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/>. */ | |
920e54ef | 25 | |
4ee9c684 | 26 | #include "libgfortran.h" |
27 | ||
28 | ||
9dc59fdf | 29 | |
30 | #define MATHFUNC(funcname) funcname | |
31 | ||
920e54ef | 32 | #if defined (HAVE_GFC_REAL_8) && defined (HAVE_COPYSIGN) && defined (HAVE_NEXTAFTER) |
33 | ||
7b6cb5bd | 34 | extern GFC_REAL_8 nearest_r8 (GFC_REAL_8 s, GFC_REAL_8 dir); |
35 | export_proto(nearest_r8); | |
36 | ||
4ee9c684 | 37 | GFC_REAL_8 |
7b6cb5bd | 38 | nearest_r8 (GFC_REAL_8 s, GFC_REAL_8 dir) |
4ee9c684 | 39 | { |
9dc59fdf | 40 | dir = MATHFUNC(copysign) (MATHFUNC(__builtin_inf) (), dir); |
4ee9c684 | 41 | if (FLT_EVAL_METHOD != 0) |
42 | { | |
43 | /* ??? Work around glibc bug on x86. */ | |
9dc59fdf | 44 | volatile GFC_REAL_8 r = MATHFUNC(nextafter) (s, dir); |
4ee9c684 | 45 | return r; |
46 | } | |
47 | else | |
9dc59fdf | 48 | return MATHFUNC(nextafter) (s, dir); |
4ee9c684 | 49 | } |
920e54ef | 50 | |
51 | #endif |