comments unified: 'End of' comment: added where omitted

BLAS/SRC/caxpy.f

@@ -133,4 +133,7 @@ SUBROUTINE CAXPY(N,CA,CX,INCX,CY,INCY)
       END IF
 *
       RETURN
+*
+*     End of CAXPY
+*
       END

BLAS/SRC/ccopy.f

@@ -119,4 +119,7 @@ SUBROUTINE CCOPY(N,CX,INCX,CY,INCY)
          END DO
       END IF
       RETURN
+*
+*     End of CCOPY
+*
       END

BLAS/SRC/cdotc.f

@@ -128,4 +128,7 @@ COMPLEX FUNCTION CDOTC(N,CX,INCX,CY,INCY)
       END IF
       CDOTC = CTEMP
       RETURN
+*
+*     End of CDOTC
+*
       END

BLAS/SRC/cdotu.f

@@ -125,4 +125,7 @@ COMPLEX FUNCTION CDOTU(N,CX,INCX,CY,INCY)
       END IF
       CDOTU = CTEMP
       RETURN
+*
+*     End of CDOTU
+*
       END

BLAS/SRC/cscal.f

@@ -111,4 +111,7 @@ SUBROUTINE CSCAL(N,CA,CX,INCX)
          END DO
       END IF
       RETURN
+*
+*     End of CSCAL
+*
       END

BLAS/SRC/csrot.f

@@ -147,4 +147,7 @@ SUBROUTINE CSROT( N, CX, INCX, CY, INCY, C, S )
          END DO
       END IF
       RETURN
+*
+*     End of CSROT
+*
       END

BLAS/SRC/csscal.f

@@ -114,4 +114,7 @@ SUBROUTINE CSSCAL(N,SA,CX,INCX)
          END DO
       END IF
       RETURN
+*
+*     End of CSSCAL
+*
       END

BLAS/SRC/cswap.f

@@ -123,4 +123,7 @@ SUBROUTINE CSWAP(N,CX,INCX,CY,INCY)
          END DO
       END IF
       RETURN
+*
+*     End of CSWAP
+*
       END

BLAS/SRC/dasum.f

@@ -125,4 +125,7 @@ DOUBLE PRECISION FUNCTION DASUM(N,DX,INCX)
       END IF
       DASUM = DTEMP
       RETURN
+*
+*     End of DASUM
+*
       END

BLAS/SRC/daxpy.f

@@ -146,4 +146,7 @@ SUBROUTINE DAXPY(N,DA,DX,INCX,DY,INCY)
          END DO
       END IF
       RETURN
+*
+*     End of DAXPY
+*
       END

BLAS/SRC/dcabs1.f

@@ -60,4 +60,7 @@ DOUBLE PRECISION FUNCTION DCABS1(Z)
 *
       DCABS1 = ABS(DBLE(Z)) + ABS(DIMAG(Z))
       RETURN
+*
+*     End of DCABS1
+*
       END

BLAS/SRC/dcopy.f

@@ -140,4 +140,7 @@ SUBROUTINE DCOPY(N,DX,INCX,DY,INCY)
          END DO
       END IF
       RETURN
+*
+*     End of DCOPY
+*
       END

BLAS/SRC/ddot.f

@@ -142,4 +142,7 @@ DOUBLE PRECISION FUNCTION DDOT(N,DX,INCX,DY,INCY)
       END IF
       DDOT = DTEMP
       RETURN
+*
+*     End of DDOT
+*
       END

BLAS/SRC/drot.f

@@ -136,4 +136,7 @@ SUBROUTINE DROT(N,DX,INCX,DY,INCY,C,S)
          END DO
       END IF
       RETURN
+*
+*     End of DROT
+*
       END

BLAS/SRC/drotm.f

@@ -194,4 +194,7 @@ SUBROUTINE DROTM(N,DX,INCX,DY,INCY,DPARAM)
          END IF
       END IF
       RETURN
+*
+*     End of DROTM
+*
       END

