Fix error if alpha=0

grisuthedragon · web-flow · commit f21a6aefd297 · 2024-08-28T22:46:32.000+02:00
diff --git a/BLAS/SRC/caxpby.f b/BLAS/SRC/caxpby.f
@@ -113,13 +113,8 @@ SUBROUTINE CAXPBY(N,CA,CX,INCX,CB,CY,INCY)
 *
 *     .. Local Scalars ..
       INTEGER I,IX,IY
-*     ..
-*     .. External Functions ..
-      REAL SCABS1
-      EXTERNAL SCABS1
 *     ..
       IF (N.LE.0) RETURN
-      IF (SCABS1(CA).EQ.0.0E+0) RETURN
       IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN
 *
 *        code for both increments equal to 1
diff --git a/BLAS/SRC/daxpby.f b/BLAS/SRC/daxpby.f
@@ -119,7 +119,6 @@ SUBROUTINE DAXPBY(N,DA,DX,INCX,DB,DY,INCY)
       INTRINSIC MOD
 *     ..
       IF (N.LE.0) RETURN
-      IF (DA.EQ.0.0d0) RETURN
       IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN
 *
 *        code for both increments equal to 1
diff --git a/BLAS/SRC/saxpby.f b/BLAS/SRC/saxpby.f
@@ -119,7 +119,6 @@ SUBROUTINE SAXPBY(N,SA,SX,INCX,SB,SY,INCY)
       INTRINSIC MOD
 *     ..
       IF (N.LE.0) RETURN
-      IF (SA.EQ.0.0) RETURN
       IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN
 *
 *        code for both increments equal to 1
diff --git a/BLAS/SRC/zaxpby.f b/BLAS/SRC/zaxpby.f
@@ -113,13 +113,8 @@ SUBROUTINE ZAXPBY(N,ZA,ZX,INCX,ZB,ZY,INCY)
 *
 *     .. Local Scalars ..
       INTEGER I,IX,IY
-*     ..
-*     .. External Functions ..
-      DOUBLE PRECISION DCABS1
-      EXTERNAL DCABS1
 *     ..
       IF (N.LE.0) RETURN
-      IF (DCABS1(ZA).EQ.0.0d0) RETURN
       IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN
 *
 *        code for both increments equal to 1

Original file line number	Diff line number	Diff line change
`@@ -113,13 +113,8 @@ SUBROUTINE CAXPBY(N,CA,CX,INCX,CB,CY,INCY)`
`113`	`113`	`*`
`114`	`114`	`* .. Local Scalars ..`
`115`	`115`	`INTEGER I,IX,IY`
`116`		`-* ..`
`117`		`-* .. External Functions ..`
`118`		`- REAL SCABS1`
`119`		`- EXTERNAL SCABS1`
`120`	`116`	`* ..`
`121`	`117`	`IF (N.LE.0) RETURN`
`122`		`- IF (SCABS1(CA).EQ.0.0E+0) RETURN`
`123`	`118`	`IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN`
`124`	`119`	`*`
`125`	`120`	`* code for both increments equal to 1`
Original file line number	Diff line number	Diff line change
`@@ -119,7 +119,6 @@ SUBROUTINE DAXPBY(N,DA,DX,INCX,DB,DY,INCY)`
`119`	`119`	`INTRINSIC MOD`
`120`	`120`	`* ..`
`121`	`121`	`IF (N.LE.0) RETURN`
`122`		`- IF (DA.EQ.0.0d0) RETURN`
`123`	`122`	`IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN`
`124`	`123`	`*`
`125`	`124`	`* code for both increments equal to 1`
Original file line number	Diff line number	Diff line change
`@@ -119,7 +119,6 @@ SUBROUTINE SAXPBY(N,SA,SX,INCX,SB,SY,INCY)`
`119`	`119`	`INTRINSIC MOD`
`120`	`120`	`* ..`
`121`	`121`	`IF (N.LE.0) RETURN`
`122`		`- IF (SA.EQ.0.0) RETURN`
`123`	`122`	`IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN`
`124`	`123`	`*`
`125`	`124`	`* code for both increments equal to 1`
Original file line number	Diff line number	Diff line change
`@@ -113,13 +113,8 @@ SUBROUTINE ZAXPBY(N,ZA,ZX,INCX,ZB,ZY,INCY)`
`113`	`113`	`*`
`114`	`114`	`* .. Local Scalars ..`
`115`	`115`	`INTEGER I,IX,IY`
`116`		`-* ..`
`117`		`-* .. External Functions ..`
`118`		`- DOUBLE PRECISION DCABS1`
`119`		`- EXTERNAL DCABS1`
`120`	`116`	`* ..`
`121`	`117`	`IF (N.LE.0) RETURN`
`122`		`- IF (DCABS1(ZA).EQ.0.0d0) RETURN`
`123`	`118`	`IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN`
`124`	`119`	`*`
`125`	`120`	`* code for both increments equal to 1`