1 /* Copyright (c) 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. */
9 #include "xbt/sysdep.h"
11 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(smpi_bench, smpi,
12 "Logging specific to SMPI (benchmarking)");
14 xbt_dict_t allocs = NULL; /* Allocated on first use */
15 xbt_dict_t samples = NULL; /* Allocated on first use */
29 void smpi_bench_destroy(void)
32 xbt_dict_free(&allocs);
35 xbt_dict_free(&samples);
39 static void smpi_execute_flops(double flops)
45 e_surf_action_state_t state;
47 host = SIMIX_host_self();
48 mutex = SIMIX_mutex_init();
49 cond = SIMIX_cond_init();
50 DEBUG1("Handle real computation time: %f flops", flops);
51 action = SIMIX_action_execute(host, "computation", flops);
52 SIMIX_mutex_lock(mutex);
53 SIMIX_register_action_to_condition(action, cond);
54 for (state = SIMIX_action_get_state(action);
55 state == SURF_ACTION_READY ||
56 state == SURF_ACTION_RUNNING;
57 state = SIMIX_action_get_state(action)) {
58 SIMIX_cond_wait(cond, mutex);
60 SIMIX_unregister_action_to_condition(action, cond);
61 SIMIX_mutex_unlock(mutex);
62 SIMIX_action_destroy(action);
63 SIMIX_cond_destroy(cond);
64 SIMIX_mutex_destroy(mutex);
67 static void smpi_execute(double duration)
69 if (duration >= xbt_cfg_get_double(_surf_cfg_set, "smpi/cpu_threshold")) {
70 DEBUG1("Sleep for %f to handle real computation time", duration);
71 smpi_execute_flops(duration *
72 xbt_cfg_get_double(_surf_cfg_set,
73 "smpi/running_power"));
77 void smpi_bench_begin(void)
79 xbt_os_timer_start(smpi_process_timer());
82 void smpi_bench_end(void)
84 xbt_os_timer_t timer = smpi_process_timer();
86 xbt_os_timer_stop(timer);
87 smpi_execute(xbt_os_timer_elapsed(timer));
90 unsigned int smpi_sleep(unsigned int secs)
92 smpi_execute((double) secs);
96 int smpi_gettimeofday(struct timeval *tv, struct timezone *tz)
98 double now = SIMIX_get_clock();
101 tv->tv_sec = (time_t) now;
102 tv->tv_usec = (suseconds_t) (now * 1e6);
107 static char *sample_location(int global, const char *file, int line)
110 return bprintf("%s:%d", file, line);
112 return bprintf("%s:%d:%d", file, line, smpi_process_index());
116 void smpi_sample_1(int global, const char *file, int line, int max)
118 char *loc = sample_location(global, file, line);
121 smpi_bench_end(); /* Take time from previous MPI call into account */
123 samples = xbt_dict_new();
125 data = xbt_dict_get_or_null(samples, loc);
127 data = (local_data_t *) xbt_new(local_data_t, 1);
132 xbt_dict_set(samples, loc, data, &free);
137 int smpi_sample_2(int global, const char *file, int line)
139 char *loc = sample_location(global, file, line);
142 xbt_assert0(samples, "You did something very inconsistent, didn't you?");
143 data = xbt_dict_get_or_null(samples, loc);
145 xbt_assert0(data, "Please, do thing in order");
147 if (!data->started) {
148 if (data->count < data->max) {
152 DEBUG1("Perform some wait of %f", data->time / (double) data->count);
153 smpi_execute(data->time / (double) data->count);
160 smpi_process_simulated_start();
161 return data->started;
164 void smpi_sample_3(int global, const char *file, int line)
166 char *loc = sample_location(global, file, line);
169 xbt_assert0(samples, "You did something very inconsistent, didn't you?");
170 data = xbt_dict_get_or_null(samples, loc);
171 if (!data || !data->started || data->count >= data->max) {
172 xbt_assert0(data, "Please, do thing in order");
175 data->time += smpi_process_simulated_elapsed();
176 DEBUG2("Average mean after %d steps is %f", data->count,
177 data->time / (double) data->count);
180 void smpi_sample_flops(double flops)
182 smpi_execute_flops(flops);
185 void *smpi_shared_malloc(size_t size, const char *file, int line)
187 char *loc = bprintf("%s:%d:%zu", file, line, size);
191 allocs = xbt_dict_new();
193 data = xbt_dict_get_or_null(allocs, loc);
195 data = (shared_data_t *) xbt_malloc0(sizeof(int) + size);
197 xbt_dict_set(allocs, loc, data, &free);
205 void smpi_shared_free(void *ptr)
207 shared_data_t *data = (shared_data_t *) ((int *) ptr - 1);
211 WARN0("Cannot free: nothing was allocated");
214 loc = xbt_dict_get_key(allocs, data);
216 WARN1("Cannot free: %p was not shared-allocated by SMPI", ptr);
220 if (data->count <= 0) {
221 xbt_dict_remove(allocs, loc);