Improve the ways to evaluate `FemField`s

pyccel / psydac

Python 3 library for isogeometric analysis

https://pyccel.github.io/psydac/

MIT License

53 stars 18 forks source link

Open anlavandier opened 2 years ago

anlavandier commented 2 years ago

Currently, there are only two ways of evaluating FemFields, point-wise, and one the grid defined in TensorFemSpace.preprocess.

Ideally the interface should allow for those 4 options:

[x] 1. scalar coordinates
[x] 2. 1D arrays of coordinates -> grid is tensor-product, but refinement is not the same in each cell
[x] 3. 2D arrays of coordinates -> grid is tensor-product, and each cell has the same number of cells
[ ] 4. N-D arrays of coordinates -> unstructured grid

anlavandier commented 2 years ago

1, 2 and 3 are done through TensorFemSpace.eval_fields once #214 is merged. 4 is still a WIP.