/* -*- 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"

int main(int argc, char *argv[])
{
    char str[10];
    int err = 0, errcodes[256], rank, rank1, nprocs;
    MPI_Comm intercomm1, intercomm2, intracomm;

    MPI_Init(&argc, &argv);

    MPI_Comm_size(MPI_COMM_WORLD, &nprocs);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank1);

    if (nprocs != 2) {
        printf("Run this program with 2 processes\n");
        MPI_Abort(MPI_COMM_WORLD, 1);
    }

    if (rank1 == 0) {
        printf("Parents spawning 2 children...\n");
        fflush(stdout);
    }

    err = MPI_Comm_spawn("spawn_merge_child1", MPI_ARGV_NULL, 2,
                         MPI_INFO_NULL, 1, MPI_COMM_WORLD, &intercomm1, errcodes);
    if (err)
        printf("Error in MPI_Comm_spawn\n");

    if (rank1 == 0) {
        printf("Parents and children merging to form new intracommunicator...\n");
        fflush(stdout);
    }

    MPI_Intercomm_merge(intercomm1, 0, &intracomm);

    if (rank1 == 0) {
        printf("Merged parents spawning 2 more children and communicating with them...\n");
        fflush(stdout);
    }

    err = MPI_Comm_spawn("spawn_merge_child2", MPI_ARGV_NULL, 2,
                         MPI_INFO_NULL, 2, intracomm, &intercomm2, errcodes);
    if (err)
        printf("Error in MPI_Comm_spawn\n");

    MPI_Comm_rank(intercomm2, &rank);

    if (rank == 2) {
        err = MPI_Recv(str, 3, MPI_CHAR, 1, 0, intercomm2, MPI_STATUS_IGNORE);
        printf("Parent received from child: %s\n", str);
        fflush(stdout);

        err = MPI_Send("bye", 4, MPI_CHAR, 1, 0, intercomm2);
    }

    MPI_Finalize();

    return 0;
}