Re: Viewing a reordered rectangular matrix...

[Barry Smith <bsmith@xxxxxxxxxxx>]

  Please run with -mat_view_info -mat_no_inode -mat_view_info_detailed
and send to petsc-maint@xxxxxxxxxxx

   Barry


On Sun, 8 Apr 2007, yaron@xxxxxxxxxxxxxxxx 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
>