Linear Algebra and the C Language/a0b1
QR Decomposition
- orth_mR(U,Q,YES); Orthogonalizes and normalizes the U matrix.
/* ------------------------------------ */
/* ------------------------------------ */
void QR_mR(
double **U,
double **Q,
double **R)
{
double **Q_T = i_mR(csize_R(Q),rsize_R(Q));
orth_mR(U,Q,YES);
transpose_mR(Q,Q_T);
mul_mR(Q_T,U, R);
f_mR(Q_T);
}
/* ------------------------------------ */
/* ------------------------------------ */