Dear PETSc Developers: I have checked the differerence between 2.3.3 and 3.0.0 in using external package.
Now, we don't need to use Matconvert() to convert the matrix to the type suitable for external package, just to use MatGetFactor(). The following codes are the realization of MatSolveXXX() after calling MatGetFactor(). 53: MatFactorInfoInitialize(&info); 54: MatGetFactor(C,MAT_SOLVER_PETSC,MAT_FACTOR_LU,&A); 55: MatLUFactorSymbolic(A,C,row,col,&info); 56: MatLUFactorNumeric(A,C,&info); 57: MatSolveTranspose(A,b,x); I have checked the description to MatSolveXXX() or MatMatSolve(). "Notes 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()." 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. Regards, Yujie -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090127/45250db3/attachment.htm>
