void matvec(const Matrix &A, const Vector &x, Vector &y)
void matmul(const Matrix &A, const Matrix &B, Matrix &C, idx block_size)
constexpr T ipow(T x) noexcept
Compute x^N at compile time via repeated squaring.
BasicVector< real > Vector
Real-valued dense vector with full backend dispatch (CPU + GPU)