band_lu_solve Subroutine

private subroutine band_lu_solve(ab_loc, n, kl_loc, ku_loc, b, x)

Banded Gaussian elimination without pivoting. Band storage: AB(ku+1+i-j, j) = A(i,j). On exit x = A^{-1} b. The matrix ab_loc is overwritten with LU factors.

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(inout) :: ab_loc(kl_loc+ku_loc+1,n)
integer, intent(in) :: n
integer, intent(in) :: kl_loc
integer, intent(in) :: ku_loc
real(kind=wp), intent(in) :: b(n)
real(kind=wp), intent(out) :: x(n)