I wonder if there is way to handle sparse nd-arrays ?
I am looking for 3D and 4D arrays to store coeficients for geometric transformation 2D -> 1D or 2D -> 2D,
with typical shape (4000x4000)x(2000x500).
My idea is to contruct the array from cython (what is the time consuming part) than re-use it many times (with something similar to a scalar product.
scipy.sparse seems limited to matrices (i.e 2D-array) is this true ? Flattening each array would be possible but there is probably a better solution... Thanks