2
Fortunately, LAPACK provides routines to deal with the $\mathbf Q$ factor from the $\mathbf A = \mathbf Q \mathbf R$ decomposition, [dgeqrf]. To find the projection of an arbitrary $\mathbf B$ onto the space orthogonal to $ \mathrm {range}(\mathbf A)$, you want to form $\mathbf C = \left(\mathbf I - \mathbf Q \mathbf Q^T\right) \mathbf B$. Here are two ...
Only top voted, non community-wiki answers of a minimum length are eligible
Related Tags
lapack × 118linear-algebra × 53
blas × 26
eigenvalues × 24
matrix × 20
fortran × 12
linear-solver × 10
c++ × 10
eigensystem × 10
matrix-factorization × 8
scipy × 6
matlab × 5
numerics × 5
c × 5
eigen × 5
parallel-computing × 4
dense-matrix × 4
banded-matrix × 4
python × 3
sparse-matrix × 3
performance × 3
intel-mkl × 3
scalapack × 3
libraries × 2
mpi × 2