Pseudospectral Hartree-Fock gradient calculations