segfault on loki

[wd15]

Segfault on loki with

from fipy import Grid1D, CellVariable
from fipy import TransientTerm, ConvectionTerm, ImplicitSourceTerm

m = Grid1D(nx=2, dx=0.5)
v = CellVariable(mesh=m)

v.constrain(1, m.facesLeft)
v.constrain(0, m.facesRight)
v.faceGrad.constrain([0], m.facesLeft)

(-TransientTerm(coeff=1 / m.x) h2. 
 ConvectionTerm(coeff=[[1]])
 + ImplicitSourceTerm(coeff=m.x)).solve(v, dt=1.)

print v

v[:] = 0
(TransientTerm(coeff=1 / m.x)
 - ConvectionTerm(coeff=[[1]])
 - ImplicitSourceTerm(coeff=m.x)).solve(v, dt=1.)

print v
raw_input('stopped')

Imported from trac ticket #467, created by wd15 on 09-20-2012 at 16:38, last modified: 09-21-2012 at 16:16

[wd15]

This is an issue with pysparse LU being unable to factorize 2x2 matrices with the upper right entry being zero. Both Trilinos and Scipy's LU work. I'm not fixing this issue.

Trac comment by wd15 on 09-21-2012 at 16:16