geometry - LAPACK: Solving large periodic banded system of equations -


i have solve large number of simultaneous equations (~1000s) solve @ every time step general mean curvature flow problem. problem defined on closed manifolds boundary condition periodic.

i using successive-over-relaxation algorithm right solve this, slow. tried dgbtrf -> dgbtrs (without periodicity condition), , quite faster.

the coefficient matrix looks this

     ⎛c₁   d₁ e₁        a₁   b₁⎞     ^      ⎢b₂   c₂ d₂ e₂  0       a₂⎥     |      ⎢a₃   b₃ c₃ d₃  .  0      ⎥     |  ← ⎢     a₄ b₄ c₄  .  .      ⎥   ~1000      ⎢      0  .  .  .  .  en₋₂⎥     |      ⎢en₋₁     0  .  .  .  dn₋₁⎥     |      ⎝dn   en        bn  cn ⎠     v 

i need solve pentadiagonal systems, not symmetric , not known positive definite.

is there way solve cyclic/periodic banded systems in lapack?

or have use general solvers, such dgetrs?


Comments

Popular posts from this blog

Load Balancing in Bluemix using custom domain and DNS SRV records -

oracle - pls-00402 alias required in select list of cursor to avoid duplicate column names -

python - Consider setting $PYTHONHOME to <prefix>[:<exec_prefix>] error -