1# create linear solver 2ksp = PETSc.KSP() 3ksp.create(PETSc.COMM_WORLD) 4# use conjugate gradients 5ksp.setType('cg') 6# and incomplete Cholesky 7ksp.getPC().setType('icc') 8# obtain sol & rhs vectors 9x, b = A.createVecs() 10x.set(0) 11b.set(1) 12# and next solve 13ksp.setOperators(A) 14ksp.setFromOptions() 15ksp.solve(b, x) 16