Update copyright years.

[thirdparty/gcc.git] / libgfortran / generated / norm2_r10.c

diff --git a/libgfortran/generated/norm2_r10.c b/libgfortran/generated/norm2_r10.c
index 64cf404e0a653d094183fa4c01d9e0ff7e836c64..4079df263fd4c8f0f2fd4ffe6b6cbcc4909d2c1a 100644 (file)
--- a/libgfortran/generated/norm2_r10.c
+++ b/libgfortran/generated/norm2_r10.c
@@ -1,5 +1,5 @@
 /* Implementation of the NORM2 intrinsic
-   Copyright 2010 Free Software Foundation, Inc.
+   Copyright (C) 2010-2020 Free Software Foundation, Inc.
    Contributed by Tobias Burnus  <burnus@net-b.de>
 
 This file is part of the GNU Fortran runtime library (libgfortran).
@@ -24,12 +24,12 @@ 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>
 
 
-#if defined (HAVE_GFC_REAL_10) && defined (HAVE_GFC_REAL_10)
+
+#if defined (HAVE_GFC_REAL_10) && defined (HAVE_GFC_REAL_10) && defined (HAVE_SQRTL) && defined (HAVE_FABSL)
+
+#define MATHFUNC(funcname) funcname ## l
 
 
 extern void norm2_r10 (gfc_array_r10 * const restrict, 
@@ -55,8 +55,15 @@ norm2_r10 (gfc_array_r10 * const restrict retarray,
   int continue_loop;
 
   /* Make dim zero based to avoid confusion.  */
-  dim = (*pdim) - 1;
   rank = GFC_DESCRIPTOR_RANK (array) - 1;
+  dim = (*pdim) - 1;
+
+  if (unlikely (dim < 0 || dim > rank))
+    {
+      runtime_error ("Dim argument incorrect in NORM intrinsic: "
+                    "is %ld, should be between 1 and %ld",
+                    (long int) dim + 1, (long int) rank + 1);
+    }
 
   len = GFC_DESCRIPTOR_EXTENT(array,dim);
   if (len < 0)
@@ -80,7 +87,7 @@ norm2_r10 (gfc_array_r10 * const restrict retarray,
        extent[n] = 0;
     }
 
-  if (retarray->data == NULL)
+  if (retarray->base_addr == NULL)
     {
       size_t alloc_size, str;
 
@@ -96,11 +103,11 @@ norm2_r10 (gfc_array_r10 * const restrict retarray,
        }
 
       retarray->offset = 0;
-      retarray->dtype = (array->dtype & ~GFC_DTYPE_RANK_MASK) | rank;
+      retarray->dtype.rank = rank;
 
-      alloc_size = sizeof (GFC_REAL_10) * GFC_DESCRIPTOR_STRIDE(retarray,rank-1)
-                  * extent[rank-1];
+      alloc_size = GFC_DESCRIPTOR_STRIDE(retarray,rank-1) * extent[rank-1];
 
+      retarray->base_addr = xmallocarray (alloc_size, sizeof (GFC_REAL_10));
       if (alloc_size == 0)
        {
          /* Make sure we have a zero-sized array.  */
@@ -108,8 +115,6 @@ norm2_r10 (gfc_array_r10 * const restrict retarray,
          return;
 
        }
-      else
-       retarray->data = internal_malloc_size (alloc_size);
     }
   else
     {
@@ -129,11 +134,11 @@ norm2_r10 (gfc_array_r10 * const restrict retarray,
       count[n] = 0;
       dstride[n] = GFC_DESCRIPTOR_STRIDE(retarray,n);
       if (extent[n] <= 0)
-       len = 0;
+       return;
     }
 
-  base = array->data;
-  dest = retarray->data;
+  base = array->base_addr;
+  dest = retarray->base_addr;
 
   continue_loop = 1;
   while (continue_loop)
@@ -144,23 +149,25 @@ norm2_r10 (gfc_array_r10 * const restrict retarray,
       {
 
        GFC_REAL_10 scale;
-       result = 0.0L;
-       scale = 1.0L;
+       result = 0;
+       scale = 1;
        if (len <= 0)
-         *dest = 0.0L;
+         *dest = 0;
        else
          {
+#if ! defined HAVE_BACK_ARG
            for (n = 0; n < len; n++, src += delta)
              {
+#endif
 
-         if (*src != 0.0L)
+         if (*src != 0)
            {
              GFC_REAL_10 absX, val;
-             absX = fabsl (*src);
+             absX = MATHFUNC(fabs) (*src);
              if (scale < absX)
                {
                  val = scale / absX;
-                 result = 1.0L + result * val * val;
+                 result = 1 + result * val * val;
                  scale = absX;
                }
              else
@@ -170,7 +177,7 @@ norm2_r10 (gfc_array_r10 * const restrict retarray,
                }
            }
              }
-           result = scale * sqrtl (result);
+           result = scale * MATHFUNC(sqrt) (result);
            *dest = result;
          }
       }
@@ -189,9 +196,9 @@ norm2_r10 (gfc_array_r10 * const restrict retarray,
          base -= sstride[n] * extent[n];
          dest -= dstride[n] * extent[n];
          n++;
-         if (n == rank)
+         if (n >= rank)
            {
-             /* Break out of the look.  */
+             /* Break out of the loop.  */
              continue_loop = 0;
              break;
            }