dolfin team mailing list archive

[Chaffra <754078@xxxxxxxxxxxxxxxxxx>]

Public bug reported:

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.

** Affects: dolfin
     Importance: Undecided
         Status: New

-- 
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:
  New

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.