packed_index Module

tools for packing 2-index and 4-index quantities by exploitinh index symmetry


Interfaces

public interface index_packed

index for the upper triangular matrix in a packed storage scheme SEE: http://www.netlib.org/lapack/lug/node123.html optimization: use precalculated array sq_half(i)=

  • private elemental function index_packed4(i, j) result(ij)

    Arguments

    Type IntentOptional Attributes Name
    integer, intent(in) :: i
    integer, intent(in) :: j

    Return Value integer(kind=8)

  • private elemental function index_packed8(i, j) result(ij)

    Arguments

    Type IntentOptional Attributes Name
    integer(kind=8), intent(in) :: i
    integer(kind=8), intent(in) :: j

    Return Value integer(kind=8)


Subroutines

public subroutine init_index_packed()

initialize index_packed: precompute index values

Arguments

None