<div><div>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.</div><div><br></div><div>Matrix Object: 1 MPI processes</div>
<div>  type: seqaij</div><div>row 0: (2, 1) </div><div>row 1: (2, 1) </div><div>row 2: (0, 1)  (1, 1)</div></div><div><br></div><div>#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</div>
<div>#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</div>
<div>#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</div>
<div>#3  0x00007ffff5b2318d in MatGetOrdering_RCM (mat=0x1ee9a00, type=0xa875b0 "rcm", row=0x7fffffffc1d0, col=0x7fffffffc1c8) at /home/jed/petsc/src/mat/order/sprcm.c:18</div><div>#4  0x00007ffff5b260fb in MatGetOrdering (mat=0x1ee9a00, type=0xa875b0 "rcm", rperm=0x7fffffffc1d0, cperm=0x7fffffffc1c8) at /home/jed/petsc/src/mat/order/sorder.c:229</div>
<div>#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</div>
<div>#6  0x00007ffff7b84f60 in dFSBuildSpace_Cont (fs=0xa864e0) at /home/jed/dohp/src/fs/impls/cont/cont.c:288</div><div>#7  0x00007ffff7b533f9 in dFSBuildSpace (fs=0xa864e0) at /home/jed/dohp/src/fs/interface/fs.c:351</div>
<div>#8  0x00007ffff7b61d00 in dFSSetFromOptions (fs=0xa864e0) at /home/jed/dohp/src/fs/interface/fsreg.c:127</div><div>#9  0x0000000000409572 in EllipSetFromOptions (elp=0x9fdc80) at /home/jed/dohp/src/fs/tests/ellip.c:256</div>
<div>#10 0x000000000040e95f in main (argc=15, argv=0x7fffffffd658) at /home/jed/dohp/src/fs/tests/ellip.c:644</div>