Xiangdong <[email protected]> writes: > I am not clear about this. In this example, where is the format declaration > (e.g.,mpiaij) of Mat jac? Why is this mat jac always distributed in > a compatible way as the DM vectors?
src/ksp/ksp/examples/tutorials/ex25.c calls KSPSetDM and sets functions to assemble the matrices and right hand side. Those objects are created internally (you can use -dm_mat_type sbaij if you like) and the user doesn't have to see them.
pgpvoWI02Z5OW.pgp
Description: PGP signature
