David A. Ham
This PR implements Dual function spaces, cofunctions in dual function spaces, arguments, and all the consequential objects needed to support dual operations and operators in UFL. This is documented in...
suppose in a complex valued inverse problem we have: ```python x = SpatialCoordinate(mesh) functional = f(x)*dx ``` where `f(x)` is some UFL expression which might be a function of x....
`ufl.constantvalue.as_ufl` checks values against the intrinsic types `float`, `complex`, and `int`. This breaks for other scalar types, for example entries in numpy arrays. The check should instead be for `numbers.Real`,...
Consider the following Firedrake MFE: ```python from firedrake import * from firedrake.__future__ import interpolate mesh = UnitSquareMesh(2, 2) Q = FunctionSpace(mesh, "CG", 1) u = Function(Q) vom = VertexOnlyMesh(mesh, ([0.5,...
When taking derivatives with respect to a `Cofunction`, ufl ends up creating `ExprList` objects containing cofunctions. Equality (`==`) then fails on the `ExprList` because the `Cofunction` is not identified as...