有关使 pardiso 更准确的任何提示？
对于牛顿线搜索，我有一个（对称）Hessian 的（条件不好的）线性系统（附加的线性约束设置为拉格朗日乘数，所以它不是正定的）。求解
- 使用 matlab 的 mldivide 需要 6 秒；
- 使用 pardiso (pardiso-project.org) 需要 1.7 秒。
但是，pardiso 的解决方案（使用默认设置时）非常不准确，并且行搜索停止。
设置以下内容显着改善了解决方案，使其与 mldivide 一样好：
info.iparm(13) = 2
我想知道是否还有其他相关参数。例如，我试图设置
pinfo.iparm(4) = 31
和其他值，但 pardiso 报告了 -4 类型的错误。
Ofir 指出了附录 D
(19) 一种快速局部内射谐波映射的子空间方法