#!/usr/bin/perl -w -I/accounts/darin/current/blib/arch -I/accounts/darin/current/blib/lib

print "1..2\n" if not @ARGV;
use PDL;
use PDL::Parallel::MPI;
mpirun(4);

MPI_Init();
$rank = get_rank();
$a=($rank+1) * sequence(4);


$b = $a->gather(root=>1);
print "rank = $rank, b=$b" if $rank == 1;
print "ok 1\n"  if $rank ==1 && $b->at(1,3) == 4;
print "nok 1\n" if $rank ==1 && $b->at(1,3) != 4;

$c = zeroes(16);
$a->gather( dest => $c );
print "rank = $rank, c=$c\n" if $rank == 0;
print "ok 2\n"  if $rank ==0 && $c->at(14) == 8;
print "nok 2\n" if $rank ==0 && $c->at(14) != 8;

MPI_Finalize();