<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Benoit :</div><div class="gmail_quote"><br></div><div class="gmail_quote">mumps options can be called from PETSc C-code as</div><div class="gmail_quote">ierr = PetscOptionsInsertString("-mat_mumps_icntl_13 1 -mat_mumps_icntl_24 1 -mat_mumps_cntl_3 1e-12");CHKERRQ(ierr); <br></div><div class="gmail_quote"><br></div><div class="gmail_quote">or see petsc/src/ksp/ksp/examples/tutorials/ex52.c</div><div class="gmail_quote"><br></div><div class="gmail_quote">I do not know how to translate these calls into python API though.</div><div class="gmail_quote"><br></div><div class="gmail_quote">Hong</div><div class="gmail_quote"><br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Dear petsc-users,<br>
<br>
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)<br>
<br>
I have written the following piece of code<br>
<br>
    ksp = PETSc.KSP()<br>
    ksp.create(PETSc.COMM_WORLD)<br>
    ksp.setOperators(A)<br>
    ksp.setType('preonly')<br>
    pc=ksp.getPC()<br>
    pc.setType('lu')<br>
    pc.setFactorSolverPackage('mumps')<br>
<br>
It seems to works, however I have 2 questions<br>
1/ how to set the Icntl parameters of mumps ? I didn't find the function MatMumpsSetIcntl.<br>
2/ Is this piece of code right or written in good petsc style ?<br>
<br>
thanks a lot,<br>
<br>
Benoit<br>
<br>
</blockquote></div><br></div></div>