/* -*- Mode: C; c-basic-offset:4 ; indent-tabs-mode:nil ; -*- */
/*
 *  (C) 2001 by Argonne National Laboratory.
 *      See COPYRIGHT in top-level directory.
 */

#include <stdio.h>
#include "mpi.h"
#include "mpitest.h"

#define TABLE_SIZE 2

int main(int argc, char **argv)
{
    int rank, size;
    double a[TABLE_SIZE];
    struct {
        double a;
        int b;
    } in[TABLE_SIZE], out[TABLE_SIZE];
    int i;
    int errors = 0, toterrors;

    /* Initialize the environment and some variables */
    MTest_Init(&argc, &argv);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
    MPI_Comm_size(MPI_COMM_WORLD, &size);

    /* Initialize the maxloc data */
    for (i = 0; i < TABLE_SIZE; i++)
        a[i] = 0;
    for (i = rank; i < TABLE_SIZE; i++)
        a[i] = (double) rank + 1.0;

    /* Copy data to the "in" buffer */
    for (i = 0; i < TABLE_SIZE; i++) {
        in[i].a = a[i];
        in[i].b = rank;
    }

    /* Reduce it! */
    MPI_Reduce(in, out, TABLE_SIZE, MPI_DOUBLE_INT, MPI_MAXLOC, 0, MPI_COMM_WORLD);
    MPI_Bcast(out, TABLE_SIZE, MPI_DOUBLE_INT, 0, MPI_COMM_WORLD);

    /* Check to see that we got the right answers */
    for (i = 0; i < TABLE_SIZE; i++)
        if (i % size == rank)
            if (out[i].b != rank) {
                printf("MAX (ranks[%d] = %d != %d\n", i, out[i].b, rank);
                errors++;
            }

    /* Initialize the minloc data */
    for (i = 0; i < TABLE_SIZE; i++)
        a[i] = 0;
    for (i = rank; i < TABLE_SIZE; i++)
        a[i] = -(double) rank - 1.0;

    /* Copy data to the "in" buffer */
    for (i = 0; i < TABLE_SIZE; i++) {
        in[i].a = a[i];
        in[i].b = rank;
    }

    /* Reduce it! */
    MPI_Allreduce(in, out, TABLE_SIZE, MPI_DOUBLE_INT, MPI_MINLOC, MPI_COMM_WORLD);

    /* Check to see that we got the right answers */
    for (i = 0; i < TABLE_SIZE; i++)
        if (i % size == rank)
            if (out[i].b != rank) {
                printf("MIN (ranks[%d] = %d != %d\n", i, out[i].b, rank);
                errors++;
            }

    /* Finish up! */
    MTest_Finalize(errors);
    MPI_Finalize();
    return MTestReturnValue(errors);
}
