|
111
|
1 `/* Implementation of the NORM2 intrinsic
|
|
145
|
2 Copyright (C) 2010-2020 Free Software Foundation, Inc.
|
|
111
|
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
|
