dolfin team mailing list archive

[Chaffra <754078@xxxxxxxxxxxxxxxxxx>]

Thanks for the reply Garth.
I guess I  do not understand how the projection works. What I am trying to do here is to setup subdomains on the mesh. This did not use to cause problems with dolfin 0.9.9. What am I doing wrong?

Thanks,
Chaffra

-- 
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.