lu_implementation Submodule


Uses

  • module~~lu_implementation~~UsesGraph module~lu_implementation lu_implementation module~tools_math tools_math module~lu_implementation->module~tools_math module~class_vector class_vector module~tools_math->module~class_vector module~class_psblas class_psblas module~tools_math->module~class_psblas module~class_vector->module~class_psblas module~class_stopwatch class_stopwatch module~class_psblas->module~class_stopwatch module~tools_psblas tools_psblas module~class_psblas->module~tools_psblas module~class_stopwatch->module~tools_psblas psb_base_mod psb_base_mod module~class_stopwatch->psb_base_mod psb_prec_mod psb_prec_mod module~tools_psblas->psb_prec_mod module~tools_psblas->psb_base_mod psb_krylov_mod psb_krylov_mod module~tools_psblas->psb_krylov_mod

Contents


Subroutines

subroutine lu_fact_2(a, ipiv, info)

Arguments

Type IntentOptional AttributesName
real(kind=psb_dpk_), intent(inout) :: a(2,2)
integer, intent(out) :: ipiv(2)
integer, intent(out) :: info

subroutine lu_fact_3(a, ipiv, info)

Arguments

Type IntentOptional AttributesName
real(kind=psb_dpk_), intent(inout) :: a(3,3)
integer, intent(out) :: ipiv(3)
integer, intent(out) :: info

subroutine lu_fact_4(a, ipiv, info)

Arguments

Type IntentOptional AttributesName
real(kind=psb_dpk_), intent(inout) :: a(4,4)
integer, intent(out) :: ipiv(4)
integer, intent(out) :: info

Module Procedures

module procedure lu_fact module subroutine lu_fact(n, a, ipiv, info)

Description: Adapted from LAPACK and unrolled for small sizes.

Arguments

Type IntentOptional AttributesName
integer, intent(in) :: n
real(kind=psb_dpk_), intent(inout) :: a(n,n)
integer, intent(out) :: ipiv(n)
integer, intent(out) :: info

module procedure lu_solve module subroutine lu_solve(n, a, ipiv, b, info)

Solve a system Ax=b where A has already been factored. Adapted from LAPACK, this is supposed to be called with small sizes (4 or thereabout)

Arguments

Type IntentOptional AttributesName
integer, intent(in) :: n
real(kind=psb_dpk_), intent(in) :: a(n,n)
integer, intent(in) :: ipiv(n)
real(kind=psb_dpk_), intent(inout) :: b(n)
integer, intent(out) :: info