I have a non-square&nbsp; SEQAIJ matrix, which I'd like to reorder. The reordering (Nested Dissection, in case that matters) goes through silently, but when I want to display the nonzero structure of the new matrix (with -mat_view_ordering_draw), I get a runtime error:<br><br>======================================<br>0]PETSC ERROR: MatSetValues_SeqAIJ() line 141 in<br>src/mat/impls/aij/seq/aij.c [0]PETSC ERROR: Argument out of range!<br>[0]PETSC ERROR: Column too large: col 905218 max 790999!<br>[0]PETSC ERROR: MatPermute_SeqAIJ() line 1824 in<br>src/mat/impls/aij/seq/aij.c [0]PETSC ERROR: MatPermute() line 3158 in<br>src/mat/interface/matrix.c [0]PETSC ERROR: MatGetOrdering() line 251 in<br>==============================================<br><br>What gives?<br><br>Yaron<br>