search

FEniCS / ffcx

Next generation FEniCS Form Compiler for finite element forms
https://fenicsproject.org
Other
144 stars 38 forks source link

Quadrature-points symbol not stored as table #650

Open jorgensd opened 7 months ago

jorgensd commented 7 months ago

MWE:

import basix.ufl
import ufl.geometry
cell = "interval"
el =basix.ufl.element("Lagrange", cell, 1, shape=(1, ))
mesh = ufl.Mesh(el)

x = ufl.geometry.CellCoordinate(mesh)
J = x[0]*ufl.dx

FFCx doesn't current create a table for the reference quadrature points, which means that the CellCoordinate and FacetCoordinate from ufl doesn't work.