Committing Sébastien Villemot patches (5) sent on Jan 6th 2016 to LAPACK Mailing list

From Sébastien Villemot:
I attach to this message 5 patches that are currently applied to the
Debian package for LAPACK, and that I think could be usefully merged in
the LAPACK SVN tree.

The patches include fixes for old outstanding bugs in the CBLAS test
programs, fixes for the build system, and fixes for various typos.

Each patch file contains a more detailed description of its purpose.

40 files changed, 50 insertions, 50 deletions

diff --git a/SRC/cgejsv.f b/SRC/cgejsv.f
index c003b5f5..e2f6cea7 100644
--- a/SRC/cgejsv.f
+++ b/SRC/cgejsv.f
@@ -317,7 +317,7 @@
 *>               the minimal requirement is LWORK >= 5*N+2*N*N. 
 *>            4.2. if JOBV.EQ.'J' the minimal requirement is 
 *>               LWORK >= 4*N+N*N.
-*>            In both cases, the allocated CWORK can accomodate blocked runs
+*>            In both cases, the allocated CWORK can accommodate blocked runs
 *>            of CGEQP3, CGEQRF, CGELQF, SUNMQR, CUNMLQ.
 *> \endverbatim
 *>
diff --git a/SRC/claed7.f b/SRC/claed7.f
index c1441393..02e0cab2 100644
--- a/SRC/claed7.f
+++ b/SRC/claed7.f
@@ -57,7 +57,7 @@
 *>
 *>       The first stage consists of deflating the size of the problem
 *>       when there are multiple eigenvalues or if there is a zero in
-*>       the Z vector.  For each such occurence the dimension of the
+*>       the Z vector.  For each such occurrence the dimension of the
 *>       secular equation problem is reduced by one.  This stage is
 *>       performed by the routine SLAED2.
 *>
diff --git a/SRC/clarrv.f b/SRC/clarrv.f
index ecedfa4d..895ad7b1 100644
--- a/SRC/clarrv.f
+++ b/SRC/clarrv.f
@@ -81,7 +81,7 @@
 *>          L is REAL array, dimension (N)
 *>          On entry, the (N-1) subdiagonal elements of the unit
 *>          bidiagonal matrix L are in elements 1 to N-1 of L
-*>          (if the matrix is not splitted.) At the end of each block
+*>          (if the matrix is not split.) At the end of each block
 *>          is stored the corresponding shift as given by SLARRE.
 *>          On exit, L is overwritten.
 *> \endverbatim
@@ -236,7 +236,7 @@
 *>          INFO is INTEGER
 *>          = 0:  successful exit
 *>
-*>          > 0:  A problem occured in CLARRV.
+*>          > 0:  A problem occurred in CLARRV.
 *>          < 0:  One of the called subroutines signaled an internal problem.
 *>                Needs inspection of the corresponding parameter IINFO
 *>                for further information.
diff --git a/SRC/clatdf.f b/SRC/clatdf.f
index 11f0dfd9..441cb317 100644
--- a/SRC/clatdf.f
+++ b/SRC/clatdf.f
@@ -58,7 +58,7 @@
 *>              Zx = +-e - f with the sign giving the greater value of
 *>              2-norm(x).  About 5 times as expensive as Default.
 *>          IJOB .ne. 2: Local look ahead strategy where
-*>              all entries of the r.h.s. b is choosen as either +1 or
+*>              all entries of the r.h.s. b is chosen as either +1 or
 *>              -1.  Default.
 *> \endverbatim
 *>
diff --git a/SRC/cstegr.f b/SRC/cstegr.f
index 6e1eae05..900e8437 100644
--- a/SRC/cstegr.f
+++ b/SRC/cstegr.f
@@ -48,7 +48,7 @@
 *> either an interval (VL,VU] or a range of indices IL:IU for the desired
 *> eigenvalues.
 *>
-*> CSTEGR is a compatability wrapper around the improved CSTEMR routine.
+*> CSTEGR is a compatibility wrapper around the improved CSTEMR routine.
 *> See SSTEMR for further details.
 *>
 *> One important change is that the ABSTOL parameter no longer provides any
