Suggesting a minor enhancement in mpdist.py:: _mpdist_vect #1082
Description
I think we can replace
with
P_ABBA = np.empty(2 * j, dtype=np.float64)
if k is None:
percentage = np.clip(percentage, 0.0, 1.0)
k = math.ceil(percentage * (2 * Q.shape[0])) # BECAUSE `P_ABBA.shape[0] == 2 * j`, see the first line
k = min(int(k), P_ABBA.shape[0] - 1)