Skip to content

num::lu_factor_blocked

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

template <std::floating_point T, class Index>
bool lu_factor_blocked(T *LU, Index *piv, idx n, idx block_size = 64) noexcept;

Blocked variant of lu_factor: same packed \(PA = LU\) contract, panel-at-a-time.