20 int main(
int argc,
char* argv[])
24 matrix<double>
A (3,3);
34 cout <<
"A=" << A << endl;
36 cout <<
"QR decomposition using Householder" << endl;
37 matrix<double> Q(3,3), R(3,3);
39 matrix<double> Z = prod(Q,R) -
A;
41 cout <<
"Q=" << Q <<endl;
42 cout <<
"R=" << R << endl;
43 cout <<
"|Q*R - A|=" << f << endl;
void HouseholderQR(const ublas::matrix< T > &M, ublas::matrix< T > &Q, ublas::matrix< T > &R)
int main(int argc, char *argv[])