Dear petsc-users, I am a beginner in petsc and I have some question with the python interface. I am trying to solve problem with mumps (or other direct sparse solver)
I have written the following piece of code ksp = PETSc.KSP() ksp.create(PETSc.COMM_WORLD) ksp.setOperators(A) ksp.setType('preonly') pc=ksp.getPC() pc.setType('lu') pc.setFactorSolverPackage('mumps') It seems to works, however I have 2 questions 1/ how to set the Icntl parameters of mumps ? I didn't find the function MatMumpsSetIcntl. 2/ Is this piece of code right or written in good petsc style ? thanks a lot, Benoit