scipy.optimize.minimize 返回零维数组？

Question

Calling the Powell minimizer for a single-dimensional problem is creating an OptimizeResult with an inaccessible value. For example:

from scipy.optimize import minimize
test = minimize(lambda x: 1.0, np.array([1.0]), method="Powell")

If I then ask for test.x I get:

array(3.58792896)

Something is wrong with that "array": I can't get the value out of it. For example, test.x[0] returns IndexError: too many indices for array . It's like it's a zero-dimensional array, or there's some other reference problem.

(A well-formed ndarray would display like array([3.58792896]) .)

What am I doing wrong?

Answer 1

That's a 0-dimensional array, but it shouldn't be. While 0-dimensional arrays are a supported concept in NumPy, that minimize call shouldn't be creating one. It looks like the devs are worried about breaking backward compatibility if they fix this , so a fix is unlikely for now.

I would recommend using numpy.atleast_1d to handle this case consistently with the cases that return a 1D array, and be forward compatible if they do eventually change the return value:

test = minimize(...)
if not test.success:
    handle_that()
result = np.atleast_1d(test.x)

For cases where you expect a 0D array and want to retrieve the stored value, index it with a tuple of 0 indices:

value = zero_d_array[()]