scale_linalgebra::linalgebra_solvelineq Interface Reference

Public Member Functions
subroutine	linalgebra_solvelineq_b1d (a, b, x)
subroutine	linalgebra_solvelineq_b2d (a, b, x)

Detailed Description

Definition at line 32 of file scale_linalgebra.F90.

Member Function/Subroutine Documentation

linalgebra_solvelineq_b1d()

subroutine scale_linalgebra::linalgebra_solvelineq::linalgebra_solvelineq_b1d	(	real(rp), dimension(:,:), intent(in)	a,
		real(rp), dimension(:), intent(in)	b,
		real(rp), dimension(size(b)), intent(out)	x )

Definition at line 96 of file scale_linalgebra.F90.

 
    real(RP), intent(in) :: A(:,:)
    real(RP), intent(in) :: b(:)
    real(RP), intent(out) :: x(size(b))
 
    real(RP) :: A_lu(size(A,1),size(A,2))
    integer :: ipiv(size(A,1))
    integer :: n, info
 
    !--------------------------------------------------------------------------- 
 
    a_lu = a
    n = size(a,1)
 
    call dgetrf(n, n, a_lu, n, ipiv, info)
    if (info /=0 ) then
      log_error("linalgebra_SolveLinEq",*) "Matrix is singular"
      call prc_abort
    end if
 
    x(:) = b
    call dgetrs('N', n, 1, a_lu, n, ipiv, x, n, info)
    if (info /=0 ) then
      log_error("linalgebra_SolveLinEq",*)  "Matrix inversion is failed"
      call prc_abort
    end if
 
    return

linalgebra_solvelineq_b2d()

subroutine scale_linalgebra::linalgebra_solvelineq::linalgebra_solvelineq_b2d	(	real(rp), dimension(:,:), intent(in)	a,
		real(rp), dimension(:,:), intent(in)	b,
		real(rp), dimension(size(b,1),size(b,2)), intent(out)	x )

Definition at line 128 of file scale_linalgebra.F90.

 
    real(RP), intent(in) :: A(:,:)
    real(RP), intent(in) :: b(:,:)
    real(RP), intent(out) :: x(size(b,1),size(b,2))
 
    real(RP) :: A_lu(size(A,1),size(A,2))
    integer :: ipiv(size(A,1))
    integer :: n, info
 
    !--------------------------------------------------------------------------- 
 
    a_lu = a
    n = size(a,1)
 
    call dgetrf(n, n, a_lu, n, ipiv, info)
    if (info /=0 ) then
      log_error("linalgebra_SolveLinEq",*) "Matrix is singular"
      call prc_abort
    end if
 
    x(:,:) = b
    call dgetrs('N', n, size(b,2), a_lu, n, ipiv, x, n, info)
    if (info /=0 ) then
      log_error("linalgebra_SolveLinEq",*)  "Matrix inversion is failed"
      call prc_abort
    end if
 
    return

---

The documentation for this interface was generated from the following file:

• scale_linalgebra.F90