[petsc-dev] What should MatGetRowIJ(..., symmetric=TRUE, ...) do with empty diagonal?

[Barry Smith]

   Jed,

   Could you perhaps save the matrix as a binary file and send it to petsc-maint at mcs.anl.gov so I can track down why it crashes here and fix it. 

  Please append this email so I know the bug I am tracking.


    Thanks 

    Barry

On Apr 23, 2011, at 12:02 PM, Jed Brown wrote:

> MatGetOrdering is failing with the following matrix. Should it fail (i.e. diagonal is required) or not? It seems to work for many larger sizes even without setting a diagonal.
> 
> Matrix Object: 1 MPI processes
>   type: seqaij
> row 0: (2, 1) 
> row 1: (2, 1) 
> row 2: (0, 1)  (1, 1)
> 
> #0  0x00007ffff5a70ffc in MatGetRowIJ_SeqAIJ_Inode_Symmetric (A=0x1ee9a00, iia=0x7fffffffbe68, jja=0x7fffffffbe60, ishift=0, oshift=1) at /home/jed/petsc/src/mat/impls/aij/seq/inode.c:108
> #1  0x00007ffff5a72808 in MatGetRowIJ_SeqAIJ_Inode (A=0x1ee9a00, oshift=1, symmetric=PETSC_TRUE, blockcompressed=PETSC_TRUE, n=0x7fffffffbe74, ia=0x7fffffffbe68, ja=0x7fffffffbe60, done=0x7fffffffbe54) at /home/jed/petsc/src/mat/impls/aij/seq/inode.c:244
> #2  0x00007ffff5705629 in MatGetRowIJ (mat=0x1ee9a00, shift=1, symmetric=PETSC_TRUE, inodecompressed=PETSC_TRUE, n=0x7fffffffbe74, ia=0x7fffffffbe68, ja=0x7fffffffbe60, done=0x7fffffffbe54) at /home/jed/petsc/src/mat/interface/matrix.c:6650
> #3  0x00007ffff5b2318d in MatGetOrdering_RCM (mat=0x1ee9a00, type=0xa875b0 "rcm", row=0x7fffffffc1d0, col=0x7fffffffc1c8) at /home/jed/petsc/src/mat/order/sprcm.c:18
> #4  0x00007ffff5b260fb in MatGetOrdering (mat=0x1ee9a00, type=0xa875b0 "rcm", rperm=0x7fffffffc1d0, cperm=0x7fffffffc1c8) at /home/jed/petsc/src/mat/order/sorder.c:229
> #5  0x00007ffff7b8312d in dMeshPopulateOrderedSet_Private (mesh=0xa01920, orderedSet=0xb000000000000018, explicitSet=0xb000000000000019, dirichletSet=0xb00000000000001a, ghostSet=0xb00000000000001b, orderingtype=0xa875b0 "rcm") at /home/jed/dohp/src/fs/impls/cont/cont.c:159
> #6  0x00007ffff7b84f60 in dFSBuildSpace_Cont (fs=0xa864e0) at /home/jed/dohp/src/fs/impls/cont/cont.c:288
> #7  0x00007ffff7b533f9 in dFSBuildSpace (fs=0xa864e0) at /home/jed/dohp/src/fs/interface/fs.c:351
> #8  0x00007ffff7b61d00 in dFSSetFromOptions (fs=0xa864e0) at /home/jed/dohp/src/fs/interface/fsreg.c:127
> #9  0x0000000000409572 in EllipSetFromOptions (elp=0x9fdc80) at /home/jed/dohp/src/fs/tests/ellip.c:256
> #10 0x000000000040e95f in main (argc=15, argv=0x7fffffffd658) at /home/jed/dohp/src/fs/tests/ellip.c:644