dolfin team mailing list archive

[Martin Sandve Alnæs <martinal@xxxxxxxxx>]

The following code:


from dolfin import *
mesh = UnitSquare(10, 10)
x = mesh.ufl_cell().x
V = FunctionSpace(mesh, "CG", 1)
M = V*V
f = project(x[0], M.sub(0))
plot(f)


fails with this message:


Traceback (most recent call last):
  File "subspaceproject.py", line 6, in <module>
    f = project(x[0], M.sub(0))
  File "/opt/fenics/dorsal-unstable/lib/python2.7/site-packages/dolfin/fem/projection.py",
line 71, in project
    A = assemble(a, form_compiler_parameters=form_compiler_parameters)
  File "/opt/fenics/dorsal-unstable/lib/python2.7/site-packages/dolfin/fem/assembling.py",
line 146, in assemble
    add_values)
  File "/opt/fenics/dorsal-unstable/lib/python2.7/site-packages/dolfin/cpp.py",
line 16735, in assemble
    return _cpp.assemble(*args)
RuntimeError: *** Error: Cannot determine ownership range for sub dofmaps.


How can I go about this? My current goal is to define a boundary
condition on a subspace from a ufl expression.

Martin