BUG: scipy.spatial.Delaunay, scipy.interpolate.LinearNDInterpolator segfault on invalid shapes about scipy

[andfoy]

This issue is not related to #20619, since it is not related to concurrency. However, I can reproduce it, which implies that some other kind of error is occurring, I can further investigate this one.

[AtsushiSakai]

In the interpolate case, the LinearNDInterpolator is using the qhull.Delaunay in its initializer.

scipy/scipy/interpolate/interpnd.pyx

Lines 329 to 330 in 4a537b7

	def _calculate_triangulation(self, points):
	self.tri = qhull.Delaunay(points)

In the latest code, I have confirmed that using the code reported in this issue does not cause a segmentation fault, thanks to the validation introduced in PR #20633:

In [1]: import numpy as np

In [2]: from scipy.interpolate import LinearNDInterpolator

In [3]: LinearNDInterpolator(np.random.randn(4, 10, 3), np.random.randn(10))
---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
Cell In[3], line 1
----> 1 LinearNDInterpolator(np.random.randn(4, 10, 3), np.random.randn(10))

File interpnd.pyx:326, in scipy.interpolate.interpnd.LinearNDInterpolator.__init__()

File interpnd.pyx:92, in scipy.interpolate.interpnd.NDInterpolatorBase.__init__()

File interpnd.pyx:330, in scipy.interpolate.interpnd.LinearNDInterpolator._calculate_triangulation()

File _qhull.pyx:1808, in scipy.spatial._qhull.Delaunay.__init__()

ValueError: Input points array must have 2 dimensions.

So, I think we can close this issue.

[lucascolley]

I see Qhull and segfault - @andfoy related to gh-20619, maybe?

[tylerjereddy]

As the cross-linked PR indicates, we just need to error out for the spatial cases since those are invalid inputs. I didn't touch the interpolate case, though it may or may not benefit from my patch if it passes through the same code path (didn't check).