dolfin team mailing list archive

Thread
Date

[Question #152123]: problem with project and cell_domains

To: dolfin@xxxxxxxxxxxxxxxxxxx
From: Garth Wells <question152123@xxxxxxxxxxxxxxxxxxxxx>
Date: Fri, 08 Apr 2011 16:18:07 -0000
Reply-to: question152123@xxxxxxxxxxxxxxxxxxxxx
Sender: bounces@xxxxxxxxxxxxx

New question #152123 on DOLFIN:
https://answers.launchpad.net/dolfin/+question/152123

The following code breaks when cell_domains is defined on the mesh

m = Interval(101, 0, 10.0)
m.data().create_mesh_function('cell_domains',m.topology().dim())
m.data().mesh_function('cell_domains').array()[:] = 0
m.data().mesh_function('cell_domains').array()[50:60] = 1

Q = FunctionSpace(m, 'CG',1)
project(Constant(10.0),Q)

___ERROR______

6678   return _cpp.solve(*args)
   6679 
   6680 def normalize(*args):

RuntimeError: *** Error: Zero pivot detected in uBLAS ILU preconditioner in row 51.

-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.