diff --git a/BLAS/SRC/ctrsm.f b/BLAS/SRC/ctrsm.f
index 6275d4a59..39dececbb 100644
--- a/BLAS/SRC/ctrsm.f
+++ b/BLAS/SRC/ctrsm.f
@@ -283,29 +283,25 @@ SUBROUTINE CTRSM(SIDE,UPLO,TRANSA,DIAG,M,N,ALPHA,A,LDA,B,LDB)
                               B(I,J) = ALPHA*B(I,J)
    30                     CONTINUE
                       END IF
-                      DO 50 K = M,1,-1
-                          IF (B(K,J).NE.ZERO) THEN
-                              IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
-                              DO 40 I = 1,K - 1
-                                  B(I,J) = B(I,J) - B(K,J)*A(I,K)
-   40                         CONTINUE
-                          END IF
+                       DO 50 K = M,1,-1
+                           IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
+                           DO 40 I = 1,K - 1
+                               B(I,J) = B(I,J) - B(K,J)*A(I,K)
+   40                     CONTINUE
    50                 CONTINUE
    60             CONTINUE
-              ELSE
-                  DO 100 J = 1,N
-                      IF (ALPHA.NE.ONE) THEN
-                          DO 70 I = 1,M
-                              B(I,J) = ALPHA*B(I,J)
+               ELSE
+                   DO 100 J = 1,N
+                       IF (ALPHA.NE.ONE) THEN
+                           DO 70 I = 1,M
+                               B(I,J) = ALPHA*B(I,J)
    70                     CONTINUE
-                      END IF
-                      DO 90 K = 1,M
-                          IF (B(K,J).NE.ZERO) THEN
-                              IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
-                              DO 80 I = K + 1,M
-                                  B(I,J) = B(I,J) - B(K,J)*A(I,K)
-   80                         CONTINUE
-                          END IF
+                       END IF
+                       DO 90 K = 1,M
+                           IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
+                           DO 80 I = K + 1,M
+                               B(I,J) = B(I,J) - B(K,J)*A(I,K)
+   80                     CONTINUE
    90                 CONTINUE
   100             CONTINUE
               END IF
diff --git a/BLAS/SRC/dtrsm.f b/BLAS/SRC/dtrsm.f
index 0cd4a91a2..4fb64bf66 100644
--- a/BLAS/SRC/dtrsm.f
+++ b/BLAS/SRC/dtrsm.f
@@ -281,13 +281,11 @@ SUBROUTINE DTRSM(SIDE,UPLO,TRANSA,DIAG,M,N,ALPHA,A,LDA,B,LDB)
                               B(I,J) = ALPHA*B(I,J)
    30                     CONTINUE
                       END IF
-                      DO 50 K = M,1,-1
-                          IF (B(K,J).NE.ZERO) THEN
-                              IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
-                              DO 40 I = 1,K - 1
-                                  B(I,J) = B(I,J) - B(K,J)*A(I,K)
-   40                         CONTINUE
-                          END IF
+                       DO 50 K = M,1,-1
+                           IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
+                           DO 40 I = 1,K - 1
+                               B(I,J) = B(I,J) - B(K,J)*A(I,K)
+   40                     CONTINUE
    50                 CONTINUE
    60             CONTINUE
               ELSE
@@ -297,13 +295,11 @@ SUBROUTINE DTRSM(SIDE,UPLO,TRANSA,DIAG,M,N,ALPHA,A,LDA,B,LDB)
                               B(I,J) = ALPHA*B(I,J)
    70                     CONTINUE
                       END IF
-                      DO 90 K = 1,M
-                          IF (B(K,J).NE.ZERO) THEN
-                              IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
-                              DO 80 I = K + 1,M
-                                  B(I,J) = B(I,J) - B(K,J)*A(I,K)
-   80                         CONTINUE
-                          END IF
+                       DO 90 K = 1,M
+                           IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
+                           DO 80 I = K + 1,M
+                               B(I,J) = B(I,J) - B(K,J)*A(I,K)
+   80                     CONTINUE
    90                 CONTINUE
   100             CONTINUE
               END IF
