I do not actually understand your point """`griddata` expects `xi` to be of shape `(M, D)`. My intention was to allow 1D `xi` of shape `(D, )` as well""". You mean you want to evaluate at a single point?

However, I think uniform interface with griddata is of high priority here. The `xi` parameter has exactly the same meaning.

The problem with atleast_2d is that it does not broadcast. If you want to evaluate it on a grid, you need to do a meshgrid+ravel dance on the user side.