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);    
}
/* ------------------------------------ */
/* ------------------------------------ */