我已经开始与JBLAS合作,但我遇到了一个问题,
doubleMatrix x = Solve.solve(A,b);
这给结果很好。但是如果我想重新计算以找到 Ax = b1 ,那么它将再次重新分解以获得 A = LU 或者 A = QR 。但这是多余的,因为 A 矩阵保持不变。
我想用它来分解A:
Decompose.LUDecomposition<DoubleMatrix> out = Decompose.lu(A);
然后将生成的对象传递给solve函数以获得所需的输出。使用这种方式,如果仅 b 值发生变化,则不需要再次进行分解。
有可能吗?请帮忙。或者我是从头开始编写该功能..