SUBROUTINE LA_GELSY( A, B, RANK=rank, &
JPVT= jpvt, RCOND= rcond, INFO= info )
type(wp), INTENT(INOUT) :: A(:,:), rhs
INTEGER, INTENT(OUT), OPTIONAL :: RANK
INTEGER, INTENT(INOUT), OPTIONAL :: JPVT(:)
REAL(wp), INTENT(IN), OPTIONAL :: RCOND
INTEGER, INTENT(OUT), OPTIONAL :: INFO
where
type ::= REAL COMPLEX
wp ::= KIND(1.0) KIND(1.0D0)
rhs ::= B(:,:) B(:)