1`/* Implementation of the SPACING intrinsic 2 Copyright (C) 2006-2016 Free Software Foundation, Inc. 3 Contributed by Steven G. Kargl <kargl@gcc.gnu.org> 4 5This file is part of the GNU Fortran 95 runtime library (libgfortran). 6 7Libgfortran is free software; you can redistribute it and/or 8modify it under the terms of the GNU General Public 9License as published by the Free Software Foundation; either 10version 3 of the License, or (at your option) any later version. 11 12Libgfortran is distributed in the hope that it will be useful, 13but WITHOUT ANY WARRANTY; without even the implied warranty of 14MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 15GNU General Public License for more details. 16 17Under Section 7 of GPL version 3, you are granted additional 18permissions described in the GCC Runtime Library Exception, version 193.1, as published by the Free Software Foundation. 20 21You should have received a copy of the GNU General Public License and 22a copy of the GCC Runtime Library Exception along with this program; 23see the files COPYING3 and COPYING.RUNTIME respectively. If not, see 24<http://www.gnu.org/licenses/>. */ 25 26#include "libgfortran.h"' 27 28include(`mtype.m4')dnl 29 30mathfunc_macro 31 32`#if defined (HAVE_'real_type`) && 'hasmathfunc(frexp)` 33 34extern 'real_type` spacing_r'kind` ('real_type` s, int p, int emin, 'real_type` tiny); 35export_proto(spacing_r'kind`); 36 37'real_type` 38spacing_r'kind` ('real_type` s, int p, int emin, 'real_type` tiny) 39{ 40 int e; 41 if (s == 0.) 42 return tiny; 43 MATHFUNC(frexp) (s, &e); 44 e = e - p; 45 e = e > emin ? e : emin; 46#if 'hasmathfunc(ldexp)` 47 return MATHFUNC(ldexp) (1., e); 48#else 49 return MATHFUNC(scalbn) (1., e); 50#endif 51} 52 53#endif' 54