#include "header.h"
int ifsym(gsl_matrix *A)// This function is ued to check if the matrix is symmetric or not
{
	int i,j;
	int n=A->size1;
	float tmp1;
	float tmp2;
	printf("n====%d\n",n);
	for(i=1;i<n;i++)
	{
	   for(j=i+1;j<n;j++)
	   {
		tmp1=gsl_matrix_get(A,i,j);
		tmp2=gsl_matrix_get(A,j,i);
		if(tmp1!=tmp2) return 0;
	   }
	}
	return 1;
}

void vsub(gsl_vector *a,gsl_vector *b,gsl_vector *ans)
{
	gsl_vector_memcpy(ans,a);
	gsl_vector_sub(ans,b);
}