Proposed new feature or change:

I use the functions np.polynomial.polynomial.polyval2d and np.polynomial.polynomial.polygrid2d on pint arrays. Pint makes use of the array_function protocol to make the standard numpy function unit-aware. For that purpose I used to add the array_function protocol to the mentioned functions myself by applying array_function_dispatch from numpy.core.overrides. But the entire core module was recently made private so I suspect that this kind of monkey-patching is not supported any longer.