diff --git a/BLAS/SRC/strsm.f b/BLAS/SRC/strsm.f
index 45fbb315b..495ee6fa3 100644
--- a/BLAS/SRC/strsm.f
+++ b/BLAS/SRC/strsm.f
@@ -281,29 +281,25 @@ SUBROUTINE STRSM(SIDE,UPLO,TRANSA,DIAG,M,N,ALPHA,A,LDA,B,LDB)
                               B(I,J) = ALPHA*B(I,J)
    30                     CONTINUE
                       END IF
-                      DO 50 K = M,1,-1
-                          IF (B(K,J).NE.ZERO) THEN
-                              IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
-                              DO 40 I = 1,K - 1
-                                  B(I,J) = B(I,J) - B(K,J)*A(I,K)
-   40                         CONTINUE
-                          END IF
+                       DO 50 K = M,1,-1
+                           IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
+                           DO 40 I = 1,K - 1
+                               B(I,J) = B(I,J) - B(K,J)*A(I,K)
+   40                     CONTINUE
    50                 CONTINUE
    60             CONTINUE
-              ELSE
-                  DO 100 J = 1,N
-                      IF (ALPHA.NE.ONE) THEN
-                          DO 70 I = 1,M
-                              B(I,J) = ALPHA*B(I,J)
+               ELSE
+                   DO 100 J = 1,N
+                       IF (ALPHA.NE.ONE) THEN
+                           DO 70 I = 1,M
+                               B(I,J) = ALPHA*B(I,J)
    70                     CONTINUE
-                      END IF
-                      DO 90 K = 1,M
-                          IF (B(K,J).NE.ZERO) THEN
-                              IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
-                              DO 80 I = K + 1,M
-                                  B(I,J) = B(I,J) - B(K,J)*A(I,K)
-   80                         CONTINUE
-                          END IF
+                       END IF
+                       DO 90 K = 1,M
+                           IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
+                           DO 80 I = K + 1,M
+                               B(I,J) = B(I,J) - B(K,J)*A(I,K)
+   80                     CONTINUE
    90                 CONTINUE
   100             CONTINUE
               END IF
diff --git a/BLAS/SRC/ztrsm.f b/BLAS/SRC/ztrsm.f
index 192190f6c..241fbcb30 100644
--- a/BLAS/SRC/ztrsm.f
+++ b/BLAS/SRC/ztrsm.f
@@ -283,29 +283,25 @@ SUBROUTINE ZTRSM(SIDE,UPLO,TRANSA,DIAG,M,N,ALPHA,A,LDA,B,LDB)
                               B(I,J) = ALPHA*B(I,J)
    30                     CONTINUE
                       END IF
-                      DO 50 K = M,1,-1
-                          IF (B(K,J).NE.ZERO) THEN
-                              IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
-                              DO 40 I = 1,K - 1
-                                  B(I,J) = B(I,J) - B(K,J)*A(I,K)
-   40                         CONTINUE
-                          END IF
+                       DO 50 K = M,1,-1
+                           IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
+                           DO 40 I = 1,K - 1
+                               B(I,J) = B(I,J) - B(K,J)*A(I,K)
+   40                     CONTINUE
    50                 CONTINUE
    60             CONTINUE
-              ELSE
-                  DO 100 J = 1,N
-                      IF (ALPHA.NE.ONE) THEN
-                          DO 70 I = 1,M
-                              B(I,J) = ALPHA*B(I,J)
+               ELSE
+                   DO 100 J = 1,N
+                       IF (ALPHA.NE.ONE) THEN
+                           DO 70 I = 1,M
+                               B(I,J) = ALPHA*B(I,J)
    70                     CONTINUE
-                      END IF
-                      DO 90 K = 1,M
-                          IF (B(K,J).NE.ZERO) THEN
-                              IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
-                              DO 80 I = K + 1,M
-                                  B(I,J) = B(I,J) - B(K,J)*A(I,K)
-   80                         CONTINUE
-                          END IF
+                       END IF
+                       DO 90 K = 1,M
+                           IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
+                           DO 80 I = K + 1,M
+                               B(I,J) = B(I,J) - B(K,J)*A(I,K)
+   80                     CONTINUE
    90                 CONTINUE
   100             CONTINUE
               END IF