丁老师 <[email protected]> writes: > which petsc function perform the reordering of the matrix
MatGetOrdering() gives you the permutation vectors. You should be doing this reordering with the *mesh*. You can represent mesh connectivity in a matrix and use MatPartitioning and MatGetOrdering, but you should use the result to redistribute the mesh before assembling the matrices containing your physics.
pgpOAsusr0Hkj.pgp
Description: PGP signature
