[petsc-dev] MatLUFactor() on MATSEQAIJ with no-inodes

[Lisandro Dalcin]

A petsc4py testcase (TestMatFactorB2 in tests/test_mat_fact.py) is
failing with petsc-dev.

Basically, I'm creating a SEQAIJ matrix, assembling a Laplace 1D with
finite direrences, and calling MatLUFactor() ... As you can see, the
matrix is not actually using inodes (as expected), but anyway it seems
that the virtual table is still using MatLUFactorNumeric_SeqAIJ_Inode.

Interestingly a similar test with a diagonal matrix (I mean, SEQAIJ
where A[i,i] = 1.0/(i+1))



testFactorILU (__main__.TestMatFactorB2) ... ok
testFactorLU (__main__.TestMatFactorB2) ...
Matrix Object:
  type=seqaij, rows=10, cols=10
  total: nonzeros=28, allocated nonzeros=50
    not using I-node routines
ERROR

======================================================================
ERROR: testFactorLU (__main__.TestMatFactorB2)
----------------------------------------------------------------------
Traceback (most recent call last):
  File "test/test_mat_fact.py", line 76, in testFactorLU
    self.A.factorLU(r,c,{})
  File "Mat.pyx", line 838, in petsc4py.PETSc.Mat.factorLU
(src/petsc4py.PETSc.c:68161)
Error: error code 77
[0] MatLUFactor() line 2391 in src/mat/interface/matrix.c
[0] MatLUFactor_SeqAIJ() line 948 in src/mat/impls/aij/seq/aijfact.c
[0] MatLUFactorNumeric() line 2561 in src/mat/interface/matrix.c
[0] MatLUFactorNumeric_SeqAIJ_Inode() line 1238 in src/mat/impls/aij/seq/inode.c
[0] Petsc has generated inconsistent data
[0] Matrix without inode information


-- 
Lisandro Dalcin
---------------
Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC)
Instituto de Desarrollo Tecnológico para la Industria Química (INTEC)
Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)
PTLC - Güemes 3450, (3000) Santa Fe, Argentina
Tel/Fax: +54-(0)342-451.1594