我正在尝试使用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

09-25 21:35