用Ojalgo语言求解线性系统

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