BLAS: inline SCABS1/DCABS1 into CAXPY/ZAXPY as statement function CABS1

Replace external function references to SCABS1 and DCABS1 with a local
statement function CABS1 in CAXPY and ZAXPY respectively, eliminating
the dependency on those auxiliary routines.

Author: nakatamaho

BLAS/SRC/caxpy.f

@@ -103,13 +103,16 @@ SUBROUTINE CAXPY(N,CA,CX,INCX,CY,INCY)
 *
 *     .. Local Scalars ..
       INTEGER I,IX,IY
+      COMPLEX CDUM
 *     ..
-*     .. External Functions ..
-      REAL SCABS1
-      EXTERNAL SCABS1
+*     .. Statement Functions ..
+      REAL CABS1
+*     ..
+*     .. Statement Function definitions ..
+      CABS1(CDUM) = ABS(REAL(CDUM)) + ABS(AIMAG(CDUM))
 *     ..
       IF (N.LE.0) RETURN
-      IF (SCABS1(CA).EQ.0.0E+0) RETURN
+      IF (CABS1(CA).EQ.0.0E+0) RETURN
       IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN
 *
 *        code for both increments equal to 1

BLAS/SRC/zaxpy.f

@@ -103,13 +103,16 @@ SUBROUTINE ZAXPY(N,ZA,ZX,INCX,ZY,INCY)
 *
 *     .. Local Scalars ..
       INTEGER I,IX,IY
+      COMPLEX*16 ZDUM
 *     ..
-*     .. External Functions ..
-      DOUBLE PRECISION DCABS1
-      EXTERNAL DCABS1
+*     .. Statement Functions ..
+      DOUBLE PRECISION CABS1
+*     ..
+*     .. Statement Function definitions ..
+      CABS1(ZDUM) = ABS(DBLE(ZDUM)) + ABS(DIMAG(ZDUM))
 *     ..
       IF (N.LE.0) RETURN
-      IF (DCABS1(ZA).EQ.0.0d0) RETURN
+      IF (CABS1(ZA).EQ.0.0D+0) RETURN
       IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN
 *
 *        code for both increments equal to 1