5-d data in tensor basis not interpolateable #163
Description
For 5-dimensional data, blob/main/src/postgkyl/data/computeInterpolationMatrices.py lacks an implementation to interpolate 5-dimensional data in tensor basis. For <= 4-d data, such an implementation is present. Trying to interpolate the data gives the following error message:
NameError: interpMatrix: Basis tensor is not supported!
Supported basis are currently 'nodal Serendipity', 'modal Serendipity', and 'modal maximal order'
This is of particular concern since neutral distributions in 2x2v simulations are output as 5-d data in this basis. The lack of such an implementation makes it difficult to check whether this distribution is correct, which is relevant given the recent change in definition of neutral velocity coordinates.
As a workaround, pretending the data is in ms basis works so long as p=1 with pgkyl sn7-Ar0_0.gkyl interp -b ms -p1, but it would be nice if a more structural solution could be found.
A minimal example can be found here: https://drive.google.com/file/d/1AhOt_Q3iXKPryrYBeBuzfqnbpr3q6dfc/view?usp=sharing