num::mgs_columns

Defined in header <linear/subspace.hpp>

template <std::floating_point T>
void mgs_columns(T *v, const T *basis, idx ldb, idx rows, idx columns,
                 T *coefficients = nullptr) noexcept;

Modified Gram--Schmidt against row-major basis columns.

Keeps the sequential projection order; project_columns then combine_columns is the faster classical variant, with weaker stability.