num::banded_solve

Defined in header <linear/banded/banded.hpp>

template <std::floating_point T, class Index>
void banded_solve(T *x, const T *ab, idx ldab, idx n, idx kl, idx ku,
                  const Index *ipiv) noexcept;

Solve \(Ax = b\) in place from banded_factor, as LAPACK gbtrs.