File.lshift of a Function from a tensorial FunctionSpace with symmetry=True throws error

3 months ago by
So I was trying to write a second order tensor expression to a file in 2017.2.0 and get thrown the error:
RuntimeError: interpolate_vertex_values: Failed to set physical value size correctly for subelements.​

Minimal not working example:

from fenics import *

mesh = UnitCubeMesh(2,2,2)

FE = FiniteElement('Lagrange', mesh.ufl_cell(), 1)
TFE = TensorElement(FE, symmetry = True)

TensorSpace = FunctionSpace(mesh, TFE)

File('bla.pvd') << project(Identity(3), TensorSpace)

Am I missing something here? Is this a bug?

Community: FEniCS Project

1 Answer

3 months ago by
Yes, that's a bug. As a workaround you can project/interpolate to an non-symmetric space on the same mesh.
Please login to add an answer/comment or follow this question.

Similar posts:
Search »