丁老师 <[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.

Attachment: pgpOAsusr0Hkj.pgp
Description: PGP signature

Reply via email to