-void smpi_sample_3(int global, const char* file, int line) {
- char* loc = sample_location(global, file, line);
- local_data_t* data;
-
- xbt_assert0(samples, "You did something very inconsistent, didn't you?");
- data = xbt_dict_get_or_null(samples, loc);
- if (!data || !data->started || data->count >= data->max) {
- xbt_assert0(data, "Please, do thing in order");
- }
- smpi_bench_end(-1, NULL);
- data->time += smpi_process_simulated_elapsed();
- DEBUG2("Average mean after %d steps is %f", data->count, data->time / (double)data->count);
+void smpi_sample_3(int global, const char *file, int line)
+{
+ char *loc = sample_location(global, file, line);
+ local_data_t *data;
+ double sample, n;
+
+ xbt_assert0(samples, "You did something very inconsistent, didn't you?");
+ data = xbt_dict_get_or_null(samples, loc);
+ if (!data || !data->started || data->count >= data->iters) {
+ xbt_assert0(data, "Please, do thing in order");
+ }
+ smpi_bench_end();
+ sample = smpi_process_simulated_elapsed();
+ data->sum += sample;
+ data->sum_pow2 += sample * sample;
+ n = (double)data->count;
+ data->mean = data->sum / n;
+ data->relstderr = sqrt((data->sum_pow2 / n - data->mean * data->mean) / n) / data->mean;
+ DEBUG4("Average mean after %d steps is %f, relative standard error is %f (sample was %f)", data->count,
+ data->mean, data->relstderr, sample);