Как повторно использовать разложенную матрицу LU для решения Ax = b в JBLAS
Я начал работать с JBLAS, но столкнулся с проблемой,
doubleMatrix x = Solve.solve(A,b);
Это дает результаты просто отлично. Но если я захочу сделать перерасчет, чтобы найти Ax = b1, то он снова разложится, чтобы получить A = LU или, возможно, A = QR. но это избыточно, так как матрица A остается неизменной.
Я хочу, чтобы это разложить A с помощью этого:
Decompose.LUDecomposition<DoubleMatrix> out = Decompose.lu(A);
а затем передать сгенерированный объект в функцию решения, чтобы получить желаемый результат. Использование этого способа не потребует повторной декомпозиции, если изменяется только значение b.
Является ли это возможным?? Пожалуйста помоги. Или я написать эту функцию с нуля..