if np.ndim(radius) == 0 and result.size == 1:
            return result[0]  # scalar outout