Functions solve_cyc_tridiag_sym¶
-
template<class
vec_t
, classvec2_t
, classvec3_t
, classvec4_t
, classmem_t
, classmem_vec_t
>
voido2scl_linalg
::
solve_cyc_tridiag_sym
(const vec_t &diag, const vec2_t &offdiag, const vec3_t &b, vec4_t &x, size_t N, mem_t &m)¶ Solve a symmetric cyclic tridiagonal linear system with user specified memory.
This function solves the system \( A x = b \) where \( A \) is a matrix of the form
* * diag[0] offdiag[0] 0 ..... offdiag[N-1] * offdiag[0] diag[1] offdiag[1] ..... * 0 offdiag[1] diag[2] * 0 0 offdiag[2] ..... * ... ... * offdiag[N-1] ...
See EngelnMullges96 .