#include "header.h"

// check symmetry

int check(gsl_matrix *A)
{
	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;	// A is not symmetric
	   }
	}
	return 1;
}

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