5903 Intro Scientific Prog:
Conjugate Gradients: Specification

1. The program should read in a matrix (A) and a RHS vector (b) from a file
• Use the redirect from standard-in. 
• Use the redirect from standard-out
2. The output should be the solution vector (x satisfying A x = b) printed to standard out.  
3. The code should use the conjugate gradients algorithm.
4. The code should be able to be run as "MyProg < in > out".  Where:
• in is a file containing MatrixFileName VectorFileName. You can assume th file names are less than 80 characters long.  Remember no magic numbers!
• out will be a file containing the solution vector in the matrix format
5. The matrix and vector storage and file formats should be as defined in class. 
6. The i/o functions should be in one file.  
7. The conjugate gradient algorithm functions should be in another.
8. You need to use either implicit rules or automatic variables in your makefile.
9. Main should be brief and well documented.  
10. We will be exchanging this code.  
11. As always ask in class if you do not understand what the specifications say.
12. You do not need (but it is a good idea) to check that the matrix and vector conform.
13. You can either write your own linear algebra or use a blas.