Second derivatives on nonaffine domains break in element table extraction

Issue #109 new
Martin Sandve Alnæs created an issue

Reproduce with:

from dolfin import *
parameters["form_compiler"]["representation"] = "uflacs"
mesh = UnitDiscMesh(mpi_comm_world(), 3, 2, 2)
V = FunctionSpace(mesh, FiniteElement("CG", mesh.ufl_cell(), 2))
f = Function(V)
print assemble(div(grad(f))*dx)

Initial incomplete attempts at fixing in

https://bitbucket.org/fenics-project/ffc/branch/martinal/fix-x-element-derivatives

https://bitbucket.org/fenics-project/ufl/branch/martinal/fix-extract-elements

Comments (7)

  1. Log in to comment