#include "Header.h"

//checks if Matrix A is symmetric
//  
//  A is symmetric if A is a square matrix (m by m) &&
//  A is symmetric if (aij == aji) of all a's in A
//
//returns 0 for false  (NOT symmetric)
//returns 1 for true   (symmetric)
int isMatSymm(const Mat * A)
{
	int i, j;
	
	//check if A is square
	if(A->n != A->m)
		return 0;
	
	//check if aij == aji for all a's in A
	for(i =0; i < A->n; i++)
	{
		for(j =0; j < A->n; j++)
		{
			//if (aij != aji)    NOT symmetric
			if(A->p[i*A->n+j] != A->p[j*A->n+i])
				return 0;
		}
	}
	
	//return symmetric if not failed above
	return 1;
}