#include "header.h"

int main(void)
{
	gsl_matrix_float A;	//We are trying to solve x in the Equation Ax=b, where A is a matrix and x and b are vectors
	gsl_vector_float b;
	gsl_vector_float x;
	gsl_vector_float ans;
	int VectorLength;	//Length of the read vector that can later be used to create the vector


	// Read in matrix A and the vector b
	A=ReadMatrix();
	b=ReadVector();

	VectorLength=A.size2;		//The vector x has the size equal to the number of columns of A

	// Create vector x
	x=CreateVector(VectorLength);  

	ConjugateGradient(A,x,b,TOLERANCE,MAXITERATIONS);	// call Conjugate Gradient function
	gsl_vector_float_fprintf(stdout,&x,"%lf");		//print the output

	return 0;
}