subroutine run(spec, rc)
type(ComboSpec), intent(in) :: spec
integer, optional, intent(out) :: rc
integer :: status
real :: tot_time
real :: tot_time_write
real :: tot_time_gather
real :: avg_time
real :: avg_time_write
real :: avg_time_gather
type(GathervKernel) :: kernel
type(BW_Benchmark) :: benchmark
integer :: writer_comm
integer :: gather_comm
integer :: i
real :: ta, tb
integer :: color, rank, npes
call MPI_Comm_rank(MPI_COMM_WORLD, rank, status)
_VERIFY(status)
call MPI_Comm_size(MPI_COMM_WORLD, npes, status)
_VERIFY(status)
color = (rank*spec%n_writers) / npes
call MPI_Comm_split(MPI_COMM_WORLD, color, 0, gather_comm, status)
_VERIFY(status)
call MPI_Comm_rank(gather_comm, rank, status)
_VERIFY(status)
call MPI_Comm_split(MPI_COMM_WORLD, rank, 0, writer_comm, status)
_VERIFY(status)
if (rank /= 0) then
writer_comm = MPI_COMM_NULL
end if
kernel = make_GathervKernel(spec, gather_comm, _RC)
if (rank == 0) then
benchmark = make_BW_Benchmark(spec, writer_comm, _RC)
end if
call write_header(MPI_COMM_WORLD, _RC)
tot_time = 0
tot_time_gather = 0
tot_time_write = 0
associate (n => spec%n_tries)
do i = 1, n
ta = time(kernel, gather_comm, _RC)
if (writer_comm /= MPI_COMM_NULL) then
tb = time(benchmark, writer_comm, _RC)
end if
tot_time_gather = tot_time_gather + ta
tot_time_write = tot_time_write + tb
tot_time = tot_time + ta + tb
end do
avg_time = tot_time / n
avg_time_gather = tot_time_gather / n
avg_time_write = tot_time_write / n
end associate
call report(spec, avg_time, avg_time_gather, avg_time_write, MPI_COMM_WORLD, _RC)
_RETURN(_SUCCESS)
end subroutine run