diff --git a/SRC/dlaed1.f b/SRC/dlaed1.f
index c37c1d21..ffee218d 100644
--- a/SRC/dlaed1.f
+++ b/SRC/dlaed1.f
@@ -54,7 +54,7 @@
 *>
 *>       The first stage consists of deflating the size of the problem
 *>       when there are multiple eigenvalues or if there is a zero in
-*>       the Z vector.  For each such occurence the dimension of the
+*>       the Z vector.  For each such occurrence the dimension of the
 *>       secular equation problem is reduced by one.  This stage is
 *>       performed by the routine DLAED2.
 *>
diff --git a/SRC/dlaed7.f b/SRC/dlaed7.f
index 658ece9a..31999b7c 100644
--- a/SRC/dlaed7.f
+++ b/SRC/dlaed7.f
@@ -59,7 +59,7 @@
 *>
 *>       The first stage consists of deflating the size of the problem
 *>       when there are multiple eigenvalues or if there is a zero in
-*>       the Z vector.  For each such occurence the dimension of the
+*>       the Z vector.  For each such occurrence the dimension of the
 *>       secular equation problem is reduced by one.  This stage is
 *>       performed by the routine DLAED8.
 *>
diff --git a/SRC/dlag2.f b/SRC/dlag2.f
index a941b940..791161e0 100644
--- a/SRC/dlag2.f
+++ b/SRC/dlag2.f
@@ -99,7 +99,7 @@
 *>          will always be positive.  If the eigenvalues are real, then
 *>          the first (real) eigenvalue is  WR1 / SCALE1 , but this may
 *>          overflow or underflow, and in fact, SCALE1 may be zero or
-*>          less than the underflow threshhold if the exact eigenvalue
+*>          less than the underflow threshold if the exact eigenvalue
 *>          is sufficiently large.
 *> \endverbatim
 *>
@@ -112,7 +112,7 @@
 *>          eigenvalues are real, then the second (real) eigenvalue is
 *>          WR2 / SCALE2 , but this may overflow or underflow, and in
 *>          fact, SCALE2 may be zero or less than the underflow
-*>          threshhold if the exact eigenvalue is sufficiently large.
+*>          threshold if the exact eigenvalue is sufficiently large.
 *> \endverbatim
 *>
 *> \param[out] WR1
@@ -266,8 +266,8 @@
 *     Note: the test of R in the following IF is to cover the case when
 *           DISCR is small and negative and is flushed to zero during
 *           the calculation of R.  On machines which have a consistent
-*           flush-to-zero threshhold and handle numbers above that
-*           threshhold correctly, it would not be necessary.
+*           flush-to-zero threshold and handle numbers above that
+*           threshold correctly, it would not be necessary.
 *
       IF( DISCR.GE.ZERO .OR. R.EQ.ZERO ) THEN
          SUM = PP + SIGN( R, PP )
diff --git a/SRC/dlamrg.f b/SRC/dlamrg.f
index 7126053e..c31e9e14 100644
--- a/SRC/dlamrg.f
+++ b/SRC/dlamrg.f
@@ -50,7 +50,7 @@
 *> \param[in] N2
 *> \verbatim
 *>          N2 is INTEGER
-*>         These arguements contain the respective lengths of the two
+*>         These arguments contain the respective lengths of the two
 *>         sorted lists to be merged.
 *> \endverbatim
 *>
diff --git a/SRC/dlarre.f b/SRC/dlarre.f
index e7eea10c..2212c574 100644
--- a/SRC/dlarre.f
+++ b/SRC/dlarre.f
@@ -244,7 +244,7 @@
 *> \verbatim
 *>          INFO is INTEGER
 *>          = 0:  successful exit
-*>          > 0:  A problem occured in DLARRE.
+*>          > 0:  A problem occurred in DLARRE.
 *>          < 0:  One of the called subroutines signaled an internal problem.
 *>                Needs inspection of the corresponding parameter IINFO
 *>                for further information.
