Ojalgoを使ってJavaの線形システムを解く

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のwikiを見ましたか？

https://github.com/optimatika/ojAlgo/wiki