#include "SparseMatrices.h"

void ReadSparseMatrix(char FileName[], sMat *A){
	FILE 	*FilePointer;
	int 	i, j, k, m, n, *Counts, *jVals;
	float 	f, *fVals;
	sRow  	*rp;

	// File OPS
	FilePointer = fopen(FileName,"r");
	fscanf(FilePointer,"5903SparseMatrix\n%i\t%i\nfloat",&m,&n);
	
	//Mem Assignments
	#define MaxNZ 20
	Counts = (int   *) malloc(m*sizeof(int));
	jVals  = (int   *) malloc(m*MaxNZ*sizeof(int));
	fVals  = (float *) malloc(m*MaxNZ*sizeof(float));
	rp     = (sRow  *) malloc(m*sizeof(sRow));
	
	// Inits
	for(i=0;i<m;i++) Counts[i]=0;
	
	// Reading the Data and closing file
	while(fscanf(FilePointer,"%i\t%i\t%f",&i,&j,&f)==3){
		jVals[i*MaxNZ + Counts[i]] = j;
		fVals[i*MaxNZ + Counts[i]] = f;
		Counts[i]++;
		};
	close(FilePointer);
	
	// Assigning values for matrix Dimensions
	A->m = m;
	A->n = n;
	A->sRow = rp;
	
	// Making rows
	for(i=0;i<m;i++){
	 MakeSparseRow(n, Counts[i], jVals+i*MaxNZ, fVals+i*MaxNZ, rp+i);
	 };
}