/* -*- Mode: C; c-basic-offset:4 ; indent-tabs-mode:nil ; -*- */
/*
 *  (C) 2001 by Argonne National Laboratory.
 *      See COPYRIGHT in top-level directory.
 */
#include "mpi.h"
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include "mpitest.h"

#define ROOT      0
#define NUM_REPS  5
#define NUM_SIZES 4

int main(int argc, char **argv)
{
    int *buf;
    int i, rank, reps, n;
    int bVerify = 1;
    int sizes[NUM_SIZES] = { 100, 64 * 1024, 128 * 1024, 1024 * 1024 };
    int num_errors = 0;

    MTest_Init(&argc, &argv);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);

    if (argc > 1) {
        if (strcmp(argv[1], "-novalidate") == 0 || strcmp(argv[1], "-noverify") == 0)
            bVerify = 0;
    }

    buf = (int *) malloc(sizes[NUM_SIZES - 1] * sizeof(int));
    memset(buf, 0, sizes[NUM_SIZES - 1] * sizeof(int));

    for (n = 0; n < NUM_SIZES; n++) {
#ifdef DEBUG
        if (rank == ROOT) {
            printf("bcasting %d MPI_INTs %d times\n", sizes[n], NUM_REPS);
            fflush(stdout);
        }
#endif
        for (reps = 0; reps < NUM_REPS; reps++) {
            if (bVerify) {
                if (rank == ROOT) {
                    for (i = 0; i < sizes[n]; i++) {
                        buf[i] = 1000000 * (n * NUM_REPS + reps) + i;
                    }
                }
                else {
                    for (i = 0; i < sizes[n]; i++) {
                        buf[i] = -1 - (n * NUM_REPS + reps);
                    }
                }
            }

#	    ifdef DEBUG
            {
                printf("rank=%d, n=%d, reps=%d\n", rank, n, reps);
            }
#           endif

            MPI_Bcast(buf, sizes[n], MPI_INT, ROOT, MPI_COMM_WORLD);

            if (bVerify) {
                num_errors = 0;
                for (i = 0; i < sizes[n]; i++) {
                    if (buf[i] != 1000000 * (n * NUM_REPS + reps) + i) {
                        num_errors++;
                        if (num_errors < 10) {
                            printf("Error: Rank=%d, n=%d, reps=%d, i=%d, buf[i]=%d expected=%d\n",
                                   rank, n, reps, i, buf[i], 1000000 * (n * NUM_REPS + reps) + i);
                            fflush(stdout);
                        }
                    }
                }
                if (num_errors >= 10) {
                    printf("Error: Rank=%d, num_errors = %d\n", rank, num_errors);
                    fflush(stdout);
                }
            }
        }
    }

    free(buf);

    MTest_Finalize(num_errors);
    MPI_Finalize();
    return 0;
}