Generalization of analytic configuration interaction (CI) gradient techniques for potential energy hypersurfaces, including a solution to the coupled perturbed Hartree–Fock equations for multiconfiguration SCF molecular wave functions