`basis` argument for calculations

[PrzeChoj]

For now, everything works on a natural orthonormal basis $I_p$.

Users may be interested in getting the decomposition for different orthonormal basis.