BEGIN_PROVIDER [integer, rank]
&BEGIN_PROVIDER [integer, rank_16]
    BEGIN_DOC
    ! calculates the rank of matrix
    END_DOC

    implicit none

    rank=nt1*nt2
    if(MOD(rank,16).eq.0)then
        rank_16=rank
    else
        rank_16=rank+16-MOD(rank,16)
    endif
END_PROVIDER