Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
Do connect all log channel manually to parent using XBT_LOG_CONNECT() too, so that...
[simgrid.git] / src / smpi / smpi_base.c
index 7f0fd9d..c119790 100644 (file)
 #include "private.h"
 
-SMPI_MPI_Global_t smpi_mpi_global = NULL;
+XBT_LOG_NEW_DEFAULT_SUBCATEGORY(smpi_base, smpi, "Logging specific to SMPI (base)");
 
-void smpi_mpi_land_func(void *x, void *y, void *z)
-{
-       *(int *)z = *(int *)x && *(int *)y;
-}
+smpi_mpi_global_t smpi_mpi_global = NULL;
 
-void smpi_mpi_sum_func(void *x, void *y, void *z)
-{
-       *(int *)z = *(int *)x + *(int *)y;
-}
+void smpi_mpi_land_func(void *a, void *b, int *length, MPI_Datatype *datatype);
 
-int inline smpi_mpi_comm_size(smpi_mpi_communicator_t comm)
+void smpi_mpi_land_func(void *a, void *b, int *length, MPI_Datatype *datatype)
 {
-       return comm->size;
+       int i;
+       if (*datatype == smpi_mpi_global->mpi_int) {
+               int *x = a, *y = b;
+               for (i = 0; i < *length; i++) {
+                       y[i] = x[i] && y[i];
+               }
+       }
 }
 
-// FIXME: smarter algorithm?
-int smpi_mpi_comm_rank(smpi_mpi_communicator_t comm, smx_host_t host)
+void smpi_mpi_sum_func(void *a, void *b, int *length, MPI_Datatype *datatype);
+
+void smpi_mpi_sum_func(void *a, void *b, int *length, MPI_Datatype *datatype)
 {
        int i;
-
-       for(i = comm->size - 1; i > 0 && host != comm->hosts[i]; i--);
-
-       return i;
+       if (*datatype == smpi_mpi_global->mpi_int) {
+               int *x = a, *y = b;
+               for (i = 0; i < *length; i++) {
+                       y[i] = x[i] + y[i];
+               }
+       }
 }
 
