Added the E->control.use_petsc flag, set to off by default
More additions to Solve_ahat_p_fhat_PETSc method