pardiso をより正確にするためのヒントはありますか?
ニュートンライン検索用に、（対称）ヘッシアンの（条件の悪い）線形システム（追加の線形制約がラグランジュ乗数として設定されているため、正定値ではありません）があります。解く
- matlab の mldivide を使用すると 6 秒かかります。
- pardiso (pardiso-project.org) を使用すると、1.7 秒かかります。
ただし、pardiso のソリューション (デフォルト設定を使用する場合) は非常に不正確であり、行検索は失速します。
以下を設定すると、ソリューションが劇的に改善され、mldivide と同程度になりました。
info.iparm(13) = 2
他に関連するパラメータがあるかどうか疑問に思っていました。たとえば、私は設定しようとしました
pinfo.iparm(4) = 31
およびその他の値ですが、pardiso はタイプ -4 のエラーを報告しました。
Ofir は付録 D を指摘した。
(19) 高速局所単射調和写像のための部分空間法