[petsc-dev] are petscfe tests functions dual or primal?

Jed Brown jed at jedbrown.org
Fri Jan 17 18:18:49 CST 2014


Geoffrey Irving <irving at naml.us> 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.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20140117/b51965c9/attachment.sig>


More information about the petsc-dev mailing list