我正在尝试使用 Ojalgo 解决 M (NxN) 线性系统 (Ax = B, B = [b1,b2,...bM])。最有效的方法是什么?我还想知道 A 是否是单数(A,B 是 PrimitiveMatrix 类型的对象)。
任何帮助将不胜感激。谢谢!
我正在尝试使用 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 维基吗?