Describe the issue:

We hit a new test failure in astropy's test suite against numpy 2.3.0's nightlies. Namely

xref astropy/astropy#18045

We were able to bisect the regression to 461d1c0 (#28624)

I apologies for not providing a more minimized and self-contained reproducer but the error we get is a bit perplexing to us, and seeing that numpy 2.3.0rc1 looks imminent, we figured it might be worth reporting now anyway, in case you guys would consider reverting the breaking PR.
I should hopefully be able to reduce the example by a lot on next business day.

Reproduce the code example:

❯ pytest astropy/cosmology/_src/tests/funcs/test_funcs.py::test_z_at_value_numpyvectorize

Error message:

Python and NumPy Versions:

Python 3.13 + numpy 2.3.0dev

Runtime Environment:

No response

Context for the issue:

For context, we currently have this test disabled for numpy >= 2.3.0dev already, so even if more inquiry is needed, it's not extremely pressing on our end.