diff --git a/SRC/dlarrf.f b/SRC/dlarrf.f
index f054caa8..de062817 100644
--- a/SRC/dlarrf.f
+++ b/SRC/dlarrf.f
@@ -51,7 +51,7 @@
 *> \param[in] N
 *> \verbatim
 *>          N is INTEGER
-*>          The order of the matrix (subblock, if the matrix splitted).
+*>          The order of the matrix (subblock, if the matrix split).
 *> \endverbatim
 *>
 *> \param[in] D
diff --git a/SRC/dlarrv.f b/SRC/dlarrv.f
index 828661f2..8f79f330 100644
--- a/SRC/dlarrv.f
+++ b/SRC/dlarrv.f
@@ -81,7 +81,7 @@
 *>          L is DOUBLE PRECISION array, dimension (N)
 *>          On entry, the (N-1) subdiagonal elements of the unit
 *>          bidiagonal matrix L are in elements 1 to N-1 of L
-*>          (if the matrix is not splitted.) At the end of each block
+*>          (if the matrix is not split.) At the end of each block
 *>          is stored the corresponding shift as given by DLARRE.
 *>          On exit, L is overwritten.
 *> \endverbatim
@@ -236,7 +236,7 @@
 *>          INFO is INTEGER
 *>          = 0:  successful exit
 *>
-*>          > 0:  A problem occured in DLARRV.
+*>          > 0:  A problem occurred in DLARRV.
 *>          < 0:  One of the called subroutines signaled an internal problem.
 *>                Needs inspection of the corresponding parameter IINFO
 *>                for further information.
diff --git a/SRC/dlasd1.f b/SRC/dlasd1.f
index 7b66d90b..c457d1ab 100644
--- a/SRC/dlasd1.f
+++ b/SRC/dlasd1.f
@@ -60,7 +60,7 @@
 *>
 *>    The first stage consists of deflating the size of the problem
 *>    when there are multiple singular values or when there are zeros in
-*>    the Z vector.  For each such occurence the dimension of the
+*>    the Z vector.  For each such occurrence the dimension of the
 *>    secular equation problem is reduced by one.  This stage is
 *>    performed by the routine DLASD2.
 *>
diff --git a/SRC/dlasd6.f b/SRC/dlasd6.f
index a5238b91..f8b7dea0 100644
--- a/SRC/dlasd6.f
+++ b/SRC/dlasd6.f
@@ -74,7 +74,7 @@
 *>
 *>       The first stage consists of deflating the size of the problem
 *>       when there are multiple singular values or if there is a zero
-*>       in the Z vector. For each such occurence the dimension of the
+*>       in the Z vector. For each such occurrence the dimension of the
 *>       secular equation problem is reduced by one. This stage is
 *>       performed by the routine DLASD7.
 *>
diff --git a/SRC/dlasdq.f b/SRC/dlasdq.f
index 6beef32a..cbfdb4a5 100644
--- a/SRC/dlasdq.f
+++ b/SRC/dlasdq.f
@@ -59,7 +59,7 @@
 *> \verbatim
 *>          UPLO is CHARACTER*1
 *>        On entry, UPLO specifies whether the input bidiagonal matrix
-*>        is upper or lower bidiagonal, and wether it is square are
+*>        is upper or lower bidiagonal, and whether it is square are
 *>        not.
 *>           UPLO = 'U' or 'u'   B is upper bidiagonal.
 *>           UPLO = 'L' or 'l'   B is lower bidiagonal.
diff --git a/SRC/dlasrt.f b/SRC/dlasrt.f
index f5d0e6cd..6d7cbb13 100644
--- a/SRC/dlasrt.f
+++ b/SRC/dlasrt.f
@@ -123,7 +123,7 @@
 *     ..
 *     .. Executable Statements ..
 *
-*     Test the input paramters.
+*     Test the input parameters.
 *
       INFO = 0
       DIR = -1
