Problem in new function ?SYTRI2 was found: the part of A below the diagonal is changed in the case UPLO='U' .

But in the description of arguments

If UPLO = 'U', the upper triangular part of the

inverse is formed and the part of A below the diagonal is not

referenced; if UPLO = 'L' the lower triangular part of the

inverse is formed and the part of A above the diagonal is

not referenced.

These elements zeroized after calling ?GEMM function in ?SYTRI2X.

CALL SGEMM('T','N',NNB,NNB,CUT,ONE,A(1,CUT+1),LDA,

$ WORK,N+NB+1, ZERO, A(CUT+1,CUT+1), LDA)

Also the error in the description of ?ORCSD/UNCSD functions

IWORK (workspace) INTEGER array, dimension (M-Q)

The dimension of IWORK is M actually.