#include "header.h"

typedef struct 
{ 
 double *  p; 
 double ** pp; 
 int       m; 
 int       n; } matrix;

matrix Matrix_Create(int m, int n){
	matrix A;
	double *p;
	double **pp;
	int i; 
	
	//Assigning storage for Matrix
	p  = (double *)  malloc(sizeof(double)  *m * n);
	pp = (double **) malloc(sizeof(double *)*m    );
	
	//Configuring pointers for rows
	for(i=0;i<m;i++){
		pp[i] = p + i*n;
		}
	
	//Setting the appropriate values
	A.m = m; A.n = n; A.p = p; A.pp = pp;
	
	return A;	
	}
	
void Matrix_Destroy(matrix A){
	free(A.p);
	free(A.pp);	
	}
	
void Matrix_Zero(matrix A){
	int i, j;
	for(i=0;i<A.m;i++){
		for(j=0;j<A.n;j++){
			A.pp[i][j]=0.0;
			}
		}
	}

void Matrix_Print(matrix A){
	int i, j;
	
	printf("\n");
	for(i=0;i<A.m;i++){
		for(j=0;j<A.n;j++){
			printf("%e ",A.pp[i][j]);
			}
		printf("\n");
		}
	printf("\n");
	}
	
void Matrix_Read(char * FileName, matrix A){
	FILE *fp;
	int i;
	
	fp = fopen(FileName,"r");
	
	for(i=0;i<A.m*A.n;i++){
		fscanf(fp, "%lf", A.p + i);
		}
	
	fclose(fp);
	}