Unexpected output when using np.meshgrid with indexing='ij' in NumPy 1.24
I'm not sure how to approach I'm having trouble with I'm upgrading from an older version and I'm working with `np.meshgrid` in NumPy 1.24 and trying to generate a grid of coordinates for a 3D plot..... However, I am getting unexpected shapes in the output when I set the `indexing` parameter to 'ij'. I expected the output arrays to match the shape of the input arrays, which are of size (4,) and (3,) respectively. Hereโs what Iโve done: ```python import numpy as np x = np.array([1, 2, 3, 4]) y = np.array([5, 6, 7]) X, Y = np.meshgrid(x, y, indexing='ij') print(X.shape) print(Y.shape) ``` I anticipated that `X` would have the shape (4, 3) and `Y` would also have the same shape (4, 3), but instead, I'm getting shapes (4, 4) for `X` and (3, 4) for `Y`. The output is: ``` (4, 4) (3, 4) ``` Iโve checked the documentation, and I think Iโm following the usage correctly. I also tried switching to the default `indexing='xy'`, and it behaved as expected, producing shapes (4, 3) for `X` and (4, 3) for `Y`. Is there something I'm missing when using 'ij'? I would appreciate any insights into why the shapes are not what I expected when using `indexing='ij'` with `np.meshgrid`. Thanks for taking the time to read this! I'm on Windows 11 using the latest version of Python. Any feedback is welcome!