<p>Dear PETSc Developers:</p><p>I have checked the differerence between 2.3.3 and 3.0.0 in using external package.&nbsp;</p><p>Now, we don&#39;t need to use Matconvert() to convert the matrix to the type suitable for external package, just to use MatGetFactor().&nbsp;</p>
<p>The following codes are the realization of MatSolveXXX() after calling MatGetFactor().</p><p>53:   MatFactorInfoInitialize(&amp;info);</p><p>54:   MatGetFactor(C,MAT_SOLVER_PETSC,MAT_FACTOR_LU,&amp;A);<br></p>&nbsp;55:   MatLUFactorSymbolic(A,C,row,col,&amp;info);<br>
&nbsp;56:   MatLUFactorNumeric(A,C,&amp;info);<br><p>&nbsp;57:   MatSolveTranspose(A,b,x);</p><p>I have checked the description to MatSolveXXX() or MatMatSolve().&nbsp;<br></p><p>&quot;Notes<br>Most users should employ the simplified KSP interface for linear solvers instead of working directly with matrix algebra routines such as this. See, e.g., KSPCreate().&quot;</p>
<p>You advise to use KSP interface for calling MatSolveXXX(). I am wondering if it is ok to directly call MatSolveXX(), why need KSP interface? thanks a lot.</p><p>Regards,</p><p>Yujie<br><br></p><br>