num::cholesky_solve

Defined in header <linear/factorization/cholesky.hpp>

template <std::floating_point T>
void cholesky_solve(T *x, const T *L, const T *b, idx n) noexcept;

Solve \(A x = b\) from a Cholesky factor, via \(L y = b\) then \(L^T x = y\).

x and b may alias.