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

[Jed Brown]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20110423/ff8651fc/attachment.html>