matlu(natom)%(nsppol,nspinor,nspinor,ndim,ndim) :: input quantity to diagonalize
natom=number of atoms
prtopt= option to define level of printing
nsppol_imp= if 1, one can diagonalize with the same matrix the Up
and Dn matlu matrix. It is convenient because one can thus have the
same interaction matrix for up and dn spins.

m_matlu/init_matlu [ Functions ]

natom = number of atoms
nspinor = number of spinorial components
nsppol = number of polarisation components
lpawu_natom(natom) = value of lpawu for all atoms
maltu <type(matlu_type)>= density matrix in the local orbital basis and related variables
OUTPUTS
maltu <type(matlu_type)>= density matrix in the local orbital basis and related variables