Ошибка с возможностью scipy.linprog? (A_ub @ x0 <= b_ub).all() имеет значение True ---но--- linprog(np.zeros_like(x0), A ⇐ Python

[Гость]

With numpy, scipy at versions

Код: Выделить всё

numpy 1.25.0
scipy 1.11.0

the following

Код: Выделить всё

scipy.optimize.linprog

call,

Код: Выделить всё

import numpy as np
from scipy.optimize import linprog

A_ub = np.array(
[[-0.15729144,  0.29943807,  0.29311432],
[-1.32475528, -2.1125364 , -1.55138585],
[ 1.00861965,  0.53283629, -0.14939833],
[ 1.07581479,  0.164022  , -1.19889684]])

b_ub = -np.ones(4)

print(linprog(np.zeros(3),
A_ub=A_ub,
b_ub=b_ub))

return unfeasible status,

Код: Выделить всё

       message: The problem is infeasible. (HiGHS Status 8: model_status is Infeasible; primal_status is At lower/fixed bound)
success: False
status: 2
fun: None
x: None
nit: 0

But the problem is actually feasible, since

Код: Выделить всё

x0 = np.array([ 229.1748166 , -507.05266751,  512.14005547])
print('x0 is feasible?', (A_ub @ x0 

Источник: [url]https://stackoverflow.com/questions/78133436/bug-with-scipy-linprog-feasibility-a-ub-x0-b-ub-all-is-true-but[/url]