Использование Ojalgo для решения линейных систем в Java
Я пытаюсь решить M (NxN) линейных систем (Ax = B, B = [b1,b2,...bM]), используя Ojalgo. Какой самый эффективный способ сделать это? Я также хотел бы знать, является ли A единственным (A, B являются объектами типа PrimitiveMatrix).
Любая помощь будет принята с благодарностью. Спасибо!
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 вики?