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

#include "mpi.h"

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

int main(int argc, char *argv[])
{
	int num_procs, myrank;
	mpf_t a, b[128];
	void *buf_a, *buf_b;
	int tag = 0, i;
	MPI_Status status;

	MPI_Init(&argc, &argv);

	_mpi_set_bnc_default_prec_decimal(50, MPI_COMM_WORLD);
	commit_mpf(&(MPI_MPF), ceil(50/log10(2.0)), MPI_COMM_WORLD);

	MPI_Comm_size(MPI_COMM_WORLD, &num_procs);
	MPI_Comm_rank(MPI_COMM_WORLD, &myrank);

	mpf_init_set_ui(a, (unsigned long)myrank);
	buf_a = allocbuf_mpf(mpf_get_prec(a), 1);
	buf_b = allocbuf_mpf(mpf_get_prec(a), num_procs);

	pack_mpf(a, 1, buf_a);
	MPI_Gather(buf_a, 1, MPI_MPF, buf_b, 1, MPI_MPF, 0, MPI_COMM_WORLD);

	if(myrank == 0)
	{
		for(i = 0; i < num_procs; i++)
			mpf_init(b[i]);
		unpack_mpf(buf_b, b[0], num_procs);
	}

	printf("Process %d: a = ", myrank);
	mpf_out_str(stdout, 10, 0, a);
	printf("\n");
	if(myrank == 0)
	{
		for(i = 0; i < num_procs; i++)
		{
			printf("b[%d] = ", i);
			mpf_out_str(stdout, 10, 0, b[i]);
			printf("\n");
		}
	}

	free(buf_a);
	free(buf_b);

	mpf_clear(a);
	if(myrank == 0)
	{
		for(i = 0; i < num_procs; i++)
			mpf_clear(b[i]);
	}

	free_mpf(&(MPI_MPF));

	MPI_Finalize();

	return EXIT_SUCCESS;
}