Mercurial > hg > CbC > CbC_gcc
comparison libgfortran/m4/norm2.m4 @ 111:04ced10e8804
gcc 7
author | kono |
---|---|
date | Fri, 27 Oct 2017 22:46:09 +0900 |
parents | |
children | 84e7813d76e9 |
comparison
equal
deleted
inserted
replaced
68:561a7518be6b | 111:04ced10e8804 |
---|---|
1 `/* Implementation of the NORM2 intrinsic | |
2 Copyright (C) 2010-2017 Free Software Foundation, Inc. | |
3 Contributed by Tobias Burnus <burnus@net-b.de> | |
4 | |
5 This file is part of the GNU Fortran 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 | |
10 version 3 of the License, or (at your option) any later version. | |
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 | |
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/>. */ | |
25 | |
26 #include "libgfortran.h"' | |
27 | |
28 include(iparm.m4)dnl | |
29 include(ifunction.m4)dnl | |
30 include(`mtype.m4')dnl | |
31 | |
32 `#if defined (HAVE_'atype_name`) && defined (HAVE_'rtype_name`) && 'hasmathfunc(sqrt) && hasmathfunc(fabs) | |
33 | |
34 mathfunc_macro | |
35 | |
36 ARRAY_FUNCTION(0, | |
37 ` 'rtype_name` scale; | |
38 result = 0; | |
39 scale = 1;', | |
40 ` if (*src != 0) | |
41 { | |
42 'rtype_name` absX, val; | |
43 absX = MATHFUNC(fabs) (*src); | |
44 if (scale < absX) | |
45 { | |
46 val = scale / absX; | |
47 result = 1 + result * val * val; | |
48 scale = absX; | |
49 } | |
50 else | |
51 { | |
52 val = absX / scale; | |
53 result += val * val; | |
54 } | |
55 }', | |
56 ` result = scale * MATHFUNC(sqrt) (result);') | |
57 | |
58 #endif |