Instead of diagonalizing the hamiltonian, perform CEPA calculation, CEPA type ncepa. This is currently available only for single configuration reference functions.