对于一般平均曲率流问题，我必须在每个时间步求解大量联立方程（约 1000 秒）。该问题是在封闭流形上定义的，因此边界条件是周期性的。
我现在正在使用连续过度松弛算法来解决这个问题，但速度很慢。我试过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？