Implements P-space discovery, and preconditioning via diagonal elements, within the solver and handlers.
