1 /* Copyright (c) 2004, 2005, 2007, 2009, 2010. The SimGrid Team.
2 * All rights reserved. */
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"
10 #include <sys/syscall.h>
14 #include <linux/futex.h>
17 #include "xbt/parmap.h"
19 #include "xbt/function_types.h"
20 #include "xbt/dynar.h"
21 #include "xbt/xbt_os_thread.h"
22 #include "xbt/sysdep.h"
24 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(xbt_parmap, xbt, "parmap: parallel map");
25 XBT_LOG_NEW_SUBCATEGORY(xbt_parmap_unit, xbt_parmap, "parmap unit testing");
30 } e_xbt_parmap_flag_t;
32 static void xbt_parmap_set_mode(xbt_parmap_t parmap, e_xbt_parmap_mode_t mode);
33 static void *xbt_parmap_worker_main(void *parmap);
35 static void xbt_parmap_posix_start(xbt_parmap_t parmap);
36 static void xbt_parmap_posix_end(xbt_parmap_t parmap);
37 static void xbt_parmap_posix_signal(xbt_parmap_t parmap);
38 static void xbt_parmap_posix_wait(xbt_parmap_t parmap);
41 static void xbt_parmap_futex_start(xbt_parmap_t parmap);
42 static void xbt_parmap_futex_end(xbt_parmap_t parmap);
43 static void xbt_parmap_futex_signal(xbt_parmap_t parmap);
44 static void xbt_parmap_futex_wait(xbt_parmap_t parmap);
45 static void futex_wait(int *uaddr, int val);
46 static void futex_wake(int *uaddr, int val);
49 static void xbt_parmap_busy_start(xbt_parmap_t parmap);
50 static void xbt_parmap_busy_end(xbt_parmap_t parmap);
51 static void xbt_parmap_busy_signal(xbt_parmap_t parmap);
52 static void xbt_parmap_busy_wait(xbt_parmap_t parmap);
56 * \brief Parallel map structure
58 typedef struct s_xbt_parmap {
59 e_xbt_parmap_flag_t status; /**< is the parmap active or being destroyed? */
60 int work; /**< index of the current round (1 is the first) */
61 int done; /**< number of rounds already done (futexes only) */
62 unsigned int thread_counter; /**< number of threads currently working */
63 unsigned int num_workers; /**< total number of worker threads including the controller */
64 void_f_pvoid_t fun; /**< function to run in parallel on each element of data */
65 xbt_dynar_t data; /**< parameters to pass to fun in parallel */
66 unsigned int index; /**< index of the next element of data to pick */
68 /* fields that depend on the synchronization mode */
69 e_xbt_parmap_mode_t mode; /**< synchronization mode */
70 void (*start_f)(xbt_parmap_t); /**< initializes the worker threads */
71 void (*end_f)(xbt_parmap_t); /**< finalizes the worker threads */
72 void (*signal_f)(xbt_parmap_t); /**< wakes the workers threads to process tasks */
73 void (*wait_f)(xbt_parmap_t); /**< waits for more work */
77 * \brief Creates a parallel map object
78 * \param num_workers number of worker threads to create
79 * \param mode how to synchronize the worker threads
80 * \return the parmap created
82 xbt_parmap_t xbt_parmap_new(unsigned int num_workers, e_xbt_parmap_mode_t mode)
85 xbt_os_thread_t worker = NULL;
87 XBT_DEBUG("Create new parmap (%u workers)", num_workers);
89 /* Initialize the thread pool data structure */
90 xbt_parmap_t parmap = xbt_new0(s_xbt_parmap_t, 1);
92 parmap->num_workers = num_workers;
93 parmap->status = XBT_PARMAP_WORK;
94 xbt_parmap_set_mode(parmap, mode);
96 /* Create the pool of worker threads */
97 for (i = 0; i < num_workers - 1; i++) {
98 worker = xbt_os_thread_create(NULL, xbt_parmap_worker_main, parmap, NULL);
99 xbt_os_thread_detach(worker);
101 parmap->start_f(parmap);
106 * \brief Destroys a parmap
107 * \param parmap the parmap to destroy
109 void xbt_parmap_destroy(xbt_parmap_t parmap)
115 parmap->status = XBT_PARMAP_DESTROY;
116 parmap->signal_f(parmap);
121 * \brief Sets the synchronization mode of a parmap.
122 * \param parmap a parallel map object
123 * \param mode the synchronization mode
125 static void xbt_parmap_set_mode(xbt_parmap_t parmap, e_xbt_parmap_mode_t mode)
127 if (mode == XBT_PARMAP_DEFAULT) {
129 mode = XBT_PARMAP_FUTEX;
131 //For now use busy wait because posix is unimplemented
132 mode = XBT_PARMAP_BUSY_WAIT;
139 case XBT_PARMAP_POSIX:
140 parmap->start_f = xbt_parmap_posix_start;
141 parmap->end_f = xbt_parmap_posix_end;
142 parmap->signal_f = xbt_parmap_posix_signal;
143 parmap->wait_f = xbt_parmap_posix_wait;
147 case XBT_PARMAP_FUTEX:
149 parmap->start_f = xbt_parmap_futex_start;
150 parmap->end_f = xbt_parmap_futex_end;
151 parmap->signal_f = xbt_parmap_futex_signal;
152 parmap->wait_f = xbt_parmap_futex_wait;
155 xbt_die("Futex is not available on this OS (maybe you are on a Mac).");
158 case XBT_PARMAP_BUSY_WAIT:
159 parmap->start_f = xbt_parmap_busy_start;
160 parmap->end_f = xbt_parmap_busy_end;
161 parmap->signal_f = xbt_parmap_busy_signal;
162 parmap->wait_f = xbt_parmap_busy_wait;
165 case XBT_PARMAP_DEFAULT:
172 * \brief Applies a list of tasks in parallel.
173 * \param parmap a parallel map object
174 * \param fun the function to call in parallel
175 * \param data each element of this dynar will be passed as an argument to fun
177 void xbt_parmap_apply(xbt_parmap_t parmap, void_f_pvoid_t fun, xbt_dynar_t data)
179 /* Assign resources to worker threads */
183 parmap->signal_f(parmap);
184 XBT_DEBUG("Job done");
188 * \brief Returns a next task to process.
190 * Worker threads call this function to get more work.
192 * \return the next task to process, or NULL if there is no more work
194 void* xbt_parmap_next(xbt_parmap_t parmap)
196 unsigned int index = __sync_fetch_and_add(&parmap->index, 1);
197 if (index < xbt_dynar_length(parmap->data)) {
198 return xbt_dynar_get_as(parmap->data, index, void*);
204 * \brief Main function of a worker thread.
205 * \param arg the parmap
207 static void *xbt_parmap_worker_main(void *arg)
209 xbt_parmap_t parmap = (xbt_parmap_t) arg;
211 XBT_DEBUG("New worker thread created");
213 /* Worker's main loop */
215 parmap->wait_f(parmap);
216 if (parmap->status == XBT_PARMAP_WORK) {
218 XBT_DEBUG("Worker got a job");
220 void* work = xbt_parmap_next(parmap);
221 while (work != NULL) {
223 work = xbt_parmap_next(parmap);
226 XBT_DEBUG("Worker has finished");
228 /* We are destroying the parmap */
230 parmap->end_f(parmap);
231 XBT_DEBUG("Shutting down worker");
238 static void futex_wait(int *uaddr, int val)
240 XBT_VERB("Waiting on futex %p", uaddr);
241 syscall(SYS_futex, uaddr, FUTEX_WAIT_PRIVATE, val, NULL, NULL, 0);
244 static void futex_wake(int *uaddr, int val)
246 XBT_VERB("Waking futex %p", uaddr);
247 syscall(SYS_futex, uaddr, FUTEX_WAKE_PRIVATE, val, NULL, NULL, 0);
251 static void xbt_parmap_posix_start(xbt_parmap_t parmap)
256 static void xbt_parmap_posix_end(xbt_parmap_t parmap)
261 static void xbt_parmap_posix_signal(xbt_parmap_t parmap)
266 static void xbt_parmap_posix_wait(xbt_parmap_t parmap)
273 * \brief Starts the parmap: waits for all workers to be ready and returns.
275 * This function is called by the controller thread.
277 * \param parmap a parmap
279 static void xbt_parmap_futex_start(xbt_parmap_t parmap)
281 int myflag = parmap->done;
282 __sync_fetch_and_add(&parmap->thread_counter, 1);
283 if (parmap->thread_counter < parmap->num_workers) {
284 /* wait for all workers to be ready */
285 futex_wait(&parmap->done, myflag);
290 * \brief Ends the parmap: wakes the controller thread when all workers terminate.
292 * This function is called by all worker threads when they end (not including
295 * \param parmap a parmap
297 static void xbt_parmap_futex_end(xbt_parmap_t parmap)
299 unsigned int mycount;
301 mycount = __sync_add_and_fetch(&parmap->thread_counter, 1);
302 if (mycount == parmap->num_workers) {
303 /* all workers have finished, wake the controller */
305 futex_wake(&parmap->done, 1);
310 * \brief Wakes all workers and waits for them to finish the tasks.
312 * This function is called by the controller thread.
314 * \param parmap a parmap
316 static void xbt_parmap_futex_signal(xbt_parmap_t parmap)
318 int myflag = parmap->done;
319 parmap->thread_counter = 0;
322 /* wake all workers */
323 futex_wake(&parmap->work, parmap->num_workers);
325 if (parmap->status == XBT_PARMAP_WORK) {
326 /* also work myself */
327 void* work = xbt_parmap_next(parmap);
328 while (work != NULL) {
330 work = xbt_parmap_next(parmap);
334 unsigned int mycount = __sync_add_and_fetch(&parmap->thread_counter, 1);
335 if (mycount < parmap->num_workers) {
336 /* some workers have not finished yet */
337 futex_wait(&parmap->done, myflag);
342 * \brief Waits for some work to process.
344 * This function is called by each worker thread (not including the controller)
345 * when it has no more work to do.
347 * \param parmap a parmap
349 static void xbt_parmap_futex_wait(xbt_parmap_t parmap)
352 unsigned int mycount;
354 myflag = parmap->work;
355 mycount = __sync_add_and_fetch(&parmap->thread_counter, 1);
356 if (mycount == parmap->num_workers) {
357 /* all workers have finished, wake the controller */
359 futex_wake(&parmap->done, 1);
362 /* wait for more work */
363 futex_wait(&parmap->work, myflag);
368 * \brief Starts the parmap: waits for all workers to be ready and returns.
370 * This function is called by the controller thread.
372 * \param parmap a parmap
374 static void xbt_parmap_busy_start(xbt_parmap_t parmap)
376 __sync_fetch_and_add(&parmap->thread_counter, 1);
377 while (parmap->thread_counter < parmap->num_workers) {
378 xbt_os_thread_yield();
383 * \brief Ends the parmap: wakes the controller thread when all workers terminate.
385 * This function is called by all worker threads when they end.
387 * \param parmap a parmap
389 static void xbt_parmap_busy_end(xbt_parmap_t parmap)
391 __sync_add_and_fetch(&parmap->thread_counter, 1);
395 * \brief Wakes all workers and waits for them to finish the tasks.
397 * This function is called by the controller thread.
399 * \param parmap a parmap
401 static void xbt_parmap_busy_signal(xbt_parmap_t parmap)
403 parmap->thread_counter = 0;
406 if (parmap->status == XBT_PARMAP_WORK) {
407 /* also work myself */
408 void* work = xbt_parmap_next(parmap);
409 while (work != NULL) {
411 work = xbt_parmap_next(parmap);
415 /* I have finished, wait for the others */
416 __sync_add_and_fetch(&parmap->thread_counter, 1);
417 while (parmap->thread_counter < parmap->num_workers) {
418 xbt_os_thread_yield();
423 * \brief Waits for some work to process.
425 * This function is called by each worker thread (not including the controller)
426 * when it has no more work to do.
428 * \param parmap a parmap
430 static void xbt_parmap_busy_wait(xbt_parmap_t parmap)
432 int work = parmap->work;
433 __sync_add_and_fetch(&parmap->thread_counter, 1);
435 /* wait for more work */
436 while (parmap->work == work) {
437 xbt_os_thread_yield();
445 XBT_TEST_SUITE("parmap", "Parallel Map");
446 XBT_LOG_EXTERNAL_DEFAULT_CATEGORY(xbt_parmap_unit);
454 //XBT_INFO("I'm job %lu", (unsigned long)arg);
457 XBT_TEST_UNIT("basic", test_parmap_basic, "Basic usage")
459 xbt_test_add("Create the parmap");
462 xbt_dynar_t data = xbt_dynar_new(sizeof(void *), NULL);
464 /* Create the parallel map */
466 parmap = xbt_parmap_new(10, XBT_PARMAP_FUTEX);
468 parmap = xbt_parmap_new(10, XBT_PARMAP_BUSY_WAIT);
470 for (j = 0; j < 100; j++) {
471 xbt_dynar_push_as(data, void *, (void *)j);
474 for (i = 0; i < 5; i++) {
475 xbt_parmap_apply(parmap, fun, data);
478 /* Destroy the parmap */
479 xbt_parmap_destroy(parmap);
480 xbt_dynar_free(&data);
483 #endif /* SIMGRID_TEST */