我必须解决大量的联立方程（~1000s），以便在每个时间步长求解一般平均曲率流问题。问题是在闭合流形上定义的，因此边界条件是周期性的。
我现在正在使用连续过度放松算法来解决这个问题，但速度很慢。我尝试dgbtrf -> dgbtrs（没有周期性条件），并且速度非常快。
系数矩阵看起来像这样
⎛c₁ d₁ e₁ a₁ b₁⎞ ^
⎢b₂ c₂ d₂ e₂ 0 a₂⎥ |
⎢a₃ b₃ c₃ d₃ . 0 ⎥ |
A ← ⎢ a₄ b₄ c₄ . . ⎥ ~1000
⎢ 0 . . . . en₋₂⎥ |
⎢en₋₁ 0 . . . dn₋₁⎥ |
⎝dn en an bn cn ⎠ v
我需要解决五边形系统，它们不对称且不为正定。
有没有办法在LAPACK中解决循环/周期性带状系统？
或者我是否必须使用常规解算器，例如dgetrs？