<div class="gmail_quote">On Mon, Nov 21, 2011 at 04:00, Thomas Witkowski <span dir="ltr"><<a href="mailto:Thomas.Witkowski@tu-dresden.de">Thomas.Witkowski@tu-dresden.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Some technical question on this point: How can I explicitly factorize a sequential matrix? Is MatLUFactor the correct function to do it?</blockquote><div><br></div><div>In general, I recommend using the KSP. You can KSPSetType(ksp,KSPPREONLY), PCSetType(pc,PCLU). Call KSPSolve() multiple times, the factorization will be reused. Or you can use MatLUFactor() directly if you really want.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"> If so, how can I provide the package (i.e. umfpack) that should be used for factorization?</blockquote></div>
<div><br></div><div>See the MatSolverPackage argument in one of these.</div><div><br></div><div><a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/PC/PCFactorSetMatSolverPackage.html">http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/PC/PCFactorSetMatSolverPackage.html</a></div>
<div><br></div><div><a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Mat/MatGetFactor.html">http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Mat/MatGetFactor.html</a></div>