-int inline smpi_mpi_comm_rank_self(smpi_mpi_communicator_t comm)
+int smpi_mpi_comm_rank(smpi_mpi_communicator_t comm)
 {
-       return smpi_mpi_comm_rank(comm, SIMIX_host_self());
+       return comm->index_to_rank_map[smpi_host_index()];
 }
 
 void smpi_mpi_init()
 {
-       smx_process_t process;
        smx_host_t host;
        smx_host_t *hosts;
-       int size;
-
+       int host_count;
+       int i;
+       smpi_host_data_t hdata;
+
+   /* Connect our log channels: that must be done manually under windows */
+   /* (should be done only once, not for each process) */
+   XBT_LOG_CONNECT(smpi_base, smpi);
+   XBT_LOG_CONNECT(smpi_bench, smpi);
+   XBT_LOG_CONNECT(smpi_kernel, smpi);
+   XBT_LOG_CONNECT(smpi_mpi, smpi);
+   XBT_LOG_CONNECT(smpi_receiver, smpi);
+   XBT_LOG_CONNECT(smpi_sender, smpi);
+   XBT_LOG_CONNECT(smpi_util, smpi);
+   
        SIMIX_mutex_lock(smpi_global->running_hosts_count_mutex);
        smpi_global->running_hosts_count++;
        SIMIX_mutex_unlock(smpi_global->running_hosts_count_mutex);
 
        // initialize some local variables
-       process = SIMIX_process_self();
-       host    = SIMIX_host_self();
-       hosts   = SIMIX_host_get_table();
-       size    = SIMIX_host_get_number();
+       host       = SIMIX_host_self();
+       hosts      = SIMIX_host_get_table();
+       host_count = SIMIX_host_get_number();
+
+       hdata      = xbt_new(s_smpi_host_data_t, 1);
+
+       for (i = 0; i < host_count && host != hosts[i]; i ++);
+
+       hdata->index = i;
+       hdata->mutex = SIMIX_mutex_init();
+       hdata->cond  = SIMIX_cond_init();
+
+       SIMIX_host_set_data(host, hdata);
 
        // node 0 sets the globals
-       if (host == hosts[0]) {
+       if (0 == i) {
+
+               smpi_global->hosts                              = hosts;
+               smpi_global->host_count                         = host_count;
 
-               smpi_mpi_global                                = xbt_new(s_SMPI_MPI_Global_t, 1);
+               smpi_mpi_global                                 = xbt_new(s_smpi_mpi_global_t, 1);
 
                // global communicator
-               smpi_mpi_global->mpi_comm_world                = xbt_new(s_smpi_mpi_communicator_t, 1);
-               smpi_mpi_global->mpi_comm_world->size          = size;
-               smpi_mpi_global->mpi_comm_world->hosts         = hosts;
-               smpi_mpi_global->mpi_comm_world->processes     = xbt_new(smx_process_t, size);
-               smpi_mpi_global->mpi_comm_world->processes[0]  = process;
-               smpi_mpi_global->mpi_comm_world->barrier_count = 0;
-               smpi_mpi_global->mpi_comm_world->barrier_mutex = SIMIX_mutex_init();
-               smpi_mpi_global->mpi_comm_world->barrier_cond  = SIMIX_cond_init();
+               smpi_mpi_global->mpi_comm_world                 = xbt_new(s_smpi_mpi_communicator_t, 1);
+               smpi_mpi_global->mpi_comm_world->size           = host_count;
+               smpi_mpi_global->mpi_comm_world->barrier_count  = 0;
+               smpi_mpi_global->mpi_comm_world->barrier_mutex  = SIMIX_mutex_init();
+               smpi_mpi_global->mpi_comm_world->barrier_cond   = SIMIX_cond_init();
+               smpi_mpi_global->mpi_comm_world->rank_to_index_map = xbt_new(int, host_count);
+               smpi_mpi_global->mpi_comm_world->index_to_rank_map = xbt_new(int, host_count);
+               for (i = 0; i < host_count; i++) {
+                       smpi_mpi_global->mpi_comm_world->rank_to_index_map[i] = i;
+                       smpi_mpi_global->mpi_comm_world->index_to_rank_map[i] = i;
+               }
 
                // mpi datatypes
-               smpi_mpi_global->mpi_byte                      = xbt_new(s_smpi_mpi_datatype_t, 1);
-               smpi_mpi_global->mpi_byte->size                = (size_t)1;
-               smpi_mpi_global->mpi_int                       = xbt_new(s_smpi_mpi_datatype_t, 1);
-               smpi_mpi_global->mpi_int->size                 = sizeof(int);
-               smpi_mpi_global->mpi_double                    = xbt_new(s_smpi_mpi_datatype_t, 1);
-               smpi_mpi_global->mpi_double->size              = sizeof(double);
+               smpi_mpi_global->mpi_byte                       = xbt_new(s_smpi_mpi_datatype_t, 1);
+               smpi_mpi_global->mpi_byte->size                 = (size_t)1;
+               smpi_mpi_global->mpi_int                        = xbt_new(s_smpi_mpi_datatype_t, 1);
+               smpi_mpi_global->mpi_int->size                  = sizeof(int);
+               smpi_mpi_global->mpi_double                     = xbt_new(s_smpi_mpi_datatype_t, 1);
+               smpi_mpi_global->mpi_double->size               = sizeof(double);
 
                // mpi operations
-               smpi_mpi_global->mpi_land                      = xbt_new(s_smpi_mpi_op_t, 1);
-               smpi_mpi_global->mpi_land->func                = smpi_mpi_land_func;
-               smpi_mpi_global->mpi_sum                       = xbt_new(s_smpi_mpi_op_t, 1);
-               smpi_mpi_global->mpi_sum->func                 = smpi_mpi_sum_func;
+               smpi_mpi_global->mpi_land                       = xbt_new(s_smpi_mpi_op_t, 1);
+               smpi_mpi_global->mpi_land->func                 = smpi_mpi_land_func;
+               smpi_mpi_global->mpi_sum                        = xbt_new(s_smpi_mpi_op_t, 1);
+               smpi_mpi_global->mpi_sum->func                  = smpi_mpi_sum_func;
 
                // signal all nodes to perform initialization
                SIMIX_mutex_lock(smpi_global->start_stop_mutex);
@@ -88,22 +117,22 @@ void smpi_mpi_init()
 
                // make sure root is done before own initialization
                SIMIX_mutex_lock(smpi_global->start_stop_mutex);
-               if (!smpi_global->root_ready) {
+               while (!smpi_global->root_ready) {
                        SIMIX_cond_wait(smpi_global->start_stop_cond, smpi_global->start_stop_mutex);
                }
                SIMIX_mutex_unlock(smpi_global->start_stop_mutex);
 
-               smpi_mpi_global->mpi_comm_world->processes[smpi_mpi_comm_rank_self(smpi_mpi_global->mpi_comm_world)] = process;
        }
 
        // wait for all nodes to signal initializatin complete
        SIMIX_mutex_lock(smpi_global->start_stop_mutex);
        smpi_global->ready_process_count++;
-       if (smpi_global->ready_process_count < 3 * size) {
-               SIMIX_cond_wait(smpi_global->start_stop_cond, smpi_global->start_stop_mutex);
-       } else {
+       if (smpi_global->ready_process_count >= 3 * host_count) {
                SIMIX_cond_broadcast(smpi_global->start_stop_cond);
        }
+       while (smpi_global->ready_process_count < 3 * host_count) {
+               SIMIX_cond_wait(smpi_global->start_stop_cond, smpi_global->start_stop_mutex);
+       }
        SIMIX_mutex_unlock(smpi_global->start_stop_mutex);
 
        return;
@@ -121,10 +150,13 @@ void smpi_mpi_finalize()
        smpi_global->ready_process_count--;
        SIMIX_mutex_unlock(smpi_global->start_stop_mutex);
 
+       SIMIX_mutex_destroy(smpi_host_mutex());
+       SIMIX_cond_destroy(smpi_host_cond());
+
        if (0 >= i) {
 
                // wake up senders/receivers
-               for (i = 0; i < smpi_mpi_global->mpi_comm_world->size; i++) {
+               for (i = 0; i < smpi_global->host_count; i++) {
                        if (SIMIX_process_is_suspended(smpi_global->sender_processes[i])) {
                                SIMIX_process_resume(smpi_global->sender_processes[i]);
                        }
@@ -135,7 +167,6 @@ void smpi_mpi_finalize()
 
                SIMIX_mutex_destroy(smpi_mpi_global->mpi_comm_world->barrier_mutex);
                SIMIX_cond_destroy(smpi_mpi_global->mpi_comm_world->barrier_cond);
-               xbt_free(smpi_mpi_global->mpi_comm_world->processes);
                xbt_free(smpi_mpi_global->mpi_comm_world);
 
                xbt_free(smpi_mpi_global->mpi_byte);
@@ -151,103 +182,82 @@ void smpi_mpi_finalize()
 
 }
 
-void smpi_barrier(smpi_mpi_communicator_t comm)
+int smpi_mpi_barrier(smpi_mpi_communicator_t comm)
 {
 
        SIMIX_mutex_lock(comm->barrier_mutex);
-       if(++comm->barrier_count < comm->size) {
-               SIMIX_cond_wait(comm->barrier_cond, comm->barrier_mutex);
-       } else {
+       ++comm->barrier_count;
+       if (comm->barrier_count > comm->size) { // only happens on second barrier...
                comm->barrier_count = 0;
+       } else if (comm->barrier_count == comm->size) {
                SIMIX_cond_broadcast(comm->barrier_cond);
        }
+       while (comm->barrier_count < comm->size) {
+               SIMIX_cond_wait(comm->barrier_cond, comm->barrier_mutex);
+       }
        SIMIX_mutex_unlock(comm->barrier_mutex);
 
-       return;
+       return MPI_SUCCESS;
 }
 
-int smpi_create_request(void *buf, int count, smpi_mpi_datatype_t datatype,
-       int src, int dst, int tag, smpi_mpi_communicator_t comm, smpi_mpi_request_t *request)
+int smpi_mpi_isend(smpi_mpi_request_t request)
 {
        int retval = MPI_SUCCESS;
+       int index  = smpi_host_index();
 
-       *request = NULL;
-
-       if (0 > count) {
-               retval = MPI_ERR_COUNT;
-       } else if (NULL == buf) {
+       if (NULL == request) {
                retval = MPI_ERR_INTERN;
-       } else if (NULL == datatype) {
-               retval = MPI_ERR_TYPE;
-       } else if (NULL == comm) {
-               retval = MPI_ERR_COMM;
-       } else if (MPI_ANY_SOURCE != src && (0 > src || comm->size <= src)) {
-               retval = MPI_ERR_RANK;
-       } else if (0 > dst || comm->size <= dst) {
-               retval = MPI_ERR_RANK;
-       } else if (0 > tag) {
-               retval = MPI_ERR_TAG;
        } else {
-               *request               = xbt_mallocator_get(smpi_global->request_mallocator);
-               (*request)->comm       = comm;
-               (*request)->src        = src;
-               (*request)->dst        = dst;
-               (*request)->tag        = tag;
-               (*request)->buf        = buf;
-               (*request)->count      = count;
-               (*request)->datatype   = datatype;
-       }
-       return retval;
-}
-
-int smpi_isend(smpi_mpi_request_t request)
-{
-       int retval = MPI_SUCCESS;
-       int rank   = smpi_mpi_comm_rank_self(smpi_mpi_global->mpi_comm_world);
+               SIMIX_mutex_lock(smpi_global->pending_send_request_queues_mutexes[index]);
+               xbt_fifo_push(smpi_global->pending_send_request_queues[index], request);
+               SIMIX_mutex_unlock(smpi_global->pending_send_request_queues_mutexes[index]);
 
-       if (NULL != request) {
-               SIMIX_mutex_lock(smpi_global->pending_send_request_queues_mutexes[rank]);
-               xbt_fifo_push(smpi_global->pending_send_request_queues[rank], request);
-               SIMIX_mutex_unlock(smpi_global->pending_send_request_queues_mutexes[rank]);
-       }
-
-       if (SIMIX_process_is_suspended(smpi_global->sender_processes[rank])) {
-               SIMIX_process_resume(smpi_global->sender_processes[rank]);
+               if (SIMIX_process_is_suspended(smpi_global->sender_processes[index])) {
+                       SIMIX_process_resume(smpi_global->sender_processes[index]);
+               }
        }
 
        return retval;
 }
 
-int smpi_irecv(smpi_mpi_request_t request)
+int smpi_mpi_irecv(smpi_mpi_request_t request)
 {
        int retval = MPI_SUCCESS;
-       int rank = smpi_mpi_comm_rank_self(smpi_mpi_global->mpi_comm_world);
+       int index = smpi_host_index();
 
-       if (NULL != request) {
-               SIMIX_mutex_lock(smpi_global->pending_recv_request_queues_mutexes[rank]);
-               xbt_fifo_push(smpi_global->pending_recv_request_queues[rank], request);
-               SIMIX_mutex_unlock(smpi_global->pending_recv_request_queues_mutexes[rank]);
-       }
+       if (NULL == request) {
+               retval = MPI_ERR_INTERN;
+       } else {
+               SIMIX_mutex_lock(smpi_global->pending_recv_request_queues_mutexes[index]);
+               xbt_fifo_push(smpi_global->pending_recv_request_queues[index], request);
+               SIMIX_mutex_unlock(smpi_global->pending_recv_request_queues_mutexes[index]);
 
-       if (SIMIX_process_is_suspended(smpi_global->receiver_processes[rank])) {
-               SIMIX_process_resume(smpi_global->receiver_processes[rank]);
+               if (SIMIX_process_is_suspended(smpi_global->receiver_processes[index])) {
+                       SIMIX_process_resume(smpi_global->receiver_processes[index]);
+               }
        }
 
        return retval;
 }
 
-int smpi_wait(smpi_mpi_request_t request, smpi_mpi_status_t *status)
+int smpi_mpi_wait(smpi_mpi_request_t request, smpi_mpi_status_t *status)
 {
-       if (NULL != request) {
+       int retval = MPI_SUCCESS;
+
+       if (NULL == request) {
+               retval = MPI_ERR_INTERN;
+       } else {
                SIMIX_mutex_lock(request->mutex);
-               if (!request->completed) {
+               while (!request->completed) {
                        SIMIX_cond_wait(request->cond, request->mutex);
                }
                if (NULL != status) {
                        status->MPI_SOURCE = request->src;
+                       status->MPI_TAG    = request->tag;
+                       status->MPI_ERROR  = MPI_SUCCESS;
                }
                SIMIX_mutex_unlock(request->mutex);
        }
 
-       return MPI_SUCCESS;
+       return retval;
 }