cannot store symmetric tensor values function to file
Issue #1065
new
from dolfin import *
mesh = UnitCubeMesh(1,1,1)
#V=TensorFunctionSpace(mesh, "CG", 1, shape=(3,3))
V=TensorFunctionSpace(mesh, "CG", 1, shape=(3,3), symmetry=True)
f = Function(V)
File("data/test.pvd") << f
gives the following error when a symmetric functionspace is used:
Traceback (most recent call last):
File "test_tensorfunction.py", line 10, in <module>
File("data/test.pvd") << f
RuntimeError: interpolate_vertex_values: Failed to set physical value size correctly for subelements.
the code works when using symmetric=False
.
thanks Florian