SUBROUTINE LA_TEST_CGELSD( M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK,& & WORK, LWORK, RWORK, IWORK, INFO) ! ! -- LAPACK95 interface driver routine (version 1.1) -- ! UNI-C, Denmark, ! September 10, 1999 ! ! .. Use Statements .. USE LA_PRECISION, ONLY: WP => SP USE F95_LAPACK, ONLY: LA_GELSD ! .. Implicit Statement .. IMPLICIT NONE ! .. Scalar Arguments .. INTEGER, INTENT(IN) :: M, N, NRHS, LDA, LDB INTEGER, INTENT(INOUT) :: INFO INTEGER, INTENT(OUT) :: RANK INTEGER, INTENT(IN) :: LWORK REAL(WP), INTENT(IN) :: RCOND ! .. Array Arguments .. COMPLEX(WP), INTENT(INOUT) :: A(1:LDA,1:N), B(1:LDB,1:NRHS) COMPLEX(WP), INTENT(OUT) :: WORK(1:LWORK) INTEGER :: IWORK(1:3 * MIN(M,N) * MAX(M,N)*MAX(M,N) + 11 * MIN(M,N)) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! REAL(WP) :: RWORK(1: 8*N+2*MAX(M,N)*MAX(M,N)+8*MAX(M,N)*MAX(M,N)+N*NRHS) REAL(WP), INTENT(OUT) :: S(1:MIN(M,N)) ! .. Parameters .. CHARACTER(LEN=8), PARAMETER :: SRNAME = 'LA_GELSD' CHARACTER(LEN=14), PARAMETER :: SRNAMT = 'LA_TEST_CGELSD' ! .. Common blocks .. INTEGER :: INFOTC COMMON /LINFO95/ INFOTC ! .. Local Scalars .. INTEGER :: I, J, IA1, IA2, IB1, IB2, IS ! .. Local Arrays .. LOGICAL, SAVE :: CTEST = .TRUE., ETEST = .TRUE. ! .. Executable Statements .. IA1 = M; IA2 = N; IB1 = MAX(1,M,N); IB2 = NRHS; IS = MIN(M,N) I = INFO / 100; J = INFO - I*100 SELECT CASE(I) CASE(0) IF (NRHS==1) THEN CALL LA_GELSD( A(1:IA1,1:IA2), B(1:IB1,1), RANK,& & S(1:IS), RCOND, INFO ) INFO = INFOTC ELSE CALL LA_GELSD( A(1:IA1,1:IA2), B(1:IB1,1:IB2), RANK,& & S(1:IS), RCOND, INFO ) INFO = INFOTC ENDIF CASE (2) IB1 = MAX(1, N, M) - 1 CASE(4) IS = MIN (M, N) - 1 CASE(:-1,1,3, 5:) CALL UESTOP(SRNAMT) END SELECT IF (I/=0) THEN SELECT CASE (NRHS) CASE (2:) CALL LA_GELSD( A(1:IA1,1:IA2), B(1:IB1,1:IB2), RANK,& & S(1:IS), RCOND, INFO ) CASE(1) CALL LA_GELSD( A(1:IA1,1:IA2), B(1:IB1,1), RANK,& & S(1:IS), RCOND, INFO ) CASE(:-1) CALL UESTOP(SRNAMT) END SELECT ENDIF CALL LA_AUX_AA01( I, CTEST, ETEST, SRNAMT ) END SUBROUTINE LA_TEST_CGELSD