#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#include "bnc.h"

#define DIM 5 

int main()
{
	long int i, j;
	DMatrix c, x, b;

	x = init_dmatrix(DIM, DIM);
	b = init_dmatrix(DIM, DIM);
	c = init_dmatrix(DIM, DIM);

	for(i = 0; i < DIM; i++)
	{
		for(j = 0; j < DIM; j++)
		{
			set_dmatrix_ij(x, i, j, (double)(i * DIM + j + 1));
			set_dmatrix_ij(b, i, j, (double)(DIM * DIM - (i * DIM + j)));
		}
	}

	add_dmatrix(c, x, b);

	print_dmatrix(c);

	free_dmatrix(x);
	free_dmatrix(b);
	free_dmatrix(c);

	return EXIT_SUCCESS;
}