diff --git a/SRC/dlatdf.f b/SRC/dlatdf.f
index be70313b..2a8fabe0 100644
--- a/SRC/dlatdf.f
+++ b/SRC/dlatdf.f
@@ -58,7 +58,7 @@
 *>              Zx = +-e - f with the sign giving the greater value
 *>              of 2-norm(x). About 5 times as expensive as Default.
 *>          IJOB .ne. 2: Local look ahead strategy where all entries of
-*>              the r.h.s. b is choosen as either +1 or -1 (Default).
+*>              the r.h.s. b is chosen as either +1 or -1 (Default).
 *> \endverbatim
 *>
 *> \param[in] N
diff --git a/SRC/dsgesv.f b/SRC/dsgesv.f
index 64e2c168..9d42c329 100644
--- a/SRC/dsgesv.f
+++ b/SRC/dsgesv.f
@@ -164,7 +164,7 @@
 *>               -3 : failure of SGETRF
 *>               -31: stop the iterative refinement after the 30th
 *>                    iterations
-*>          > 0: iterative refinement has been sucessfully used.
+*>          > 0: iterative refinement has been successfully used.
 *>               Returns the number of iterations
 *> \endverbatim
 *>
diff --git a/SRC/dsposv.f b/SRC/dsposv.f
index bb72199b..49896a31 100644
--- a/SRC/dsposv.f
+++ b/SRC/dsposv.f
@@ -168,7 +168,7 @@
 *>               -3 : failure of SPOTRF
 *>               -31: stop the iterative refinement after the 30th
 *>                    iterations
-*>          > 0: iterative refinement has been sucessfully used.
+*>          > 0: iterative refinement has been successfully used.
 *>               Returns the number of iterations
 *> \endverbatim
 *>
diff --git a/SRC/dstegr.f b/SRC/dstegr.f
index 298e1c76..541d4655 100644
--- a/SRC/dstegr.f
+++ b/SRC/dstegr.f
@@ -48,7 +48,7 @@
 *> either an interval (VL,VU] or a range of indices IL:IU for the desired
 *> eigenvalues.
 *>
-*> DSTEGR is a compatability wrapper around the improved DSTEMR routine.
+*> DSTEGR is a compatibility wrapper around the improved DSTEMR routine.
 *> See DSTEMR for further details.
 *>
 *> One important change is that the ABSTOL parameter no longer provides any
diff --git a/SRC/slaed1.f b/SRC/slaed1.f
index 74eeb633..bef873a1 100644
--- a/SRC/slaed1.f
+++ b/SRC/slaed1.f
@@ -54,7 +54,7 @@
 *>
 *>       The first stage consists of deflating the size of the problem
 *>       when there are multiple eigenvalues or if there is a zero in
-*>       the Z vector.  For each such occurence the dimension of the
+*>       the Z vector.  For each such occurrence the dimension of the
 *>       secular equation problem is reduced by one.  This stage is
 *>       performed by the routine SLAED2.
 *>
diff --git a/SRC/slaed7.f b/SRC/slaed7.f
index 3d3d6292..8e3e8d93 100644
--- a/SRC/slaed7.f
+++ b/SRC/slaed7.f
@@ -59,7 +59,7 @@
 *>
 *>       The first stage consists of deflating the size of the problem
 *>       when there are multiple eigenvalues or if there is a zero in
-*>       the Z vector.  For each such occurence the dimension of the
+*>       the Z vector.  For each such occurrence the dimension of the
 *>       secular equation problem is reduced by one.  This stage is
 *>       performed by the routine SLAED8.
 *>
diff --git a/SRC/slag2.f b/SRC/slag2.f
index ad04333a..03baa0c5 100644
--- a/SRC/slag2.f
+++ b/SRC/slag2.f
@@ -99,7 +99,7 @@
 *>          will always be positive.  If the eigenvalues are real, then
 *>          the first (real) eigenvalue is  WR1 / SCALE1 , but this may
 *>          overflow or underflow, and in fact, SCALE1 may be zero or
-*>          less than the underflow threshhold if the exact eigenvalue
+*>          less than the underflow threshold if the exact eigenvalue
 *>          is sufficiently large.
 *> \endverbatim
 *>
