Joe Dean
Joe Dean
Hello, The convergence rate of an H(div)-conforming field computed from an H(curl)-conforming field with the discrete curl operator appears to be suboptimal. Consider the following minimal failing example, where a...
This PR adds support for mixed domain assembly. Forms with test functions, trial functions, and coefficients defined over different meshes can be assembled. The user must provide a map to...
Currently, `Form` computes which entities to integrate over from a `MeshTags` object. This is fine for many common use cases, but in more specialised applications, it is useful to be...
Interpolation on a mesh where a process has no cells currently causes a segfault. Returning an empty array of interpolation coordinates on processes with no cells (see https://github.com/FEniCS/dolfinx/tree/jpdean/submesh_function_fix) fixes this....
This PR allows the creation of kernels for mixed dimensional integrals of codimension 1. It is a step towards https://github.com/FEniCS/dolfinx/projects/9 The implementation is short and simple but is a bit...
Minimal failing example: ``` import ffc.codegeneration.jit import ufl cell = ufl.quadrilateral element = ufl.FiniteElement("HDiv Trace", cell, 1) compiled_element, module = ffc.codegeneration.jit.compile_elements([element]) ``` The above fails at `A = e.element.A` in...
This PR adds a divergence conforming discontinuous Galerkin demo for the Navier-Stokes equations. It demonstrates: - how to implement DG methods - how to use predefined and custom restriction operators...
Small fix for https://github.com/FEniCS/ufl/pull/120
Currently, UFL checks that all integrals of a particular type in a form have the same subdomain data (using the UFL ID). I'm not sure this should be UFL's job....
A call to NBX was added in [this](https://github.com/FEniCS/dolfinx/pull/3099) PR to fix a bug in sub map creation, however, it should be possible to compute the sub map destination ranks without...