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

#define USE_GMP
#define USE_MPFR
#include "bnc.h"

main()
{
	MPFVector x, y;
	MPFMatrix a;

	set_bnc_default_prec_decimal(100);

	a = init_mpfmatrix(2, 2);
	x = init_mpfvector(2);
	y = init_mpfvector(2);

	set_mpfmatrix_ij_d(a, 0, 0, 1.0);
	set_mpfmatrix_ij_d(a, 0, 1, 2.0);
	set_mpfmatrix_ij_d(a, 1, 0, 3.0);
	set_mpfmatrix_ij_d(a, 1, 1, 4.0);

	set_mpfvector_i_d(x, 0, 1.0);
	set_mpfvector_i_d(x, 1, 2.0);

	print_mpfmatrix(a);
	print_mpfvector(x);

	mul_mpfmatrix_mpfvec(y, a, x);

	print_mpfvector(y);

	free_mpfmatrix(a);
	free_mpfvector(x);
	free_mpfvector(y);

	return EXIT_SUCCESS;
}