@@ -112,7 +112,7 @@
 *>          eigenvalues are real, then the second (real) eigenvalue is
 *>          WR2 / SCALE2 , but this may overflow or underflow, and in
 *>          fact, SCALE2 may be zero or less than the underflow
-*>          threshhold if the exact eigenvalue is sufficiently large.
+*>          threshold if the exact eigenvalue is sufficiently large.
 *> \endverbatim
 *>
 *> \param[out] WR1
@@ -266,8 +266,8 @@
 *     Note: the test of R in the following IF is to cover the case when
 *           DISCR is small and negative and is flushed to zero during
 *           the calculation of R.  On machines which have a consistent
-*           flush-to-zero threshhold and handle numbers above that
-*           threshhold correctly, it would not be necessary.
+*           flush-to-zero threshold and handle numbers above that
+*           threshold correctly, it would not be necessary.
 *
       IF( DISCR.GE.ZERO .OR. R.EQ.ZERO ) THEN
          SUM = PP + SIGN( R, PP )
diff --git a/SRC/slamrg.f b/SRC/slamrg.f
index 6229abd6..815c6bff 100644
--- a/SRC/slamrg.f
+++ b/SRC/slamrg.f
@@ -50,7 +50,7 @@
 *> \param[in] N2
 *> \verbatim
 *>          N2 is INTEGER
-*>         These arguements contain the respective lengths of the two
+*>         These arguments contain the respective lengths of the two
 *>         sorted lists to be merged.
 *> \endverbatim
 *>
diff --git a/SRC/slarre.f b/SRC/slarre.f
index a5b9f2fd..2df98c7d 100644
--- a/SRC/slarre.f
+++ b/SRC/slarre.f
@@ -244,7 +244,7 @@
 *> \verbatim
 *>          INFO is INTEGER
 *>          = 0:  successful exit
-*>          > 0:  A problem occured in SLARRE.
+*>          > 0:  A problem occurred in SLARRE.
 *>          < 0:  One of the called subroutines signaled an internal problem.
 *>                Needs inspection of the corresponding parameter IINFO
 *>                for further information.
diff --git a/SRC/slarrf.f b/SRC/slarrf.f
index 058e5027..b136e61e 100644
--- a/SRC/slarrf.f
+++ b/SRC/slarrf.f
@@ -51,7 +51,7 @@
 *> \param[in] N
 *> \verbatim
 *>          N is INTEGER
-*>          The order of the matrix (subblock, if the matrix splitted).
+*>          The order of the matrix (subblock, if the matrix split).
 *> \endverbatim
 *>
 *> \param[in] D
diff --git a/SRC/slarrv.f b/SRC/slarrv.f
index 73847f39..bd7d5282 100644
--- a/SRC/slarrv.f
+++ b/SRC/slarrv.f
@@ -81,7 +81,7 @@
 *>          L is REAL array, dimension (N)
 *>          On entry, the (N-1) subdiagonal elements of the unit
 *>          bidiagonal matrix L are in elements 1 to N-1 of L
-*>          (if the matrix is not splitted.) At the end of each block
+*>          (if the matrix is not split.) At the end of each block
 *>          is stored the corresponding shift as given by SLARRE.
 *>          On exit, L is overwritten.
 *> \endverbatim
@@ -236,7 +236,7 @@
 *>          INFO is INTEGER
 *>          = 0:  successful exit
 *>
-*>          > 0:  A problem occured in SLARRV.
+*>          > 0:  A problem occurred in SLARRV.
 *>          < 0:  One of the called subroutines signaled an internal problem.
 *>                Needs inspection of the corresponding parameter IINFO
 *>                for further information.
diff --git a/SRC/slasd1.f b/SRC/slasd1.f
index ae076a0f..4157180c 100644
--- a/SRC/slasd1.f
+++ b/SRC/slasd1.f
@@ -60,7 +60,7 @@
 *>
 *>    The first stage consists of deflating the size of the problem
 *>    when there are multiple singular values or when there are zeros in
