Function rgsl::linear_algebra::solve_symm_cyc_tridiag

pub fn solve_symm_cyc_tridiag(diag: &VectorF64, e: &VectorF64, b: &VectorF64, x: &VectorF64) -> Value

This function solves the general N-by-N system A x = b where A is symmetric cyclic tridiagonal (N >= 3). The cyclic off-diagonal vector e must have the same number of elements as the diagonal vector diag. The form of A for the 4-by-4 case is shown below,

A = ( d_0 e_0 0 e_3 ) ( e_0 d_1 e_1 0 ) ( 0 e_1 d_2 e_2 ) ( e_3 0 e_2 d_3 )