1627f7eb2Smrg`/* Implementation of the NORM2 intrinsic 2*4c3eb207Smrg Copyright (C) 2010-2020 Free Software Foundation, Inc. 3627f7eb2Smrg Contributed by Tobias Burnus <burnus@net-b.de> 4627f7eb2Smrg 5627f7eb2SmrgThis file is part of the GNU Fortran runtime library (libgfortran). 6627f7eb2Smrg 7627f7eb2SmrgLibgfortran is free software; you can redistribute it and/or 8627f7eb2Smrgmodify it under the terms of the GNU General Public 9627f7eb2SmrgLicense as published by the Free Software Foundation; either 10627f7eb2Smrgversion 3 of the License, or (at your option) any later version. 11627f7eb2Smrg 12627f7eb2SmrgLibgfortran is distributed in the hope that it will be useful, 13627f7eb2Smrgbut WITHOUT ANY WARRANTY; without even the implied warranty of 14627f7eb2SmrgMERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 15627f7eb2SmrgGNU General Public License for more details. 16627f7eb2Smrg 17627f7eb2SmrgUnder Section 7 of GPL version 3, you are granted additional 18627f7eb2Smrgpermissions described in the GCC Runtime Library Exception, version 19627f7eb2Smrg3.1, as published by the Free Software Foundation. 20627f7eb2Smrg 21627f7eb2SmrgYou should have received a copy of the GNU General Public License and 22627f7eb2Smrga copy of the GCC Runtime Library Exception along with this program; 23627f7eb2Smrgsee the files COPYING3 and COPYING.RUNTIME respectively. If not, see 24627f7eb2Smrg<http://www.gnu.org/licenses/>. */ 25627f7eb2Smrg 26627f7eb2Smrg#include "libgfortran.h"' 27627f7eb2Smrg 28627f7eb2Smrginclude(iparm.m4)dnl 29627f7eb2Smrginclude(ifunction.m4)dnl 30627f7eb2Smrginclude(`mtype.m4')dnl 31627f7eb2Smrg 32627f7eb2Smrg`#if defined (HAVE_'atype_name`) && defined (HAVE_'rtype_name`) && 'hasmathfunc(sqrt) && hasmathfunc(fabs) 33627f7eb2Smrg 34627f7eb2Smrgmathfunc_macro 35627f7eb2Smrg 36627f7eb2SmrgARRAY_FUNCTION(0, 37627f7eb2Smrg` 'rtype_name` scale; 38627f7eb2Smrg result = 0; 39627f7eb2Smrg scale = 1;', 40627f7eb2Smrg` if (*src != 0) 41627f7eb2Smrg { 42627f7eb2Smrg 'rtype_name` absX, val; 43627f7eb2Smrg absX = MATHFUNC(fabs) (*src); 44627f7eb2Smrg if (scale < absX) 45627f7eb2Smrg { 46627f7eb2Smrg val = scale / absX; 47627f7eb2Smrg result = 1 + result * val * val; 48627f7eb2Smrg scale = absX; 49627f7eb2Smrg } 50627f7eb2Smrg else 51627f7eb2Smrg { 52627f7eb2Smrg val = absX / scale; 53627f7eb2Smrg result += val * val; 54627f7eb2Smrg } 55627f7eb2Smrg }', 56627f7eb2Smrg` result = scale * MATHFUNC(sqrt) (result);') 57627f7eb2Smrg 58627f7eb2Smrg#endif 59