Cyclic pentadiagonal compact solve via Sherman-Morrison-Woodbury rank-4 correction. A_cyc = A_np + U*W^T where U = [e_1,e_2,e_{n-1},e_n] and W^T encodes the four corner coupling pairs; A_np is the non-periodic Thomas factorisation in (ffr,...).
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=dp), | intent(out), | dimension(:, :) | :: | du | ||
| real(kind=dp), | intent(in), | dimension(:, :) | :: | u | ||
| real(kind=dp), | intent(in), | dimension(:, :) | :: | u_s | ||
| real(kind=dp), | intent(in), | dimension(:, :) | :: | u_e | ||
| integer, | intent(in) | :: | n_tds | |||
| integer, | intent(in) | :: | n_rhs | |||
| real(kind=dp), | intent(in), | dimension(:, :) | :: | coeffs_s | ||
| real(kind=dp), | intent(in), | dimension(:, :) | :: | coeffs_e | ||
| real(kind=dp), | intent(in), | dimension(:) | :: | coeffs | ||
| real(kind=dp), | intent(in), | dimension(:) | :: | ffr | ||
| real(kind=dp), | intent(in), | dimension(:) | :: | faf | ||
| real(kind=dp), | intent(in), | dimension(:) | :: | fsa | ||
| real(kind=dp), | intent(in), | dimension(:) | :: | fbw | ||
| real(kind=dp), | intent(in) | :: | beta_lhs | |||
| real(kind=dp), | intent(in) | :: | beta_lhs_s | |||
| real(kind=dp), | intent(in) | :: | alp | |||
| real(kind=dp), | intent(in) | :: | bet |