subroutine init_get_helement_tj_guga
call init_tmat(lat)
call init_spin_free_tmat(lat)
call setup_exchange_matrix(lat)
call setup_spin_free_exchange(lat)
get_umat_el => get_umat_heisenberg_spin_free
end subroutine init_get_helement_tj_guga