about MatSolve() and KSP

[Yujie]

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>