dense_to_banded Interface

public interface dense_to_banded

An interface to the dense to banded matrix conversion routines.


Module Procedures

private subroutine dense_to_banded_dbl(a, kl, ku, x, err)

Converts a banded matrix stored in dense format to a compressed form.

The banded matrix is stored in a compressed form supplied column by column. The following code segment transfers between a full matrix to the bonded matrix storage scheme. \code{fortran} do j = 1, n k = ku + 1 - j do i = max(1, j - ku), min(n, j + kl) a(k + i, j) = matrix(i, j) end do end do \endcode

Arguments

Type IntentOptional Attributes Name
real(kind=real64), intent(in), dimension(:,:) :: a

The matrix to convert.

integer(kind=int32), intent(in) :: kl

The number of subdiagonals. Must be at least 0.

integer(kind=int32), intent(in) :: ku

The number of superdiagonals. Must be at least 0.

real(kind=real64), intent(out), dimension(:,:) :: x

The (KL+KU+1)-by-N banded matrix.

class(errors), intent(inout), optional, target :: err

An error object to report any errors that occur.

private subroutine dense_to_banded_cmplx(a, kl, ku, x, err)

Converts a banded matrix stored in dense format to a compressed form.

The banded matrix is stored in a compressed form supplied column by column. The following code segment transfers between a full matrix to the bonded matrix storage scheme. \code{fortran} do j = 1, n k = ku + 1 - j do i = max(1, j - ku), min(n, j + kl) a(k + i, j) = matrix(i, j) end do end do \endcode

Arguments

Type IntentOptional Attributes Name
complex(kind=real64), intent(in), dimension(:,:) :: a

The matrix to convert.

integer(kind=int32), intent(in) :: kl

The number of subdiagonals. Must be at least 0.

integer(kind=int32), intent(in) :: ku

The number of superdiagonals. Must be at least 0.

complex(kind=real64), intent(out), dimension(:,:) :: x

The (KL+KU+1)-by-N banded matrix.

class(errors), intent(inout), optional, target :: err

An error object to report any errors that occur.