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"
9 #include <sys/syscall.h>
13 #include <linux/futex.h>
15 #include "xbt/xbt_os_thread.h"
18 #include "parmap_private.h"
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");
23 static void *_xbt_parmap_worker_main(void *parmap);
25 static void futex_wait(int *uaddr, int val);
26 static void futex_wake(int *uaddr, int val);
28 xbt_parmap_t xbt_parmap_new(unsigned int num_workers)
31 xbt_os_thread_t worker = NULL;
33 XBT_DEBUG("Create new parmap (%u workers)", num_workers);
35 /* Initialize the thread pool data structure */
36 xbt_parmap_t parmap = xbt_new0(s_xbt_parmap_t, 1);
38 parmap->sync_event = xbt_new0(s_xbt_event_t, 1);
40 parmap->num_workers = num_workers;
41 parmap->status = PARMAP_WORK;
43 parmap->sync_event->threads_to_wait = num_workers;
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);
51 xbt_event_init(parmap->sync_event);
56 void xbt_parmap_destroy(xbt_parmap_t parmap)
58 XBT_DEBUG("Destroy parmap %p", parmap);
59 parmap->status = PARMAP_DESTROY;
61 xbt_event_signal(parmap->sync_event);
62 xbt_free(parmap->sync_event);
67 void xbt_parmap_apply(xbt_parmap_t parmap, void_f_pvoid_t fun, xbt_dynar_t data)
69 /* Assign resources to worker threads*/
73 xbt_event_signal(parmap->sync_event);
75 XBT_DEBUG("Job done");
78 static void *_xbt_parmap_worker_main(void *arg)
80 unsigned int data_start, data_end, data_size, worker_id;
81 xbt_parmap_t parmap = (xbt_parmap_t)arg;
83 /* Fetch a worker id */
84 worker_id = __sync_fetch_and_add(&parmap->workers_max_id, 1);
85 xbt_os_thread_set_extra_data((void *)(unsigned long)worker_id);
87 XBT_DEBUG("New worker thread created (%u)", worker_id);
89 /* Worker's main loop */
92 xbt_event_wait(parmap->sync_event);
94 if(parmap->status == PARMAP_WORK){
95 XBT_DEBUG("Worker %u got a job", worker_id);
97 /* Compute how much data does every worker gets */
98 data_size = (xbt_dynar_length(parmap->data) / parmap->num_workers)
99 + ((xbt_dynar_length(parmap->data) % parmap->num_workers) ? 1 : 0);
101 /* Each worker data segment starts in a position associated with its id*/
102 data_start = data_size * worker_id;
104 /* The end of the worker data segment must be bounded by the end of the data vector */
105 data_end = MIN(data_start + data_size, xbt_dynar_length(parmap->data));
107 XBT_DEBUG("Worker %u: data_start=%u data_end=%u (data_size=%u)",
108 worker_id, data_start, data_end, data_size);
110 /* While the worker don't pass the end of it data segment apply the function */
111 while(data_start < data_end){
112 parmap->fun(*(void **)xbt_dynar_get_ptr(parmap->data, data_start));
116 /* We are destroying the parmap */
119 xbt_event_end(parmap->sync_event);
121 XBT_DEBUG("Shutting down worker %u", worker_id);
128 static void futex_wait(int *uaddr, int val)
130 XBT_VERB("Waiting on futex %p", uaddr);
131 syscall(SYS_futex, uaddr, FUTEX_WAIT_PRIVATE, val, NULL, NULL, 0);
134 static void futex_wake(int *uaddr, int val)
136 XBT_VERB("Waking futex %p", uaddr);
137 syscall(SYS_futex, uaddr, FUTEX_WAKE_PRIVATE, val, NULL, NULL, 0);
140 void xbt_event_init(xbt_event_t event)
142 int myflag = event->done;
143 if(event->thread_counter < event->threads_to_wait)
144 futex_wait(&event->done, myflag);
147 void xbt_event_signal(xbt_event_t event)
149 int myflag = event->done;
150 event->thread_counter = 0;
152 futex_wake(&event->work, event->threads_to_wait);
153 futex_wait(&event->done, myflag);
156 void xbt_event_wait(xbt_event_t event)
159 unsigned int mycount;
161 myflag = event->work;
162 mycount = __sync_add_and_fetch(&event->thread_counter, 1);
163 if(mycount == event->threads_to_wait){
165 futex_wake(&event->done, 1);
168 futex_wait(&event->work, myflag);
171 void xbt_event_end(xbt_event_t event)
174 unsigned int mycount;
176 myflag = event->work;
177 mycount = __sync_add_and_fetch(&event->thread_counter, 1);
178 if(mycount == event->threads_to_wait){
180 futex_wake(&event->done, 1);
189 XBT_TEST_SUITE("parmap", "Parallel Map");
190 XBT_LOG_EXTERNAL_DEFAULT_CATEGORY(xbt_parmap_unit);
200 //XBT_INFO("I'm job %lu", (unsigned long)arg);
203 XBT_TEST_UNIT("basic", test_parmap_basic, "Basic usage")
205 xbt_test_add("Create the parmap");
208 xbt_dynar_t data = xbt_dynar_new(sizeof(void *), NULL);
210 /* Create the parallel map */
211 parmap = xbt_parmap_new(10);
213 for(j=0; j < 100; j++){
214 xbt_dynar_push_as(data, void *, (void *)j);
218 xbt_parmap_apply(parmap, fun, data);
220 /* Destroy the parmap */
221 xbt_parmap_destroy(parmap);
222 xbt_dynar_free(&data);
225 #endif /* SIMGRID_TEST */