Please run with -mat_view_info -mat_no_inode -mat_view_info_detailed and send to petsc-maint at mcs.anl.gov
Barry On Sun, 8 Apr 2007, yaron at oak-research.com wrote: > 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? > > Yaron >
