subroutine init_get_helement_tj
get_helement_lattice_ex_mat => get_helement_tJ_ex_mat
get_helement_lattice_general => get_helement_tJ_general
call init_tmat(lat)
call setup_exchange_matrix(lat)
end subroutine init_get_helement_tj