582
583
spEDM=spEDM, spuEDM=spuEDM, E_dim2=N_mu)
584
585
call weight_DM( N_Elec, Elecs, N_mu, mus, na_u, lasto, &
585
spDM, spDMneq, spEDM, n_s, sc_off)
586
sp_dist, sparse_pattern, Ss, &
587
spDM, spDMneq, spEDM, n_s, sc_off, DE_NEGF)
587
589
call update_DM(sp_dist,sparse_pattern, n_nzs, &
588
590
DM(:,ispin), spDM, Ef=Ef, &