Our push-forward operators take a callable function in the logical domain and an evaluation grid in the logical coordinates, and return the values of the push-forwarded function on the corresponding physical grid. To this end they can take advantage of several optimizations as shown in file psydac/core/field_evaluation_kernels.py.
These kernels are used by the high-level class Pushforward in module psydac.feec.pushforward: the __init__ constructor takes a fairly generic evaluation grid, and the __call__ method evaluates the correct push-forward (H1, H-curl, H-div, or L2) of multiple fields based on the kind parameter of their spaces.
Unfortunately, the functions push_* in module psydac.feec.pull_push provide an alternative interface which is less efficient, but extensively used across Psydac. In particular, these functions are based on the __call__ method of FemField objects, which in turn calls the method eval_field of FemSpace.
In general eval_field is implemented in the subclasses of FemSpace for a single-point evaluation, and it is mostly in pure Python.
Our push-forward operators take a callable function in the logical domain and an evaluation grid in the logical coordinates, and return the values of the push-forwarded function on the corresponding physical grid. To this end they can take advantage of several optimizations as shown in file
psydac/core/field_evaluation_kernels.py.These kernels are used by the high-level class
Pushforwardin modulepsydac.feec.pushforward: the__init__constructor takes a fairly generic evaluation grid, and the__call__method evaluates the correct push-forward (H1, H-curl, H-div, or L2) of multiple fields based on thekindparameter of their spaces.Unfortunately, the functions
push_*in modulepsydac.feec.pull_pushprovide an alternative interface which is less efficient, but extensively used across Psydac. In particular, these functions are based on the__call__method ofFemFieldobjects, which in turn calls the methodeval_fieldofFemSpace.In general
eval_fieldis implemented in the subclasses ofFemSpacefor a single-point evaluation, and it is mostly in pure Python.