1e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat/*
2e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat * blktrace output analysis: generate a timeline & gather statistics
3e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *
4e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat * Copyright (C) 2006 Alan D. Brunelle <Alan.Brunelle@hp.com>
5e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *
6e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *  This program is free software; you can redistribute it and/or modify
7e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *  it under the terms of the GNU General Public License as published by
8e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *  the Free Software Foundation; either version 2 of the License, or
9e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *  (at your option) any later version.
10e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *
11e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *  This program is distributed in the hope that it will be useful,
12e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *  but WITHOUT ANY WARRANTY; without even the implied warranty of
13e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *  GNU General Public License for more details.
15e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *
16e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *  You should have received a copy of the GNU General Public License
17e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *  along with this program; if not, write to the Free Software
18e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
19e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat *
20e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat */
21e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat#include <float.h>
22e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat#include "globals.h"
23e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
24e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatstruct seek_bkt {
25e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct rb_node rb_node;
26e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	long long sectors;
27e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	int nseeks;
28e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat};
29e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
30e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat/* Seeks per second */
31e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatstruct sps_bkt {
32e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	double t_start, t_last;
33e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	unsigned long nseeks;
34e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat};
35e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
36e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatstruct seeki {
37e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	FILE *rfp, *wfp, *cfp, *sps_fp;
38e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct rb_root root;
39e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct sps_bkt sps;
40e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	long long tot_seeks;
41e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	double total_sectors;
42e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	long long last_start, last_end;
43e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat};
44e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
45e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatstatic FILE *seek_open(char *str, char rw)
46e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
47e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	FILE *fp;
48e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	char *oname;
49e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
50e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (seek_name == NULL) return NULL;
51e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
52e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	oname = malloc(strlen(seek_name) + strlen(str) + 32);
53e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sprintf(oname, "%s_%s_%c.dat", seek_name, str, rw);
54e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if ((fp = my_fopen(oname, "w")) == NULL)
55e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		perror(oname);
56e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	else
57e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		add_file(fp, oname);
58e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
59e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	return fp;
60e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
61e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
62e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatstatic void __insert(struct rb_root *root, long long sectors)
63e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
64e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct seek_bkt *sbp;
65e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct rb_node *parent = NULL;
66e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct rb_node **p = &root->rb_node;
67e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
68e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	while (*p) {
69e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		parent = *p;
70e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		sbp = rb_entry(parent, struct seek_bkt, rb_node);
71e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		if (sectors < sbp->sectors)
72e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			p = &(*p)->rb_left;
73e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		else if (sectors > sbp->sectors)
74e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			p = &(*p)->rb_right;
75e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		else {
76e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			sbp->nseeks++;
77e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			return;
78e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		}
79e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	}
80e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
81e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sbp = malloc(sizeof(struct seek_bkt));
82e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sbp->nseeks = 1;
83e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sbp->sectors = sectors;
84e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
85e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	rb_link_node(&sbp->rb_node, parent, p);
86e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	rb_insert_color(&sbp->rb_node, root);
87e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
88e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
89e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatstatic void __destroy(struct rb_node *n)
90e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
91e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (n) {
92e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		struct seek_bkt *sbp = rb_entry(n, struct seek_bkt, rb_node);
93e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
94e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		__destroy(n->rb_left);
95e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		__destroy(n->rb_right);
96e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		free(sbp);
97e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	}
98e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
99e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
100e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatstatic void sps_emit(struct seeki *sip)
101e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
102e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	double tstamp, s_p_s;
103e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct sps_bkt *sps = &sip->sps;
104e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	double delta = sps->t_last - sps->t_start;
105e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
106e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if ((sps->nseeks == 1) || (delta < DBL_EPSILON)) {
107e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		s_p_s = (double)(sps->nseeks);
108e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		tstamp = sps->t_start;
109e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	} else {
110e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
111e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		s_p_s = (double)(sps->nseeks) / delta;
112e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		tstamp = sps->t_start + (delta / 2);
113e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	}
114e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
115e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	fprintf(sip->sps_fp, "%15.9lf %.2lf\n", sps->t_start, s_p_s);
116e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
117e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sps->t_start = 0;
118e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sps->nseeks = 0;
119e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
120e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
121e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatstatic void sps_add(struct seeki *sip, double t)
122e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
123e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (sip->sps_fp) {
124e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		struct sps_bkt *sps = &sip->sps;
125e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
126e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		if (sps->nseeks != 0 && ((t - sps->t_start) >= 1.0))
127e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			sps_emit(sip);
128e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
129e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		sps->t_last = t;
130e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		if (sps->nseeks == 0) {
131e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			sps->t_start = t;
132e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			sps->nseeks = 1;
133e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		} else
134e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			sps->nseeks++;
135e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	}
136e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
137e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
138e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatstatic int __median(struct rb_node *n, long long sofar, long long target,
139e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		    long long *rvp)
140e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
141e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct seek_bkt *sbp;
142e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
143e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sbp = rb_entry(n, struct seek_bkt, rb_node);
144e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if ((sofar + sbp->nseeks) >= target) {
145e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		*rvp = sbp->sectors;
146e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		return 1;
147e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	}
148e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
149e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (n->rb_left && __median(n->rb_left, sofar, target, rvp))
150e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		return 1;
151e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
152e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (n->rb_right && __median(n->rb_right, sofar, target, rvp))
153e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		return 1;
154e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
155e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	return 0;
156e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
157e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
158e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
159e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatstatic void __mode(struct rb_node *n, struct mode *mp)
160e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
161e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct seek_bkt *sbp;
162e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
163e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (n->rb_left)
164e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		__mode(n->rb_left, mp);
165e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (n->rb_right)
166e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		__mode(n->rb_right, mp);
167e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
168e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sbp = rb_entry(n, struct seek_bkt, rb_node);
169e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (mp->modes == NULL) {
170e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		mp->modes = malloc(sizeof(long long));
171e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		mp->nmds = 0;
172e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	} else if (sbp->nseeks > mp->most_seeks)
173e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		mp->nmds = 0;
174e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	else if (sbp->nseeks == mp->most_seeks)
175e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		mp->modes = realloc(mp->modes, (mp->nmds + 1) *
176e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat							sizeof(long long));
177e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	else
178e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		return;
179e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
180e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	mp->most_seeks = sbp->nseeks;
181e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	mp->modes[mp->nmds++] = sbp->sectors;
182e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
183e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
184e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatlong long seek_dist(struct seeki *sip, struct io *iop)
185e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
186e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	long long dist;
187e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	long long start = BIT_START(iop), end = BIT_END(iop);
188e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
189e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (seek_absolute)
190e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		dist = start - sip->last_end;
191e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	else {
192e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		/* Some overlap means no seek */
193e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		if (((sip->last_start <= start) && (start <= sip->last_end)) ||
194e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		    ((sip->last_start <= end) && (end <= sip->last_end)))
195e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			dist = 0;
196e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		else if (start > sip->last_end)
197e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			dist = start - sip->last_end;
198e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		else
199e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			dist = start - sip->last_start;
200e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
201e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	}
202e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
203e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sip->last_start = start;
204e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sip->last_end = end;
205e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	return dist;
206e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
207e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
208e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatvoid *seeki_alloc(char *str)
209e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
210e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct seeki *sip = malloc(sizeof(struct seeki));
211e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
212e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sip->rfp = seek_open(str, 'r');
213e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sip->wfp = seek_open(str, 'w');
214e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sip->cfp = seek_open(str, 'c');
215e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sip->tot_seeks = 0;
216e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sip->total_sectors = 0.0;
217e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sip->last_start = sip->last_end = 0;
218e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	memset(&sip->root, 0, sizeof(sip->root));
219e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
220e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (sps_name) {
221e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		char *oname;
222e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
223e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		memset(&sip->sps, 0, sizeof(sip->sps));
224e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
225e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		oname = malloc(strlen(sps_name) + strlen(str) + 32);
226e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		sprintf(oname, "%s_%s.dat", sps_name, str);
227e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		if ((sip->sps_fp = my_fopen(oname, "w")) == NULL)
228e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			perror(oname);
229e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		else
230e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			add_file(sip->sps_fp, oname);
231e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	} else
232e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		sip->sps_fp = NULL;
233e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
234e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	return sip;
235e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
236e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
237e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatvoid seeki_free(void *param)
238e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
239e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct seeki *sip = param;
240e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
241e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (sip->sps_fp && sip->sps.nseeks != 0)
242e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		sps_emit(sip);
243e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
244e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	/*
245e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	 * Associated files are cleaned up by seek_clean
246e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	 */
247e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	__destroy(sip->root.rb_node);
248e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	free(sip);
249e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
250e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
251e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatvoid seeki_add(void *handle, struct io *iop)
252e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
253e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct seeki *sip = handle;
254e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	char rw = IOP_READ(iop) ? 'r' : 'w';
255e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	long long dist = seek_dist(sip, iop);
256e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	double tstamp = BIT_TIME(iop->t.time);
257e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	FILE *fp = IOP_READ(iop) ? sip->rfp : sip->wfp;
258e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
259e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (fp)
260e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		fprintf(fp, "%15.9lf %13lld %c\n", tstamp, dist, rw);
261e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (sip->cfp)
262e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		fprintf(sip->cfp, "%15.9lf %13lld %c\n", tstamp, dist, rw);
263e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
264e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	dist = llabs(dist);
265e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sip->tot_seeks++;
266e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sip->total_sectors += dist;
267e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	__insert(&sip->root, dist);
268e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
269e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	sps_add(sip, tstamp);
270e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
271e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
272e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatlong long seeki_nseeks(void *handle)
273e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
274e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	return ((struct seeki *)handle)->tot_seeks;
275e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
276e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
277e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatdouble seeki_mean(void *handle)
278e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
279e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct seeki *sip = handle;
280e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	return sip->total_sectors / sip->tot_seeks;
281e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
282e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
283e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatlong long seeki_median(void *handle)
284e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
285e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	long long rval = 0LL;
286e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct seeki *sip = handle;
287e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
288e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (sip->root.rb_node)
289e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		(void)__median(sip->root.rb_node, 0LL, sip->tot_seeks / 2,
290e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat			       &rval);
291e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
292e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	return rval;
293e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
294e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
295e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehatint seeki_mode(void *handle, struct mode *mp)
296e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat{
297e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct seeki *sip = handle;
298e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	struct rb_root *root = &sip->root;
299e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
300e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	memset(mp, 0, sizeof(struct mode));
301e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	if (root->rb_node)
302e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat		__mode(root->rb_node, mp);
303e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat
304e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat	return mp->nmds;
305e20e1347b9914aa05e30548c15d7cd5e412cc0e2San Mehat}
306