Viewing a reordered rectangular matrix...

yaron at yaron at
Sun Apr 8 18:22:27 CDT 2007

I have a non-square  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:

0]PETSC ERROR: MatSetValues_SeqAIJ() line 141 in
src/mat/impls/aij/seq/aij.c [0]PETSC ERROR: Argument out of range!
[0]PETSC ERROR: Column too large: col 905218 max 790999!
[0]PETSC ERROR: MatPermute_SeqAIJ() line 1824 in
src/mat/impls/aij/seq/aij.c [0]PETSC ERROR: MatPermute() line 3158 in
src/mat/interface/matrix.c [0]PETSC ERROR: MatGetOrdering() line 251 in

What gives?

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <>

More information about the petsc-users mailing list