solvels()

subroutine, private ompi_kronlininv::solvels ( real(dp), dimension(:,:), intent(in)  A, real(dp), dimension(:,:), intent(in)  B, real(dp), dimension(:,:), intent(out)  sol  )

private

Solves a linear system AX = B, real numbers

Parameters

[in]	A	coefficients matrix (linear operator)
[in]	B	matrix
[out]	sol	solution to the linear system (matrix)

Date

5/8/2016 - Initial Version

Definition at line 1207 of file ompi_kronlininv.f08.

    real(dp),intent(in) :: a(:,:),b(:,:)
    real(dp),intent(out) :: sol(:,:)
    integer,allocatable :: ipiv(:)
    integer :: n,nrhs,lda,ldb
    integer :: info
    real(dp),allocatable :: a2(:,:)
    ! External procedures defined in LAPACK
    external dgesv
    n = size(a,1)
    nrhs = size(b,2)
    allocate(a2(n,n))
    a2 = a
    lda = size(a,1)
    allocate(ipiv(n))
    ldb = size(b,1)
    sol = b
    call dgesv(n,nrhs,a2,lda,ipiv,sol,ldb,info)
    if (info/=0) then
       write(*,*) "linear system solver failed...'"
       print*,'info: ',info
       stop
    end if

Here is the caller graph for this function: