num::lu_solve

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

template <std::floating_point T, class Index>
void lu_solve(T *x, const T *LU, const Index *piv, const T *b, idx n) noexcept;

Solve \(A x = b\) from a packed \(PA = LU\) factorization.