-*>    the Z vector.  For each such occurence the dimension of the
+*>    the Z vector.  For each such occurrence the dimension of the
 *>    secular equation problem is reduced by one.  This stage is
 *>    performed by the routine SLASD2.
 *>
diff --git a/SRC/slasd6.f b/SRC/slasd6.f
index f79a06d0..e0cb1da9 100644
--- a/SRC/slasd6.f
+++ b/SRC/slasd6.f
@@ -74,7 +74,7 @@
 *>
 *>       The first stage consists of deflating the size of the problem
 *>       when there are multiple singular values or if there is a zero
-*>       in the Z vector. For each such occurence the dimension of the
+*>       in the Z vector. For each such occurrence the dimension of the
 *>       secular equation problem is reduced by one. This stage is
 *>       performed by the routine SLASD7.
 *>
diff --git a/SRC/slasdq.f b/SRC/slasdq.f
index 289ed855..defbe352 100644
--- a/SRC/slasdq.f
+++ b/SRC/slasdq.f
@@ -59,7 +59,7 @@
 *> \verbatim
 *>          UPLO is CHARACTER*1
 *>        On entry, UPLO specifies whether the input bidiagonal matrix
-*>        is upper or lower bidiagonal, and wether it is square are
+*>        is upper or lower bidiagonal, and whether it is square are
 *>        not.
 *>           UPLO = 'U' or 'u'   B is upper bidiagonal.
 *>           UPLO = 'L' or 'l'   B is lower bidiagonal.
diff --git a/SRC/slasrt.f b/SRC/slasrt.f
index e93c0d6d..05a492ad 100644
--- a/SRC/slasrt.f
+++ b/SRC/slasrt.f
@@ -123,7 +123,7 @@
 *     ..
 *     .. Executable Statements ..
 *
-*     Test the input paramters.
+*     Test the input parameters.
 *
       INFO = 0
       DIR = -1
diff --git a/SRC/slatdf.f b/SRC/slatdf.f
index 51773d4e..4f8c0391 100644
--- a/SRC/slatdf.f
+++ b/SRC/slatdf.f
@@ -58,7 +58,7 @@
 *>              Zx = +-e - f with the sign giving the greater value
 *>              of 2-norm(x). About 5 times as expensive as Default.
 *>          IJOB .ne. 2: Local look ahead strategy where all entries of
-*>              the r.h.s. b is choosen as either +1 or -1 (Default).
+*>              the r.h.s. b is chosen as either +1 or -1 (Default).
 *> \endverbatim
 *>
 *> \param[in] N
diff --git a/SRC/sstegr.f b/SRC/sstegr.f
index d98c451f..f943f2c8 100644
--- a/SRC/sstegr.f
+++ b/SRC/sstegr.f
@@ -48,7 +48,7 @@
 *> either an interval (VL,VU] or a range of indices IL:IU for the desired
 *> eigenvalues.
 *>
-*> SSTEGR is a compatability wrapper around the improved SSTEMR routine.
+*> SSTEGR is a compatibility wrapper around the improved SSTEMR routine.
 *> See SSTEMR for further details.
 *>
 *> One important change is that the ABSTOL parameter no longer provides any
diff --git a/SRC/zcgesv.f b/SRC/zcgesv.f
index d7d0a9d2..a27ce0d1 100644
--- a/SRC/zcgesv.f
+++ b/SRC/zcgesv.f
@@ -170,7 +170,7 @@
 *>               -3 : failure of CGETRF
 *>               -31: stop the iterative refinement after the 30th
 *>                    iterations
-*>          > 0: iterative refinement has been sucessfully used.
+*>          > 0: iterative refinement has been successfully used.
 *>               Returns the number of iterations
 *> \endverbatim
 *>
diff --git a/SRC/zcposv.f b/SRC/zcposv.f
index dfa114d9..983b6ef9 100644
--- a/SRC/zcposv.f
+++ b/SRC/zcposv.f
@@ -178,7 +178,7 @@
 *>               -3 : failure of CPOTRF
 *>               -31: stop the iterative refinement after the 30th
 *>                    iterations
