Second elasticity tensor does not show minor symmetries

[dokken]

Projection onto a quadrature space would only evaluate the data at the relevant quadrature points.

See for instance

Projecting into a quadrature space

The following works for me: import dolfin as df df.parameters["form_compiler"]["representation"] = "quadrature" mesh = df.UnitSquareMesh(10, 10) # Step 1: Create a Quadrature element specifying the quadrature scheme quad_element = df.VectorElement( "Quadrature", mesh.ufl_cell(), degree=4, quad_scheme="default" ) V = df.FunctionSpace(mesh, quad_element) u = df.TrialFunction(V) uf = df.Function(V) v = df.TestFunction(V) f = df.Expression(("sin(x[0])*cos(x[1])", "cos(x[0])*sin(x[1])"),…