La función MinimumResidualsSolve
fracasa en su implementación para matrices (no virtuales).
Matrix mB = Diag(3, 1);
Matrix mC = Col(1,2,3);
Matrix mX1 = VMat2Mat(MinimumResidualsSolve(Mat2VMat(mB), Mat2VMat(mC)));
//> (1, 2, 3)
Matrix mX2 = MinimumResidualsSolve(mB, mC);
//> (?, ?, ?)
Matrix mB2 = SetMat([[
[[2500, 1264, 1231]],
[[1264, 1264, 1623]],
[[1231, 1623, 1231]]
]]);
Matrix mC2 = Col(2021,1922,546);
Matrix mX3 = VMat2Mat(MinimumResidualsSolve(Mat2VMat(mB2), Mat2VMat(mC2)));
//> (0.69, -1.66, 1.93)
Matrix mX4 = MinimumResidualsSolve(mB2, mC2);
//> (222448.05, -390853.09, 111433.23)
Recuerdo que corregí algo cuando hice la versión de
VMatrix
y seguro que olvidé trasladarlo a Matrix.