I'm writing code that requires 2D polynomial interpolation. I realized while working that, contrary to the documentation, functions like polygrid2d or leggrid2d do NOT compute cartesian product of x ...