subroutine init_get_helement_heisenberg
get_helement_lattice_ex_mat => get_helement_heisenberg_ex_mat
get_helement_lattice_general => get_helement_heisenberg_general
call setup_exchange_matrix(lat)
end subroutine init_get_helement_heisenberg