BLAS/SRC/drotmg.f

@@ -254,6 +254,9 @@ SUBROUTINE DROTMG(DD1,DD2,DX1,DY1,DPARAM)
 
       DPARAM(1) = DFLAG
       RETURN
+*
+*     End of DROTMG
+*
       END
 
 

BLAS/SRC/dscal.f

@@ -130,4 +130,7 @@ SUBROUTINE DSCAL(N,DA,DX,INCX)
          END DO
       END IF
       RETURN
+*
+*     End of DSCAL
+*
       END

BLAS/SRC/dsdot.f

@@ -166,4 +166,7 @@ DOUBLE PRECISION FUNCTION DSDOT(N,SX,INCX,SY,INCY)
          END DO
       END IF
       RETURN
+*
+*     End of DSDOT
+*
       END

BLAS/SRC/dswap.f

@@ -147,4 +147,7 @@ SUBROUTINE DSWAP(N,DX,INCX,DY,INCY)
          END DO
       END IF
       RETURN
+*
+*     End of DSWAP
+*
       END

BLAS/SRC/dzasum.f

@@ -112,4 +112,7 @@ DOUBLE PRECISION FUNCTION DZASUM(N,ZX,INCX)
       END IF
       DZASUM = STEMP
       RETURN
+*
+*     End of DZASUM
+*
       END

BLAS/SRC/icamax.f

@@ -121,4 +121,7 @@ INTEGER FUNCTION ICAMAX(N,CX,INCX)
          END DO
       END IF
       RETURN
+*
+*     End of ICAMAX
+*
       END

BLAS/SRC/idamax.f

@@ -120,4 +120,7 @@ INTEGER FUNCTION IDAMAX(N,DX,INCX)
          END DO
       END IF
       RETURN
+*
+*     End of IDAMAX
+*
       END

BLAS/SRC/isamax.f

@@ -120,4 +120,7 @@ INTEGER FUNCTION ISAMAX(N,SX,INCX)
          END DO
       END IF
       RETURN
+*
+*     End of ISAMAX
+*
       END

BLAS/SRC/izamax.f

@@ -121,4 +121,7 @@ INTEGER FUNCTION IZAMAX(N,ZX,INCX)
          END DO
       END IF
       RETURN
+*
+*     End of IZAMAX
+*
       END

BLAS/SRC/sasum.f

@@ -126,4 +126,7 @@ REAL FUNCTION SASUM(N,SX,INCX)
       END IF
       SASUM = STEMP
       RETURN
+*
+*     End of SASUM
+*
       END

BLAS/SRC/saxpy.f

@@ -146,4 +146,7 @@ SUBROUTINE SAXPY(N,SA,SX,INCX,SY,INCY)
          END DO
       END IF
       RETURN
+*
+*     End of SAXPY
+*
       END

BLAS/SRC/scabs1.f

@@ -59,4 +59,7 @@ REAL FUNCTION SCABS1(Z)
 *     ..
       SCABS1 = ABS(REAL(Z)) + ABS(AIMAG(Z))
       RETURN
+*
+*     End of SCABS1
+*
       END

BLAS/SRC/scasum.f

@@ -111,4 +111,7 @@ REAL FUNCTION SCASUM(N,CX,INCX)
       END IF
       SCASUM = STEMP
       RETURN
+*
+*     End of SCASUM
+*
       END

BLAS/SRC/scopy.f

@@ -140,4 +140,7 @@ SUBROUTINE SCOPY(N,SX,INCX,SY,INCY)
          END DO
       END IF
       RETURN
+*
+*     End of SCOPY
+*
       END

BLAS/SRC/sdot.f

@@ -142,4 +142,7 @@ REAL FUNCTION SDOT(N,SX,INCX,SY,INCY)
       END IF
       SDOT = STEMP
       RETURN
+*
+*     End of SDOT
+*
       END

BLAS/SRC/sdsdot.f

