hermitian_banded_matrix_t Derived Type

private pure function get_element(this, row, col) result(element)

Retrieves the element at position (row, col) of the original matrix. See the LAPACK documentation, element $a_{ij}$ of the original matrix is stored at position $(kd + 1 + i - j, j)$ if uplo = "U" and at position $(1 + i - j, j)$ if uplo = "L" in the banded storage.

Arguments

Return Value complex(kind=dp)

the element at the original position (row, col)

private pure subroutine set_element(this, row, col, element)

Sets the element $a_{ij}$ of the original array into the banded structure. The row and col arguments refer to the row and column indices of the element in the original array. This routine has no effect if the location falls outside of the banded structure.

Arguments

private pure function get_total_nb_elements(this)

Returns the total number of elements inside the banded matrix

Arguments

Return Value integer

private pure function is_compatible_with(this, other)

Checks if a Hermitian band matrix is compatible with another Hermitian band matrix. This implies that the following attributes should be equal: - number of rows/columns - number of sub/superdiagonals - storage of upper or lower triangular part - dimensions of the banded matrices themselves Returns .true. if all criteria are satisfied, .false. otherwise.

Arguments

Return Value logical

private pure subroutine destroy(this)

Destructor, deallocates the datastructure.

Arguments