This branch of the code was never tested and does not run. One simple way to do this is to replace the in-house `CGSolverSteihaug` with a `petsc4py.PETSc.KSP` of type `"stcg"` See https://petsc.org/release/manualpages/KSP/KSPSTCG/