numpy
diff --git a/‎numpy/lib/function_base.py
Lines changed: 6 additions & 19 deletions b/‎numpy/lib/function_base.py
Lines changed: 6 additions & 19 deletions
@@ -3943,42 +3943,29 @@ def _quantile_ureduce_func(a, q, axis=None, out=None, overwrite_input=False,
         if np.issubdtype(a.dtype, np.inexact):
             indices_above = concatenate((indices_above, [-1]))
 
-        weights_above = indices - indices_below
-        weights_below = 1 - weights_above
-
-        weights_shape = [1, ] * ap.ndim
-        weights_shape[axis] = len(indices)
-        weights_below.shape = weights_shape
-        weights_above.shape = weights_shape
-
         ap.partition(concatenate((indices_below, indices_above)), axis=axis)
 
         ap = np.moveaxis(ap, axis, 0)
-        weights_below = np.moveaxis(weights_below, axis, 0)
-        weights_above = np.moveaxis(weights_above, axis, 0)
         axis = 0
 
+        weights_shape = [1] * ap.ndim
+        weights_shape[axis] = len(indices)
+        weights_above = (indices - indices_below).reshape(weights_shape)
+
         # Check if the array contains any nan's
         if np.issubdtype(a.dtype, np.inexact):
             indices_above = indices_above[:-1]
             n = np.isnan(ap[-1:, ...])
 
-        x1 = take(ap, indices_below, axis=axis) * weights_below
+        x1 = take(ap, indices_below, axis=axis) * (1 - weights_above)
         x2 = take(ap, indices_above, axis=axis) * weights_above
 
-        # ensure axis with q-th is first
-        x1 = np.moveaxis(x1, axis, 0)
-        x2 = np.moveaxis(x2, axis, 0)
-
         if zerod:
             x1 = x1.squeeze(0)
             x2 = x2.squeeze(0)
 
-        if out is not None:
-            r = add(x1, x2, out=out)
-        else:
-            r = add(x1, x2)
+        r = add(x1, x2, out=out)
 
     if np.any(n):
         if zerod: