[petsc-users] Field split degree of freedom ordering

[Alexander Lindsay]

In the block matrices documentation, it's stated: "Note that for interlaced
storage the number of rows/columns of each block must be the same size" Is
interlacing defined in a global sense, or a process-local sense? So
explicitly, if I don't want the same size restriction, do I need to ensure
that globally all of my block 1 dofs are numbered after my block 0 dofs? Or
do I need to follow that on a process-local level? Essentially in libMesh
we always follow rank-major ordering. I'm asking whether for unequal row
sizes, in order to split, would we need to strictly follow variable-major
ordering (splitting here meaning splitting by variable)?

Alex
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221101/ddaa4f4c/attachment.html>