URL
https://opencores.org/ocsvn/openrisc/openrisc/trunk
Subversion Repositories openrisc
[/] [openrisc/] [trunk/] [gnu-dev/] [or1k-gcc/] [libgfortran/] [m4/] [norm2.m4] - Rev 734
Go to most recent revision | Compare with Previous | Blame | View Log
`/* Implementation of the NORM2 intrinsic Copyright 2010 Free Software Foundation, Inc. Contributed by Tobias Burnus <burnus@net-b.de> This file is part of the GNU Fortran runtime library (libgfortran). Libgfortran is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 3 of the License, or (at your option) any later version. Libgfortran is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. Under Section 7 of GPL version 3, you are granted additional permissions described in the GCC Runtime Library Exception, version 3.1, as published by the Free Software Foundation. You should have received a copy of the GNU General Public License and a copy of the GCC Runtime Library Exception along with this program; see the files COPYING3 and COPYING.RUNTIME respectively. If not, see <http://www.gnu.org/licenses/>. */ #include "libgfortran.h" #include <stdlib.h> #include <math.h> #include <assert.h>' include(iparm.m4)dnl include(ifunction.m4)dnl include(`mtype.m4')dnl `#if defined (HAVE_'atype_name`) && defined (HAVE_'rtype_name`) && 'hasmathfunc(sqrt) && hasmathfunc(fabs) mathfunc_macro ARRAY_FUNCTION(0, ` 'rtype_name` scale; result = 0; scale = 1;', ` if (*src != 0) { 'rtype_name` absX, val; absX = MATHFUNC(fabs) (*src); if (scale < absX) { val = scale / absX; result = 1 + result * val * val; scale = absX; } else { val = absX / scale; result += val * val; } }', ` result = scale * MATHFUNC(sqrt) (result);') #endif
Go to most recent revision | Compare with Previous | Blame | View Log