Bug in convexify eigen-reflect?

[jgillis]

from casadi import *
A = blockcat([[4,-4,0],[-4,10,0],[0,0,2]])
print(A,np.linalg.eig(A))
print(evalf(convexify(A,{"strategy":"eigen-reflect"})))
print(evalf(convexify(sparsify(A),{"strategy":"eigen-reflect"})))

[jgillis]

Reproduces in tri-diagonalize ->un-tri-diagonalize roundtrip, disabling all other algorithms.

[jgillis]

Bug dissappears for any eps!=0 for A = blockcat([[4,-4,eps],[-4,10,0],[eps,0,2]]).
Caused by sigma==0 clause in casadi_cvx_house

[victorfors]

There seems to be a typo in this fix in casadi_cvx.hpp

// Note: third edition has *beta = 0
// Note: fourth edition has *beta = -2*(v0<0)
*beta = 2*(v0<0);

Missing the minus sign *beta = -2*(v0<0);