c//@implements LapackFlow

      SUBROUTINE SOLVE(N, NRHS, A, B, X)
      INTEGER N
      INTEGER NRHS
      REAL A(N,N)
      REAL B(N,NRHS)
      REAL X(N,NRHS)
      INTEGER INFO
      REAL AF(N,N)
      REAL R(N)
      REAL C(N)
      INTEGER IPIV(N)
      CHARACTER EQUED
      REAL RCOND
      REAL FERR(NRHS)
      REAL BERR(NRHS)
      REAL WORK(4*N)
      INTEGER IWORK(N)

      call SGESVX('E', 'T', N, NRHS, A, N, AF, N, IPIV, EQUED, R, C, B,
     *N, X, N, RCOND, FERR, BERR, WORK, IWORK, INFO)
      RETURN
      END



      SUBROUTINE SOLVELS(M, N, NRHS, A, BX)
      INTEGER M
      INTEGER N
      INTEGER NRHS
      REAL A(N,N)
      REAL BX(N,NRHS)
      INTEGER INFO
      REAL SV(N)
      REAL RANK
      INTEGER LWORK
      REAL WORK(10000)
      INTEGER IWORK(10000)
      LWORK = 10000

      call SGELSD(M, N, NRHS, A, M, BX, M, SV, -1, RANK, WORK, LWORK, 
     *IWORK, INFO)
      RETURN
      END
      


      SUBROUTINE EIG(N, A, D, V)
      REAL D(N)
      REAL E(N)
      REAL TAU(N-1)
      REAL WORK(28*N)
      INTEGER INFO
      REAL FDUMMY
      INTEGER IDUMMY
      INTEGER IDUMMY2
      REAL ISUPPZ(2*N)
      INTEGER IWORK(10*N)
c      CALL SSYEV('V', 'U', N, A, N, D, WORK, 18*N, INFO)
      CALL SSYEVR('V', 'A', 'U', N, A, N, FDUMMY, FDUMMY, IDUMMY, 
     *IDUMMY, DLAMCH, IDUMMY2, D, V, N, ISUPPZ, WORK, 28*N, IWORK, 10*N,
     *INFO)

      RETURN
      END



      SUBROUTINE SVD(M, N, A, U, SIGMA, VT)
      INTEGER LWORK
      REAL WORK(15*N*N)
      INTEGER IWORK(8*N)
      INTEGER INFO

      CALL SGESDD('A', M, N, A, M, SIGMA, U, M, VT, N, WORK, 15*N*N, 
     *IWORK, INFO)

c      CALL SGESVD('A', 'A', M, N, A, M, SIGMA, U, M, VT, N, WORK, 
c     *15*N*N, INFO)

      RETURN
      END