cannot store symmetric tensor values function to file

Issue #1065 new
Florian Bruckner created an issue
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

Comments (0)

  1. Log in to comment