dolfin team mailing list archive

[Garth Wells <754078@xxxxxxxxxxxxxxxxxx>]

** Converted to question:
   https://answers.launchpad.net/dolfin/+question/152123

-- 
You received this bug notification because you are a member of DOLFIN
Team, which is subscribed to DOLFIN.
https://bugs.launchpad.net/bugs/754078

Title:
  problem with project and cell_domains

Status in DOLFIN:
  Invalid

Bug description:
  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.