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
These elements zeroized after calling ?GEMM function in ?SYTRI2X.
$ 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.