Geoffrey Irving <[email protected]> writes: > Does that mean the *Project* functions do not produce optimal answers?
It is a "nodal" projection rather than an L^2 projection (which requires solving with a mass matrix). The dual space from the Ciarlet definition typically consists of point values at nodes, plus perhaps some averages, gradients at points, etc. Point values are the typical case. The dual space is used to define any continuity between elements. PetscFE does the change of basis from a "prime basis" (something that can be stably evaluated) to the finite element basis defined in terms of the dual space. See the Brenner and Scott or the FIAT paper for further details on this.
pgpXma112CgqD.pgp
Description: PGP signature
