Jed, Could you perhaps save the matrix as a binary file and send it to [email protected] 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