stat.c

   1 #include <stdio.h>
   2 #include <string.h>
   3 #include <sys/time.h>
   4 #include <sys/types.h>
   5 #include <sys/stat.h>
   6 #include <dirent.h>
   7 #include <libgen.h>
   8 #include <math.h>
   9
  10 #include "fio.h"
  11
  12 static struct itimerval itimer;
  13 static struct list_head disk_list = LIST_HEAD_INIT(disk_list);
  14
  15 /*
  16  * Cheasy number->string conversion, complete with carry rounding error.
  17  */
  18 static char *num2str(unsigned long num, int maxlen, int base)
  19 {
  20         /*
  21          * could be passed in for 10^3 base, but every caller expects
  22          * 2^10 base right now.
  23          */
  24         const unsigned int thousand = 1024;
  25         char postfix[] = { 'K', 'M', 'G', 'P', 'E' };
  26         char *buf;
  27         int i;
  28
  29         buf = malloc(128);
  30
  31         for (i = 0; base > 1; i++)
  32                 base /= thousand;
  33
  34         do {
  35                 int len, carry = 0;
  36
  37                 len = sprintf(buf, "%'lu", num);
  38                 if (len <= maxlen) {
  39                         buf[len] = postfix[i];
  40                         buf[len + 1] = '\0';
  41                         return buf;
  42                 }
  43
  44                 if ((num % thousand) >= (thousand / 2))
  45                         carry = 1;
  46
  47                 num /= thousand;
  48                 num += carry;
  49                 i++;
  50         } while (i <= 5);
  51
  52         return buf;
  53 }
  54
  55 static int get_io_ticks(struct disk_util *du, struct disk_util_stat *dus)
  56 {
  57         unsigned in_flight;
  58         char line[256];
  59         FILE *f;
  60         char *p;
  61
  62         f = fopen(du->path, "r");
  63         if (!f)
  64                 return 1;
  65
  66         p = fgets(line, sizeof(line), f);
  67         if (!p) {
  68                 fclose(f);
  69                 return 1;
  70         }
  71
  72         if (sscanf(p, "%u %u %llu %u %u %u %llu %u %u %u %u\n", &dus->ios[0], &dus->merges[0], &dus->sectors[0], &dus->ticks[0], &dus->ios[1], &dus->merges[1], &dus->sectors[1], &dus->ticks[1], &in_flight, &dus->io_ticks, &dus->time_in_queue) != 11) {
  73                 fclose(f);
  74                 return 1;
  75         }
  76
  77         fclose(f);
  78         return 0;
  79 }
  80
  81 static void update_io_tick_disk(struct disk_util *du)
  82 {
  83         struct disk_util_stat __dus, *dus, *ldus;
  84         struct timeval t;
  85
  86         if (get_io_ticks(du, &__dus))
  87                 return;
  88
  89         dus = &du->dus;
  90         ldus = &du->last_dus;
  91
  92         dus->sectors[0] += (__dus.sectors[0] - ldus->sectors[0]);
  93         dus->sectors[1] += (__dus.sectors[1] - ldus->sectors[1]);
  94         dus->ios[0] += (__dus.ios[0] - ldus->ios[0]);
  95         dus->ios[1] += (__dus.ios[1] - ldus->ios[1]);
  96         dus->merges[0] += (__dus.merges[0] - ldus->merges[0]);
  97         dus->merges[1] += (__dus.merges[1] - ldus->merges[1]);
  98         dus->ticks[0] += (__dus.ticks[0] - ldus->ticks[0]);
  99         dus->ticks[1] += (__dus.ticks[1] - ldus->ticks[1]);
 100         dus->io_ticks += (__dus.io_ticks - ldus->io_ticks);
 101         dus->time_in_queue += (__dus.time_in_queue - ldus->time_in_queue);
 102
 103         fio_gettime(&t, NULL);
 104         du->msec += mtime_since(&du->time, &t);
 105         memcpy(&du->time, &t, sizeof(t));
 106         memcpy(ldus, &__dus, sizeof(__dus));
 107 }
 108
 109 void update_io_ticks(void)
 110 {
 111         struct list_head *entry;
 112         struct disk_util *du;
 113
 114         list_for_each(entry, &disk_list) {
 115                 du = list_entry(entry, struct disk_util, list);
 116                 update_io_tick_disk(du);
 117         }
 118 }
 119
 120 static int disk_util_exists(dev_t dev)
 121 {
 122         struct list_head *entry;
 123         struct disk_util *du;
 124
 125         list_for_each(entry, &disk_list) {
 126                 du = list_entry(entry, struct disk_util, list);
 127
 128                 if (du->dev == dev)
 129                         return 1;
 130         }
 131
 132         return 0;
 133 }
 134
 135 static void disk_util_add(dev_t dev, char *path)
 136 {
 137         struct disk_util *du, *__du;
 138         struct list_head *entry;
 139
 140         du = malloc(sizeof(*du));
 141         memset(du, 0, sizeof(*du));
 142         INIT_LIST_HEAD(&du->list);
 143         sprintf(du->path, "%s/stat", path);
 144         du->name = strdup(basename(path));
 145         du->dev = dev;
 146
 147         list_for_each(entry, &disk_list) {
 148                 __du = list_entry(entry, struct disk_util, list);
 149
 150                 if (!strcmp(du->name, __du->name)) {
 151                         free(du->name);
 152                         free(du);
 153                         return;
 154                 }
 155         }
 156
 157         fio_gettime(&du->time, NULL);
 158         get_io_ticks(du, &du->last_dus);
 159
 160         list_add_tail(&du->list, &disk_list);
 161 }
 162
 163 static int check_dev_match(dev_t dev, char *path)
 164 {
 165         unsigned int major, minor;
 166         char line[256], *p;
 167         FILE *f;
 168
 169         f = fopen(path, "r");
 170         if (!f) {
 171                 perror("open path");
 172                 return 1;
 173         }
 174
 175         p = fgets(line, sizeof(line), f);
 176         if (!p) {
 177                 fclose(f);
 178                 return 1;
 179         }
 180
 181         if (sscanf(p, "%u:%u", &major, &minor) != 2) {
 182                 fclose(f);
 183                 return 1;
 184         }
 185
 186         if (((major << 8) | minor) == dev) {
 187                 fclose(f);
 188                 return 0;
 189         }
 190
 191         fclose(f);
 192         return 1;
 193 }
 194
 195 static int find_block_dir(dev_t dev, char *path)
 196 {
 197         struct dirent *dir;
 198         struct stat st;
 199         int found = 0;
 200         DIR *D;
 201
 202         D = opendir(path);
 203         if (!D)
 204                 return 0;
 205
 206         while ((dir = readdir(D)) != NULL) {
 207                 char full_path[256];
 208
 209                 if (!strcmp(dir->d_name, ".") || !strcmp(dir->d_name, ".."))
 210                         continue;
 211
 212                 sprintf(full_path, "%s/%s", path, dir->d_name);
 213
 214                 if (!strcmp(dir->d_name, "dev")) {
 215                         if (!check_dev_match(dev, full_path)) {
 216                                 found = 1;
 217                                 break;
 218                         }
 219                 }
 220
 221                 if (lstat(full_path, &st) == -1) {
 222                         perror("stat");
 223                         break;
 224                 }
 225
 226                 if (!S_ISDIR(st.st_mode) || S_ISLNK(st.st_mode))
 227                         continue;
 228
 229                 found = find_block_dir(dev, full_path);
 230                 if (found) {
 231                         strcpy(path, full_path);
 232                         break;
 233                 }
 234         }
 235
 236         closedir(D);
 237         return found;
 238 }
 239
 240 void init_disk_util(struct thread_data *td)
 241 {
 242         struct fio_file *f;
 243         struct stat st;
 244         char foo[PATH_MAX], tmp[PATH_MAX];
 245         dev_t dev;
 246         char *p;
 247
 248         if (!td->do_disk_util || (td->io_ops->flags & (FIO_NETIO | FIO_NULLIO)))
 249                 return;
 250
 251         /*
 252          * Just use the same file, they are on the same device.
 253          */
 254         f = &td->files[0];
 255         if (!stat(f->file_name, &st)) {
 256                 if (S_ISBLK(st.st_mode))
 257                         dev = st.st_rdev;
 258                 else
 259                         dev = st.st_dev;
 260         } else {
 261                 /*
 262                  * must be a file, open "." in that path
 263                  */
 264                 strncpy(foo, f->file_name, PATH_MAX - 1);
 265                 p = dirname(foo);
 266                 if (stat(p, &st)) {
 267                         perror("disk util stat");
 268                         return;
 269                 }
 270
 271                 dev = st.st_dev;
 272         }
 273
 274         if (disk_util_exists(dev))
 275                 return;
 276
 277         sprintf(foo, "/sys/block");
 278         if (!find_block_dir(dev, foo))
 279                 return;
 280
 281         /*
 282          * If there's a ../queue/ directory there, we are inside a partition.
 283          * Check if that is the case and jump back. For loop/md/dm etc we
 284          * are already in the right spot.
 285          */
 286         sprintf(tmp, "%s/../queue", foo);
 287         if (!stat(tmp, &st)) {
 288                 p = dirname(foo);
 289                 sprintf(tmp, "%s/queue", p);
 290                 if (stat(tmp, &st)) {
 291                         log_err("unknown sysfs layout\n");
 292                         return;
 293                 }
 294                 strncpy(tmp, p, PATH_MAX - 1);
 295                 sprintf(foo, "%s", tmp);
 296         }
 297
 298         if (td->ioscheduler)
 299                 td->sysfs_root = strdup(foo);
 300
 301         disk_util_add(dev, foo);
 302 }
 303
 304 void disk_util_timer_arm(void)
 305 {
 306         itimer.it_value.tv_sec = 0;
 307         itimer.it_value.tv_usec = DISK_UTIL_MSEC * 1000;
 308         setitimer(ITIMER_REAL, &itimer, NULL);
 309 }
 310
 311 void update_rusage_stat(struct thread_data *td)
 312 {
 313         struct thread_stat *ts = &td->ts;
 314
 315         getrusage(RUSAGE_SELF, &ts->ru_end);
 316
 317         ts->usr_time += mtime_since(&ts->ru_start.ru_utime, &ts->ru_end.ru_utime);
 318         ts->sys_time += mtime_since(&ts->ru_start.ru_stime, &ts->ru_end.ru_stime);
 319         ts->ctx += ts->ru_end.ru_nvcsw + ts->ru_end.ru_nivcsw - (ts->ru_start.ru_nvcsw + ts->ru_start.ru_nivcsw);
 320
 321         memcpy(&ts->ru_start, &ts->ru_end, sizeof(ts->ru_end));
 322 }
 323
 324 static int calc_lat(struct io_stat *is, unsigned long *min, unsigned long *max,
 325                     double *mean, double *dev)
 326 {
 327         double n = is->samples;
 328
 329         if (is->samples == 0)
 330                 return 0;
 331
 332         *min = is->min_val;
 333         *max = is->max_val;
 334
 335         n = (double) is->samples;
 336         *mean = is->mean;
 337
 338         if (n > 1.0)
 339                 *dev = sqrt(is->S / (n - 1.0));
 340         else
 341                 *dev = -1.0;
 342
 343         return 1;
 344 }
 345
 346 static void show_group_stats(struct group_run_stats *rs, int id)
 347 {
 348         char *p1, *p2, *p3, *p4;
 349         const char *ddir_str[] = { "   READ", "  WRITE" };
 350         int i;
 351
 352         fprintf(f_out, "\nRun status group %d (all jobs):\n", id);
 353
 354         for (i = 0; i <= DDIR_WRITE; i++) {
 355                 if (!rs->max_run[i])
 356                         continue;
 357
 358                 p1 = num2str(rs->io_kb[i], 6, 1);
 359                 p2 = num2str(rs->agg[i], 6, 1);
 360                 p3 = num2str(rs->min_bw[i], 6, 1);
 361                 p4 = num2str(rs->max_bw[i], 6, 1);
 362
 363                 fprintf(f_out, "%s: io=%siB, aggrb=%siB/s, minb=%siB/s, maxb=%siB/s, mint=%llumsec, maxt=%llumsec\n", ddir_str[i], p1, p2, p3, p4, rs->min_run[0], rs->max_run[0]);
 364
 365                 free(p1);
 366                 free(p2);
 367                 free(p3);
 368                 free(p4);
 369         }
 370 }
 371
 372 static void show_disk_util(void)
 373 {
 374         struct disk_util_stat *dus;
 375         struct list_head *entry, *next;
 376         struct disk_util *du;
 377         double util;
 378
 379         fprintf(f_out, "\nDisk stats (read/write):\n");
 380
 381         list_for_each(entry, &disk_list) {
 382                 du = list_entry(entry, struct disk_util, list);
 383                 dus = &du->dus;
 384
 385                 util = (double) 100 * du->dus.io_ticks / (double) du->msec;
 386                 if (util > 100.0)
 387                         util = 100.0;
 388
 389                 fprintf(f_out, "  %s: ios=%u/%u, merge=%u/%u, ticks=%u/%u, in_queue=%u, util=%3.2f%%\n", du->name, dus->ios[0], dus->ios[1], dus->merges[0], dus->merges[1], dus->ticks[0], dus->ticks[1], dus->time_in_queue, util);
 390         }
 391
 392         /*
 393          * now free the list
 394          */
 395         list_for_each_safe(entry, next, &disk_list) {
 396                 list_del(entry);
 397                 du = list_entry(entry, struct disk_util, list);
 398                 free(du->name);
 399                 free(du);
 400         }
 401 }
 402
 403 static void show_ddir_status(struct thread_data *td, struct group_run_stats *rs,
 404                              int ddir)
 405 {
 406         const char *ddir_str[] = { "read ", "write" };
 407         struct thread_stat *ts;
 408         unsigned long min, max;
 409         unsigned long long bw;
 410         double mean, dev;
 411         char *io_p, *bw_p;
 412
 413         if (!td->runtime[ddir])
 414                 return;
 415
 416         bw = td->io_bytes[ddir] / td->runtime[ddir];
 417         io_p = num2str(td->io_bytes[ddir] >> 10, 6, 1);
 418         bw_p = num2str(bw, 6, 1);
 419
 420         fprintf(f_out, "  %s: io=%siB, bw=%siB/s, runt=%6lumsec\n", ddir_str[ddir], io_p, bw_p, td->runtime[ddir]);
 421
 422         free(io_p);
 423         free(bw_p);
 424
 425         ts = &td->ts;
 426         if (calc_lat(&ts->slat_stat[ddir], &min, &max, &mean, &dev))
 427                 fprintf(f_out, "    slat (msec): min=%5lu, max=%5lu, avg=%5.02f, stdev=%5.02f\n", min, max, mean, dev);
 428
 429         if (calc_lat(&ts->clat_stat[ddir], &min, &max, &mean, &dev))
 430                 fprintf(f_out, "    clat (msec): min=%5lu, max=%5lu, avg=%5.02f, stdev=%5.02f\n", min, max, mean, dev);
 431
 432         if (calc_lat(&ts->bw_stat[ddir], &min, &max, &mean, &dev)) {
 433                 double p_of_agg;
 434
 435                 p_of_agg = mean * 100 / (double) rs->agg[ddir];
 436                 fprintf(f_out, "    bw (KiB/s) : min=%5lu, max=%5lu, per=%3.2f%%, avg=%5.02f, stdev=%5.02f\n", min, max, p_of_agg, mean, dev);
 437         }
 438 }
 439
 440 static void show_thread_status(struct thread_data *td,
 441                                struct group_run_stats *rs)
 442 {
 443         double usr_cpu, sys_cpu;
 444         unsigned long runtime;
 445         double io_u_dist[FIO_IO_U_MAP_NR];
 446         double io_u_lat[FIO_IO_U_LAT_NR];
 447         int i;
 448
 449         if (!(td->io_bytes[0] + td->io_bytes[1]))
 450                 return;
 451
 452         if (!td->error)
 453                 fprintf(f_out, "%s: (groupid=%d): err=%2d: pid=%d\n",td->name, td->groupid, td->error, td->pid);
 454         else
 455                 fprintf(f_out, "%s: (groupid=%d): err=%2d (%s): pid=%d\n",td->name, td->groupid, td->error, td->verror, td->pid);
 456
 457         show_ddir_status(td, rs, td->ddir);
 458         if (td->io_bytes[td->ddir ^ 1])
 459                 show_ddir_status(td, rs, td->ddir ^ 1);
 460
 461         runtime = mtime_since(&td->epoch, &td->end_time);
 462         if (runtime) {
 463                 double runt = (double) runtime;
 464
 465                 usr_cpu = (double) td->ts.usr_time * 100 / runt;
 466                 sys_cpu = (double) td->ts.sys_time * 100 / runt;
 467         } else {
 468                 usr_cpu = 0;
 469                 sys_cpu = 0;
 470         }
 471
 472         fprintf(f_out, "  cpu          : usr=%3.2f%%, sys=%3.2f%%, ctx=%lu\n", usr_cpu, sys_cpu, td->ts.ctx);
 473
 474         /*
 475          * Do depth distribution calculations
 476          */
 477         for (i = 0; i < FIO_IO_U_MAP_NR; i++) {
 478                 io_u_dist[i] = (double) td->io_u_map[i] / (double) td->total_io_u;
 479                 io_u_dist[i] *= 100.0;
 480         }
 481
 482         fprintf(f_out, "  IO depths    : 1=%3.1f%%, 2=%3.1f%%, 4=%3.1f%%, 8=%3.1f%%, 16=%3.1f%%, 32=%3.1f%%, >=64=%3.1f%%\n", io_u_dist[0], io_u_dist[1], io_u_dist[2], io_u_dist[3], io_u_dist[4], io_u_dist[5], io_u_dist[6]);
 483
 484         /*
 485          * Do latency distribution calculations
 486          */
 487         for (i = 0; i < FIO_IO_U_LAT_NR; i++) {
 488                 io_u_lat[i] = (double) td->io_u_lat[i] / (double) td->total_io_u;
 489                 io_u_lat[i] *= 100.0;
 490         }
 491
 492         fprintf(f_out, "     lat (msec): 2=%3.1f%%, 4=%3.1f%%, 10=%3.1f%%, 20=%3.1f%%, 50=%3.1f%%, 100=%3.1f%%\n", io_u_lat[0], io_u_lat[1], io_u_lat[2], io_u_lat[3], io_u_lat[4], io_u_lat[5]);
 493         fprintf(f_out, "     lat (msec): 250=%3.1f%%, 500=%3.1f%%, 750=%3.1f%%, 1000=%3.1f%%, >=2000=%3.1f%%\n", io_u_lat[6], io_u_lat[7], io_u_lat[8], io_u_lat[9], io_u_lat[10]);
 494
 495         if (td->description)
 496                 fprintf(f_out, "%s\n", td->description);
 497 }
 498
 499 static void show_ddir_status_terse(struct thread_data *td,
 500                                    struct group_run_stats *rs, int ddir)
 501 {
 502         struct thread_stat *ts = &td->ts;
 503         unsigned long min, max;
 504         unsigned long long bw;
 505         double mean, dev;
 506
 507         bw = 0;
 508         if (td->runtime[ddir])
 509                 bw = td->io_bytes[ddir] / td->runtime[ddir];
 510
 511         fprintf(f_out, ",%llu,%llu,%lu", td->io_bytes[ddir] >> 10, bw, td->runtime[ddir]);
 512
 513         if (calc_lat(&ts->slat_stat[ddir], &min, &max, &mean, &dev))
 514                 fprintf(f_out, ",%lu,%lu,%f,%f", min, max, mean, dev);
 515         else
 516                 fprintf(f_out, ",%lu,%lu,%f,%f", 0UL, 0UL, 0.0, 0.0);
 517
 518         if (calc_lat(&ts->clat_stat[ddir], &min, &max, &mean, &dev))
 519                 fprintf(f_out, ",%lu,%lu,%f,%f", min, max, mean, dev);
 520         else
 521                 fprintf(f_out, ",%lu,%lu,%f,%f", 0UL, 0UL, 0.0, 0.0);
 522
 523         if (calc_lat(&ts->bw_stat[ddir], &min, &max, &mean, &dev)) {
 524                 double p_of_agg;
 525
 526                 p_of_agg = mean * 100 / (double) rs->agg[ddir];
 527                 fprintf(f_out, ",%lu,%lu,%f%%,%f,%f", min, max, p_of_agg, mean, dev);
 528         } else
 529                 fprintf(f_out, ",%lu,%lu,%f%%,%f,%f", 0UL, 0UL, 0.0, 0.0, 0.0);
 530
 531 }
 532
 533
 534 static void show_thread_status_terse(struct thread_data *td,
 535                                      struct group_run_stats *rs)
 536 {
 537         double usr_cpu, sys_cpu;
 538
 539         fprintf(f_out, "%s,%d,%d",td->name, td->groupid, td->error);
 540
 541         show_ddir_status_terse(td, rs, 0);
 542         show_ddir_status_terse(td, rs, 1);
 543
 544         if (td->runtime[0] + td->runtime[1]) {
 545                 double runt = (double) (td->runtime[0] + td->runtime[1]);
 546
 547                 usr_cpu = (double) td->ts.usr_time * 100 / runt;
 548                 sys_cpu = (double) td->ts.sys_time * 100 / runt;
 549         } else {
 550                 usr_cpu = 0;
 551                 sys_cpu = 0;
 552         }
 553
 554         fprintf(f_out, ",%f%%,%f%%,%lu\n", usr_cpu, sys_cpu, td->ts.ctx);
 555 }
 556
 557 void show_run_stats(void)
 558 {
 559         struct group_run_stats *runstats, *rs;
 560         struct thread_data *td;
 561         int i;
 562
 563         runstats = malloc(sizeof(struct group_run_stats) * (groupid + 1));
 564
 565         for (i = 0; i < groupid + 1; i++) {
 566                 rs = &runstats[i];
 567
 568                 memset(rs, 0, sizeof(*rs));
 569                 rs->min_bw[0] = rs->min_run[0] = ~0UL;
 570                 rs->min_bw[1] = rs->min_run[1] = ~0UL;
 571         }
 572
 573         for_each_td(td, i) {
 574                 unsigned long long rbw, wbw;
 575
 576                 rs = &runstats[td->groupid];
 577
 578                 if (td->runtime[0] < rs->min_run[0] || !rs->min_run[0])
 579                         rs->min_run[0] = td->runtime[0];
 580                 if (td->runtime[0] > rs->max_run[0])
 581                         rs->max_run[0] = td->runtime[0];
 582                 if (td->runtime[1] < rs->min_run[1] || !rs->min_run[1])
 583                         rs->min_run[1] = td->runtime[1];
 584                 if (td->runtime[1] > rs->max_run[1])
 585                         rs->max_run[1] = td->runtime[1];
 586
 587                 rbw = wbw = 0;
 588                 if (td->runtime[0])
 589                         rbw = td->io_bytes[0] / (unsigned long long) td->runtime[0];
 590                 if (td->runtime[1])
 591                         wbw = td->io_bytes[1] / (unsigned long long) td->runtime[1];
 592
 593                 if (rbw < rs->min_bw[0])
 594                         rs->min_bw[0] = rbw;
 595                 if (wbw < rs->min_bw[1])
 596                         rs->min_bw[1] = wbw;
 597                 if (rbw > rs->max_bw[0])
 598                         rs->max_bw[0] = rbw;
 599                 if (wbw > rs->max_bw[1])
 600                         rs->max_bw[1] = wbw;
 601
 602                 rs->io_kb[0] += td->io_bytes[0] >> 10;
 603                 rs->io_kb[1] += td->io_bytes[1] >> 10;
 604         }
 605
 606         for (i = 0; i < groupid + 1; i++) {
 607                 rs = &runstats[i];
 608
 609                 if (rs->max_run[0])
 610                         rs->agg[0] = (rs->io_kb[0]*1024) / rs->max_run[0];
 611                 if (rs->max_run[1])
 612                         rs->agg[1] = (rs->io_kb[1]*1024) / rs->max_run[1];
 613         }
 614
 615         /*
 616          * don't overwrite last signal output
 617          */
 618         if (!terse_output)
 619                 printf("\n");
 620
 621         for_each_td(td, i) {
 622                 rs = &runstats[td->groupid];
 623
 624                 if (terse_output)
 625                         show_thread_status_terse(td, rs);
 626                 else
 627                         show_thread_status(td, rs);
 628         }
 629
 630         if (!terse_output) {
 631                 for (i = 0; i < groupid + 1; i++)
 632                         show_group_stats(&runstats[i], i);
 633
 634                 show_disk_util();
 635         }
 636
 637         free(runstats);
 638 }
 639
 640 static inline void add_stat_sample(struct io_stat *is, unsigned long data)
 641 {
 642         double val = data;
 643         double delta, n;
 644
 645         if (data > is->max_val)
 646                 is->max_val = data;
 647         if (data < is->min_val)
 648                 is->min_val = data;
 649
 650         delta = val - is->mean;
 651         n = is->samples + 1.0;
 652         is->mean += delta / n;
 653         is->S += delta * (val - is->mean);
 654
 655         is->samples++;
 656 }
 657
 658 static void __add_log_sample(struct io_log *iolog, unsigned long val,
 659                              enum fio_ddir ddir, unsigned long time)
 660 {
 661         if (iolog->nr_samples == iolog->max_samples) {
 662                 int new_size = sizeof(struct io_sample) * iolog->max_samples*2;
 663
 664                 iolog->log = realloc(iolog->log, new_size);
 665                 iolog->max_samples <<= 1;
 666         }
 667
 668         iolog->log[iolog->nr_samples].val = val;
 669         iolog->log[iolog->nr_samples].time = time;
 670         iolog->log[iolog->nr_samples].ddir = ddir;
 671         iolog->nr_samples++;
 672 }
 673
 674 static void add_log_sample(struct thread_data *td, struct io_log *iolog,
 675                            unsigned long val, enum fio_ddir ddir)
 676 {
 677         __add_log_sample(iolog, val, ddir, mtime_since_now(&td->epoch));
 678 }
 679
 680 void add_agg_sample(unsigned long val, enum fio_ddir ddir)
 681 {
 682         struct io_log *iolog = agg_io_log[ddir];
 683
 684         __add_log_sample(iolog, val, ddir, mtime_since_genesis());
 685 }
 686
 687 void add_clat_sample(struct thread_data *td, enum fio_ddir ddir,
 688                      unsigned long msec)
 689 {
 690         struct thread_stat *ts = &td->ts;
 691
 692         add_stat_sample(&ts->clat_stat[ddir], msec);
 693
 694         if (ts->clat_log)
 695                 add_log_sample(td, ts->clat_log, msec, ddir);
 696 }
 697
 698 void add_slat_sample(struct thread_data *td, enum fio_ddir ddir,
 699                      unsigned long msec)
 700 {
 701         struct thread_stat *ts = &td->ts;
 702
 703         add_stat_sample(&ts->slat_stat[ddir], msec);
 704
 705         if (ts->slat_log)
 706                 add_log_sample(td, ts->slat_log, msec, ddir);
 707 }
 708
 709 void add_bw_sample(struct thread_data *td, enum fio_ddir ddir,
 710                    struct timeval *t)
 711 {
 712         struct thread_stat *ts = &td->ts;
 713         unsigned long spent = mtime_since(&ts->stat_sample_time[ddir], t);
 714         unsigned long rate;
 715
 716         if (spent < td->bw_avg_time)
 717                 return;
 718
 719         rate = (td->this_io_bytes[ddir] - ts->stat_io_bytes[ddir]) / spent;
 720         add_stat_sample(&ts->bw_stat[ddir], rate);
 721
 722         if (ts->bw_log)
 723                 add_log_sample(td, ts->bw_log, rate, ddir);
 724
 725         fio_gettime(&ts->stat_sample_time[ddir], NULL);
 726         ts->stat_io_bytes[ddir] = td->this_io_bytes[ddir];
 727 }