@@ -157,4 +157,7 @@ REAL FUNCTION SDSDOT(N,SB,SX,INCX,SY,INCY)
       END IF
       SDSDOT = DSDOT
       RETURN
+*
+*     End of SDSDOT
+*
       END

BLAS/SRC/srot.f

@@ -136,4 +136,7 @@ SUBROUTINE SROT(N,SX,INCX,SY,INCY,C,S)
          END DO
       END IF
       RETURN
+*
+*     End of SROT
+*
       END

BLAS/SRC/srotm.f

@@ -195,4 +195,7 @@ SUBROUTINE SROTM(N,SX,INCX,SY,INCY,SPARAM)
          END IF
       END IF
       RETURN
+*
+*     End of SROTM
+*
       END

BLAS/SRC/srotmg.f

@@ -254,8 +254,7 @@ SUBROUTINE SROTMG(SD1,SD2,SX1,SY1,SPARAM)
 
       SPARAM(1) = SFLAG
       RETURN
+*
+*     End of SROTMG
+*
       END
-
-
-
-

BLAS/SRC/sscal.f

@@ -130,4 +130,7 @@ SUBROUTINE SSCAL(N,SA,SX,INCX)
          END DO
       END IF
       RETURN
+*
+*     End of SSCAL
+*
       END

BLAS/SRC/sswap.f

@@ -147,4 +147,7 @@ SUBROUTINE SSWAP(N,SX,INCX,SY,INCY)
          END DO
       END IF
       RETURN
+*
+*     End of SSWAP
+*
       END

BLAS/SRC/xerbla_array.f

@@ -113,4 +113,7 @@ SUBROUTINE XERBLA_ARRAY(SRNAME_ARRAY, SRNAME_LEN, INFO)
       CALL XERBLA( SRNAME, INFO )
 
       RETURN
+*
+*     End of XERBLA_ARRAY
+*
       END

BLAS/SRC/zaxpy.f

@@ -133,4 +133,7 @@ SUBROUTINE ZAXPY(N,ZA,ZX,INCX,ZY,INCY)
       END IF
 *
       RETURN
+*
+*     End of ZAXPY
+*
       END

BLAS/SRC/zcopy.f

@@ -119,4 +119,7 @@ SUBROUTINE ZCOPY(N,ZX,INCX,ZY,INCY)
          END DO
       END IF
       RETURN
+*
+*     End of ZCOPY
+*
       END

BLAS/SRC/zdotc.f

@@ -128,4 +128,7 @@ COMPLEX*16 FUNCTION ZDOTC(N,ZX,INCX,ZY,INCY)
       END IF
       ZDOTC = ZTEMP
       RETURN
+*
+*     End of ZDOTC
+*
       END

BLAS/SRC/zdotu.f

@@ -125,4 +125,7 @@ COMPLEX*16 FUNCTION ZDOTU(N,ZX,INCX,ZY,INCY)
       END IF
       ZDOTU = ZTEMP
       RETURN
+*
+*     End of ZDOTU
+*
       END

BLAS/SRC/zdrot.f

@@ -147,4 +147,7 @@ SUBROUTINE ZDROT( N, ZX, INCX, ZY, INCY, C, S )
          END DO
       END IF
       RETURN
+*
+*     End of ZDROT
+*
       END

BLAS/SRC/zdscal.f

@@ -114,4 +114,7 @@ SUBROUTINE ZDSCAL(N,DA,ZX,INCX)
          END DO
       END IF
       RETURN
+*
+*     End of ZDSCAL
+*
       END

BLAS/SRC/zscal.f

@@ -111,4 +111,7 @@ SUBROUTINE ZSCAL(N,ZA,ZX,INCX)
          END DO
       END IF
       RETURN
+*
+*     End of ZSCAL
+*
       END

BLAS/SRC/zswap.f

@@ -123,4 +123,7 @@ SUBROUTINE ZSWAP(N,ZX,INCX,ZY,INCY)
          END DO
       END IF
       RETURN
+*
+*     End of ZSWAP
+*
       END