num::csr_lu_solve

Defined in header <linear/solvers/ilu.hpp>

template <std::floating_point T, std::integral Index>
void csr_lu_solve(T *x, const T *val, const Index *row_ptr, const Index *col_idx,
                  const Index *diagonal, const T *b,
                  std::type_identity_t<Index> n) noexcept;

Solve \(LUx = b\) for factors packed by ilu0_factor.

Forward substitution against the implicit unit-diagonal L, then backward substitution against U. x may alias b.