The gradient and the laplacian of the density used to compute the xc potential are now computed directly from the wavefunctions.