It computes the average of the diagonal elements of the LDA+U potential for each spin For details see notes on LDA+U implementation