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

#include "bnc.h"

#define DIM 5

int main(int argc, char *argv[])
{
    DMatrix dmat_ans, dmat1, dmat2;
    long int i, j;
	double start_wtime, end_wtime;

    dmat_ans = init_dmatrix(DIM, DIM);
    dmat1 = init_dmatrix(DIM, DIM);
    dmat2 = init_dmatrix(DIM, DIM);
    for(i = 0; i < DIM; i++)
	{
        for(j = 0; j < DIM; j++)
        {
            set_dmatrix_ij(dmat1, i, j, (double)(i*DIM + j + 1));
            set_dmatrix_ij(dmat2, i, j, (double)(DIM * DIM - (i*DIM +j)));
        }
    }

    mul_dmatrix(dmat_ans, dmat1, dmat2);

    printf("BNC:\n"); print_dmatrix(dmat_ans);

    free_dmatrix(dmat_ans);
    free_dmatrix(dmat1);
    free_dmatrix(dmat2);

	return EXIT_SUCCESS;
}