-*>          > 0: iterative refinement has been sucessfully used.
+*>          > 0: iterative refinement has been successfully used.
 *>               Returns the number of iterations
 *> \endverbatim
 *>
diff --git a/SRC/zgejsv.f b/SRC/zgejsv.f
index 644e67a3..89e19e63 100644
--- a/SRC/zgejsv.f
+++ b/SRC/zgejsv.f
@@ -318,7 +318,7 @@
 *>               the minimal requirement is LWORK >= 5*N+2*N*N. 
 *>            4.2. if JOBV.EQ.'J' the minimal requirement is 
 *>               LWORK >= 4*N+N*N.
-*>            In both cases, the allocated CWORK can accomodate blocked runs
+*>            In both cases, the allocated CWORK can accommodate blocked runs
 *>            of ZGEQP3, ZGEQRF, ZGELQF, SUNMQR, CUNMLQ.
 *> \endverbatim
 *>
diff --git a/SRC/zlaed7.f b/SRC/zlaed7.f
index ae6e9a36..9a4b4162 100644
--- a/SRC/zlaed7.f
+++ b/SRC/zlaed7.f
@@ -57,7 +57,7 @@
 *>
 *>       The first stage consists of deflating the size of the problem
 *>       when there are multiple eigenvalues or if there is a zero in
-*>       the Z vector.  For each such occurence the dimension of the
+*>       the Z vector.  For each such occurrence the dimension of the
 *>       secular equation problem is reduced by one.  This stage is
 *>       performed by the routine DLAED2.
 *>
diff --git a/SRC/zlarrv.f b/SRC/zlarrv.f
index 3992f14d..acc6c7f1 100644
--- a/SRC/zlarrv.f
+++ b/SRC/zlarrv.f
@@ -81,7 +81,7 @@
 *>          L is DOUBLE PRECISION array, dimension (N)
 *>          On entry, the (N-1) subdiagonal elements of the unit
 *>          bidiagonal matrix L are in elements 1 to N-1 of L
-*>          (if the matrix is not splitted.) At the end of each block
+*>          (if the matrix is not split.) At the end of each block
 *>          is stored the corresponding shift as given by DLARRE.
 *>          On exit, L is overwritten.
 *> \endverbatim
@@ -236,7 +236,7 @@
 *>          INFO is INTEGER
 *>          = 0:  successful exit
 *>
-*>          > 0:  A problem occured in ZLARRV.
+*>          > 0:  A problem occurred in ZLARRV.
 *>          < 0:  One of the called subroutines signaled an internal problem.
 *>                Needs inspection of the corresponding parameter IINFO
 *>                for further information.
diff --git a/SRC/zlatdf.f b/SRC/zlatdf.f
index e90bfede..23b4713c 100644
--- a/SRC/zlatdf.f
+++ b/SRC/zlatdf.f
@@ -58,7 +58,7 @@
 *>              Zx = +-e - f with the sign giving the greater value of
 *>              2-norm(x).  About 5 times as expensive as Default.
 *>          IJOB .ne. 2: Local look ahead strategy where
-*>              all entries of the r.h.s. b is choosen as either +1 or
+*>              all entries of the r.h.s. b is chosen as either +1 or
 *>              -1.  Default.
 *> \endverbatim
 *>
diff --git a/SRC/zstegr.f b/SRC/zstegr.f
index 16a4e789..5213e493 100644
--- a/SRC/zstegr.f
+++ b/SRC/zstegr.f
@@ -48,7 +48,7 @@
 *> either an interval (VL,VU] or a range of indices IL:IU for the desired
 *> eigenvalues.
 *>
-*> ZSTEGR is a compatability wrapper around the improved ZSTEMR routine.
+*> ZSTEGR is a compatibility wrapper around the improved ZSTEMR routine.
 *> See DSTEMR for further details.
 *>
 *> One important change is that the ABSTOL parameter no longer provides any