HouseholderQR.h

using namespace boost::numeric;
/* Ist  ublas::matrix<T> Basisklasse v. matrix_sparse? Nein, aber es scheint, dass
dass boost-Funktionen spars u. dense koennen u. auch mischen koennen. Vermutlich
sind dort die Matrixargumente als Templates implementiert.*/
template< class v, class M>//class T>
 int backwardSubstitution ( v vec, M R);//ublas::vector<T>& vec,   const ublas::matrix<T>& R)

//template<class T>
//int backwardSubstitution ( ublas::vector<T>& vec, 
//                  const ublas::matrix<T>& R);

 template<class T>
 int householderQR (const ublas::matrix<T>& M, 
           ublas::matrix<T>& Q, 
           ublas::matrix<T>& R);

 template<class T>
 int householderQR (const ublas::mapped_matrix<T>& M, 
           ublas::mapped_matrix<T>& Q, 
           ublas::mapped_matrix<T>& R);