2 #include "steadystate.h"
3 #include "helper_thread.h"
5 void steadystate_setup()
8 struct thread_data *td, *prev_td;
14 * if group reporting is enabled, identify the last td
15 * for each group and use it for storing steady state
21 if (!td->o.group_reporting)
24 if (prev_groupid != td->groupid) {
26 prev_td->ss.last_in_group = 1;
27 prev_groupid = td->groupid;
32 if (prev_td != NULL && prev_td->o.group_reporting)
33 prev_td->ss.last_in_group = 1;
36 void steadystate_check()
38 int i, j, ddir, prev_groupid, group_ramp_time_over = 0;
39 unsigned long rate_time;
40 struct thread_data *td, *td2;
42 unsigned long group_bw = 0, group_iops = 0;
43 unsigned long long td_iops;
44 unsigned long long td_bytes;
48 struct steadystate_data *ss = &td->ss;
50 if (!ss->dur || td->runstate <= TD_SETTING_UP || td->runstate >= TD_EXITED || ss->attained)
55 if (!td->o.group_reporting ||
56 (td->o.group_reporting && td->groupid != prev_groupid)) {
59 group_ramp_time_over = 0;
61 prev_groupid = td->groupid;
63 fio_gettime(&now, NULL);
64 if (ss->ramp_time && !ss->ramp_time_over)
66 * Begin recording data one second after ss->ramp_time
69 if (utime_since(&td->epoch, &now) >= (ss->ramp_time + 1000000L))
70 ss->ramp_time_over = 1;
72 for (ddir = DDIR_READ; ddir < DDIR_RWDIR_CNT; ddir++) {
73 td_iops += td->io_blocks[ddir];
74 td_bytes += td->io_bytes[ddir];
77 rate_time = mtime_since(&ss->prev_time, &now);
78 memcpy(&ss->prev_time, &now, sizeof(now));
81 * Begin monitoring when job starts but don't actually use
82 * data in checking stopping criterion until ss->ramp_time is
83 * over. This ensures that we will have a sane value in
84 * prev_iops/bw the first time through after ss->ramp_time
87 if (ss->ramp_time_over) {
88 group_bw += 1000 * (td_bytes - ss->prev_bytes) / rate_time;
89 group_iops += 1000 * (td_iops - ss->prev_iops) / rate_time;
90 ++group_ramp_time_over;
92 ss->prev_iops = td_iops;
93 ss->prev_bytes = td_bytes;
95 if (td->o.group_reporting && !ss->last_in_group)
98 /* don't begin checking criterion until ss->ramp_time is over for at least one thread in group */
99 if (!group_ramp_time_over)
102 dprint(FD_STEADYSTATE, "steadystate_check() thread: %d, groupid: %u, rate_msec: %ld, iops: %lu, bw: %lu, head: %d, tail: %d\n",
103 i, td->groupid, rate_time, group_iops, group_bw, ss->head, ss->tail);
105 if (ss->evaluate(group_iops, group_bw, td))
107 if (td->o.group_reporting)
108 for_each_td(td2, j) {
109 if (td2->groupid == td->groupid) {
110 td2->ss.attained = 1;
111 fio_mark_td_terminate(td2);
116 fio_mark_td_terminate(td);
122 bool steadystate_slope(unsigned long iops, unsigned long bw, struct thread_data *td)
127 struct steadystate_data *ss = &td->ss;
129 ss->cache[ss->tail] = ss->check_iops ? iops : bw;
131 if (ss->tail < ss->head || (ss->tail - ss->head == ss->dur - 1))
133 if (ss->sum_y == 0) /* first time through */
135 for(i = 0; i < ss->dur; i++)
137 ss->sum_y += ss->cache[i];
141 ss->sum_xy += ss->cache[x] * i;
143 } else { /* easy to update the sums */
144 ss->sum_y -= ss->oldest_y;
145 ss->sum_y += ss->cache[ss->tail];
146 ss->sum_xy = ss->sum_xy - ss->sum_y + ss->dur * ss->cache[ss->tail];
149 ss->oldest_y = ss->cache[ss->head];
152 * calculate slope as (sum_xy - sum_x * sum_y / n) / (sum_(x^2) - (sum_x)^2 / n)
153 * This code assumes that all x values are equally spaced when they are often
154 * off by a few milliseconds. This assumption greatly simplifies the
157 slope = (ss->sum_xy - (double) ss->sum_x * ss->sum_y / ss->dur) / (ss->sum_x_sq - (double) ss->sum_x * ss->sum_x / ss->dur);
158 ss->criterion = ss->pct ? slope / (ss->sum_y / ss->dur) * 100.0: slope;
160 dprint(FD_STEADYSTATE, "sum_y: %llu, sum_xy: %llu, slope: %f, criterion: %f, limit: %f\n",
161 ss->sum_y, ss->sum_xy, slope, ss->criterion, ss->limit);
163 result = ss->criterion * (ss->criterion < 0.0 ? -1 : 1);
164 if (result < ss->limit)
168 ss->tail = (ss->tail + 1) % ss->dur;
169 if (ss->tail <= ss->head)
170 ss->head = (ss->head + 1) % ss->dur;
174 bool steadystate_deviation(unsigned long iops, unsigned long bw, struct thread_data *td)
181 struct steadystate_data *ss = &td->ss;
183 ss->cache[ss->tail] = ss->check_iops ? iops : bw;
185 if (ss->tail < ss->head || (ss->tail - ss->head == ss->dur - 1))
187 if (ss->sum_y == 0) /* first time through */
189 for(i = 0; i < ss->dur; i++)
190 ss->sum_y += ss->cache[i];
191 } else { /* easy to update the sum */
192 ss->sum_y -= ss->oldest_y;
193 ss->sum_y += ss->cache[ss->tail];
196 ss->oldest_y = ss->cache[ss->head];
197 mean = (double) ss->sum_y / ss->dur;
200 for (i = 0; i < ss->dur; i++)
202 diff = (double) ss->cache[i] - mean;
203 deviation = max(deviation, diff * (diff < 0.0 ? -1 : 1));
206 ss->criterion = ss->pct ? deviation / mean * 100.0 : deviation;
208 dprint(FD_STEADYSTATE, "sum_y: %llu, mean: %f, max diff: %f, objective: %f, limit: %f\n", ss->sum_y, mean, deviation, ss->criterion, ss->limit);
210 if (ss->criterion < ss->limit)
214 ss->tail = (ss->tail + 1) % ss->dur;
215 if (ss->tail <= ss->head)
216 ss->head = (ss->head + 1) % ss->dur;