65 INTEGER,
INTENT(IN) :: nspins, natoms
67 INTEGER :: i, iatom, iblock, j, jatom
72 admm_dm%method = admm_control%method
73 admm_dm%eps_filter = admm_control%eps_filter
75 ALLOCATE (admm_dm%mcweeny_history(nspins))
79 ALLOCATE (admm_dm%block_map(natoms, natoms))
80 admm_dm%block_map(:, :) = 0
81 DO iblock = 1,
SIZE(admm_control%blocks)
82 DO i = 1,
SIZE(admm_control%blocks(iblock)%list)
83 iatom = admm_control%blocks(iblock)%list(i)
84 DO j = 1,
SIZE(admm_control%blocks(iblock)%list)
85 jatom = admm_control%blocks(iblock)%list(j)
86 admm_dm%block_map(iatom, jatom) = 1