HW12: Conjugate Gradients

The purpose is to solve Ax=b by using the conjugate gradients algorithm. A is a known n by m matrix. b is a known n by 1 vector. x is an unknown m by 1 vector. We want to solve for x. 

The program will standard-out the solution vector x. The file called 'filename' contains two file names: Matfile (for A) and Vecfile (for b). The file io.c let us read a matrix A and a vector b, and create a zero vector x. The conjugate gradient algorithm functions are in fcg.c.

The program uses gsl library. The gsl/gsl_blas.h is included in the header file.

./rebecca_hw12 < filename > out