src/xbt/parmap.c

   1 /* Copyright (c) 2004, 2005, 2007, 2009, 2010. The SimGrid Team.
   2  * All rights reserved.                                                     */
   3
   4 /* This program is free software; you can redistribute it and/or modify it
   5  * under the terms of the license (GNU LGPL) which comes with this package. */
   6 #include "gras_config.h"
   7 #include <unistd.h>
   8 #ifndef _XBT_WIN32
   9 #include <sys/syscall.h>
  10 #endif
  11
  12 #ifdef HAVE_FUTEX_H
  13         #include <linux/futex.h>
  14 #else
  15         #include "xbt/xbt_os_thread.h"
  16 #endif
  17 #include <errno.h>
  18 #include "parmap_private.h"
  19
  20 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(xbt_parmap, xbt, "parmap: parallel map");
  21 XBT_LOG_NEW_SUBCATEGORY(xbt_parmap_unit, xbt_parmap, "parmap unit testing");
  22
  23 static void *_xbt_parmap_worker_main(void *parmap);
  24 #ifdef HAVE_FUTEX_H
  25         static void futex_wait(int *uaddr, int val);
  26         static void futex_wake(int *uaddr, int val);
  27 #endif
  28 xbt_parmap_t xbt_parmap_new(unsigned int num_workers)
  29 {
  30   unsigned int i;
  31   xbt_os_thread_t worker = NULL;
  32
  33   XBT_DEBUG("Create new parmap (%u workers)", num_workers);
  34
  35   /* Initialize the thread pool data structure */
  36   xbt_parmap_t parmap = xbt_new0(s_xbt_parmap_t, 1);
  37 #ifdef HAVE_FUTEX_H
  38   parmap->sync_event = xbt_new0(s_xbt_event_t, 1);
  39 #endif
  40   parmap->num_workers = num_workers;
  41   parmap->status = PARMAP_WORK;
  42 #ifdef HAVE_FUTEX_H
  43   parmap->sync_event->threads_to_wait = num_workers;
  44 #endif
  45   /* Create the pool of worker threads */
  46   for(i=0; i < num_workers; i++){
  47     worker = xbt_os_thread_create(NULL, _xbt_parmap_worker_main, parmap, NULL);
  48     xbt_os_thread_detach(worker);
  49   }
  50 #ifdef HAVE_FUTEX_H
  51   xbt_event_init(parmap->sync_event);
  52 #endif
  53   return parmap;
  54 }
  55
  56 void xbt_parmap_destroy(xbt_parmap_t parmap)
  57 {
  58   XBT_DEBUG("Destroy parmap %p", parmap);
  59   parmap->status = PARMAP_DESTROY;
  60 #ifdef HAVE_FUTEX_H
  61   xbt_event_signal(parmap->sync_event);
  62   xbt_free(parmap->sync_event);
  63 #endif
  64   xbt_free(parmap);
  65 }
  66
  67  void xbt_parmap_apply(xbt_parmap_t parmap, void_f_pvoid_t fun, xbt_dynar_t data)
  68 {
  69   /* Assign resources to worker threads*/
  70   parmap->fun = fun;
  71   parmap->data = data;
  72   parmap->index = 0;
  73 #ifdef HAVE_FUTEX_H
  74   xbt_event_signal(parmap->sync_event);
  75 #endif
  76   XBT_DEBUG("Job done");
  77 }
  78
  79 static void *_xbt_parmap_worker_main(void *arg)
  80 {
  81   unsigned int worker_id;
  82   xbt_parmap_t parmap = (xbt_parmap_t)arg;
  83
  84   /* Fetch a worker id */
  85   worker_id = __sync_fetch_and_add(&parmap->workers_max_id, 1);
  86   xbt_os_thread_set_extra_data((void *)(unsigned long)worker_id);
  87
  88   XBT_DEBUG("New worker thread created (%u)", worker_id);
  89
  90   /* Worker's main loop */
  91   while(1){
  92 #ifdef HAVE_FUTEX_H
  93     xbt_event_wait(parmap->sync_event);
  94 #endif
  95     if(parmap->status == PARMAP_WORK){
  96       unsigned int i;
  97       unsigned int n = 0;
  98
  99       XBT_DEBUG("Worker %u got a job", worker_id);
 100
 101       while ((i = __sync_fetch_and_add(&parmap->index, 1))
 102              < xbt_dynar_length(parmap->data)) {
 103         parmap->fun(xbt_dynar_get_as(parmap->data, i, void*));
 104         n++;
 105       }
 106
 107       XBT_DEBUG("Worker %u processed %u tasks", worker_id, n);
 108
 109     /* We are destroying the parmap */
 110     }else{
 111 #ifdef HAVE_FUTEX_H
 112       xbt_event_end(parmap->sync_event);
 113 #endif
 114       XBT_DEBUG("Shutting down worker %u", worker_id);
 115       return NULL;
 116     }
 117   }
 118 }
 119
 120 #ifdef HAVE_FUTEX_H
 121 static void futex_wait(int *uaddr, int val)
 122 {
 123   XBT_VERB("Waiting on futex %p", uaddr);
 124   syscall(SYS_futex, uaddr, FUTEX_WAIT_PRIVATE, val, NULL, NULL, 0);
 125 }
 126
 127 static void futex_wake(int *uaddr, int val)
 128 {
 129   XBT_VERB("Waking futex %p", uaddr);
 130   syscall(SYS_futex, uaddr, FUTEX_WAKE_PRIVATE, val, NULL, NULL, 0);
 131 }
 132
 133 void xbt_event_init(xbt_event_t event)
 134 {
 135   int myflag = event->done;
 136   if(event->thread_counter < event->threads_to_wait)
 137     futex_wait(&event->done, myflag);
 138 }
 139
 140 void xbt_event_signal(xbt_event_t event)
 141 {
 142   int myflag = event->done;
 143   event->thread_counter = 0;
 144   event->work++;
 145   futex_wake(&event->work, event->threads_to_wait);
 146   futex_wait(&event->done, myflag);
 147 }
 148
 149 void xbt_event_wait(xbt_event_t event)
 150 {
 151   int myflag;
 152   unsigned int mycount;
 153
 154   myflag = event->work;
 155   mycount = __sync_add_and_fetch(&event->thread_counter, 1);
 156   if(mycount == event->threads_to_wait){
 157     event->done++;
 158     futex_wake(&event->done, 1);
 159   }
 160
 161   futex_wait(&event->work, myflag);
 162 }
 163
 164 void xbt_event_end(xbt_event_t event)
 165 {
 166   int myflag;
 167   unsigned int mycount;
 168
 169   myflag = event->work;
 170   mycount = __sync_add_and_fetch(&event->thread_counter, 1);
 171   if(mycount == event->threads_to_wait){
 172     event->done++;
 173     futex_wake(&event->done, 1);
 174   }
 175 }
 176 #endif
 177
 178 #ifdef SIMGRID_TEST
 179 #include "xbt.h"
 180 #include "xbt/ex.h"
 181
 182 XBT_TEST_SUITE("parmap", "Parallel Map");
 183 XBT_LOG_EXTERNAL_DEFAULT_CATEGORY(xbt_parmap_unit);
 184
 185
 186
 187 xbt_parmap_t parmap;
 188
 189 void fun(void *arg);
 190
 191 void fun(void *arg)
 192 {
 193   //XBT_INFO("I'm job %lu", (unsigned long)arg);
 194 }
 195
 196 XBT_TEST_UNIT("basic", test_parmap_basic, "Basic usage")
 197 {
 198   xbt_test_add("Create the parmap");
 199
 200   unsigned long i,j;
 201   xbt_dynar_t data = xbt_dynar_new(sizeof(void *), NULL);
 202
 203   /* Create the parallel map */
 204   parmap = xbt_parmap_new(10);
 205
 206   for(j=0; j < 100; j++){
 207     xbt_dynar_push_as(data, void *, (void *)j);
 208   }
 209
 210   for(i=0; i < 5; i++)
 211     xbt_parmap_apply(parmap, fun, data);
 212
 213   /* Destroy the parmap */
 214   xbt_parmap_destroy(parmap);
 215   xbt_dynar_free(&data);
 216 }
 217
 218 #endif /* SIMGRID_TEST */