#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[128], b;
	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);

	if(myrank == 0)
	{
		for(i = 0; i < num_procs; i++)
			mpf_init_set_ui(a[i], i);
		buf_a = allocbuf_mpf(mpf_get_prec(a[0]), num_procs);
		pack_mpf(a[0], num_procs, buf_a);
	}
	mpf_init(b);
	buf_b = allocbuf_mpf(mpf_get_prec(b), 1);

	MPI_Scatter(buf_a, 1, MPI_MPF, buf_b, 1, MPI_MPF, 0, MPI_COMM_WORLD);
	unpack_mpf(buf_b, b, 1);

	printf("Process %d: b =", myrank);
	mpf_out_str(stdout, 10, 0, b);
	printf("\n");

	MPI_Finalize();

	return EXIT_SUCCESS;
}