我正在尝试使用Ojalgo解决M(NxN)个线性系统(Ax = B,B = [b1,b2,... bM])。最有效的方法是什么?我还想知道A是否为奇数(A,B是PrimitiveMatrix类型的对象)。
任何帮助将不胜感激。谢谢!
最佳答案
PrimitiveMatrix具有固定/受限功能集,您也无法控制操作方式。如果需要选项和控制,请切换为使用PrimitiveDenseStore(或任何MatrixStore实现)。然后,您需要做的就是:
final LU<Double> tmpLU = LU.PRIMITIVE.make();
tmpLU.decompose(A);
if (tmpLU.isSquareAndNotSingular()) {
x = tmpLU.solve(b);
} else {
// Do something else...
}
您是否看过ojAlgo Wiki?
https://github.com/optimatika/ojAlgo/wiki