1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
|
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html><head><meta http-equiv="Content-Type" content="text/html;charset=iso-8859-1">
<title>Qwt User's Guide: qwt_scldiv.cpp Source File</title>
<link href="doxygen.css" rel="stylesheet" type="text/css">
</head><body>
<!-- Generated by Doxygen 1.3.8 -->
<div class="qindex"><a class="qindex" href="index.html">Main Page</a> | <a class="qindex" href="hierarchy.html">Class Hierarchy</a> | <a class="qindex" href="classes.html">Alphabetical List</a> | <a class="qindex" href="annotated.html">Class List</a> | <a class="qindex" href="files.html">File List</a> | <a class="qindex" href="functions.html">Class Members</a> | <a class="qindex" href="globals.html">File Members</a> | <a class="qindex" href="pages.html">Related Pages</a></div>
<h1>qwt_scldiv.cpp</h1><pre class="fragment"><div>00001 <span class="comment">/* -*- mode: C++ ; c-file-style: "stroustrup" -*- *****************************</span>
00002 <span class="comment"> * Qwt Widget Library</span>
00003 <span class="comment"> * Copyright (C) 1997 Josef Wilgen</span>
00004 <span class="comment"> * Copyright (C) 2002 Uwe Rathmann</span>
00005 <span class="comment"> * </span>
00006 <span class="comment"> * This library is free software; you can redistribute it and/or</span>
00007 <span class="comment"> * modify it under the terms of the Qwt License, Version 1.0</span>
00008 <span class="comment"> *****************************************************************************/</span>
00009
00010 <span class="preprocessor">#include "qwt_scldiv.h"</span>
00011 <span class="preprocessor">#include "<a class="code" href="qwt__math_8h.html">qwt_math.h</a>"</span>
00012
00013 <span class="keyword">static</span> <span class="keyword">const</span> <span class="keywordtype">double</span> step_eps = 1.0e-3;
00014 <span class="keyword">static</span> <span class="keyword">const</span> <span class="keywordtype">double</span> border_eps = 1.0e-10;
00015
00016 <span class="keyword">static</span> <span class="keywordtype">bool</span> qwtLimRange(<span class="keywordtype">double</span> &val, <span class="keywordtype">double</span> v1, <span class="keywordtype">double</span> v2,
00017 <span class="keywordtype">double</span> eps_rel = 0.0, <span class="keywordtype">double</span> eps_abs = 0.0)
00018 {
00019 <span class="keywordtype">bool</span> rv = TRUE;
00020 <span class="keywordtype">double</span> vmin = qwtMin(v1, v2);
00021 <span class="keywordtype">double</span> vmax = qwtMax(v1, v2);
00022 <span class="keywordtype">double</span> delta_min = qwtMax(qwtAbs(eps_rel * vmin), qwtAbs(eps_abs));
00023 <span class="keywordtype">double</span> delta_max = qwtMax(qwtAbs(eps_rel * vmax), qwtAbs(eps_abs));
00024
00025 <span class="keywordflow">if</span> (val < vmin)
00026 {
00027 <span class="keywordflow">if</span> (val < vmin - delta_min)
00028 rv = FALSE;
00029 val = vmin;
00030 }
00031 <span class="keywordflow">else</span> <span class="keywordflow">if</span> (val > vmax)
00032 {
00033 <span class="keywordflow">if</span> (val > vmax + delta_max)
00034 rv = FALSE;
00035 val = vmax;
00036 }
00037 <span class="keywordflow">return</span> rv;
00038 }
00039
00040
<a name="l00042"></a><a class="code" href="class_qwt_scale_div.html#a0">00042</a> <a class="code" href="class_qwt_scale_div.html#a0">QwtScaleDiv::QwtScaleDiv</a>()
00043 {
00044 d_lBound = 0.0;
00045 d_hBound = 0.0;
00046 d_majStep = 0.0;
00047 d_log = FALSE;
00048 }
00049
<a name="l00086"></a><a class="code" href="class_qwt_scale_div.html#a12">00086</a> <span class="keywordtype">bool</span> <a class="code" href="class_qwt_scale_div.html#a12">QwtScaleDiv::rebuild</a>(<span class="keywordtype">double</span> x1, <span class="keywordtype">double</span> x2,
00087 <span class="keywordtype">int</span> maxMajSteps, <span class="keywordtype">int</span> maxMinSteps, <span class="keywordtype">bool</span> log, <span class="keywordtype">double</span> step, <span class="keywordtype">bool</span> ascend)
00088 {
00089 <span class="keywordtype">int</span> rv;
00090
00091 d_lBound = qwtMin(x1, x2);
00092 d_hBound = qwtMax(x1, x2);
00093 d_log = log;
00094
00095 <span class="keywordflow">if</span> (d_log)
00096 rv = buildLogDiv(maxMajSteps,maxMinSteps,step);
00097 <span class="keywordflow">else</span>
00098 rv = buildLinDiv(maxMajSteps, maxMinSteps, step);
00099
00100 <span class="keywordflow">if</span> ((!ascend) && (x2 < x1))
00101 {
00102 d_lBound = x1;
00103 d_hBound = x2;
00104 qwtTwistArray(d_majMarks.data(), d_majMarks.size());
00105 qwtTwistArray(d_minMarks.data(), d_minMarks.size());
00106 }
00107
00108 <span class="keywordflow">return</span> rv;
00109 }
00110
00129 <span class="keywordtype">bool</span> QwtScaleDiv::buildLinDiv(<span class="keywordtype">int</span> maxMajSteps, <span class="keywordtype">int</span> maxMinSteps, <span class="keywordtype">double</span> step)
00130 {
00131 <span class="keywordtype">int</span> nMaj, nMin, minSize, i0,i,k;
00132 <span class="keywordtype">double</span> val, mval;
00133 <span class="keywordtype">double</span> firstTick, lastTick;
00134 <span class="keywordtype">double</span> minStep;
00135 QwtArray<double> buffer;
00136 <span class="keywordtype">bool</span> rv = TRUE;
00137
00138 <span class="comment">// parameter range check</span>
00139 maxMajSteps = qwtMax(1, maxMajSteps);
00140 maxMinSteps = qwtMax(0, maxMinSteps);
00141 step = qwtAbs(step);
00142
00143 <span class="comment">// detach arrays</span>
00144 d_majMarks.duplicate(0,0);
00145 d_minMarks.duplicate(0,0);
00146
00147 <span class="keywordflow">if</span> (d_lBound == d_hBound) <span class="keywordflow">return</span> TRUE;
00148
00149 <span class="comment">//</span>
00150 <span class="comment">// Set up major divisions</span>
00151 <span class="comment">//</span>
00152 <span class="keywordflow">if</span> (step == 0.0)
00153 d_majStep = qwtCeil125(qwtAbs(d_hBound - d_lBound) * 0.999999
00154 / <span class="keywordtype">double</span>(maxMajSteps));
00155 <span class="keywordflow">else</span>
00156 d_majStep = step;
00157
00158 <span class="keywordflow">if</span> (d_majStep == 0.0) <span class="keywordflow">return</span> TRUE;
00159
00160 firstTick = ceil( (d_lBound - step_eps * d_majStep) / d_majStep) * d_majStep;
00161 lastTick = floor( (d_hBound + step_eps * d_majStep) / d_majStep) * d_majStep;
00162
00163 nMaj = qwtMin(10000, <span class="keywordtype">int</span>(floor ((lastTick - firstTick) / d_majStep + 0.5)) + 1);
00164
00165 <span class="keywordflow">if</span> ((rv = d_majMarks.resize(nMaj)))
00166 qwtLinSpace(d_majMarks.data(), d_majMarks.size(), firstTick, lastTick);
00167 <span class="keywordflow">else</span>
00168 <span class="keywordflow">return</span> FALSE;
00169
00170 <span class="comment">//</span>
00171 <span class="comment">// Set up minor divisions</span>
00172 <span class="comment">//</span>
00173 <span class="keywordflow">if</span> (maxMinSteps < 1) <span class="comment">// no minor divs</span>
00174 <span class="keywordflow">return</span> TRUE;
00175
00176 minStep = qwtCeil125( d_majStep / <span class="keywordtype">double</span>(maxMinSteps) );
00177
00178 <span class="keywordflow">if</span> (minStep == 0.0) <span class="keywordflow">return</span> TRUE;
00179
00180 nMin = qwtAbs(<span class="keywordtype">int</span>(floor(d_majStep / minStep + 0.5))) - 1; <span class="comment">// # minor steps per interval</span>
00181
00182 <span class="comment">// Do the minor steps fit into the interval?</span>
00183 <span class="keywordflow">if</span> ( qwtAbs(<span class="keywordtype">double</span>(nMin + 1) * minStep - d_majStep)
00184 > step_eps * d_majStep)
00185 {
00186 nMin = 1;
00187 minStep = d_majStep * 0.5;
00188 }
00189
00190 <span class="comment">// Are there minor ticks below the first major tick?</span>
00191 <span class="keywordflow">if</span> (d_majMarks[0] > d_lBound )
00192 i0 = -1;
00193 <span class="keywordflow">else</span>
00194 i0 = 0;
00195
00196 <span class="comment">// resize buffer to the maximum possible number of minor ticks</span>
00197 rv = buffer.resize(nMin * (nMaj + 1));
00198
00199 <span class="comment">// calculate minor ticks</span>
00200 <span class="keywordflow">if</span> (rv)
00201 {
00202 minSize = 0;
00203 <span class="keywordflow">for</span> (i = i0; i < (<span class="keywordtype">int</span>)d_majMarks.size(); i++)
00204 {
00205 <span class="keywordflow">if</span> (i >= 0)
00206 val = d_majMarks[i];
00207 <span class="keywordflow">else</span>
00208 val = d_majMarks[0] - d_majStep;
00209
00210 <span class="keywordflow">for</span> (k=0; k< nMin; k++)
00211 {
00212 mval = (val += minStep);
00213 <span class="keywordflow">if</span> (qwtLimRange(mval, d_lBound, d_hBound, border_eps))
00214 {
00215 buffer[minSize] = mval;
00216 minSize++;
00217 }
00218 }
00219 }
00220 d_minMarks.duplicate(buffer.data(), minSize);
00221 }
00222
00223 <span class="keywordflow">return</span> rv;
00224 }
00225
00232 <span class="keywordtype">bool</span> QwtScaleDiv::buildLogDiv(<span class="keywordtype">int</span> maxMajSteps, <span class="keywordtype">int</span> maxMinSteps, <span class="keywordtype">double</span> majStep)
00233 {
00234 <span class="keywordtype">double</span> firstTick, lastTick;
00235 <span class="keywordtype">double</span> lFirst, lLast;
00236 <span class="keywordtype">double</span> val, sval, minStep, minFactor;
00237 <span class="keywordtype">int</span> nMaj, nMin, minSize, i, k, k0, kstep, kmax, i0;
00238 <span class="keywordtype">int</span> rv = TRUE;
00239 <span class="keywordtype">double</span> width;
00240
00241 QwtArray<double> buffer;
00242
00243
00244 <span class="comment">// Parameter range check</span>
00245 maxMajSteps = qwtMax(1, qwtAbs(maxMajSteps));
00246 maxMinSteps = qwtMax(0, qwtAbs(maxMinSteps));
00247 majStep = qwtAbs(majStep);
00248
00249 <span class="comment">// boundary check</span>
00250 qwtLimRange(d_hBound, LOG_MIN, LOG_MAX);
00251 qwtLimRange(d_lBound, LOG_MIN, LOG_MAX);
00252
00253 <span class="comment">// detach arrays</span>
00254 d_majMarks.duplicate(0,0);
00255 d_minMarks.duplicate(0,0);
00256
00257 <span class="keywordflow">if</span> (d_lBound == d_hBound) <span class="keywordflow">return</span> TRUE;
00258
00259 <span class="comment">// scale width in decades</span>
00260 width = log10(d_hBound) - log10(d_lBound);
00261
00262 <span class="comment">// scale width is less than one decade -> build linear scale</span>
00263 <span class="keywordflow">if</span> (width < 1.0)
00264 {
00265 rv = buildLinDiv(maxMajSteps, maxMinSteps, 0.0);
00266 <span class="comment">// convert step width to decades</span>
00267 <span class="keywordflow">if</span> (d_majStep > 0)
00268 d_majStep = log10(d_majStep);
00269
00270 <span class="keywordflow">return</span> rv;
00271 }
00272
00273 <span class="comment">//</span>
00274 <span class="comment">// Set up major scale divisions</span>
00275 <span class="comment">//</span>
00276 <span class="keywordflow">if</span> (majStep == 0.0)
00277 d_majStep = qwtCeil125( width * 0.999999 / <span class="keywordtype">double</span>(maxMajSteps));
00278 <span class="keywordflow">else</span>
00279 d_majStep = majStep;
00280
00281 <span class="comment">// major step must be >= 1 decade</span>
00282 d_majStep = qwtMax(d_majStep, 1.0);
00283
00284
00285 lFirst = ceil((log10(d_lBound) - step_eps * d_majStep) / d_majStep) * d_majStep;
00286 lLast = floor((log10(d_hBound) + step_eps * d_majStep) / d_majStep) * d_majStep;
00287
00288 firstTick = pow(10.0, lFirst);
00289 lastTick = pow(10.0, lLast);
00290
00291 nMaj = qwtMin(10000, <span class="keywordtype">int</span>(floor (qwtAbs(lLast - lFirst) / d_majStep + 0.5)) + 1);
00292
00293 <span class="keywordflow">if</span> (d_majMarks.resize(nMaj))
00294 qwtLogSpace(d_majMarks.data(), d_majMarks.size(), firstTick, lastTick);
00295 <span class="keywordflow">else</span>
00296 <span class="keywordflow">return</span> FALSE;
00297
00298
00299 <span class="comment">//</span>
00300 <span class="comment">// Set up minor scale divisions</span>
00301 <span class="comment">//</span>
00302
00303 <span class="keywordflow">if</span> ((d_majMarks.size() < 1) || (maxMinSteps < 1)) <span class="keywordflow">return</span> TRUE; <span class="comment">// no minor marks</span>
00304
00305 <span class="keywordflow">if</span> (d_majStep < 1.1) <span class="comment">// major step width is one decade</span>
00306 {
00307 <span class="keywordflow">if</span> (maxMinSteps >= 8)
00308 {
00309 k0 = 2;
00310 kmax = 9;
00311 kstep = 1;
00312 minSize = (d_majMarks.size() + 1) * 8;
00313 }
00314 <span class="keywordflow">else</span> <span class="keywordflow">if</span> (maxMinSteps >= 4)
00315 {
00316 k0 = 2;
00317 kmax = 8;
00318 kstep = 2;
00319 minSize = (d_majMarks.size() + 1) * 4;
00320 }
00321 <span class="keywordflow">else</span> <span class="keywordflow">if</span> (maxMinSteps >= 2)
00322 {
00323 k0 = 2;
00324 kmax = 5;
00325 kstep = 3;
00326 minSize = (d_majMarks.size() + 1) * 2;
00327 }
00328 <span class="keywordflow">else</span>
00329 {
00330 k0 = 5;
00331 kmax = 5;
00332 kstep = 1;
00333 minSize = (d_majMarks.size() + 1);
00334 }
00335
00336 <span class="comment">// resize buffer to the max. possible number of minor marks</span>
00337 buffer.resize(minSize);
00338
00339 <span class="comment">// Are there minor ticks below the first major tick?</span>
00340 <span class="keywordflow">if</span> ( d_lBound < firstTick )
00341 i0 = -1;
00342 <span class="keywordflow">else</span>
00343 i0 = 0;
00344
00345 minSize = 0;
00346 <span class="keywordflow">for</span> (i = i0; i< (<span class="keywordtype">int</span>)d_majMarks.size(); i++)
00347 {
00348 <span class="keywordflow">if</span> (i >= 0)
00349 val = d_majMarks[i];
00350 <span class="keywordflow">else</span>
00351 val = d_majMarks[0] / pow(10.0, d_majStep);
00352
00353 <span class="keywordflow">for</span> (k=k0; k<= kmax; k+=kstep)
00354 {
00355 sval = val * double(k);
00356 <span class="keywordflow">if</span> (qwtLimRange(sval, d_lBound, d_hBound, border_eps))
00357 {
00358 buffer[minSize] = sval;
00359 minSize++;
00360 }
00361 }
00362 }
00363
00364 <span class="comment">// copy values into the minMarks array</span>
00365 d_minMarks.duplicate(buffer.data(), minSize);
00366
00367 }
00368 <span class="keywordflow">else</span> <span class="comment">// major step > one decade</span>
00369 {
00370
00371 <span class="comment">// substep width in decades, at least one decade</span>
00372 minStep = qwtCeil125( (d_majStep - step_eps * (d_majStep / <span class="keywordtype">double</span>(maxMinSteps)))
00373 / <span class="keywordtype">double</span>(maxMinSteps) );
00374 minStep = qwtMax(1.0, minStep);
00375
00376 <span class="comment">// # subticks per interval</span>
00377 nMin = int(floor (d_majStep / minStep + 0.5)) - 1;
00378
00379 <span class="comment">// Do the minor steps fit into the interval?</span>
00380 <span class="keywordflow">if</span> ( qwtAbs( <span class="keywordtype">double</span>(nMin + 1) * minStep - d_majStep) > step_eps * d_majStep)
00381 nMin = 0;
00382
00383 <span class="keywordflow">if</span> (nMin < 1) <span class="keywordflow">return</span> TRUE; <span class="comment">// no subticks</span>
00384
00385 <span class="comment">// resize buffer to max. possible number of subticks</span>
00386 buffer.resize((d_majMarks.size() + 1) * nMin );
00387
00388 <span class="comment">// substep factor = 10^substeps</span>
00389 minFactor = qwtMax(pow(10.0, minStep), 10.0);
00390
00391 <span class="comment">// Are there minor ticks below the first major tick?</span>
00392 <span class="keywordflow">if</span> ( d_lBound < firstTick )
00393 i0 = -1;
00394 <span class="keywordflow">else</span>
00395 i0 = 0;
00396
00397 minSize = 0;
00398 <span class="keywordflow">for</span> (i = i0; i< (<span class="keywordtype">int</span>)d_majMarks.size(); i++)
00399 {
00400 <span class="keywordflow">if</span> (i >= 0)
00401 val = d_majMarks[i];
00402 <span class="keywordflow">else</span>
00403 val = firstTick / pow(10.0, d_majStep);
00404
00405 <span class="keywordflow">for</span> (k=0; k< nMin; k++)
00406 {
00407 sval = (val *= minFactor);
00408 <span class="keywordflow">if</span> (qwtLimRange(sval, d_lBound, d_hBound, border_eps))
00409 {
00410 buffer[minSize] = sval;
00411 minSize++;
00412 }
00413 }
00414 }
00415 d_minMarks.duplicate(buffer.data(), minSize);
00416 }
00417
00418 <span class="keywordflow">return</span> rv;
00419 }
00420
<a name="l00425"></a><a class="code" href="class_qwt_scale_div.html#a1">00425</a> <span class="keywordtype">int</span> <a class="code" href="class_qwt_scale_div.html#a1">QwtScaleDiv::operator==</a>(<span class="keyword">const</span> <a class="code" href="class_qwt_scale_div.html">QwtScaleDiv</a> &s)<span class="keyword"> const</span>
00426 <span class="keyword"></span>{
00427 <span class="keywordflow">if</span> (d_lBound != s.<a class="code" href="class_qwt_scale_div.html#r0">d_lBound</a>) <span class="keywordflow">return</span> 0;
00428 <span class="keywordflow">if</span> (d_hBound != s.<a class="code" href="class_qwt_scale_div.html#r1">d_hBound</a>) <span class="keywordflow">return</span> 0;
00429 <span class="keywordflow">if</span> (d_log != s.<a class="code" href="class_qwt_scale_div.html#r3">d_log</a>) <span class="keywordflow">return</span> 0;
00430 <span class="keywordflow">if</span> (d_majStep != s.<a class="code" href="class_qwt_scale_div.html#r2">d_majStep</a>) <span class="keywordflow">return</span> 0;
00431 <span class="keywordflow">if</span> (d_majMarks != s.<a class="code" href="class_qwt_scale_div.html#r4">d_majMarks</a>) <span class="keywordflow">return</span> 0;
00432 <span class="keywordflow">return</span> (d_minMarks == s.<a class="code" href="class_qwt_scale_div.html#r5">d_minMarks</a>);
00433 }
00434
<a name="l00439"></a><a class="code" href="class_qwt_scale_div.html#a2">00439</a> <span class="keywordtype">int</span> <a class="code" href="class_qwt_scale_div.html#a2">QwtScaleDiv::operator!=</a>(<span class="keyword">const</span> <a class="code" href="class_qwt_scale_div.html">QwtScaleDiv</a> &s)<span class="keyword"> const</span>
00440 <span class="keyword"></span>{
00441 <span class="keywordflow">return</span> (!(*<span class="keyword">this</span> == s));
00442 }
00443
<a name="l00445"></a><a class="code" href="class_qwt_scale_div.html#a11">00445</a> <span class="keywordtype">void</span> <a class="code" href="class_qwt_scale_div.html#a11">QwtScaleDiv::reset</a>()
00446 {
00447 <span class="comment">// detach arrays</span>
00448 d_majMarks.duplicate(0,0);
00449 d_minMarks.duplicate(0,0);
00450
00451 d_lBound = 0.0;
00452 d_hBound = 0.0;
00453 d_majStep = 0.0;
00454 d_log = FALSE;
00455 }
</div></pre><hr size="1"><address style="align: right;"><small>Generated on Tue Nov 16 21:12:21 2004 for Qwt User's Guide by
<a href="http://www.doxygen.org/index.html">
<img src="doxygen.png" alt="doxygen" align="middle" border=0 ></a> 1.3.8 </small></address>
</body>
</html>
|
