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
_______________________________________________
Mailing list: https://launchpad.net/~dolfin
Post to : dolfin@lists.launchpad.net
Unsubscribe : https://launchpad.net/~dolfin
More help : https://help.launchpad.net/ListHelp
