使用Ojalgo解决Java中的线性系统

Question

我试图用Ojalgo解决M（NxN）线性系统（Ax = B，B = [b1，b2，... bM]）。最有效的方法是什么？我也想知道A是否是单数（A，B是PrimitiveMatrix类型的对象）。

非常感谢任何帮助。谢谢！

Answer 1

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维基吗？

https://github.com/optimatika/ojAlgo/wiki