Geoffrey Irving <[email protected]> writes:
> For anything other than MPI_SUM, the evaluation points would likely
> want to be different, so the structure of the code would be
> significantly changed.  For example, MPI_MAX probably wants evaluation
> points at the vertices, which can't be optimally coded using the same
> code structure used for MPI_SUM (since different elements share
> vertices).

That's fair, though extrema can also occur internally to an element.

> The main difference is that the quadrature rule has to be specified
> independent of the PetscFE objects, since the PetscFE elements for the
> different fields might have incompatible quadrature rules.

How is that different from evaluating coupling terms between fields in
different spaces?

Attachment: pgprmvZcEar5q.pgp
Description: PGP signature

Reply via email to