Hi, and thanks again for this great package. I use complex numbers and I was getting crazy outputs until I found the following.
The _MPI.DOUBLE used to inform the block size in https://github.com/scientificcomputing/fenicsx_ii/blob/main/src/fenicsx_ii/interpolation.py#L124 is not suited to dtype=np.complex128. The effect is that recv_basis_functions is a shuffled version of basis_values_on_V.
I have proposed a fix on:
https://github.com/pierricmora/fenicsx_ii/blob/complex/src/fenicsx_ii/interpolation.py#L128