Open
Description
I have the following 🍇 about the current implementation:
- Array-valued interpolations are allocating;
- Different methods for scalar and array interpolation per interpolation type yields code duplication;
- Support for array-valued interpolation is inconsistent.
As far as I'm aware, for every interpolation type the array-valued interpolation can be defined by a per-index scalar interpolation.
I propose:
- Array Interpolation objects get a pre-allocated value, derivative and integral on which computations are done in-place;
- We define a function
get_u(A.u, i)
which gets the appropriate view of a slice ofu
.
Related issues:
Metadata
Metadata
Assignees
Labels
No labels