Variational cluster approximation and exact diagonalization for calculation of Hubbard models. (March 2019)