5903 Intro Scientific Prog:
QR Factorization: Specification

  1. Write a subroutine to compute the QR factorization of a symmetric matrix A using Householder refections. You do not need to implement a new matrix class for full symmetric matrices but you can if you want. I looked at the wikipedia page http://en.wikipedia.org/wiki/QR_decomposition and it looks OK.
  2. You must organize your code into seperate files and make the "makefile" reasonably compact.  Provide some estimate of the operation count in the header of the procedure.

Additions and changes are in red.