blob: c9e163852c934e8fa0d7fc9bd54fc4f78a80368e (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
#pragma once
#include "matrix.hpp"
#include "vector.hpp"
namespace linalg {
struct QRResult {
Matrix Q;
Matrix R;
};
// Classical Gram-Schmidt.
// Mathematically natural but numerically fragile: orthogonality of Q
// degrades rapidly on ill-conditioned inputs.
// Provided for comparison — prefer modified_gs or householder in practice.
//
// Throws DimensionMismatchError if rows < cols.
// Throws SingularMatrixError if a column is (nearly) linearly dependent.
QRResult qr_classical_gs(const Matrix& A, double zero_tolerance = 1e-14);
// Modified Gram-Schmidt.
// Subtracts each projection immediately on the running vector rather than
// on the original column. Algebraically equivalent to classical GS but
// numerically much better — round-off stays local instead of accumulating.
//
// Same exceptions as classical GS.
QRResult qr_modified_gs(const Matrix& A, double zero_tolerance = 1e-14);
// Householder QR.
// Applies a sequence of orthogonal reflections to zero out below-diagonal
// entries column by column. Backward-stable and the standard choice for
// dense QR. Works correctly on rank-deficient matrices (zero pivots
// produce zero diagonal entries in R without throwing).
//
// Throws DimensionMismatchError if rows < cols.
QRResult qr_householder(const Matrix& A);
} // namespace linalg
|