ff1859d6622cedfb85b3fc0d13c2dc9811604134
1 /****************************************************************************
2 * RRDtool 1.3.2 Copyright by Tobi Oetiker, 1997-2008
3 ****************************************************************************
4 * rrd__graph.c produce graphs from data in rrdfiles
5 ****************************************************************************/
8 #include <sys/stat.h>
10 #ifdef WIN32
11 #include "strftime.h"
12 #endif
14 #include "rrd_tool.h"
16 /* for basename */
17 #ifdef HAVE_LIBGEN_H
18 # include <libgen.h>
19 #else
20 #include "plbasename.h"
21 #endif
23 #if defined(WIN32) && !defined(__CYGWIN__) && !defined(__CYGWIN32__)
24 #include <io.h>
25 #include <fcntl.h>
26 #endif
28 #ifdef HAVE_TIME_H
29 #include <time.h>
30 #endif
32 #ifdef HAVE_LOCALE_H
33 #include <locale.h>
34 #endif
36 #include "rrd_graph.h"
37 #include "rrd_client.h"
39 /* some constant definitions */
43 #ifndef RRD_DEFAULT_FONT
44 /* there is special code later to pick Cour.ttf when running on windows */
45 #define RRD_DEFAULT_FONT "DejaVu Sans Mono,Bitstream Vera Sans Mono,monospace,Courier"
46 #endif
48 text_prop_t text_prop[] = {
49 {8.0, RRD_DEFAULT_FONT,NULL}
50 , /* default */
51 {9.0, RRD_DEFAULT_FONT,NULL}
52 , /* title */
53 {7.0, RRD_DEFAULT_FONT,NULL}
54 , /* axis */
55 {8.0, RRD_DEFAULT_FONT,NULL}
56 , /* unit */
57 {8.0, RRD_DEFAULT_FONT,NULL} /* legend */
58 ,
59 {5.5, RRD_DEFAULT_FONT,NULL} /* watermark */
60 };
62 xlab_t xlab[] = {
63 {0, 0, TMT_SECOND, 30, TMT_MINUTE, 5, TMT_MINUTE, 5, 0, "%H:%M"}
64 ,
65 {2, 0, TMT_MINUTE, 1, TMT_MINUTE, 5, TMT_MINUTE, 5, 0, "%H:%M"}
66 ,
67 {5, 0, TMT_MINUTE, 2, TMT_MINUTE, 10, TMT_MINUTE, 10, 0, "%H:%M"}
68 ,
69 {10, 0, TMT_MINUTE, 5, TMT_MINUTE, 20, TMT_MINUTE, 20, 0, "%H:%M"}
70 ,
71 {30, 0, TMT_MINUTE, 10, TMT_HOUR, 1, TMT_HOUR, 1, 0, "%H:%M"}
72 ,
73 {60, 0, TMT_MINUTE, 30, TMT_HOUR, 2, TMT_HOUR, 2, 0, "%H:%M"}
74 ,
75 {60, 24 * 3600, TMT_MINUTE, 30, TMT_HOUR, 2, TMT_HOUR, 6, 0, "%a %H:%M"}
76 ,
77 {180, 0, TMT_HOUR, 1, TMT_HOUR, 6, TMT_HOUR, 6, 0, "%H:%M"}
78 ,
79 {180, 24 * 3600, TMT_HOUR, 1, TMT_HOUR, 6, TMT_HOUR, 12, 0, "%a %H:%M"}
80 ,
81 /*{300, 0, TMT_HOUR,3, TMT_HOUR,12, TMT_HOUR,12, 12*3600,"%a %p"}, this looks silly */
82 {600, 0, TMT_HOUR, 6, TMT_DAY, 1, TMT_DAY, 1, 24 * 3600, "%a"}
83 ,
84 {1200, 0, TMT_HOUR, 6, TMT_DAY, 1, TMT_DAY, 1, 24 * 3600, "%d"}
85 ,
86 {1800, 0, TMT_HOUR, 12, TMT_DAY, 1, TMT_DAY, 2, 24 * 3600, "%a %d"}
87 ,
88 {2400, 0, TMT_HOUR, 12, TMT_DAY, 1, TMT_DAY, 2, 24 * 3600, "%a"}
89 ,
90 {3600, 0, TMT_DAY, 1, TMT_WEEK, 1, TMT_WEEK, 1, 7 * 24 * 3600, "Week %V"}
91 ,
92 {3 * 3600, 0, TMT_WEEK, 1, TMT_MONTH, 1, TMT_WEEK, 2, 7 * 24 * 3600,
93 "Week %V"}
94 ,
95 {6 * 3600, 0, TMT_MONTH, 1, TMT_MONTH, 1, TMT_MONTH, 1, 30 * 24 * 3600,
96 "%b"}
97 ,
98 {48 * 3600, 0, TMT_MONTH, 1, TMT_MONTH, 3, TMT_MONTH, 3, 30 * 24 * 3600,
99 "%b"}
100 ,
101 {315360, 0, TMT_MONTH, 3, TMT_YEAR, 1, TMT_YEAR, 1, 365 * 24 * 3600, "%Y"}
102 ,
103 {10 * 24 * 3600, 0, TMT_YEAR, 1, TMT_YEAR, 1, TMT_YEAR, 1,
104 365 * 24 * 3600, "%y"}
105 ,
106 {-1, 0, TMT_MONTH, 0, TMT_MONTH, 0, TMT_MONTH, 0, 0, ""}
107 };
109 /* sensible y label intervals ...*/
111 ylab_t ylab[] = {
112 {0.1, {1, 2, 5, 10}
113 }
114 ,
115 {0.2, {1, 5, 10, 20}
116 }
117 ,
118 {0.5, {1, 2, 4, 10}
119 }
120 ,
121 {1.0, {1, 2, 5, 10}
122 }
123 ,
124 {2.0, {1, 5, 10, 20}
125 }
126 ,
127 {5.0, {1, 2, 4, 10}
128 }
129 ,
130 {10.0, {1, 2, 5, 10}
131 }
132 ,
133 {20.0, {1, 5, 10, 20}
134 }
135 ,
136 {50.0, {1, 2, 4, 10}
137 }
138 ,
139 {100.0, {1, 2, 5, 10}
140 }
141 ,
142 {200.0, {1, 5, 10, 20}
143 }
144 ,
145 {500.0, {1, 2, 4, 10}
146 }
147 ,
148 {0.0, {0, 0, 0, 0}
149 }
150 };
153 gfx_color_t graph_col[] = /* default colors */
154 {
155 {1.00, 1.00, 1.00, 1.00}, /* canvas */
156 {0.95, 0.95, 0.95, 1.00}, /* background */
157 {0.81, 0.81, 0.81, 1.00}, /* shade A */
158 {0.62, 0.62, 0.62, 1.00}, /* shade B */
159 {0.56, 0.56, 0.56, 0.75}, /* grid */
160 {0.87, 0.31, 0.31, 0.60}, /* major grid */
161 {0.00, 0.00, 0.00, 1.00}, /* font */
162 {0.50, 0.12, 0.12, 1.00}, /* arrow */
163 {0.12, 0.12, 0.12, 1.00}, /* axis */
164 {0.00, 0.00, 0.00, 1.00} /* frame */
165 };
168 /* #define DEBUG */
170 #ifdef DEBUG
171 # define DPRINT(x) (void)(printf x, printf("\n"))
172 #else
173 # define DPRINT(x)
174 #endif
177 /* initialize with xtr(im,0); */
178 int xtr(
179 image_desc_t *im,
180 time_t mytime)
181 {
182 static double pixie;
184 if (mytime == 0) {
185 pixie = (double) im->xsize / (double) (im->end - im->start);
186 return im->xorigin;
187 }
188 return (int) ((double) im->xorigin + pixie * (mytime - im->start));
189 }
191 /* translate data values into y coordinates */
192 double ytr(
193 image_desc_t *im,
194 double value)
195 {
196 static double pixie;
197 double yval;
199 if (isnan(value)) {
200 if (!im->logarithmic)
201 pixie = (double) im->ysize / (im->maxval - im->minval);
202 else
203 pixie =
204 (double) im->ysize / (log10(im->maxval) - log10(im->minval));
205 yval = im->yorigin;
206 } else if (!im->logarithmic) {
207 yval = im->yorigin - pixie * (value - im->minval);
208 } else {
209 if (value < im->minval) {
210 yval = im->yorigin;
211 } else {
212 yval = im->yorigin - pixie * (log10(value) - log10(im->minval));
213 }
214 }
215 return yval;
216 }
220 /* conversion function for symbolic entry names */
223 #define conv_if(VV,VVV) \
224 if (strcmp(#VV, string) == 0) return VVV ;
226 enum gf_en gf_conv(
227 char *string)
228 {
230 conv_if(PRINT, GF_PRINT);
231 conv_if(GPRINT, GF_GPRINT);
232 conv_if(COMMENT, GF_COMMENT);
233 conv_if(HRULE, GF_HRULE);
234 conv_if(VRULE, GF_VRULE);
235 conv_if(LINE, GF_LINE);
236 conv_if(AREA, GF_AREA);
237 conv_if(STACK, GF_STACK);
238 conv_if(TICK, GF_TICK);
239 conv_if(TEXTALIGN, GF_TEXTALIGN);
240 conv_if(DEF, GF_DEF);
241 conv_if(CDEF, GF_CDEF);
242 conv_if(VDEF, GF_VDEF);
243 conv_if(XPORT, GF_XPORT);
244 conv_if(SHIFT, GF_SHIFT);
246 return (enum gf_en)(-1);
247 }
249 enum gfx_if_en if_conv(
250 char *string)
251 {
253 conv_if(PNG, IF_PNG);
254 conv_if(SVG, IF_SVG);
255 conv_if(EPS, IF_EPS);
256 conv_if(PDF, IF_PDF);
258 return (enum gfx_if_en)(-1);
259 }
261 enum tmt_en tmt_conv(
262 char *string)
263 {
265 conv_if(SECOND, TMT_SECOND);
266 conv_if(MINUTE, TMT_MINUTE);
267 conv_if(HOUR, TMT_HOUR);
268 conv_if(DAY, TMT_DAY);
269 conv_if(WEEK, TMT_WEEK);
270 conv_if(MONTH, TMT_MONTH);
271 conv_if(YEAR, TMT_YEAR);
272 return (enum tmt_en)(-1);
273 }
275 enum grc_en grc_conv(
276 char *string)
277 {
279 conv_if(BACK, GRC_BACK);
280 conv_if(CANVAS, GRC_CANVAS);
281 conv_if(SHADEA, GRC_SHADEA);
282 conv_if(SHADEB, GRC_SHADEB);
283 conv_if(GRID, GRC_GRID);
284 conv_if(MGRID, GRC_MGRID);
285 conv_if(FONT, GRC_FONT);
286 conv_if(ARROW, GRC_ARROW);
287 conv_if(AXIS, GRC_AXIS);
288 conv_if(FRAME, GRC_FRAME);
290 return (enum grc_en)(-1);
291 }
293 enum text_prop_en text_prop_conv(
294 char *string)
295 {
297 conv_if(DEFAULT, TEXT_PROP_DEFAULT);
298 conv_if(TITLE, TEXT_PROP_TITLE);
299 conv_if(AXIS, TEXT_PROP_AXIS);
300 conv_if(UNIT, TEXT_PROP_UNIT);
301 conv_if(LEGEND, TEXT_PROP_LEGEND);
302 conv_if(WATERMARK, TEXT_PROP_WATERMARK);
303 return (enum text_prop_en)(-1);
304 }
307 #undef conv_if
309 int im_free(
310 image_desc_t *im)
311 {
312 unsigned long i, ii;
313 cairo_status_t status = (cairo_status_t) 0;
315 if (im == NULL)
316 return 0;
318 if (im->daemon_addr != NULL)
319 free(im->daemon_addr);
321 for (i = 0; i < (unsigned) im->gdes_c; i++) {
322 if (im->gdes[i].data_first) {
323 /* careful here, because a single pointer can occur several times */
324 free(im->gdes[i].data);
325 if (im->gdes[i].ds_namv) {
326 for (ii = 0; ii < im->gdes[i].ds_cnt; ii++)
327 free(im->gdes[i].ds_namv[ii]);
328 free(im->gdes[i].ds_namv);
329 }
330 }
331 /* free allocated memory used for dashed lines */
332 if (im->gdes[i].p_dashes != NULL)
333 free(im->gdes[i].p_dashes);
335 free(im->gdes[i].p_data);
336 free(im->gdes[i].rpnp);
337 }
338 free(im->gdes);
339 if (im->font_options)
340 cairo_font_options_destroy(im->font_options);
342 if (im->cr) {
343 status = cairo_status(im->cr);
344 cairo_destroy(im->cr);
345 }
346 if (im->rendered_image) {
347 free(im->rendered_image);
348 }
350 if (im->layout) {
351 g_object_unref (im->layout);
352 }
354 if (im->surface)
355 cairo_surface_destroy(im->surface);
357 if (status)
358 fprintf(stderr, "OOPS: Cairo has issues it can't even die: %s\n",
359 cairo_status_to_string(status));
361 return 0;
362 }
364 /* find SI magnitude symbol for the given number*/
365 void auto_scale(
366 image_desc_t *im, /* image description */
367 double *value,
368 char **symb_ptr,
369 double *magfact)
370 {
372 char *symbol[] = { "a", /* 10e-18 Atto */
373 "f", /* 10e-15 Femto */
374 "p", /* 10e-12 Pico */
375 "n", /* 10e-9 Nano */
376 "u", /* 10e-6 Micro */
377 "m", /* 10e-3 Milli */
378 " ", /* Base */
379 "k", /* 10e3 Kilo */
380 "M", /* 10e6 Mega */
381 "G", /* 10e9 Giga */
382 "T", /* 10e12 Tera */
383 "P", /* 10e15 Peta */
384 "E"
385 }; /* 10e18 Exa */
387 int symbcenter = 6;
388 int sindex;
390 if (*value == 0.0 || isnan(*value)) {
391 sindex = 0;
392 *magfact = 1.0;
393 } else {
394 sindex = floor(log(fabs(*value)) / log((double) im->base));
395 *magfact = pow((double) im->base, (double) sindex);
396 (*value) /= (*magfact);
397 }
398 if (sindex <= symbcenter && sindex >= -symbcenter) {
399 (*symb_ptr) = symbol[sindex + symbcenter];
400 } else {
401 (*symb_ptr) = "?";
402 }
403 }
406 static char si_symbol[] = {
407 'a', /* 10e-18 Atto */
408 'f', /* 10e-15 Femto */
409 'p', /* 10e-12 Pico */
410 'n', /* 10e-9 Nano */
411 'u', /* 10e-6 Micro */
412 'm', /* 10e-3 Milli */
413 ' ', /* Base */
414 'k', /* 10e3 Kilo */
415 'M', /* 10e6 Mega */
416 'G', /* 10e9 Giga */
417 'T', /* 10e12 Tera */
418 'P', /* 10e15 Peta */
419 'E', /* 10e18 Exa */
420 };
421 static const int si_symbcenter = 6;
423 /* find SI magnitude symbol for the numbers on the y-axis*/
424 void si_unit(
425 image_desc_t *im /* image description */
426 )
427 {
429 double digits, viewdigits = 0;
431 digits =
432 floor(log(max(fabs(im->minval), fabs(im->maxval))) /
433 log((double) im->base));
435 if (im->unitsexponent != 9999) {
436 /* unitsexponent = 9, 6, 3, 0, -3, -6, -9, etc */
437 viewdigits = floor((double)(im->unitsexponent / 3));
438 } else {
439 viewdigits = digits;
440 }
442 im->magfact = pow((double) im->base, digits);
444 #ifdef DEBUG
445 printf("digits %6.3f im->magfact %6.3f\n", digits, im->magfact);
446 #endif
448 im->viewfactor = im->magfact / pow((double) im->base, viewdigits);
450 if (((viewdigits + si_symbcenter) < sizeof(si_symbol)) &&
451 ((viewdigits + si_symbcenter) >= 0))
452 im->symbol = si_symbol[(int) viewdigits + si_symbcenter];
453 else
454 im->symbol = '?';
455 }
457 /* move min and max values around to become sensible */
459 void expand_range(
460 image_desc_t *im)
461 {
462 double sensiblevalues[] = { 1000.0, 900.0, 800.0, 750.0, 700.0,
463 600.0, 500.0, 400.0, 300.0, 250.0,
464 200.0, 125.0, 100.0, 90.0, 80.0,
465 75.0, 70.0, 60.0, 50.0, 40.0, 30.0,
466 25.0, 20.0, 10.0, 9.0, 8.0,
467 7.0, 6.0, 5.0, 4.0, 3.5, 3.0,
468 2.5, 2.0, 1.8, 1.5, 1.2, 1.0,
469 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1, 0.0, -1
470 };
472 double scaled_min, scaled_max;
473 double adj;
474 int i;
478 #ifdef DEBUG
479 printf("Min: %6.2f Max: %6.2f MagFactor: %6.2f\n",
480 im->minval, im->maxval, im->magfact);
481 #endif
483 if (isnan(im->ygridstep)) {
484 if (im->extra_flags & ALTAUTOSCALE) {
485 /* measure the amplitude of the function. Make sure that
486 graph boundaries are slightly higher then max/min vals
487 so we can see amplitude on the graph */
488 double delt, fact;
490 delt = im->maxval - im->minval;
491 adj = delt * 0.1;
492 fact = 2.0 * pow(10.0,
493 floor(log10
494 (max(fabs(im->minval), fabs(im->maxval)) /
495 im->magfact)) - 2);
496 if (delt < fact) {
497 adj = (fact - delt) * 0.55;
498 #ifdef DEBUG
499 printf
500 ("Min: %6.2f Max: %6.2f delt: %6.2f fact: %6.2f adj: %6.2f\n",
501 im->minval, im->maxval, delt, fact, adj);
502 #endif
503 }
504 im->minval -= adj;
505 im->maxval += adj;
506 } else if (im->extra_flags & ALTAUTOSCALE_MIN) {
507 /* measure the amplitude of the function. Make sure that
508 graph boundaries are slightly lower than min vals
509 so we can see amplitude on the graph */
510 adj = (im->maxval - im->minval) * 0.1;
511 im->minval -= adj;
512 } else if (im->extra_flags & ALTAUTOSCALE_MAX) {
513 /* measure the amplitude of the function. Make sure that
514 graph boundaries are slightly higher than max vals
515 so we can see amplitude on the graph */
516 adj = (im->maxval - im->minval) * 0.1;
517 im->maxval += adj;
518 } else {
519 scaled_min = im->minval / im->magfact;
520 scaled_max = im->maxval / im->magfact;
522 for (i = 1; sensiblevalues[i] > 0; i++) {
523 if (sensiblevalues[i - 1] >= scaled_min &&
524 sensiblevalues[i] <= scaled_min)
525 im->minval = sensiblevalues[i] * (im->magfact);
527 if (-sensiblevalues[i - 1] <= scaled_min &&
528 -sensiblevalues[i] >= scaled_min)
529 im->minval = -sensiblevalues[i - 1] * (im->magfact);
531 if (sensiblevalues[i - 1] >= scaled_max &&
532 sensiblevalues[i] <= scaled_max)
533 im->maxval = sensiblevalues[i - 1] * (im->magfact);
535 if (-sensiblevalues[i - 1] <= scaled_max &&
536 -sensiblevalues[i] >= scaled_max)
537 im->maxval = -sensiblevalues[i] * (im->magfact);
538 }
539 }
540 } else {
541 /* adjust min and max to the grid definition if there is one */
542 im->minval = (double) im->ylabfact * im->ygridstep *
543 floor(im->minval / ((double) im->ylabfact * im->ygridstep));
544 im->maxval = (double) im->ylabfact * im->ygridstep *
545 ceil(im->maxval / ((double) im->ylabfact * im->ygridstep));
546 }
548 #ifdef DEBUG
549 fprintf(stderr, "SCALED Min: %6.2f Max: %6.2f Factor: %6.2f\n",
550 im->minval, im->maxval, im->magfact);
551 #endif
552 }
555 void apply_gridfit(
556 image_desc_t *im)
557 {
558 if (isnan(im->minval) || isnan(im->maxval))
559 return;
560 ytr(im, DNAN);
561 if (im->logarithmic) {
562 double ya, yb, ypix, ypixfrac;
563 double log10_range = log10(im->maxval) - log10(im->minval);
565 ya = pow((double) 10, floor(log10(im->minval)));
566 while (ya < im->minval)
567 ya *= 10;
568 if (ya > im->maxval)
569 return; /* don't have y=10^x gridline */
570 yb = ya * 10;
571 if (yb <= im->maxval) {
572 /* we have at least 2 y=10^x gridlines.
573 Make sure distance between them in pixels
574 are an integer by expanding im->maxval */
575 double y_pixel_delta = ytr(im, ya) - ytr(im, yb);
576 double factor = y_pixel_delta / floor(y_pixel_delta);
577 double new_log10_range = factor * log10_range;
578 double new_ymax_log10 = log10(im->minval) + new_log10_range;
580 im->maxval = pow(10, new_ymax_log10);
581 ytr(im, DNAN); /* reset precalc */
582 log10_range = log10(im->maxval) - log10(im->minval);
583 }
584 /* make sure first y=10^x gridline is located on
585 integer pixel position by moving scale slightly
586 downwards (sub-pixel movement) */
587 ypix = ytr(im, ya) + im->ysize; /* add im->ysize so it always is positive */
588 ypixfrac = ypix - floor(ypix);
589 if (ypixfrac > 0 && ypixfrac < 1) {
590 double yfrac = ypixfrac / im->ysize;
592 im->minval = pow(10, log10(im->minval) - yfrac * log10_range);
593 im->maxval = pow(10, log10(im->maxval) - yfrac * log10_range);
594 ytr(im, DNAN); /* reset precalc */
595 }
596 } else {
597 /* Make sure we have an integer pixel distance between
598 each minor gridline */
599 double ypos1 = ytr(im, im->minval);
600 double ypos2 = ytr(im, im->minval + im->ygrid_scale.gridstep);
601 double y_pixel_delta = ypos1 - ypos2;
602 double factor = y_pixel_delta / floor(y_pixel_delta);
603 double new_range = factor * (im->maxval - im->minval);
604 double gridstep = im->ygrid_scale.gridstep;
605 double minor_y, minor_y_px, minor_y_px_frac;
607 if (im->maxval > 0.0)
608 im->maxval = im->minval + new_range;
609 else
610 im->minval = im->maxval - new_range;
611 ytr(im, DNAN); /* reset precalc */
612 /* make sure first minor gridline is on integer pixel y coord */
613 minor_y = gridstep * floor(im->minval / gridstep);
614 while (minor_y < im->minval)
615 minor_y += gridstep;
616 minor_y_px = ytr(im, minor_y) + im->ysize; /* ensure > 0 by adding ysize */
617 minor_y_px_frac = minor_y_px - floor(minor_y_px);
618 if (minor_y_px_frac > 0 && minor_y_px_frac < 1) {
619 double yfrac = minor_y_px_frac / im->ysize;
620 double range = im->maxval - im->minval;
622 im->minval = im->minval - yfrac * range;
623 im->maxval = im->maxval - yfrac * range;
624 ytr(im, DNAN); /* reset precalc */
625 }
626 calc_horizontal_grid(im); /* recalc with changed im->maxval */
627 }
628 }
630 /* reduce data reimplementation by Alex */
632 void reduce_data(
633 enum cf_en cf, /* which consolidation function ? */
634 unsigned long cur_step, /* step the data currently is in */
635 time_t *start, /* start, end and step as requested ... */
636 time_t *end, /* ... by the application will be ... */
637 unsigned long *step, /* ... adjusted to represent reality */
638 unsigned long *ds_cnt, /* number of data sources in file */
639 rrd_value_t **data)
640 { /* two dimensional array containing the data */
641 int i, reduce_factor = ceil((double) (*step) / (double) cur_step);
642 unsigned long col, dst_row, row_cnt, start_offset, end_offset, skiprows =
643 0;
644 rrd_value_t *srcptr, *dstptr;
646 (*step) = cur_step * reduce_factor; /* set new step size for reduced data */
647 dstptr = *data;
648 srcptr = *data;
649 row_cnt = ((*end) - (*start)) / cur_step;
651 #ifdef DEBUG
652 #define DEBUG_REDUCE
653 #endif
654 #ifdef DEBUG_REDUCE
655 printf("Reducing %lu rows with factor %i time %lu to %lu, step %lu\n",
656 row_cnt, reduce_factor, *start, *end, cur_step);
657 for (col = 0; col < row_cnt; col++) {
658 printf("time %10lu: ", *start + (col + 1) * cur_step);
659 for (i = 0; i < *ds_cnt; i++)
660 printf(" %8.2e", srcptr[*ds_cnt * col + i]);
661 printf("\n");
662 }
663 #endif
665 /* We have to combine [reduce_factor] rows of the source
666 ** into one row for the destination. Doing this we also
667 ** need to take care to combine the correct rows. First
668 ** alter the start and end time so that they are multiples
669 ** of the new step time. We cannot reduce the amount of
670 ** time so we have to move the end towards the future and
671 ** the start towards the past.
672 */
673 end_offset = (*end) % (*step);
674 start_offset = (*start) % (*step);
676 /* If there is a start offset (which cannot be more than
677 ** one destination row), skip the appropriate number of
678 ** source rows and one destination row. The appropriate
679 ** number is what we do know (start_offset/cur_step) of
680 ** the new interval (*step/cur_step aka reduce_factor).
681 */
682 #ifdef DEBUG_REDUCE
683 printf("start_offset: %lu end_offset: %lu\n", start_offset, end_offset);
684 printf("row_cnt before: %lu\n", row_cnt);
685 #endif
686 if (start_offset) {
687 (*start) = (*start) - start_offset;
688 skiprows = reduce_factor - start_offset / cur_step;
689 srcptr += skiprows * *ds_cnt;
690 for (col = 0; col < (*ds_cnt); col++)
691 *dstptr++ = DNAN;
692 row_cnt -= skiprows;
693 }
694 #ifdef DEBUG_REDUCE
695 printf("row_cnt between: %lu\n", row_cnt);
696 #endif
698 /* At the end we have some rows that are not going to be
699 ** used, the amount is end_offset/cur_step
700 */
701 if (end_offset) {
702 (*end) = (*end) - end_offset + (*step);
703 skiprows = end_offset / cur_step;
704 row_cnt -= skiprows;
705 }
706 #ifdef DEBUG_REDUCE
707 printf("row_cnt after: %lu\n", row_cnt);
708 #endif
710 /* Sanity check: row_cnt should be multiple of reduce_factor */
711 /* if this gets triggered, something is REALLY WRONG ... we die immediately */
713 if (row_cnt % reduce_factor) {
714 printf("SANITY CHECK: %lu rows cannot be reduced by %i \n",
715 row_cnt, reduce_factor);
716 printf("BUG in reduce_data()\n");
717 exit(1);
718 }
720 /* Now combine reduce_factor intervals at a time
721 ** into one interval for the destination.
722 */
724 for (dst_row = 0; (long int) row_cnt >= reduce_factor; dst_row++) {
725 for (col = 0; col < (*ds_cnt); col++) {
726 rrd_value_t newval = DNAN;
727 unsigned long validval = 0;
729 for (i = 0; i < reduce_factor; i++) {
730 if (isnan(srcptr[i * (*ds_cnt) + col])) {
731 continue;
732 }
733 validval++;
734 if (isnan(newval))
735 newval = srcptr[i * (*ds_cnt) + col];
736 else {
737 switch (cf) {
738 case CF_HWPREDICT:
739 case CF_MHWPREDICT:
740 case CF_DEVSEASONAL:
741 case CF_DEVPREDICT:
742 case CF_SEASONAL:
743 case CF_AVERAGE:
744 newval += srcptr[i * (*ds_cnt) + col];
745 break;
746 case CF_MINIMUM:
747 newval = min(newval, srcptr[i * (*ds_cnt) + col]);
748 break;
749 case CF_FAILURES:
750 /* an interval contains a failure if any subintervals contained a failure */
751 case CF_MAXIMUM:
752 newval = max(newval, srcptr[i * (*ds_cnt) + col]);
753 break;
754 case CF_LAST:
755 newval = srcptr[i * (*ds_cnt) + col];
756 break;
757 }
758 }
759 }
760 if (validval == 0) {
761 newval = DNAN;
762 } else {
763 switch (cf) {
764 case CF_HWPREDICT:
765 case CF_MHWPREDICT:
766 case CF_DEVSEASONAL:
767 case CF_DEVPREDICT:
768 case CF_SEASONAL:
769 case CF_AVERAGE:
770 newval /= validval;
771 break;
772 case CF_MINIMUM:
773 case CF_FAILURES:
774 case CF_MAXIMUM:
775 case CF_LAST:
776 break;
777 }
778 }
779 *dstptr++ = newval;
780 }
781 srcptr += (*ds_cnt) * reduce_factor;
782 row_cnt -= reduce_factor;
783 }
784 /* If we had to alter the endtime, we didn't have enough
785 ** source rows to fill the last row. Fill it with NaN.
786 */
787 if (end_offset)
788 for (col = 0; col < (*ds_cnt); col++)
789 *dstptr++ = DNAN;
790 #ifdef DEBUG_REDUCE
791 row_cnt = ((*end) - (*start)) / *step;
792 srcptr = *data;
793 printf("Done reducing. Currently %lu rows, time %lu to %lu, step %lu\n",
794 row_cnt, *start, *end, *step);
795 for (col = 0; col < row_cnt; col++) {
796 printf("time %10lu: ", *start + (col + 1) * (*step));
797 for (i = 0; i < *ds_cnt; i++)
798 printf(" %8.2e", srcptr[*ds_cnt * col + i]);
799 printf("\n");
800 }
801 #endif
802 }
805 /* get the data required for the graphs from the
806 relevant rrds ... */
808 int data_fetch(
809 image_desc_t *im)
810 {
811 int i, ii;
812 int skip;
814 /* pull the data from the rrd files ... */
815 for (i = 0; i < (int) im->gdes_c; i++) {
816 /* only GF_DEF elements fetch data */
817 if (im->gdes[i].gf != GF_DEF)
818 continue;
820 skip = 0;
821 /* do we have it already ? */
822 for (ii = 0; ii < i; ii++) {
823 if (im->gdes[ii].gf != GF_DEF)
824 continue;
825 if ((strcmp(im->gdes[i].rrd, im->gdes[ii].rrd) == 0)
826 && (im->gdes[i].cf == im->gdes[ii].cf)
827 && (im->gdes[i].cf_reduce == im->gdes[ii].cf_reduce)
828 && (im->gdes[i].start_orig == im->gdes[ii].start_orig)
829 && (im->gdes[i].end_orig == im->gdes[ii].end_orig)
830 && (im->gdes[i].step_orig == im->gdes[ii].step_orig)) {
831 /* OK, the data is already there.
832 ** Just copy the header portion
833 */
834 im->gdes[i].start = im->gdes[ii].start;
835 im->gdes[i].end = im->gdes[ii].end;
836 im->gdes[i].step = im->gdes[ii].step;
837 im->gdes[i].ds_cnt = im->gdes[ii].ds_cnt;
838 im->gdes[i].ds_namv = im->gdes[ii].ds_namv;
839 im->gdes[i].data = im->gdes[ii].data;
840 im->gdes[i].data_first = 0;
841 skip = 1;
842 }
843 if (skip)
844 break;
845 }
846 if (!skip) {
847 unsigned long ft_step = im->gdes[i].step; /* ft_step will record what we got from fetch */
849 /* Flush the file if
850 * - a connection to the daemon has been established
851 * - this is the first occurrence of that RRD file
852 */
853 if (rrdc_is_connected(im->daemon_addr))
854 {
855 int status;
857 status = 0;
858 for (ii = 0; ii < i; ii++)
859 {
860 if (strcmp (im->gdes[i].rrd, im->gdes[ii].rrd) == 0)
861 {
862 status = 1;
863 break;
864 }
865 }
867 if (status == 0)
868 {
869 status = rrdc_flush (im->gdes[i].rrd);
870 if (status != 0)
871 {
872 rrd_set_error ("rrdc_flush (%s) failed with status %i.",
873 im->gdes[i].rrd, status);
874 return (-1);
875 }
876 }
877 } /* if (rrdc_is_connected()) */
879 if ((rrd_fetch_fn(im->gdes[i].rrd,
880 im->gdes[i].cf,
881 &im->gdes[i].start,
882 &im->gdes[i].end,
883 &ft_step,
884 &im->gdes[i].ds_cnt,
885 &im->gdes[i].ds_namv,
886 &im->gdes[i].data)) == -1) {
887 return -1;
888 }
889 im->gdes[i].data_first = 1;
891 if (ft_step < im->gdes[i].step) {
892 reduce_data(im->gdes[i].cf_reduce,
893 ft_step,
894 &im->gdes[i].start,
895 &im->gdes[i].end,
896 &im->gdes[i].step,
897 &im->gdes[i].ds_cnt, &im->gdes[i].data);
898 } else {
899 im->gdes[i].step = ft_step;
900 }
901 }
903 /* lets see if the required data source is really there */
904 for (ii = 0; ii < (int) im->gdes[i].ds_cnt; ii++) {
905 if (strcmp(im->gdes[i].ds_namv[ii], im->gdes[i].ds_nam) == 0) {
906 im->gdes[i].ds = ii;
907 }
908 }
909 if (im->gdes[i].ds == -1) {
910 rrd_set_error("No DS called '%s' in '%s'",
911 im->gdes[i].ds_nam, im->gdes[i].rrd);
912 return -1;
913 }
915 }
916 return 0;
917 }
919 /* evaluate the expressions in the CDEF functions */
921 /*************************************************************
922 * CDEF stuff
923 *************************************************************/
925 long find_var_wrapper(
926 void *arg1,
927 char *key)
928 {
929 return find_var((image_desc_t *) arg1, key);
930 }
932 /* find gdes containing var*/
933 long find_var(
934 image_desc_t *im,
935 char *key)
936 {
937 long ii;
939 for (ii = 0; ii < im->gdes_c - 1; ii++) {
940 if ((im->gdes[ii].gf == GF_DEF
941 || im->gdes[ii].gf == GF_VDEF || im->gdes[ii].gf == GF_CDEF)
942 && (strcmp(im->gdes[ii].vname, key) == 0)) {
943 return ii;
944 }
945 }
946 return -1;
947 }
949 /* find the greatest common divisor for all the numbers
950 in the 0 terminated num array */
951 long lcd(
952 long *num)
953 {
954 long rest;
955 int i;
957 for (i = 0; num[i + 1] != 0; i++) {
958 do {
959 rest = num[i] % num[i + 1];
960 num[i] = num[i + 1];
961 num[i + 1] = rest;
962 } while (rest != 0);
963 num[i + 1] = num[i];
964 }
965 /* return i==0?num[i]:num[i-1]; */
966 return num[i];
967 }
969 /* run the rpn calculator on all the VDEF and CDEF arguments */
970 int data_calc(
971 image_desc_t *im)
972 {
974 int gdi;
975 int dataidx;
976 long *steparray, rpi;
977 int stepcnt;
978 time_t now;
979 rpnstack_t rpnstack;
981 rpnstack_init(&rpnstack);
983 for (gdi = 0; gdi < im->gdes_c; gdi++) {
984 /* Look for GF_VDEF and GF_CDEF in the same loop,
985 * so CDEFs can use VDEFs and vice versa
986 */
987 switch (im->gdes[gdi].gf) {
988 case GF_XPORT:
989 break;
990 case GF_SHIFT:{
991 graph_desc_t *vdp = &im->gdes[im->gdes[gdi].vidx];
993 /* remove current shift */
994 vdp->start -= vdp->shift;
995 vdp->end -= vdp->shift;
997 /* vdef */
998 if (im->gdes[gdi].shidx >= 0)
999 vdp->shift = im->gdes[im->gdes[gdi].shidx].vf.val;
1000 /* constant */
1001 else
1002 vdp->shift = im->gdes[gdi].shval;
1004 /* normalize shift to multiple of consolidated step */
1005 vdp->shift = (vdp->shift / (long) vdp->step) * (long) vdp->step;
1007 /* apply shift */
1008 vdp->start += vdp->shift;
1009 vdp->end += vdp->shift;
1010 break;
1011 }
1012 case GF_VDEF:
1013 /* A VDEF has no DS. This also signals other parts
1014 * of rrdtool that this is a VDEF value, not a CDEF.
1015 */
1016 im->gdes[gdi].ds_cnt = 0;
1017 if (vdef_calc(im, gdi)) {
1018 rrd_set_error("Error processing VDEF '%s'",
1019 im->gdes[gdi].vname);
1020 rpnstack_free(&rpnstack);
1021 return -1;
1022 }
1023 break;
1024 case GF_CDEF:
1025 im->gdes[gdi].ds_cnt = 1;
1026 im->gdes[gdi].ds = 0;
1027 im->gdes[gdi].data_first = 1;
1028 im->gdes[gdi].start = 0;
1029 im->gdes[gdi].end = 0;
1030 steparray = NULL;
1031 stepcnt = 0;
1032 dataidx = -1;
1034 /* Find the variables in the expression.
1035 * - VDEF variables are substituted by their values
1036 * and the opcode is changed into OP_NUMBER.
1037 * - CDEF variables are analized for their step size,
1038 * the lowest common denominator of all the step
1039 * sizes of the data sources involved is calculated
1040 * and the resulting number is the step size for the
1041 * resulting data source.
1042 */
1043 for (rpi = 0; im->gdes[gdi].rpnp[rpi].op != OP_END; rpi++) {
1044 if (im->gdes[gdi].rpnp[rpi].op == OP_VARIABLE ||
1045 im->gdes[gdi].rpnp[rpi].op == OP_PREV_OTHER) {
1046 long ptr = im->gdes[gdi].rpnp[rpi].ptr;
1048 if (im->gdes[ptr].ds_cnt == 0) { /* this is a VDEF data source */
1049 #if 0
1050 printf
1051 ("DEBUG: inside CDEF '%s' processing VDEF '%s'\n",
1052 im->gdes[gdi].vname, im->gdes[ptr].vname);
1053 printf("DEBUG: value from vdef is %f\n",
1054 im->gdes[ptr].vf.val);
1055 #endif
1056 im->gdes[gdi].rpnp[rpi].val = im->gdes[ptr].vf.val;
1057 im->gdes[gdi].rpnp[rpi].op = OP_NUMBER;
1058 } else { /* normal variables and PREF(variables) */
1060 /* add one entry to the array that keeps track of the step sizes of the
1061 * data sources going into the CDEF. */
1062 if ((steparray =
1063 (long*)rrd_realloc(steparray,
1064 (++stepcnt +
1065 1) * sizeof(*steparray))) == NULL) {
1066 rrd_set_error("realloc steparray");
1067 rpnstack_free(&rpnstack);
1068 return -1;
1069 };
1071 steparray[stepcnt - 1] = im->gdes[ptr].step;
1073 /* adjust start and end of cdef (gdi) so
1074 * that it runs from the latest start point
1075 * to the earliest endpoint of any of the
1076 * rras involved (ptr)
1077 */
1079 if (im->gdes[gdi].start < im->gdes[ptr].start)
1080 im->gdes[gdi].start = im->gdes[ptr].start;
1082 if (im->gdes[gdi].end == 0 ||
1083 im->gdes[gdi].end > im->gdes[ptr].end)
1084 im->gdes[gdi].end = im->gdes[ptr].end;
1086 /* store pointer to the first element of
1087 * the rra providing data for variable,
1088 * further save step size and data source
1089 * count of this rra
1090 */
1091 im->gdes[gdi].rpnp[rpi].data =
1092 im->gdes[ptr].data + im->gdes[ptr].ds;
1093 im->gdes[gdi].rpnp[rpi].step = im->gdes[ptr].step;
1094 im->gdes[gdi].rpnp[rpi].ds_cnt = im->gdes[ptr].ds_cnt;
1096 /* backoff the *.data ptr; this is done so
1097 * rpncalc() function doesn't have to treat
1098 * the first case differently
1099 */
1100 } /* if ds_cnt != 0 */
1101 } /* if OP_VARIABLE */
1102 } /* loop through all rpi */
1104 /* move the data pointers to the correct period */
1105 for (rpi = 0; im->gdes[gdi].rpnp[rpi].op != OP_END; rpi++) {
1106 if (im->gdes[gdi].rpnp[rpi].op == OP_VARIABLE ||
1107 im->gdes[gdi].rpnp[rpi].op == OP_PREV_OTHER) {
1108 long ptr = im->gdes[gdi].rpnp[rpi].ptr;
1109 long diff =
1110 im->gdes[gdi].start - im->gdes[ptr].start;
1112 if (diff > 0)
1113 im->gdes[gdi].rpnp[rpi].data +=
1114 (diff / im->gdes[ptr].step) *
1115 im->gdes[ptr].ds_cnt;
1116 }
1117 }
1119 if (steparray == NULL) {
1120 rrd_set_error("rpn expressions without DEF"
1121 " or CDEF variables are not supported");
1122 rpnstack_free(&rpnstack);
1123 return -1;
1124 }
1125 steparray[stepcnt] = 0;
1126 /* Now find the resulting step. All steps in all
1127 * used RRAs have to be visited
1128 */
1129 im->gdes[gdi].step = lcd(steparray);
1130 free(steparray);
1131 if ((im->gdes[gdi].data = (rrd_value_t*)malloc(((im->gdes[gdi].end -
1132 im->gdes[gdi].start)
1133 / im->gdes[gdi].step)
1134 * sizeof(double))) == NULL) {
1135 rrd_set_error("malloc im->gdes[gdi].data");
1136 rpnstack_free(&rpnstack);
1137 return -1;
1138 }
1140 /* Step through the new cdef results array and
1141 * calculate the values
1142 */
1143 for (now = im->gdes[gdi].start + im->gdes[gdi].step;
1144 now <= im->gdes[gdi].end; now += im->gdes[gdi].step) {
1145 rpnp_t *rpnp = im->gdes[gdi].rpnp;
1147 /* 3rd arg of rpn_calc is for OP_VARIABLE lookups;
1148 * in this case we are advancing by timesteps;
1149 * we use the fact that time_t is a synonym for long
1150 */
1151 if (rpn_calc(rpnp, &rpnstack, (long) now,
1152 im->gdes[gdi].data, ++dataidx) == -1) {
1153 /* rpn_calc sets the error string */
1154 rpnstack_free(&rpnstack);
1155 return -1;
1156 }
1157 } /* enumerate over time steps within a CDEF */
1158 break;
1159 default:
1160 continue;
1161 }
1162 } /* enumerate over CDEFs */
1163 rpnstack_free(&rpnstack);
1164 return 0;
1165 }
1167 /* from http://www.cygnus-software.com/papers/comparingfloats/comparingfloats.htm */
1168 /* yes we are loosing precision by doing tos with floats instead of doubles
1169 but it seems more stable this way. */
1171 static int AlmostEqual2sComplement(
1172 float A,
1173 float B,
1174 int maxUlps)
1175 {
1177 int aInt = *(int *) &A;
1178 int bInt = *(int *) &B;
1179 int intDiff;
1181 /* Make sure maxUlps is non-negative and small enough that the
1182 default NAN won't compare as equal to anything. */
1184 /* assert(maxUlps > 0 && maxUlps < 4 * 1024 * 1024); */
1186 /* Make aInt lexicographically ordered as a twos-complement int */
1188 if (aInt < 0)
1189 aInt = 0x80000000l - aInt;
1191 /* Make bInt lexicographically ordered as a twos-complement int */
1193 if (bInt < 0)
1194 bInt = 0x80000000l - bInt;
1196 intDiff = abs(aInt - bInt);
1198 if (intDiff <= maxUlps)
1199 return 1;
1201 return 0;
1202 }
1204 /* massage data so, that we get one value for each x coordinate in the graph */
1205 int data_proc(
1206 image_desc_t *im)
1207 {
1208 long i, ii;
1209 double pixstep = (double) (im->end - im->start)
1210 / (double) im->xsize; /* how much time
1211 passes in one pixel */
1212 double paintval;
1213 double minval = DNAN, maxval = DNAN;
1215 unsigned long gr_time;
1217 /* memory for the processed data */
1218 for (i = 0; i < im->gdes_c; i++) {
1219 if ((im->gdes[i].gf == GF_LINE) ||
1220 (im->gdes[i].gf == GF_AREA) || (im->gdes[i].gf == GF_TICK)) {
1221 if ((im->gdes[i].p_data = (rrd_value_t*)malloc((im->xsize + 1)
1222 * sizeof(rrd_value_t))) == NULL) {
1223 rrd_set_error("malloc data_proc");
1224 return -1;
1225 }
1226 }
1227 }
1229 for (i = 0; i < im->xsize; i++) { /* for each pixel */
1230 long vidx;
1232 gr_time = im->start + pixstep * i; /* time of the current step */
1233 paintval = 0.0;
1235 for (ii = 0; ii < im->gdes_c; ii++) {
1236 double value;
1238 switch (im->gdes[ii].gf) {
1239 case GF_LINE:
1240 case GF_AREA:
1241 case GF_TICK:
1242 if (!im->gdes[ii].stack)
1243 paintval = 0.0;
1244 value = im->gdes[ii].yrule;
1245 if (isnan(value) || (im->gdes[ii].gf == GF_TICK)) {
1246 /* The time of the data doesn't necessarily match
1247 ** the time of the graph. Beware.
1248 */
1249 vidx = im->gdes[ii].vidx;
1250 if (im->gdes[vidx].gf == GF_VDEF) {
1251 value = im->gdes[vidx].vf.val;
1252 } else
1253 if (((long int) gr_time >=
1254 (long int) im->gdes[vidx].start)
1255 && ((long int) gr_time <=
1256 (long int) im->gdes[vidx].end)) {
1257 value = im->gdes[vidx].data[(unsigned long)
1258 floor((double)
1259 (gr_time -
1260 im->gdes[vidx].
1261 start)
1262 /
1263 im->gdes[vidx].step)
1264 * im->gdes[vidx].ds_cnt +
1265 im->gdes[vidx].ds];
1266 } else {
1267 value = DNAN;
1268 }
1269 };
1271 if (!isnan(value)) {
1272 paintval += value;
1273 im->gdes[ii].p_data[i] = paintval;
1274 /* GF_TICK: the data values are not
1275 ** relevant for min and max
1276 */
1277 if (finite(paintval) && im->gdes[ii].gf != GF_TICK) {
1278 if ((isnan(minval) || paintval < minval) &&
1279 !(im->logarithmic && paintval <= 0.0))
1280 minval = paintval;
1281 if (isnan(maxval) || paintval > maxval)
1282 maxval = paintval;
1283 }
1284 } else {
1285 im->gdes[ii].p_data[i] = DNAN;
1286 }
1287 break;
1288 case GF_STACK:
1289 rrd_set_error
1290 ("STACK should already be turned into LINE or AREA here");
1291 return -1;
1292 break;
1293 default:
1294 break;
1295 }
1296 }
1297 }
1299 /* if min or max have not been asigned a value this is because
1300 there was no data in the graph ... this is not good ...
1301 lets set these to dummy values then ... */
1303 if (im->logarithmic) {
1304 if (isnan(minval) || isnan(maxval) || maxval <= 0) {
1305 minval = 0.0; /* catching this right away below */
1306 maxval = 5.1;
1307 }
1308 /* in logarithm mode, where minval is smaller or equal
1309 to 0 make the beast just way smaller than maxval */
1310 if (minval <= 0) {
1311 minval = maxval / 10e8;
1312 }
1313 } else {
1314 if (isnan(minval) || isnan(maxval)) {
1315 minval = 0.0;
1316 maxval = 1.0;
1317 }
1318 }
1320 /* adjust min and max values given by the user */
1321 /* for logscale we add something on top */
1322 if (isnan(im->minval)
1323 || ((!im->rigid) && im->minval > minval)
1324 ) {
1325 if (im->logarithmic)
1326 im->minval = minval / 2.0;
1327 else
1328 im->minval = minval;
1329 }
1330 if (isnan(im->maxval)
1331 || (!im->rigid && im->maxval < maxval)
1332 ) {
1333 if (im->logarithmic)
1334 im->maxval = maxval * 2.0;
1335 else
1336 im->maxval = maxval;
1337 }
1339 /* make sure min is smaller than max */
1340 if (im->minval > im->maxval) {
1341 if (im->minval > 0)
1342 im->minval = 0.99 * im->maxval;
1343 else
1344 im->minval = 1.01 * im->maxval;
1345 }
1347 /* make sure min and max are not equal */
1348 if (AlmostEqual2sComplement(im->minval, im->maxval, 4)) {
1349 if (im->maxval > 0)
1350 im->maxval *= 1.01;
1351 else
1352 im->maxval *= 0.99;
1354 /* make sure min and max are not both zero */
1355 if (AlmostEqual2sComplement(im->maxval, 0, 4)) {
1356 im->maxval = 1.0;
1357 }
1358 }
1359 return 0;
1360 }
1364 /* identify the point where the first gridline, label ... gets placed */
1366 time_t find_first_time(
1367 time_t start, /* what is the initial time */
1368 enum tmt_en baseint, /* what is the basic interval */
1369 long basestep /* how many if these do we jump a time */
1370 )
1371 {
1372 struct tm tm;
1374 localtime_r(&start, &tm);
1376 switch (baseint) {
1377 case TMT_SECOND:
1378 tm. tm_sec -= tm.tm_sec % basestep;
1380 break;
1381 case TMT_MINUTE:
1382 tm. tm_sec = 0;
1383 tm. tm_min -= tm.tm_min % basestep;
1385 break;
1386 case TMT_HOUR:
1387 tm. tm_sec = 0;
1388 tm. tm_min = 0;
1389 tm. tm_hour -= tm.tm_hour % basestep;
1391 break;
1392 case TMT_DAY:
1393 /* we do NOT look at the basestep for this ... */
1394 tm. tm_sec = 0;
1395 tm. tm_min = 0;
1396 tm. tm_hour = 0;
1398 break;
1399 case TMT_WEEK:
1400 /* we do NOT look at the basestep for this ... */
1401 tm. tm_sec = 0;
1402 tm. tm_min = 0;
1403 tm. tm_hour = 0;
1404 tm. tm_mday -= tm.tm_wday - 1; /* -1 because we want the monday */
1406 if (tm.tm_wday == 0)
1407 tm. tm_mday -= 7; /* we want the *previous* monday */
1409 break;
1410 case TMT_MONTH:
1411 tm. tm_sec = 0;
1412 tm. tm_min = 0;
1413 tm. tm_hour = 0;
1414 tm. tm_mday = 1;
1415 tm. tm_mon -= tm.tm_mon % basestep;
1417 break;
1419 case TMT_YEAR:
1420 tm. tm_sec = 0;
1421 tm. tm_min = 0;
1422 tm. tm_hour = 0;
1423 tm. tm_mday = 1;
1424 tm. tm_mon = 0;
1425 tm. tm_year -= (
1426 tm.tm_year + 1900) %basestep;
1428 }
1429 return mktime(&tm);
1430 }
1432 /* identify the point where the next gridline, label ... gets placed */
1433 time_t find_next_time(
1434 time_t current, /* what is the initial time */
1435 enum tmt_en baseint, /* what is the basic interval */
1436 long basestep /* how many if these do we jump a time */
1437 )
1438 {
1439 struct tm tm;
1440 time_t madetime;
1442 localtime_r(¤t, &tm);
1444 do {
1445 switch (baseint) {
1446 case TMT_SECOND:
1447 tm. tm_sec += basestep;
1449 break;
1450 case TMT_MINUTE:
1451 tm. tm_min += basestep;
1453 break;
1454 case TMT_HOUR:
1455 tm. tm_hour += basestep;
1457 break;
1458 case TMT_DAY:
1459 tm. tm_mday += basestep;
1461 break;
1462 case TMT_WEEK:
1463 tm. tm_mday += 7 * basestep;
1465 break;
1466 case TMT_MONTH:
1467 tm. tm_mon += basestep;
1469 break;
1470 case TMT_YEAR:
1471 tm. tm_year += basestep;
1472 }
1473 madetime = mktime(&tm);
1474 } while (madetime == -1); /* this is necessary to skip impssible times
1475 like the daylight saving time skips */
1476 return madetime;
1478 }
1481 /* calculate values required for PRINT and GPRINT functions */
1483 int print_calc(
1484 image_desc_t *im)
1485 {
1486 long i, ii, validsteps;
1487 double printval;
1488 struct tm tmvdef;
1489 int graphelement = 0;
1490 long vidx;
1491 int max_ii;
1492 double magfact = -1;
1493 char *si_symb = "";
1494 char *percent_s;
1495 int prline_cnt = 0;
1497 /* wow initializing tmvdef is quite a task :-) */
1498 time_t now = time(NULL);
1500 localtime_r(&now, &tmvdef);
1501 for (i = 0; i < im->gdes_c; i++) {
1502 vidx = im->gdes[i].vidx;
1503 switch (im->gdes[i].gf) {
1504 case GF_PRINT:
1505 case GF_GPRINT:
1506 /* PRINT and GPRINT can now print VDEF generated values.
1507 * There's no need to do any calculations on them as these
1508 * calculations were already made.
1509 */
1510 if (im->gdes[vidx].gf == GF_VDEF) { /* simply use vals */
1511 printval = im->gdes[vidx].vf.val;
1512 localtime_r(&im->gdes[vidx].vf.when, &tmvdef);
1513 } else { /* need to calculate max,min,avg etcetera */
1514 max_ii = ((im->gdes[vidx].end - im->gdes[vidx].start)
1515 / im->gdes[vidx].step * im->gdes[vidx].ds_cnt);
1516 printval = DNAN;
1517 validsteps = 0;
1518 for (ii = im->gdes[vidx].ds;
1519 ii < max_ii; ii += im->gdes[vidx].ds_cnt) {
1520 if (!finite(im->gdes[vidx].data[ii]))
1521 continue;
1522 if (isnan(printval)) {
1523 printval = im->gdes[vidx].data[ii];
1524 validsteps++;
1525 continue;
1526 }
1528 switch (im->gdes[i].cf) {
1529 case CF_HWPREDICT:
1530 case CF_MHWPREDICT:
1531 case CF_DEVPREDICT:
1532 case CF_DEVSEASONAL:
1533 case CF_SEASONAL:
1534 case CF_AVERAGE:
1535 validsteps++;
1536 printval += im->gdes[vidx].data[ii];
1537 break;
1538 case CF_MINIMUM:
1539 printval = min(printval, im->gdes[vidx].data[ii]);
1540 break;
1541 case CF_FAILURES:
1542 case CF_MAXIMUM:
1543 printval = max(printval, im->gdes[vidx].data[ii]);
1544 break;
1545 case CF_LAST:
1546 printval = im->gdes[vidx].data[ii];
1547 }
1548 }
1549 if (im->gdes[i].cf == CF_AVERAGE || im->gdes[i].cf > CF_LAST) {
1550 if (validsteps > 1) {
1551 printval = (printval / validsteps);
1552 }
1553 }
1554 } /* prepare printval */
1556 if ((percent_s = strstr(im->gdes[i].format, "%S")) != NULL) {
1557 /* Magfact is set to -1 upon entry to print_calc. If it
1558 * is still less than 0, then we need to run auto_scale.
1559 * Otherwise, put the value into the correct units. If
1560 * the value is 0, then do not set the symbol or magnification
1561 * so next the calculation will be performed again. */
1562 if (magfact < 0.0) {
1563 auto_scale(im, &printval, &si_symb, &magfact);
1564 if (printval == 0.0)
1565 magfact = -1.0;
1566 } else {
1567 printval /= magfact;
1568 }
1569 *(++percent_s) = 's';
1570 } else if (strstr(im->gdes[i].format, "%s") != NULL) {
1571 auto_scale(im, &printval, &si_symb, &magfact);
1572 }
1574 if (im->gdes[i].gf == GF_PRINT) {
1575 rrd_infoval_t prline;
1577 if (im->gdes[i].strftm) {
1578 prline.u_str = (char*)malloc((FMT_LEG_LEN + 2) * sizeof(char));
1579 strftime(prline.u_str,
1580 FMT_LEG_LEN, im->gdes[i].format, &tmvdef);
1581 } else if (bad_format(im->gdes[i].format)) {
1582 rrd_set_error
1583 ("bad format for PRINT in '%s'", im->gdes[i].format);
1584 return -1;
1585 } else {
1586 prline.u_str =
1587 sprintf_alloc(im->gdes[i].format, printval, si_symb);
1588 }
1589 grinfo_push(im,
1590 sprintf_alloc
1591 ("print[%ld]", prline_cnt++), RD_I_STR, prline);
1592 free(prline.u_str);
1593 } else {
1594 /* GF_GPRINT */
1596 if (im->gdes[i].strftm) {
1597 strftime(im->gdes[i].legend,
1598 FMT_LEG_LEN, im->gdes[i].format, &tmvdef);
1599 } else {
1600 if (bad_format(im->gdes[i].format)) {
1601 rrd_set_error
1602 ("bad format for GPRINT in '%s'",
1603 im->gdes[i].format);
1604 return -1;
1605 }
1606 #ifdef HAVE_SNPRINTF
1607 snprintf(im->gdes[i].legend,
1608 FMT_LEG_LEN - 2,
1609 im->gdes[i].format, printval, si_symb);
1610 #else
1611 sprintf(im->gdes[i].legend,
1612 im->gdes[i].format, printval, si_symb);
1613 #endif
1614 }
1615 graphelement = 1;
1616 }
1617 break;
1618 case GF_LINE:
1619 case GF_AREA:
1620 case GF_TICK:
1621 graphelement = 1;
1622 break;
1623 case GF_HRULE:
1624 if (isnan(im->gdes[i].yrule)) { /* we must set this here or the legend printer can not decide to print the legend */
1625 im->gdes[i].yrule = im->gdes[vidx].vf.val;
1626 };
1627 graphelement = 1;
1628 break;
1629 case GF_VRULE:
1630 if (im->gdes[i].xrule == 0) { /* again ... the legend printer needs it */
1631 im->gdes[i].xrule = im->gdes[vidx].vf.when;
1632 };
1633 graphelement = 1;
1634 break;
1635 case GF_COMMENT:
1636 case GF_TEXTALIGN:
1637 case GF_DEF:
1638 case GF_CDEF:
1639 case GF_VDEF:
1640 #ifdef WITH_PIECHART
1641 case GF_PART:
1642 #endif
1643 case GF_SHIFT:
1644 case GF_XPORT:
1645 break;
1646 case GF_STACK:
1647 rrd_set_error
1648 ("STACK should already be turned into LINE or AREA here");
1649 return -1;
1650 break;
1651 }
1652 }
1653 return graphelement;
1654 }
1658 /* place legends with color spots */
1659 int leg_place(
1660 image_desc_t *im,
1661 int calc_width)
1662 {
1663 /* graph labels */
1664 int interleg = im->text_prop[TEXT_PROP_LEGEND].size * 2.0;
1665 int border = im->text_prop[TEXT_PROP_LEGEND].size * 2.0;
1666 int fill = 0, fill_last;
1667 double legendwidth; // = im->ximg - 2 * border;
1668 int leg_c = 0;
1669 double leg_x = border;
1670 int leg_y = 0; //im->yimg;
1671 int leg_y_prev = 0; // im->yimg;
1672 int leg_cc;
1673 double glue = 0;
1674 int i, ii, mark = 0;
1675 char default_txtalign = TXA_JUSTIFIED; /*default line orientation */
1676 int *legspace;
1677 char *tab;
1678 char saved_legend[FMT_LEG_LEN + 5];
1680 if(calc_width){
1681 legendwidth = 0;
1682 }
1683 else{
1684 legendwidth = im->legendwidth - 2 * border;
1685 }
1688 if (!(im->extra_flags & NOLEGEND) && !(im->extra_flags & ONLY_GRAPH)) {
1689 if ((legspace = (int*)malloc(im->gdes_c * sizeof(int))) == NULL) {
1690 rrd_set_error("malloc for legspace");
1691 return -1;
1692 }
1694 for (i = 0; i < im->gdes_c; i++) {
1695 char prt_fctn; /*special printfunctions */
1696 if(calc_width){
1697 strcpy(saved_legend, im->gdes[i].legend);
1698 }
1700 fill_last = fill;
1701 /* hide legends for rules which are not displayed */
1702 if (im->gdes[i].gf == GF_TEXTALIGN) {
1703 default_txtalign = im->gdes[i].txtalign;
1704 }
1706 if (!(im->extra_flags & FORCE_RULES_LEGEND)) {
1707 if (im->gdes[i].gf == GF_HRULE
1708 && (im->gdes[i].yrule <
1709 im->minval || im->gdes[i].yrule > im->maxval))
1710 im->gdes[i].legend[0] = '\0';
1711 if (im->gdes[i].gf == GF_VRULE
1712 && (im->gdes[i].xrule <
1713 im->start || im->gdes[i].xrule > im->end))
1714 im->gdes[i].legend[0] = '\0';
1715 }
1717 /* turn \\t into tab */
1718 while ((tab = strstr(im->gdes[i].legend, "\\t"))) {
1719 memmove(tab, tab + 1, strlen(tab));
1720 tab[0] = (char) 9;
1721 }
1723 leg_cc = strlen(im->gdes[i].legend);
1724 /* is there a controle code at the end of the legend string ? */
1725 if (leg_cc >= 2 && im->gdes[i].legend[leg_cc - 2] == '\\') {
1726 prt_fctn = im->gdes[i].legend[leg_cc - 1];
1727 leg_cc -= 2;
1728 im->gdes[i].legend[leg_cc] = '\0';
1729 } else {
1730 prt_fctn = '\0';
1731 }
1732 /* only valid control codes */
1733 if (prt_fctn != 'l' && prt_fctn != 'n' && /* a synonym for l */
1734 prt_fctn != 'r' &&
1735 prt_fctn != 'j' &&
1736 prt_fctn != 'c' &&
1737 prt_fctn != 's' && prt_fctn != '\0' && prt_fctn != 'g') {
1738 free(legspace);
1739 rrd_set_error
1740 ("Unknown control code at the end of '%s\\%c'",
1741 im->gdes[i].legend, prt_fctn);
1742 return -1;
1743 }
1744 /* \n -> \l */
1745 if (prt_fctn == 'n') {
1746 prt_fctn = 'l';
1747 }
1749 /* remove exess space from the end of the legend for \g */
1750 while (prt_fctn == 'g' &&
1751 leg_cc > 0 && im->gdes[i].legend[leg_cc - 1] == ' ') {
1752 leg_cc--;
1753 im->gdes[i].legend[leg_cc] = '\0';
1754 }
1756 if (leg_cc != 0) {
1758 /* no interleg space if string ends in \g */
1759 legspace[i] = (prt_fctn == 'g' ? 0 : interleg);
1760 if (fill > 0) {
1761 fill += legspace[i];
1762 }
1763 fill +=
1764 gfx_get_text_width(im,
1765 fill + border,
1766 im->
1767 text_prop
1768 [TEXT_PROP_LEGEND].
1769 font_desc,
1770 im->tabwidth, im->gdes[i].legend);
1771 leg_c++;
1772 } else {
1773 legspace[i] = 0;
1774 }
1775 /* who said there was a special tag ... ? */
1776 if (prt_fctn == 'g') {
1777 prt_fctn = '\0';
1778 }
1780 if (prt_fctn == '\0') {
1781 if(calc_width && (fill > legendwidth)){
1782 legendwidth = fill;
1783 }
1784 if (i == im->gdes_c - 1 || fill > legendwidth) {
1785 /* just one legend item is left right or center */
1786 switch (default_txtalign) {
1787 case TXA_RIGHT:
1788 prt_fctn = 'r';
1789 break;
1790 case TXA_CENTER:
1791 prt_fctn = 'c';
1792 break;
1793 case TXA_JUSTIFIED:
1794 prt_fctn = 'j';
1795 break;
1796 default:
1797 prt_fctn = 'l';
1798 break;
1799 }
1800 }
1801 /* is it time to place the legends ? */
1802 if (fill > legendwidth) {
1803 if (leg_c > 1) {
1804 /* go back one */
1805 i--;
1806 fill = fill_last;
1807 leg_c--;
1808 }
1809 }
1810 if (leg_c == 1 && prt_fctn == 'j') {
1811 prt_fctn = 'l';
1812 }
1813 }
1815 if (prt_fctn != '\0') {
1816 leg_x = border;
1817 if (leg_c >= 2 && prt_fctn == 'j') {
1818 glue = (double)(legendwidth - fill) / (double)(leg_c - 1);
1819 } else {
1820 glue = 0;
1821 }
1822 if (prt_fctn == 'c')
1823 leg_x = (double)(legendwidth - fill) / 2.0;
1824 if (prt_fctn == 'r')
1825 leg_x = legendwidth - fill - border;
1826 for (ii = mark; ii <= i; ii++) {
1827 if (im->gdes[ii].legend[0] == '\0')
1828 continue; /* skip empty legends */
1829 im->gdes[ii].leg_x = leg_x;
1830 im->gdes[ii].leg_y = leg_y + border;
1831 leg_x +=
1832 (double)gfx_get_text_width(im, leg_x,
1833 im->
1834 text_prop
1835 [TEXT_PROP_LEGEND].
1836 font_desc,
1837 im->tabwidth, im->gdes[ii].legend)
1838 +(double)legspace[ii]
1839 + glue;
1840 }
1841 leg_y_prev = leg_y;
1842 if (leg_x > border || prt_fctn == 's')
1843 leg_y += im->text_prop[TEXT_PROP_LEGEND].size * 1.8;
1844 if (prt_fctn == 's')
1845 leg_y -= im->text_prop[TEXT_PROP_LEGEND].size;
1847 if(calc_width && (fill > legendwidth)){
1848 legendwidth = fill;
1849 }
1850 fill = 0;
1851 leg_c = 0;
1852 mark = ii;
1853 }
1855 if(calc_width){
1856 strcpy(im->gdes[i].legend, saved_legend);
1857 }
1858 }
1860 if(calc_width){
1861 im->legendwidth = legendwidth + 2 * border;
1862 }
1863 else{
1864 im->legendheight = leg_y + border * 0.6;
1865 }
1866 free(legspace);
1867 }
1868 return 0;
1869 }
1871 /* create a grid on the graph. it determines what to do
1872 from the values of xsize, start and end */
1874 /* the xaxis labels are determined from the number of seconds per pixel
1875 in the requested graph */
1877 int calc_horizontal_grid(
1878 image_desc_t
1879 *im)
1880 {
1881 double range;
1882 double scaledrange;
1883 int pixel, i;
1884 int gridind = 0;
1885 int decimals, fractionals;
1887 im->ygrid_scale.labfact = 2;
1888 range = im->maxval - im->minval;
1889 scaledrange = range / im->magfact;
1890 /* does the scale of this graph make it impossible to put lines
1891 on it? If so, give up. */
1892 if (isnan(scaledrange)) {
1893 return 0;
1894 }
1896 /* find grid spaceing */
1897 pixel = 1;
1898 if (isnan(im->ygridstep)) {
1899 if (im->extra_flags & ALTYGRID) {
1900 /* find the value with max number of digits. Get number of digits */
1901 decimals =
1902 ceil(log10
1903 (max(fabs(im->maxval), fabs(im->minval)) *
1904 im->viewfactor / im->magfact));
1905 if (decimals <= 0) /* everything is small. make place for zero */
1906 decimals = 1;
1907 im->ygrid_scale.gridstep =
1908 pow((double) 10,
1909 floor(log10(range * im->viewfactor / im->magfact))) /
1910 im->viewfactor * im->magfact;
1911 if (im->ygrid_scale.gridstep == 0) /* range is one -> 0.1 is reasonable scale */
1912 im->ygrid_scale.gridstep = 0.1;
1913 /* should have at least 5 lines but no more then 15 */
1914 if (range / im->ygrid_scale.gridstep < 5
1915 && im->ygrid_scale.gridstep >= 30)
1916 im->ygrid_scale.gridstep /= 10;
1917 if (range / im->ygrid_scale.gridstep > 15)
1918 im->ygrid_scale.gridstep *= 10;
1919 if (range / im->ygrid_scale.gridstep > 5) {
1920 im->ygrid_scale.labfact = 1;
1921 if (range / im->ygrid_scale.gridstep > 8
1922 || im->ygrid_scale.gridstep <
1923 1.8 * im->text_prop[TEXT_PROP_AXIS].size)
1924 im->ygrid_scale.labfact = 2;
1925 } else {
1926 im->ygrid_scale.gridstep /= 5;
1927 im->ygrid_scale.labfact = 5;
1928 }
1929 fractionals =
1930 floor(log10
1931 (im->ygrid_scale.gridstep *
1932 (double) im->ygrid_scale.labfact * im->viewfactor /
1933 im->magfact));
1934 if (fractionals < 0) { /* small amplitude. */
1935 int len = decimals - fractionals + 1;
1937 if (im->unitslength < len + 2)
1938 im->unitslength = len + 2;
1939 sprintf(im->ygrid_scale.labfmt,
1940 "%%%d.%df%s", len,
1941 -fractionals, (im->symbol != ' ' ? " %c" : ""));
1942 } else {
1943 int len = decimals + 1;
1945 if (im->unitslength < len + 2)
1946 im->unitslength = len + 2;
1947 sprintf(im->ygrid_scale.labfmt,
1948 "%%%d.0f%s", len, (im->symbol != ' ' ? " %c" : ""));
1949 }
1950 } else { /* classic rrd grid */
1951 for (i = 0; ylab[i].grid > 0; i++) {
1952 pixel = im->ysize / (scaledrange / ylab[i].grid);
1953 gridind = i;
1954 if (pixel >= 5)
1955 break;
1956 }
1958 for (i = 0; i < 4; i++) {
1959 if (pixel * ylab[gridind].lfac[i] >=
1960 1.8 * im->text_prop[TEXT_PROP_AXIS].size) {
1961 im->ygrid_scale.labfact = ylab[gridind].lfac[i];
1962 break;
1963 }
1964 }
1966 im->ygrid_scale.gridstep = ylab[gridind].grid * im->magfact;
1967 }
1968 } else {
1969 im->ygrid_scale.gridstep = im->ygridstep;
1970 im->ygrid_scale.labfact = im->ylabfact;
1971 }
1972 return 1;
1973 }
1975 int draw_horizontal_grid(
1976 image_desc_t
1977 *im)
1978 {
1979 int i;
1980 double scaledstep;
1981 char graph_label[100];
1982 int nlabels = 0;
1983 double X0 = im->xorigin;
1984 double X1 = im->xorigin + im->xsize;
1985 int sgrid = (int) (im->minval / im->ygrid_scale.gridstep - 1);
1986 int egrid = (int) (im->maxval / im->ygrid_scale.gridstep + 1);
1987 double MaxY;
1988 double second_axis_magfact = 0;
1989 char *second_axis_symb = "";
1991 scaledstep =
1992 im->ygrid_scale.gridstep /
1993 (double) im->magfact * (double) im->viewfactor;
1994 MaxY = scaledstep * (double) egrid;
1995 for (i = sgrid; i <= egrid; i++) {
1996 double Y0 = ytr(im,
1997 im->ygrid_scale.gridstep * i);
1998 double YN = ytr(im,
1999 im->ygrid_scale.gridstep * (i + 1));
2001 if (floor(Y0 + 0.5) >=
2002 im->yorigin - im->ysize && floor(Y0 + 0.5) <= im->yorigin) {
2003 /* Make sure at least 2 grid labels are shown, even if it doesn't agree
2004 with the chosen settings. Add a label if required by settings, or if
2005 there is only one label so far and the next grid line is out of bounds. */
2006 if (i % im->ygrid_scale.labfact == 0
2007 || (nlabels == 1
2008 && (YN < im->yorigin - im->ysize || YN > im->yorigin))) {
2009 if (im->symbol == ' ') {
2010 if (im->extra_flags & ALTYGRID) {
2011 sprintf(graph_label,
2012 im->ygrid_scale.labfmt,
2013 scaledstep * (double) i);
2014 } else {
2015 if (MaxY < 10) {
2016 sprintf(graph_label, "%4.1f",
2017 scaledstep * (double) i);
2018 } else {
2019 sprintf(graph_label, "%4.0f",
2020 scaledstep * (double) i);
2021 }
2022 }
2023 } else {
2024 char sisym = (i == 0 ? ' ' : im->symbol);
2026 if (im->extra_flags & ALTYGRID) {
2027 sprintf(graph_label,
2028 im->ygrid_scale.labfmt,
2029 scaledstep * (double) i, sisym);
2030 } else {
2031 if (MaxY < 10) {
2032 sprintf(graph_label, "%4.1f %c",
2033 scaledstep * (double) i, sisym);
2034 } else {
2035 sprintf(graph_label, "%4.0f %c",
2036 scaledstep * (double) i, sisym);
2037 }
2038 }
2039 }
2040 nlabels++;
2041 if (im->second_axis_scale != 0){
2042 char graph_label_right[100];
2043 double sval = im->ygrid_scale.gridstep*(double)i*im->second_axis_scale+im->second_axis_shift;
2044 if (im->second_axis_format[0] == '\0'){
2045 if (!second_axis_magfact){
2046 double dummy = im->ygrid_scale.gridstep*(double)(sgrid+egrid)/2.0*im->second_axis_scale+im->second_axis_shift;
2047 auto_scale(im,&dummy,&second_axis_symb,&second_axis_magfact);
2048 }
2049 sval /= second_axis_magfact;
2051 if(MaxY < 10) {
2052 sprintf(graph_label_right,"%5.1f %s",sval,second_axis_symb);
2053 } else {
2054 sprintf(graph_label_right,"%5.0f %s",sval,second_axis_symb);
2055 }
2056 }
2057 else {
2058 sprintf(graph_label_right,im->second_axis_format,sval);
2059 }
2060 gfx_text ( im,
2061 X1+7, Y0,
2062 im->graph_col[GRC_FONT],
2063 im->text_prop[TEXT_PROP_AXIS].font_desc,
2064 im->tabwidth,0.0, GFX_H_LEFT, GFX_V_CENTER,
2065 graph_label_right );
2066 }
2068 gfx_text(im,
2069 X0 -
2070 im->
2071 text_prop[TEXT_PROP_AXIS].
2072 size, Y0,
2073 im->graph_col[GRC_FONT],
2074 im->
2075 text_prop[TEXT_PROP_AXIS].
2076 font_desc,
2077 im->tabwidth, 0.0,
2078 GFX_H_RIGHT, GFX_V_CENTER, graph_label);
2079 gfx_line(im, X0 - 2, Y0, X0, Y0,
2080 MGRIDWIDTH, im->graph_col[GRC_MGRID]);
2081 gfx_line(im, X1, Y0, X1 + 2, Y0,
2082 MGRIDWIDTH, im->graph_col[GRC_MGRID]);
2083 gfx_dashed_line(im, X0 - 2, Y0,
2084 X1 + 2, Y0,
2085 MGRIDWIDTH,
2086 im->
2087 graph_col
2088 [GRC_MGRID],
2089 im->grid_dash_on, im->grid_dash_off);
2090 } else if (!(im->extra_flags & NOMINOR)) {
2091 gfx_line(im,
2092 X0 - 2, Y0,
2093 X0, Y0, GRIDWIDTH, im->graph_col[GRC_GRID]);
2094 gfx_line(im, X1, Y0, X1 + 2, Y0,
2095 GRIDWIDTH, im->graph_col[GRC_GRID]);
2096 gfx_dashed_line(im, X0 - 1, Y0,
2097 X1 + 1, Y0,
2098 GRIDWIDTH,
2099 im->
2100 graph_col[GRC_GRID],
2101 im->grid_dash_on, im->grid_dash_off);
2102 }
2103 }
2104 }
2105 return 1;
2106 }
2108 /* this is frexp for base 10 */
2109 double frexp10(
2110 double,
2111 double *);
2112 double frexp10(
2113 double x,
2114 double *e)
2115 {
2116 double mnt;
2117 int iexp;
2119 iexp = floor(log((double)fabs(x)) / log((double)10));
2120 mnt = x / pow(10.0, iexp);
2121 if (mnt >= 10.0) {
2122 iexp++;
2123 mnt = x / pow(10.0, iexp);
2124 }
2125 *e = iexp;
2126 return mnt;
2127 }
2130 /* logaritmic horizontal grid */
2131 int horizontal_log_grid(
2132 image_desc_t
2133 *im)
2134 {
2135 double yloglab[][10] = {
2136 {
2137 1.0, 10., 0.0, 0.0, 0.0, 0.0, 0.0,
2138 0.0, 0.0, 0.0}, {
2139 1.0, 5.0, 10., 0.0, 0.0, 0.0, 0.0,
2140 0.0, 0.0, 0.0}, {
2141 1.0, 2.0, 5.0, 7.0, 10., 0.0, 0.0,
2142 0.0, 0.0, 0.0}, {
2143 1.0, 2.0, 4.0,
2144 6.0, 8.0, 10.,
2145 0.0,
2146 0.0, 0.0, 0.0}, {
2147 1.0,
2148 2.0,
2149 3.0,
2150 4.0,
2151 5.0,
2152 6.0,
2153 7.0,
2154 8.0,
2155 9.0,
2156 10.},
2157 {
2158 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} /* last line */
2159 };
2160 int i, j, val_exp, min_exp;
2161 double nex; /* number of decades in data */
2162 double logscale; /* scale in logarithmic space */
2163 int exfrac = 1; /* decade spacing */
2164 int mid = -1; /* row in yloglab for major grid */
2165 double mspac; /* smallest major grid spacing (pixels) */
2166 int flab; /* first value in yloglab to use */
2167 double value, tmp, pre_value;
2168 double X0, X1, Y0;
2169 char graph_label[100];
2171 nex = log10(im->maxval / im->minval);
2172 logscale = im->ysize / nex;
2173 /* major spacing for data with high dynamic range */
2174 while (logscale * exfrac < 3 * im->text_prop[TEXT_PROP_LEGEND].size) {
2175 if (exfrac == 1)
2176 exfrac = 3;
2177 else
2178 exfrac += 3;
2179 }
2181 /* major spacing for less dynamic data */
2182 do {
2183 /* search best row in yloglab */
2184 mid++;
2185 for (i = 0; yloglab[mid][i + 1] < 10.0; i++);
2186 mspac = logscale * log10(10.0 / yloglab[mid][i]);
2187 }
2188 while (mspac >
2189 2 * im->text_prop[TEXT_PROP_LEGEND].size && yloglab[mid][0] > 0);
2190 if (mid)
2191 mid--;
2192 /* find first value in yloglab */
2193 for (flab = 0;
2194 yloglab[mid][flab] < 10
2195 && frexp10(im->minval, &tmp) > yloglab[mid][flab]; flab++);
2196 if (yloglab[mid][flab] == 10.0) {
2197 tmp += 1.0;
2198 flab = 0;
2199 }
2200 val_exp = tmp;
2201 if (val_exp % exfrac)
2202 val_exp += abs(-val_exp % exfrac);
2203 X0 = im->xorigin;
2204 X1 = im->xorigin + im->xsize;
2205 /* draw grid */
2206 pre_value = DNAN;
2207 while (1) {
2209 value = yloglab[mid][flab] * pow(10.0, val_exp);
2210 if (AlmostEqual2sComplement(value, pre_value, 4))
2211 break; /* it seems we are not converging */
2212 pre_value = value;
2213 Y0 = ytr(im, value);
2214 if (floor(Y0 + 0.5) <= im->yorigin - im->ysize)
2215 break;
2216 /* major grid line */
2217 gfx_line(im,
2218 X0 - 2, Y0, X0, Y0, MGRIDWIDTH, im->graph_col[GRC_MGRID]);
2219 gfx_line(im, X1, Y0, X1 + 2, Y0,
2220 MGRIDWIDTH, im->graph_col[GRC_MGRID]);
2221 gfx_dashed_line(im, X0 - 2, Y0,
2222 X1 + 2, Y0,
2223 MGRIDWIDTH,
2224 im->
2225 graph_col
2226 [GRC_MGRID], im->grid_dash_on, im->grid_dash_off);
2227 /* label */
2228 if (im->extra_flags & FORCE_UNITS_SI) {
2229 int scale;
2230 double pvalue;
2231 char symbol;
2233 scale = floor(val_exp / 3.0);
2234 if (value >= 1.0)
2235 pvalue = pow(10.0, val_exp % 3);
2236 else
2237 pvalue = pow(10.0, ((val_exp + 1) % 3) + 2);
2238 pvalue *= yloglab[mid][flab];
2239 if (((scale + si_symbcenter) < (int) sizeof(si_symbol))
2240 && ((scale + si_symbcenter) >= 0))
2241 symbol = si_symbol[scale + si_symbcenter];
2242 else
2243 symbol = '?';
2244 sprintf(graph_label, "%3.0f %c", pvalue, symbol);
2245 } else {
2246 sprintf(graph_label, "%3.0e", value);
2247 }
2248 if (im->second_axis_scale != 0){
2249 char graph_label_right[100];
2250 double sval = value*im->second_axis_scale+im->second_axis_shift;
2251 if (im->second_axis_format[0] == '\0'){
2252 if (im->extra_flags & FORCE_UNITS_SI) {
2253 double mfac = 1;
2254 char *symb = "";
2255 auto_scale(im,&sval,&symb,&mfac);
2256 sprintf(graph_label_right,"%4.0f %s", sval,symb);
2257 }
2258 else {
2259 sprintf(graph_label_right,"%3.0e", sval);
2260 }
2261 }
2262 else {
2263 sprintf(graph_label_right,im->second_axis_format,sval);
2264 }
2266 gfx_text ( im,
2267 X1+7, Y0,
2268 im->graph_col[GRC_FONT],
2269 im->text_prop[TEXT_PROP_AXIS].font_desc,
2270 im->tabwidth,0.0, GFX_H_LEFT, GFX_V_CENTER,
2271 graph_label_right );
2272 }
2274 gfx_text(im,
2275 X0 -
2276 im->
2277 text_prop[TEXT_PROP_AXIS].
2278 size, Y0,
2279 im->graph_col[GRC_FONT],
2280 im->
2281 text_prop[TEXT_PROP_AXIS].
2282 font_desc,
2283 im->tabwidth, 0.0,
2284 GFX_H_RIGHT, GFX_V_CENTER, graph_label);
2285 /* minor grid */
2286 if (mid < 4 && exfrac == 1) {
2287 /* find first and last minor line behind current major line
2288 * i is the first line and j tha last */
2289 if (flab == 0) {
2290 min_exp = val_exp - 1;
2291 for (i = 1; yloglab[mid][i] < 10.0; i++);
2292 i = yloglab[mid][i - 1] + 1;
2293 j = 10;
2294 } else {
2295 min_exp = val_exp;
2296 i = yloglab[mid][flab - 1] + 1;
2297 j = yloglab[mid][flab];
2298 }
2300 /* draw minor lines below current major line */
2301 for (; i < j; i++) {
2303 value = i * pow(10.0, min_exp);
2304 if (value < im->minval)
2305 continue;
2306 Y0 = ytr(im, value);
2307 if (floor(Y0 + 0.5) <= im->yorigin - im->ysize)
2308 break;
2309 /* draw lines */
2310 gfx_line(im,
2311 X0 - 2, Y0,
2312 X0, Y0, GRIDWIDTH, im->graph_col[GRC_GRID]);
2313 gfx_line(im, X1, Y0, X1 + 2, Y0,
2314 GRIDWIDTH, im->graph_col[GRC_GRID]);
2315 gfx_dashed_line(im, X0 - 1, Y0,
2316 X1 + 1, Y0,
2317 GRIDWIDTH,
2318 im->
2319 graph_col[GRC_GRID],
2320 im->grid_dash_on, im->grid_dash_off);
2321 }
2322 } else if (exfrac > 1) {
2323 for (i = val_exp - exfrac / 3 * 2; i < val_exp; i += exfrac / 3) {
2324 value = pow(10.0, i);
2325 if (value < im->minval)
2326 continue;
2327 Y0 = ytr(im, value);
2328 if (floor(Y0 + 0.5) <= im->yorigin - im->ysize)
2329 break;
2330 /* draw lines */
2331 gfx_line(im,
2332 X0 - 2, Y0,
2333 X0, Y0, GRIDWIDTH, im->graph_col[GRC_GRID]);
2334 gfx_line(im, X1, Y0, X1 + 2, Y0,
2335 GRIDWIDTH, im->graph_col[GRC_GRID]);
2336 gfx_dashed_line(im, X0 - 1, Y0,
2337 X1 + 1, Y0,
2338 GRIDWIDTH,
2339 im->
2340 graph_col[GRC_GRID],
2341 im->grid_dash_on, im->grid_dash_off);
2342 }
2343 }
2345 /* next decade */
2346 if (yloglab[mid][++flab] == 10.0) {
2347 flab = 0;
2348 val_exp += exfrac;
2349 }
2350 }
2352 /* draw minor lines after highest major line */
2353 if (mid < 4 && exfrac == 1) {
2354 /* find first and last minor line below current major line
2355 * i is the first line and j tha last */
2356 if (flab == 0) {
2357 min_exp = val_exp - 1;
2358 for (i = 1; yloglab[mid][i] < 10.0; i++);
2359 i = yloglab[mid][i - 1] + 1;
2360 j = 10;
2361 } else {
2362 min_exp = val_exp;
2363 i = yloglab[mid][flab - 1] + 1;
2364 j = yloglab[mid][flab];
2365 }
2367 /* draw minor lines below current major line */
2368 for (; i < j; i++) {
2370 value = i * pow(10.0, min_exp);
2371 if (value < im->minval)
2372 continue;
2373 Y0 = ytr(im, value);
2374 if (floor(Y0 + 0.5) <= im->yorigin - im->ysize)
2375 break;
2376 /* draw lines */
2377 gfx_line(im,
2378 X0 - 2, Y0, X0, Y0, GRIDWIDTH, im->graph_col[GRC_GRID]);
2379 gfx_line(im, X1, Y0, X1 + 2, Y0,
2380 GRIDWIDTH, im->graph_col[GRC_GRID]);
2381 gfx_dashed_line(im, X0 - 1, Y0,
2382 X1 + 1, Y0,
2383 GRIDWIDTH,
2384 im->
2385 graph_col[GRC_GRID],
2386 im->grid_dash_on, im->grid_dash_off);
2387 }
2388 }
2389 /* fancy minor gridlines */
2390 else if (exfrac > 1) {
2391 for (i = val_exp - exfrac / 3 * 2; i < val_exp; i += exfrac / 3) {
2392 value = pow(10.0, i);
2393 if (value < im->minval)
2394 continue;
2395 Y0 = ytr(im, value);
2396 if (floor(Y0 + 0.5) <= im->yorigin - im->ysize)
2397 break;
2398 /* draw lines */
2399 gfx_line(im,
2400 X0 - 2, Y0, X0, Y0, GRIDWIDTH, im->graph_col[GRC_GRID]);
2401 gfx_line(im, X1, Y0, X1 + 2, Y0,
2402 GRIDWIDTH, im->graph_col[GRC_GRID]);
2403 gfx_dashed_line(im, X0 - 1, Y0,
2404 X1 + 1, Y0,
2405 GRIDWIDTH,
2406 im->
2407 graph_col[GRC_GRID],
2408 im->grid_dash_on, im->grid_dash_off);
2409 }
2410 }
2412 return 1;
2413 }
2416 void vertical_grid(
2417 image_desc_t *im)
2418 {
2419 int xlab_sel; /* which sort of label and grid ? */
2420 time_t ti, tilab, timajor;
2421 long factor;
2422 char graph_label[100];
2423 double X0, Y0, Y1; /* points for filled graph and more */
2424 struct tm tm;
2426 /* the type of time grid is determined by finding
2427 the number of seconds per pixel in the graph */
2428 if (im->xlab_user.minsec == -1) {
2429 factor = (im->end - im->start) / im->xsize;
2430 xlab_sel = 0;
2431 while (xlab[xlab_sel + 1].minsec !=
2432 -1 && xlab[xlab_sel + 1].minsec <= factor) {
2433 xlab_sel++;
2434 } /* pick the last one */
2435 while (xlab[xlab_sel - 1].minsec ==
2436 xlab[xlab_sel].minsec
2437 && xlab[xlab_sel].length > (im->end - im->start)) {
2438 xlab_sel--;
2439 } /* go back to the smallest size */
2440 im->xlab_user.gridtm = xlab[xlab_sel].gridtm;
2441 im->xlab_user.gridst = xlab[xlab_sel].gridst;
2442 im->xlab_user.mgridtm = xlab[xlab_sel].mgridtm;
2443 im->xlab_user.mgridst = xlab[xlab_sel].mgridst;
2444 im->xlab_user.labtm = xlab[xlab_sel].labtm;
2445 im->xlab_user.labst = xlab[xlab_sel].labst;
2446 im->xlab_user.precis = xlab[xlab_sel].precis;
2447 im->xlab_user.stst = xlab[xlab_sel].stst;
2448 }
2450 /* y coords are the same for every line ... */
2451 Y0 = im->yorigin;
2452 Y1 = im->yorigin - im->ysize;
2453 /* paint the minor grid */
2454 if (!(im->extra_flags & NOMINOR)) {
2455 for (ti = find_first_time(im->start,
2456 im->
2457 xlab_user.
2458 gridtm,
2459 im->
2460 xlab_user.
2461 gridst),
2462 timajor =
2463 find_first_time(im->start,
2464 im->xlab_user.
2465 mgridtm,
2466 im->xlab_user.
2467 mgridst);
2468 ti < im->end;
2469 ti =
2470 find_next_time(ti, im->xlab_user.gridtm, im->xlab_user.gridst)
2471 ) {
2472 /* are we inside the graph ? */
2473 if (ti < im->start || ti > im->end)
2474 continue;
2475 while (timajor < ti) {
2476 timajor = find_next_time(timajor,
2477 im->
2478 xlab_user.
2479 mgridtm, im->xlab_user.mgridst);
2480 }
2481 if (ti == timajor)
2482 continue; /* skip as falls on major grid line */
2483 X0 = xtr(im, ti);
2484 gfx_line(im, X0, Y1 - 2, X0, Y1,
2485 GRIDWIDTH, im->graph_col[GRC_GRID]);
2486 gfx_line(im, X0, Y0, X0, Y0 + 2,
2487 GRIDWIDTH, im->graph_col[GRC_GRID]);
2488 gfx_dashed_line(im, X0, Y0 + 1, X0,
2489 Y1 - 1, GRIDWIDTH,
2490 im->
2491 graph_col[GRC_GRID],
2492 im->grid_dash_on, im->grid_dash_off);
2493 }
2494 }
2496 /* paint the major grid */
2497 for (ti = find_first_time(im->start,
2498 im->
2499 xlab_user.
2500 mgridtm,
2501 im->
2502 xlab_user.
2503 mgridst);
2504 ti < im->end;
2505 ti = find_next_time(ti, im->xlab_user.mgridtm, im->xlab_user.mgridst)
2506 ) {
2507 /* are we inside the graph ? */
2508 if (ti < im->start || ti > im->end)
2509 continue;
2510 X0 = xtr(im, ti);
2511 gfx_line(im, X0, Y1 - 2, X0, Y1,
2512 MGRIDWIDTH, im->graph_col[GRC_MGRID]);
2513 gfx_line(im, X0, Y0, X0, Y0 + 3,
2514 MGRIDWIDTH, im->graph_col[GRC_MGRID]);
2515 gfx_dashed_line(im, X0, Y0 + 3, X0,
2516 Y1 - 2, MGRIDWIDTH,
2517 im->
2518 graph_col
2519 [GRC_MGRID], im->grid_dash_on, im->grid_dash_off);
2520 }
2521 /* paint the labels below the graph */
2522 for (ti =
2523 find_first_time(im->start -
2524 im->xlab_user.
2525 precis / 2,
2526 im->xlab_user.
2527 labtm,
2528 im->xlab_user.
2529 labst);
2530 ti <=
2531 im->end -
2532 im->xlab_user.precis / 2;
2533 ti = find_next_time(ti, im->xlab_user.labtm, im->xlab_user.labst)
2534 ) {
2535 tilab = ti + im->xlab_user.precis / 2; /* correct time for the label */
2536 /* are we inside the graph ? */
2537 if (tilab < im->start || tilab > im->end)
2538 continue;
2539 #if HAVE_STRFTIME
2540 localtime_r(&tilab, &tm);
2541 strftime(graph_label, 99, im->xlab_user.stst, &tm);
2542 #else
2543 # error "your libc has no strftime I guess we'll abort the exercise here."
2544 #endif
2545 gfx_text(im,
2546 xtr(im, tilab),
2547 Y0 + 3,
2548 im->graph_col[GRC_FONT],
2549 im->
2550 text_prop[TEXT_PROP_AXIS].
2551 font_desc,
2552 im->tabwidth, 0.0,
2553 GFX_H_CENTER, GFX_V_TOP, graph_label);
2554 }
2556 }
2559 void axis_paint(
2560 image_desc_t *im)
2561 {
2562 /* draw x and y axis */
2563 /* gfx_line ( im->canvas, im->xorigin+im->xsize,im->yorigin,
2564 im->xorigin+im->xsize,im->yorigin-im->ysize,
2565 GRIDWIDTH, im->graph_col[GRC_AXIS]);
2567 gfx_line ( im->canvas, im->xorigin,im->yorigin-im->ysize,
2568 im->xorigin+im->xsize,im->yorigin-im->ysize,
2569 GRIDWIDTH, im->graph_col[GRC_AXIS]); */
2571 gfx_line(im, im->xorigin - 4,
2572 im->yorigin,
2573 im->xorigin + im->xsize +
2574 4, im->yorigin, MGRIDWIDTH, im->graph_col[GRC_AXIS]);
2575 gfx_line(im, im->xorigin,
2576 im->yorigin + 4,
2577 im->xorigin,
2578 im->yorigin - im->ysize -
2579 4, MGRIDWIDTH, im->graph_col[GRC_AXIS]);
2580 /* arrow for X and Y axis direction */
2581 gfx_new_area(im, im->xorigin + im->xsize + 2, im->yorigin - 3, im->xorigin + im->xsize + 2, im->yorigin + 3, im->xorigin + im->xsize + 7, im->yorigin, /* horyzontal */
2582 im->graph_col[GRC_ARROW]);
2583 gfx_close_path(im);
2584 gfx_new_area(im, im->xorigin - 3, im->yorigin - im->ysize - 2, im->xorigin + 3, im->yorigin - im->ysize - 2, im->xorigin, im->yorigin - im->ysize - 7, /* vertical */
2585 im->graph_col[GRC_ARROW]);
2586 gfx_close_path(im);
2587 if (im->second_axis_scale != 0){
2588 gfx_line ( im, im->xorigin+im->xsize,im->yorigin+4,
2589 im->xorigin+im->xsize,im->yorigin-im->ysize-4,
2590 MGRIDWIDTH, im->graph_col[GRC_AXIS]);
2591 gfx_new_area ( im,
2592 im->xorigin+im->xsize-2, im->yorigin-im->ysize-2,
2593 im->xorigin+im->xsize+3, im->yorigin-im->ysize-2,
2594 im->xorigin+im->xsize, im->yorigin-im->ysize-7, /* LINEOFFSET */
2595 im->graph_col[GRC_ARROW]);
2596 gfx_close_path(im);
2597 }
2599 }
2601 void grid_paint(
2602 image_desc_t *im)
2603 {
2604 long i;
2605 int res = 0;
2606 double X0, Y0; /* points for filled graph and more */
2607 struct gfx_color_t water_color;
2609 if (im->draw_3d_border > 0) {
2610 /* draw 3d border */
2611 i = im->draw_3d_border;
2612 gfx_new_area(im, 0, im->yimg,
2613 i, im->yimg - i, i, i, im->graph_col[GRC_SHADEA]);
2614 gfx_add_point(im, im->ximg - i, i);
2615 gfx_add_point(im, im->ximg, 0);
2616 gfx_add_point(im, 0, 0);
2617 gfx_close_path(im);
2618 gfx_new_area(im, i, im->yimg - i,
2619 im->ximg - i,
2620 im->yimg - i, im->ximg - i, i, im->graph_col[GRC_SHADEB]);
2621 gfx_add_point(im, im->ximg, 0);
2622 gfx_add_point(im, im->ximg, im->yimg);
2623 gfx_add_point(im, 0, im->yimg);
2624 gfx_close_path(im);
2625 }
2626 if (im->draw_x_grid == 1)
2627 vertical_grid(im);
2628 if (im->draw_y_grid == 1) {
2629 if (im->logarithmic) {
2630 res = horizontal_log_grid(im);
2631 } else {
2632 res = draw_horizontal_grid(im);
2633 }
2635 /* dont draw horizontal grid if there is no min and max val */
2636 if (!res) {
2637 char *nodata = "No Data found";
2639 gfx_text(im, im->ximg / 2,
2640 (2 * im->yorigin -
2641 im->ysize) / 2,
2642 im->graph_col[GRC_FONT],
2643 im->
2644 text_prop[TEXT_PROP_AXIS].
2645 font_desc,
2646 im->tabwidth, 0.0,
2647 GFX_H_CENTER, GFX_V_CENTER, nodata);
2648 }
2649 }
2651 /* yaxis unit description */
2652 if (im->ylegend[0] != '\0'){
2653 gfx_text(im,
2654 im->xOriginLegendY+10,
2655 im->yOriginLegendY,
2656 im->graph_col[GRC_FONT],
2657 im->
2658 text_prop[TEXT_PROP_UNIT].
2659 font_desc,
2660 im->tabwidth,
2661 RRDGRAPH_YLEGEND_ANGLE, GFX_H_CENTER, GFX_V_CENTER, im->ylegend);
2663 }
2664 if (im->second_axis_legend[0] != '\0'){
2665 gfx_text( im,
2666 im->xOriginLegendY2+10,
2667 im->yOriginLegendY2,
2668 im->graph_col[GRC_FONT],
2669 im->text_prop[TEXT_PROP_UNIT].font_desc,
2670 im->tabwidth,
2671 RRDGRAPH_YLEGEND_ANGLE,
2672 GFX_H_CENTER, GFX_V_CENTER,
2673 im->second_axis_legend);
2674 }
2676 /* graph title */
2677 gfx_text(im,
2678 im->xOriginTitle, im->yOriginTitle+6,
2679 im->graph_col[GRC_FONT],
2680 im->
2681 text_prop[TEXT_PROP_TITLE].
2682 font_desc,
2683 im->tabwidth, 0.0, GFX_H_CENTER, GFX_V_TOP, im->title);
2684 /* rrdtool 'logo' */
2685 if (!(im->extra_flags & NO_RRDTOOL_TAG)){
2686 water_color = im->graph_col[GRC_FONT];
2687 water_color.alpha = 0.3;
2688 double xpos = im->legendposition == EAST ? im->xOriginLegendY : im->ximg - 4;
2689 gfx_text(im, xpos, 5,
2690 water_color,
2691 im->
2692 text_prop[TEXT_PROP_WATERMARK].
2693 font_desc, im->tabwidth,
2694 -90, GFX_H_LEFT, GFX_V_TOP, "RRDTOOL / TOBI OETIKER");
2695 }
2696 /* graph watermark */
2697 if (im->watermark[0] != '\0') {
2698 water_color = im->graph_col[GRC_FONT];
2699 water_color.alpha = 0.3;
2700 gfx_text(im,
2701 im->ximg / 2, im->yimg - 6,
2702 water_color,
2703 im->
2704 text_prop[TEXT_PROP_WATERMARK].
2705 font_desc, im->tabwidth, 0,
2706 GFX_H_CENTER, GFX_V_BOTTOM, im->watermark);
2707 }
2709 /* graph labels */
2710 if (!(im->extra_flags & NOLEGEND) && !(im->extra_flags & ONLY_GRAPH)) {
2711 for (i = 0; i < im->gdes_c; i++) {
2712 if (im->gdes[i].legend[0] == '\0')
2713 continue;
2714 /* im->gdes[i].leg_y is the bottom of the legend */
2715 X0 = im->xOriginLegend + im->gdes[i].leg_x;
2716 Y0 = im->legenddirection == TOP_DOWN ? im->yOriginLegend + im->gdes[i].leg_y : im->yOriginLegend + im->legendheight - im->gdes[i].leg_y;
2717 gfx_text(im, X0, Y0,
2718 im->graph_col[GRC_FONT],
2719 im->
2720 text_prop
2721 [TEXT_PROP_LEGEND].font_desc,
2722 im->tabwidth, 0.0,
2723 GFX_H_LEFT, GFX_V_BOTTOM, im->gdes[i].legend);
2724 /* The legend for GRAPH items starts with "M " to have
2725 enough space for the box */
2726 if (im->gdes[i].gf != GF_PRINT &&
2727 im->gdes[i].gf != GF_GPRINT && im->gdes[i].gf != GF_COMMENT) {
2728 double boxH, boxV;
2729 double X1, Y1;
2731 boxH = gfx_get_text_width(im, 0,
2732 im->
2733 text_prop
2734 [TEXT_PROP_LEGEND].
2735 font_desc,
2736 im->tabwidth, "o") * 1.2;
2737 boxV = boxH;
2738 /* shift the box up a bit */
2739 Y0 -= boxV * 0.4;
2740 /* make sure transparent colors show up the same way as in the graph */
2741 gfx_new_area(im,
2742 X0, Y0 - boxV,
2743 X0, Y0, X0 + boxH, Y0, im->graph_col[GRC_BACK]);
2744 gfx_add_point(im, X0 + boxH, Y0 - boxV);
2745 gfx_close_path(im);
2746 gfx_new_area(im, X0, Y0 - boxV, X0,
2747 Y0, X0 + boxH, Y0, im->gdes[i].col);
2748 gfx_add_point(im, X0 + boxH, Y0 - boxV);
2749 gfx_close_path(im);
2750 cairo_save(im->cr);
2751 cairo_new_path(im->cr);
2752 cairo_set_line_width(im->cr, 1.0);
2753 X1 = X0 + boxH;
2754 Y1 = Y0 - boxV;
2755 gfx_line_fit(im, &X0, &Y0);
2756 gfx_line_fit(im, &X1, &Y1);
2757 cairo_move_to(im->cr, X0, Y0);
2758 cairo_line_to(im->cr, X1, Y0);
2759 cairo_line_to(im->cr, X1, Y1);
2760 cairo_line_to(im->cr, X0, Y1);
2761 cairo_close_path(im->cr);
2762 cairo_set_source_rgba(im->cr,
2763 im->
2764 graph_col
2765 [GRC_FRAME].
2766 red,
2767 im->
2768 graph_col
2769 [GRC_FRAME].
2770 green,
2771 im->
2772 graph_col
2773 [GRC_FRAME].
2774 blue, im->graph_col[GRC_FRAME].alpha);
2775 if (im->gdes[i].dash) {
2776 /* make box borders in legend dashed if the graph is dashed */
2777 double dashes[] = {
2778 3.0
2779 };
2780 cairo_set_dash(im->cr, dashes, 1, 0.0);
2781 }
2782 cairo_stroke(im->cr);
2783 cairo_restore(im->cr);
2784 }
2785 }
2786 }
2787 }
2790 /*****************************************************
2791 * lazy check make sure we rely need to create this graph
2792 *****************************************************/
2794 int lazy_check(
2795 image_desc_t *im)
2796 {
2797 FILE *fd = NULL;
2798 int size = 1;
2799 struct stat imgstat;
2801 if (im->lazy == 0)
2802 return 0; /* no lazy option */
2803 if (strlen(im->graphfile) == 0)
2804 return 0; /* inmemory option */
2805 if (stat(im->graphfile, &imgstat) != 0)
2806 return 0; /* can't stat */
2807 /* one pixel in the existing graph is more then what we would
2808 change here ... */
2809 if (time(NULL) - imgstat.st_mtime > (im->end - im->start) / im->xsize)
2810 return 0;
2811 if ((fd = fopen(im->graphfile, "rb")) == NULL)
2812 return 0; /* the file does not exist */
2813 switch (im->imgformat) {
2814 case IF_PNG:
2815 size = PngSize(fd, &(im->ximg), &(im->yimg));
2816 break;
2817 default:
2818 size = 1;
2819 }
2820 fclose(fd);
2821 return size;
2822 }
2825 int graph_size_location(
2826 image_desc_t
2827 *im,
2828 int elements)
2829 {
2830 /* The actual size of the image to draw is determined from
2831 ** several sources. The size given on the command line is
2832 ** the graph area but we need more as we have to draw labels
2833 ** and other things outside the graph area. If the option
2834 ** --full-size-mode is selected the size defines the total
2835 ** image size and the size available for the graph is
2836 ** calculated.
2837 */
2839 /** +---+-----------------------------------+
2840 ** | y |...............graph title.........|
2841 ** | +---+-------------------------------+
2842 ** | a | y | |
2843 ** | x | | |
2844 ** | i | a | |
2845 ** | s | x | main graph area |
2846 ** | | i | |
2847 ** | t | s | |
2848 ** | i | | |
2849 ** | t | l | |
2850 ** | l | b +-------------------------------+
2851 ** | e | l | x axis labels |
2852 ** +---+---+-------------------------------+
2853 ** |....................legends............|
2854 ** +---------------------------------------+
2855 ** | watermark |
2856 ** +---------------------------------------+
2857 */
2859 int Xvertical = 0, Xvertical2 = 0, Ytitle =
2860 0, Xylabel = 0, Xmain = 0, Ymain =
2861 0, Yxlabel = 0, Xspacing = 15, Yspacing = 15, Ywatermark = 4;
2863 // no legends and no the shall be plotted it's easy
2864 if (im->extra_flags & ONLY_GRAPH) {
2865 im->xorigin = 0;
2866 im->ximg = im->xsize;
2867 im->yimg = im->ysize;
2868 im->yorigin = im->ysize;
2869 ytr(im, DNAN);
2870 return 0;
2871 }
2873 if(im->watermark[0] != '\0') {
2874 Ywatermark = im->text_prop[TEXT_PROP_WATERMARK].size * 2;
2875 }
2877 // calculate the width of the left vertical legend
2878 if (im->ylegend[0] != '\0') {
2879 Xvertical = im->text_prop[TEXT_PROP_UNIT].size * 2;
2880 }
2882 // calculate the width of the right vertical legend
2883 if (im->second_axis_legend[0] != '\0') {
2884 Xvertical2 = im->text_prop[TEXT_PROP_UNIT].size * 2;
2885 }
2886 else{
2887 Xvertical2 = Xspacing;
2888 }
2890 if (im->title[0] != '\0') {
2891 /* The title is placed "inbetween" two text lines so it
2892 ** automatically has some vertical spacing. The horizontal
2893 ** spacing is added here, on each side.
2894 */
2895 /* if necessary, reduce the font size of the title until it fits the image width */
2896 Ytitle = im->text_prop[TEXT_PROP_TITLE].size * 2.6 + 10;
2897 }
2898 else{
2899 // we have no title; get a little clearing from the top
2900 Ytitle = 1.5 * Yspacing;
2901 }
2903 if (elements) {
2904 if (im->draw_x_grid) {
2905 // calculate the height of the horizontal labelling
2906 Yxlabel = im->text_prop[TEXT_PROP_AXIS].size * 2.5;
2907 }
2908 if (im->draw_y_grid || im->forceleftspace) {
2909 // calculate the width of the vertical labelling
2910 Xylabel =
2911 gfx_get_text_width(im, 0,
2912 im->text_prop[TEXT_PROP_AXIS].font_desc,
2913 im->tabwidth, "0") * im->unitslength;
2914 }
2915 }
2917 // add some space to the labelling
2918 Xylabel += Xspacing;
2920 /* If the legend is printed besides the graph the width has to be
2921 ** calculated first. Placing the legend north or south of the
2922 ** graph requires the width calculation first, so the legend is
2923 ** skipped for the moment.
2924 */
2925 im->legendheight = 0;
2926 im->legendwidth = 0;
2927 if (!(im->extra_flags & NOLEGEND)) {
2928 if(im->legendposition == WEST || im->legendposition == EAST){
2929 if (leg_place(im, 1) == -1){
2930 return -1;
2931 }
2932 }
2933 }
2935 if (im->extra_flags & FULL_SIZE_MODE) {
2937 /* The actual size of the image to draw has been determined by the user.
2938 ** The graph area is the space remaining after accounting for the legend,
2939 ** the watermark, the axis labels, and the title.
2940 */
2941 im->ximg = im->xsize;
2942 im->yimg = im->ysize;
2943 Xmain = im->ximg;
2944 Ymain = im->yimg;
2946 /* Now calculate the total size. Insert some spacing where
2947 desired. im->xorigin and im->yorigin need to correspond
2948 with the lower left corner of the main graph area or, if
2949 this one is not set, the imaginary box surrounding the
2950 pie chart area. */
2951 /* Initial size calculation for the main graph area */
2953 Xmain -= Xylabel;// + Xspacing;
2954 if((im->legendposition == WEST || im->legendposition == EAST) && !(im->extra_flags & NOLEGEND) ){
2955 Xmain -= im->legendwidth;// + Xspacing;
2956 }
2957 if (im->second_axis_scale != 0){
2958 Xmain -= Xylabel;
2959 }
2960 if (!(im->extra_flags & NO_RRDTOOL_TAG)){
2961 Xmain -= Xspacing;
2962 }
2964 Xmain -= Xvertical + Xvertical2;
2966 /* limit the remaining space to 0 */
2967 if(Xmain < 1){
2968 Xmain = 1;
2969 }
2970 im->xsize = Xmain;
2972 /* Putting the legend north or south, the height can now be calculated */
2973 if (!(im->extra_flags & NOLEGEND)) {
2974 if(im->legendposition == NORTH || im->legendposition == SOUTH){
2975 im->legendwidth = im->ximg;
2976 if (leg_place(im, 0) == -1){
2977 return -1;
2978 }
2979 }
2980 }
2982 if( (im->legendposition == NORTH || im->legendposition == SOUTH) && !(im->extra_flags & NOLEGEND) ){
2983 Ymain -= Yxlabel + im->legendheight;
2984 }
2985 else{
2986 Ymain -= Yxlabel;
2987 }
2989 /* reserve space for the title *or* some padding above the graph */
2990 Ymain -= Ytitle;
2992 /* reserve space for padding below the graph */
2993 if (im->extra_flags & NOLEGEND) {
2994 Ymain -= Yspacing;
2995 }
2997 if (im->watermark[0] != '\0') {
2998 Ymain -= Ywatermark;
2999 }
3000 /* limit the remaining height to 0 */
3001 if(Ymain < 1){
3002 Ymain = 1;
3003 }
3004 im->ysize = Ymain;
3005 } else { /* dimension options -width and -height refer to the dimensions of the main graph area */
3007 /* The actual size of the image to draw is determined from
3008 ** several sources. The size given on the command line is
3009 ** the graph area but we need more as we have to draw labels
3010 ** and other things outside the graph area.
3011 */
3013 if (elements) {
3014 Xmain = im->xsize; // + Xspacing;
3015 Ymain = im->ysize;
3016 }
3018 im->ximg = Xmain + Xylabel;
3019 if (!(im->extra_flags & NO_RRDTOOL_TAG)){
3020 im->ximg += Xspacing;
3021 }
3023 if( (im->legendposition == WEST || im->legendposition == EAST) && !(im->extra_flags & NOLEGEND) ){
3024 im->ximg += im->legendwidth;// + Xspacing;
3025 }
3026 if (im->second_axis_scale != 0){
3027 im->ximg += Xylabel;
3028 }
3030 im->ximg += Xvertical + Xvertical2;
3032 if (!(im->extra_flags & NOLEGEND)) {
3033 if(im->legendposition == NORTH || im->legendposition == SOUTH){
3034 im->legendwidth = im->ximg;
3035 if (leg_place(im, 0) == -1){
3036 return -1;
3037 }
3038 }
3039 }
3041 im->yimg = Ymain + Yxlabel;
3042 if( (im->legendposition == NORTH || im->legendposition == SOUTH) && !(im->extra_flags & NOLEGEND) ){
3043 im->yimg += im->legendheight;
3044 }
3046 /* reserve space for the title *or* some padding above the graph */
3047 if (Ytitle) {
3048 im->yimg += Ytitle;
3049 } else {
3050 im->yimg += 1.5 * Yspacing;
3051 }
3052 /* reserve space for padding below the graph */
3053 if (im->extra_flags & NOLEGEND) {
3054 im->yimg += Yspacing;
3055 }
3057 if (im->watermark[0] != '\0') {
3058 im->yimg += Ywatermark;
3059 }
3060 }
3063 /* In case of putting the legend in west or east position the first
3064 ** legend calculation might lead to wrong positions if some items
3065 ** are not aligned on the left hand side (e.g. centered) as the
3066 ** legendwidth wight have been increased after the item was placed.
3067 ** In this case the positions have to be recalculated.
3068 */
3069 if (!(im->extra_flags & NOLEGEND)) {
3070 if(im->legendposition == WEST || im->legendposition == EAST){
3071 if (leg_place(im, 0) == -1){
3072 return -1;
3073 }
3074 }
3075 }
3077 /* After calculating all dimensions
3078 ** it is now possible to calculate
3079 ** all offsets.
3080 */
3081 switch(im->legendposition){
3082 case NORTH:
3083 im->xOriginTitle = Xvertical + Xylabel + (im->xsize / 2);
3084 im->yOriginTitle = 0;
3086 im->xOriginLegend = 0;
3087 im->yOriginLegend = Ytitle;
3089 im->xOriginLegendY = 0;
3090 im->yOriginLegendY = Ytitle + im->legendheight + (Ymain / 2) + Yxlabel;
3092 im->xorigin = Xvertical + Xylabel;
3093 im->yorigin = Ytitle + im->legendheight + Ymain;
3095 im->xOriginLegendY2 = Xvertical + Xylabel + Xmain;
3096 if (im->second_axis_scale != 0){
3097 im->xOriginLegendY2 += Xylabel;
3098 }
3099 im->yOriginLegendY2 = Ytitle + im->legendheight + (Ymain / 2) + Yxlabel;
3101 break;
3103 case WEST:
3104 im->xOriginTitle = im->legendwidth + Xvertical + Xylabel + im->xsize / 2;
3105 im->yOriginTitle = 0;
3107 im->xOriginLegend = 0;
3108 im->yOriginLegend = Ytitle;
3110 im->xOriginLegendY = im->legendwidth;
3111 im->yOriginLegendY = Ytitle + (Ymain / 2);
3113 im->xorigin = im->legendwidth + Xvertical + Xylabel;
3114 im->yorigin = Ytitle + Ymain;
3116 im->xOriginLegendY2 = im->legendwidth + Xvertical + Xylabel + Xmain;
3117 if (im->second_axis_scale != 0){
3118 im->xOriginLegendY2 += Xylabel;
3119 }
3120 im->yOriginLegendY2 = Ytitle + (Ymain / 2);
3122 break;
3124 case SOUTH:
3125 im->xOriginTitle = Xvertical + Xylabel + im->xsize / 2;
3126 im->yOriginTitle = 0;
3128 im->xOriginLegend = 0;
3129 im->yOriginLegend = Ytitle + Ymain + Yxlabel;
3131 im->xOriginLegendY = 0;
3132 im->yOriginLegendY = Ytitle + (Ymain / 2);
3134 im->xorigin = Xvertical + Xylabel;
3135 im->yorigin = Ytitle + Ymain;
3137 im->xOriginLegendY2 = Xvertical + Xylabel + Xmain;
3138 if (im->second_axis_scale != 0){
3139 im->xOriginLegendY2 += Xylabel;
3140 }
3141 im->yOriginLegendY2 = Ytitle + (Ymain / 2);
3143 break;
3145 case EAST:
3146 im->xOriginTitle = Xvertical + Xylabel + im->xsize / 2;
3147 im->yOriginTitle = 0;
3149 im->xOriginLegend = Xvertical + Xylabel + Xmain + Xvertical2;
3150 if (im->second_axis_scale != 0){
3151 im->xOriginLegend += Xylabel;
3152 }
3153 im->yOriginLegend = Ytitle;
3155 im->xOriginLegendY = 0;
3156 im->yOriginLegendY = Ytitle + (Ymain / 2);
3158 im->xorigin = Xvertical + Xylabel;
3159 im->yorigin = Ytitle + Ymain;
3161 im->xOriginLegendY2 = Xvertical + Xylabel + Xmain;
3162 if (im->second_axis_scale != 0){
3163 im->xOriginLegendY2 += Xylabel;
3164 }
3165 im->yOriginLegendY2 = Ytitle + (Ymain / 2);
3167 if (!(im->extra_flags & NO_RRDTOOL_TAG)){
3168 im->xOriginTitle += Xspacing;
3169 im->xOriginLegend += Xspacing;
3170 im->xOriginLegendY += Xspacing;
3171 im->xorigin += Xspacing;
3172 im->xOriginLegendY2 += Xspacing;
3173 }
3174 break;
3175 }
3177 xtr(im, 0);
3178 ytr(im, DNAN);
3179 return 0;
3180 }
3182 static cairo_status_t cairo_output(
3183 void *closure,
3184 const unsigned char
3185 *data,
3186 unsigned int length)
3187 {
3188 image_desc_t *im = (image_desc_t*)closure;
3190 im->rendered_image =
3191 (unsigned char*)realloc(im->rendered_image, im->rendered_image_size + length);
3192 if (im->rendered_image == NULL)
3193 return CAIRO_STATUS_WRITE_ERROR;
3194 memcpy(im->rendered_image + im->rendered_image_size, data, length);
3195 im->rendered_image_size += length;
3196 return CAIRO_STATUS_SUCCESS;
3197 }
3199 /* draw that picture thing ... */
3200 int graph_paint(
3201 image_desc_t *im)
3202 {
3203 int i, ii;
3204 int lazy = lazy_check(im);
3205 double areazero = 0.0;
3206 graph_desc_t *lastgdes = NULL;
3207 rrd_infoval_t info;
3209 // PangoFontMap *font_map = pango_cairo_font_map_get_default();
3211 /* pull the data from the rrd files ... */
3212 if (data_fetch(im) == -1)
3213 return -1;
3214 /* evaluate VDEF and CDEF operations ... */
3215 if (data_calc(im) == -1)
3216 return -1;
3217 /* calculate and PRINT and GPRINT definitions. We have to do it at
3218 * this point because it will affect the length of the legends
3219 * if there are no graph elements (i==0) we stop here ...
3220 * if we are lazy, try to quit ...
3221 */
3222 i = print_calc(im);
3223 if (i < 0)
3224 return -1;
3226 /* if we want and can be lazy ... quit now */
3227 if (i == 0)
3228 return 0;
3230 /**************************************************************
3231 *** Calculating sizes and locations became a bit confusing ***
3232 *** so I moved this into a separate function. ***
3233 **************************************************************/
3234 if (graph_size_location(im, i) == -1)
3235 return -1;
3237 info.u_cnt = im->xorigin;
3238 grinfo_push(im, sprintf_alloc("graph_left"), RD_I_CNT, info);
3239 info.u_cnt = im->yorigin - im->ysize;
3240 grinfo_push(im, sprintf_alloc("graph_top"), RD_I_CNT, info);
3241 info.u_cnt = im->xsize;
3242 grinfo_push(im, sprintf_alloc("graph_width"), RD_I_CNT, info);
3243 info.u_cnt = im->ysize;
3244 grinfo_push(im, sprintf_alloc("graph_height"), RD_I_CNT, info);
3245 info.u_cnt = im->ximg;
3246 grinfo_push(im, sprintf_alloc("image_width"), RD_I_CNT, info);
3247 info.u_cnt = im->yimg;
3248 grinfo_push(im, sprintf_alloc("image_height"), RD_I_CNT, info);
3249 info.u_cnt = im->start;
3250 grinfo_push(im, sprintf_alloc("graph_start"), RD_I_CNT, info);
3251 info.u_cnt = im->end;
3252 grinfo_push(im, sprintf_alloc("graph_end"), RD_I_CNT, info);
3254 /* if we want and can be lazy ... quit now */
3255 if (lazy)
3256 return 0;
3258 /* get actual drawing data and find min and max values */
3259 if (data_proc(im) == -1)
3260 return -1;
3261 if (!im->logarithmic) {
3262 si_unit(im);
3263 }
3265 /* identify si magnitude Kilo, Mega Giga ? */
3266 if (!im->rigid && !im->logarithmic)
3267 expand_range(im); /* make sure the upper and lower limit are
3268 sensible values */
3270 info.u_val = im->minval;
3271 grinfo_push(im, sprintf_alloc("value_min"), RD_I_VAL, info);
3272 info.u_val = im->maxval;
3273 grinfo_push(im, sprintf_alloc("value_max"), RD_I_VAL, info);
3276 if (!calc_horizontal_grid(im))
3277 return -1;
3278 /* reset precalc */
3279 ytr(im, DNAN);
3280 /* if (im->gridfit)
3281 apply_gridfit(im); */
3282 /* the actual graph is created by going through the individual
3283 graph elements and then drawing them */
3284 cairo_surface_destroy(im->surface);
3285 switch (im->imgformat) {
3286 case IF_PNG:
3287 im->surface =
3288 cairo_image_surface_create(CAIRO_FORMAT_ARGB32,
3289 im->ximg * im->zoom,
3290 im->yimg * im->zoom);
3291 break;
3292 case IF_PDF:
3293 im->gridfit = 0;
3294 im->surface = strlen(im->graphfile)
3295 ? cairo_pdf_surface_create(im->graphfile, im->ximg * im->zoom,
3296 im->yimg * im->zoom)
3297 : cairo_pdf_surface_create_for_stream
3298 (&cairo_output, im, im->ximg * im->zoom, im->yimg * im->zoom);
3299 break;
3300 case IF_EPS:
3301 im->gridfit = 0;
3302 im->surface = strlen(im->graphfile)
3303 ?
3304 cairo_ps_surface_create(im->graphfile, im->ximg * im->zoom,
3305 im->yimg * im->zoom)
3306 : cairo_ps_surface_create_for_stream
3307 (&cairo_output, im, im->ximg * im->zoom, im->yimg * im->zoom);
3308 break;
3309 case IF_SVG:
3310 im->gridfit = 0;
3311 im->surface = strlen(im->graphfile)
3312 ?
3313 cairo_svg_surface_create(im->
3314 graphfile,
3315 im->ximg * im->zoom, im->yimg * im->zoom)
3316 : cairo_svg_surface_create_for_stream
3317 (&cairo_output, im, im->ximg * im->zoom, im->yimg * im->zoom);
3318 cairo_svg_surface_restrict_to_version
3319 (im->surface, CAIRO_SVG_VERSION_1_1);
3320 break;
3321 };
3322 cairo_destroy(im->cr);
3323 im->cr = cairo_create(im->surface);
3324 cairo_set_antialias(im->cr, im->graph_antialias);
3325 cairo_scale(im->cr, im->zoom, im->zoom);
3326 // pango_cairo_font_map_set_resolution(PANGO_CAIRO_FONT_MAP(font_map), 100);
3327 gfx_new_area(im, 0, 0, 0, im->yimg,
3328 im->ximg, im->yimg, im->graph_col[GRC_BACK]);
3329 gfx_add_point(im, im->ximg, 0);
3330 gfx_close_path(im);
3331 gfx_new_area(im, im->xorigin,
3332 im->yorigin,
3333 im->xorigin +
3334 im->xsize, im->yorigin,
3335 im->xorigin +
3336 im->xsize,
3337 im->yorigin - im->ysize, im->graph_col[GRC_CANVAS]);
3338 gfx_add_point(im, im->xorigin, im->yorigin - im->ysize);
3339 gfx_close_path(im);
3340 cairo_rectangle(im->cr, im->xorigin, im->yorigin - im->ysize - 1.0,
3341 im->xsize, im->ysize + 2.0);
3342 cairo_clip(im->cr);
3343 if (im->minval > 0.0)
3344 areazero = im->minval;
3345 if (im->maxval < 0.0)
3346 areazero = im->maxval;
3347 for (i = 0; i < im->gdes_c; i++) {
3348 switch (im->gdes[i].gf) {
3349 case GF_CDEF:
3350 case GF_VDEF:
3351 case GF_DEF:
3352 case GF_PRINT:
3353 case GF_GPRINT:
3354 case GF_COMMENT:
3355 case GF_TEXTALIGN:
3356 case GF_HRULE:
3357 case GF_VRULE:
3358 case GF_XPORT:
3359 case GF_SHIFT:
3360 break;
3361 case GF_TICK:
3362 for (ii = 0; ii < im->xsize; ii++) {
3363 if (!isnan(im->gdes[i].p_data[ii])
3364 && im->gdes[i].p_data[ii] != 0.0) {
3365 if (im->gdes[i].yrule > 0) {
3366 gfx_line(im,
3367 im->xorigin + ii,
3368 im->yorigin + 1.0,
3369 im->xorigin + ii,
3370 im->yorigin -
3371 im->gdes[i].yrule *
3372 im->ysize, 1.0, im->gdes[i].col);
3373 } else if (im->gdes[i].yrule < 0) {
3374 gfx_line(im,
3375 im->xorigin + ii,
3376 im->yorigin - im->ysize - 1.0,
3377 im->xorigin + ii,
3378 im->yorigin - im->ysize -
3379 im->gdes[i].
3380 yrule *
3381 im->ysize, 1.0, im->gdes[i].col);
3382 }
3383 }
3384 }
3385 break;
3386 case GF_LINE:
3387 case GF_AREA:
3388 /* fix data points at oo and -oo */
3389 for (ii = 0; ii < im->xsize; ii++) {
3390 if (isinf(im->gdes[i].p_data[ii])) {
3391 if (im->gdes[i].p_data[ii] > 0) {
3392 im->gdes[i].p_data[ii] = im->maxval;
3393 } else {
3394 im->gdes[i].p_data[ii] = im->minval;
3395 }
3397 }
3398 } /* for */
3400 /* *******************************************************
3401 a ___. (a,t)
3402 | | ___
3403 ____| | | |
3404 | |___|
3405 -------|--t-1--t--------------------------------
3407 if we know the value at time t was a then
3408 we draw a square from t-1 to t with the value a.
3410 ********************************************************* */
3411 if (im->gdes[i].col.alpha != 0.0) {
3412 /* GF_LINE and friend */
3413 if (im->gdes[i].gf == GF_LINE) {
3414 double last_y = 0.0;
3415 int draw_on = 0;
3417 cairo_save(im->cr);
3418 cairo_new_path(im->cr);
3419 cairo_set_line_width(im->cr, im->gdes[i].linewidth);
3420 if (im->gdes[i].dash) {
3421 cairo_set_dash(im->cr,
3422 im->gdes[i].p_dashes,
3423 im->gdes[i].ndash, im->gdes[i].offset);
3424 }
3426 for (ii = 1; ii < im->xsize; ii++) {
3427 if (isnan(im->gdes[i].p_data[ii])
3428 || (im->slopemode == 1
3429 && isnan(im->gdes[i].p_data[ii - 1]))) {
3430 draw_on = 0;
3431 continue;
3432 }
3433 if (draw_on == 0) {
3434 last_y = ytr(im, im->gdes[i].p_data[ii]);
3435 if (im->slopemode == 0) {
3436 double x = ii - 1 + im->xorigin;
3437 double y = last_y;
3439 gfx_line_fit(im, &x, &y);
3440 cairo_move_to(im->cr, x, y);
3441 x = ii + im->xorigin;
3442 y = last_y;
3443 gfx_line_fit(im, &x, &y);
3444 cairo_line_to(im->cr, x, y);
3445 } else {
3446 double x = ii - 1 + im->xorigin;
3447 double y =
3448 ytr(im, im->gdes[i].p_data[ii - 1]);
3449 gfx_line_fit(im, &x, &y);
3450 cairo_move_to(im->cr, x, y);
3451 x = ii + im->xorigin;
3452 y = last_y;
3453 gfx_line_fit(im, &x, &y);
3454 cairo_line_to(im->cr, x, y);
3455 }
3456 draw_on = 1;
3457 } else {
3458 double x1 = ii + im->xorigin;
3459 double y1 = ytr(im, im->gdes[i].p_data[ii]);
3461 if (im->slopemode == 0
3462 && !AlmostEqual2sComplement(y1, last_y, 4)) {
3463 double x = ii - 1 + im->xorigin;
3464 double y = y1;
3466 gfx_line_fit(im, &x, &y);
3467 cairo_line_to(im->cr, x, y);
3468 };
3469 last_y = y1;
3470 gfx_line_fit(im, &x1, &y1);
3471 cairo_line_to(im->cr, x1, y1);
3472 };
3473 }
3474 cairo_set_source_rgba(im->cr,
3475 im->gdes[i].
3476 col.red,
3477 im->gdes[i].
3478 col.green,
3479 im->gdes[i].
3480 col.blue, im->gdes[i].col.alpha);
3481 cairo_set_line_cap(im->cr, CAIRO_LINE_CAP_ROUND);
3482 cairo_set_line_join(im->cr, CAIRO_LINE_JOIN_ROUND);
3483 cairo_stroke(im->cr);
3484 cairo_restore(im->cr);
3485 } else {
3486 int idxI = -1;
3487 double *foreY =
3488 (double *) malloc(sizeof(double) * im->xsize * 2);
3489 double *foreX =
3490 (double *) malloc(sizeof(double) * im->xsize * 2);
3491 double *backY =
3492 (double *) malloc(sizeof(double) * im->xsize * 2);
3493 double *backX =
3494 (double *) malloc(sizeof(double) * im->xsize * 2);
3495 int drawem = 0;
3497 for (ii = 0; ii <= im->xsize; ii++) {
3498 double ybase, ytop;
3500 if (idxI > 0 && (drawem != 0 || ii == im->xsize)) {
3501 int cntI = 1;
3502 int lastI = 0;
3504 while (cntI < idxI
3505 &&
3506 AlmostEqual2sComplement(foreY
3507 [lastI],
3508 foreY[cntI], 4)
3509 &&
3510 AlmostEqual2sComplement(foreY
3511 [lastI],
3512 foreY
3513 [cntI + 1], 4)) {
3514 cntI++;
3515 }
3516 gfx_new_area(im,
3517 backX[0], backY[0],
3518 foreX[0], foreY[0],
3519 foreX[cntI],
3520 foreY[cntI], im->gdes[i].col);
3521 while (cntI < idxI) {
3522 lastI = cntI;
3523 cntI++;
3524 while (cntI < idxI
3525 &&
3526 AlmostEqual2sComplement(foreY
3527 [lastI],
3528 foreY[cntI], 4)
3529 &&
3530 AlmostEqual2sComplement(foreY
3531 [lastI],
3532 foreY
3533 [cntI
3534 + 1], 4)) {
3535 cntI++;
3536 }
3537 gfx_add_point(im, foreX[cntI], foreY[cntI]);
3538 }
3539 gfx_add_point(im, backX[idxI], backY[idxI]);
3540 while (idxI > 1) {
3541 lastI = idxI;
3542 idxI--;
3543 while (idxI > 1
3544 &&
3545 AlmostEqual2sComplement(backY
3546 [lastI],
3547 backY[idxI], 4)
3548 &&
3549 AlmostEqual2sComplement(backY
3550 [lastI],
3551 backY
3552 [idxI
3553 - 1], 4)) {
3554 idxI--;
3555 }
3556 gfx_add_point(im, backX[idxI], backY[idxI]);
3557 }
3558 idxI = -1;
3559 drawem = 0;
3560 gfx_close_path(im);
3561 }
3562 if (drawem != 0) {
3563 drawem = 0;
3564 idxI = -1;
3565 }
3566 if (ii == im->xsize)
3567 break;
3568 if (im->slopemode == 0 && ii == 0) {
3569 continue;
3570 }
3571 if (isnan(im->gdes[i].p_data[ii])) {
3572 drawem = 1;
3573 continue;
3574 }
3575 ytop = ytr(im, im->gdes[i].p_data[ii]);
3576 if (lastgdes && im->gdes[i].stack) {
3577 ybase = ytr(im, lastgdes->p_data[ii]);
3578 } else {
3579 ybase = ytr(im, areazero);
3580 }
3581 if (ybase == ytop) {
3582 drawem = 1;
3583 continue;
3584 }
3586 if (ybase > ytop) {
3587 double extra = ytop;
3589 ytop = ybase;
3590 ybase = extra;
3591 }
3592 if (im->slopemode == 0) {
3593 backY[++idxI] = ybase - 0.2;
3594 backX[idxI] = ii + im->xorigin - 1;
3595 foreY[idxI] = ytop + 0.2;
3596 foreX[idxI] = ii + im->xorigin - 1;
3597 }
3598 backY[++idxI] = ybase - 0.2;
3599 backX[idxI] = ii + im->xorigin;
3600 foreY[idxI] = ytop + 0.2;
3601 foreX[idxI] = ii + im->xorigin;
3602 }
3603 /* close up any remaining area */
3604 free(foreY);
3605 free(foreX);
3606 free(backY);
3607 free(backX);
3608 } /* else GF_LINE */
3609 }
3610 /* if color != 0x0 */
3611 /* make sure we do not run into trouble when stacking on NaN */
3612 for (ii = 0; ii < im->xsize; ii++) {
3613 if (isnan(im->gdes[i].p_data[ii])) {
3614 if (lastgdes && (im->gdes[i].stack)) {
3615 im->gdes[i].p_data[ii] = lastgdes->p_data[ii];
3616 } else {
3617 im->gdes[i].p_data[ii] = areazero;
3618 }
3619 }
3620 }
3621 lastgdes = &(im->gdes[i]);
3622 break;
3623 case GF_STACK:
3624 rrd_set_error
3625 ("STACK should already be turned into LINE or AREA here");
3626 return -1;
3627 break;
3628 } /* switch */
3629 }
3630 cairo_reset_clip(im->cr);
3632 /* grid_paint also does the text */
3633 if (!(im->extra_flags & ONLY_GRAPH))
3634 grid_paint(im);
3635 if (!(im->extra_flags & ONLY_GRAPH))
3636 axis_paint(im);
3637 /* the RULES are the last thing to paint ... */
3638 for (i = 0; i < im->gdes_c; i++) {
3640 switch (im->gdes[i].gf) {
3641 case GF_HRULE:
3642 if (im->gdes[i].yrule >= im->minval
3643 && im->gdes[i].yrule <= im->maxval) {
3644 cairo_save(im->cr);
3645 if (im->gdes[i].dash) {
3646 cairo_set_dash(im->cr,
3647 im->gdes[i].p_dashes,
3648 im->gdes[i].ndash, im->gdes[i].offset);
3649 }
3650 gfx_line(im, im->xorigin,
3651 ytr(im, im->gdes[i].yrule),
3652 im->xorigin + im->xsize,
3653 ytr(im, im->gdes[i].yrule), 1.0, im->gdes[i].col);
3654 cairo_stroke(im->cr);
3655 cairo_restore(im->cr);
3656 }
3657 break;
3658 case GF_VRULE:
3659 if (im->gdes[i].xrule >= im->start
3660 && im->gdes[i].xrule <= im->end) {
3661 cairo_save(im->cr);
3662 if (im->gdes[i].dash) {
3663 cairo_set_dash(im->cr,
3664 im->gdes[i].p_dashes,
3665 im->gdes[i].ndash, im->gdes[i].offset);
3666 }
3667 gfx_line(im,
3668 xtr(im, im->gdes[i].xrule),
3669 im->yorigin, xtr(im,
3670 im->
3671 gdes[i].
3672 xrule),
3673 im->yorigin - im->ysize, 1.0, im->gdes[i].col);
3674 cairo_stroke(im->cr);
3675 cairo_restore(im->cr);
3676 }
3677 break;
3678 default:
3679 break;
3680 }
3681 }
3684 switch (im->imgformat) {
3685 case IF_PNG:
3686 {
3687 cairo_status_t status;
3689 status = strlen(im->graphfile) ?
3690 cairo_surface_write_to_png(im->surface, im->graphfile)
3691 : cairo_surface_write_to_png_stream(im->surface, &cairo_output,
3692 im);
3694 if (status != CAIRO_STATUS_SUCCESS) {
3695 rrd_set_error("Could not save png to '%s'", im->graphfile);
3696 return 1;
3697 }
3698 break;
3699 }
3700 default:
3701 if (strlen(im->graphfile)) {
3702 cairo_show_page(im->cr);
3703 } else {
3704 cairo_surface_finish(im->surface);
3705 }
3706 break;
3707 }
3709 return 0;
3710 }
3713 /*****************************************************
3714 * graph stuff
3715 *****************************************************/
3717 int gdes_alloc(
3718 image_desc_t *im)
3719 {
3721 im->gdes_c++;
3722 if ((im->gdes = (graph_desc_t *)
3723 rrd_realloc(im->gdes, (im->gdes_c)
3724 * sizeof(graph_desc_t))) == NULL) {
3725 rrd_set_error("realloc graph_descs");
3726 return -1;
3727 }
3730 im->gdes[im->gdes_c - 1].step = im->step;
3731 im->gdes[im->gdes_c - 1].step_orig = im->step;
3732 im->gdes[im->gdes_c - 1].stack = 0;
3733 im->gdes[im->gdes_c - 1].linewidth = 0;
3734 im->gdes[im->gdes_c - 1].debug = 0;
3735 im->gdes[im->gdes_c - 1].start = im->start;
3736 im->gdes[im->gdes_c - 1].start_orig = im->start;
3737 im->gdes[im->gdes_c - 1].end = im->end;
3738 im->gdes[im->gdes_c - 1].end_orig = im->end;
3739 im->gdes[im->gdes_c - 1].vname[0] = '\0';
3740 im->gdes[im->gdes_c - 1].data = NULL;
3741 im->gdes[im->gdes_c - 1].ds_namv = NULL;
3742 im->gdes[im->gdes_c - 1].data_first = 0;
3743 im->gdes[im->gdes_c - 1].p_data = NULL;
3744 im->gdes[im->gdes_c - 1].rpnp = NULL;
3745 im->gdes[im->gdes_c - 1].p_dashes = NULL;
3746 im->gdes[im->gdes_c - 1].shift = 0.0;
3747 im->gdes[im->gdes_c - 1].dash = 0;
3748 im->gdes[im->gdes_c - 1].ndash = 0;
3749 im->gdes[im->gdes_c - 1].offset = 0;
3750 im->gdes[im->gdes_c - 1].col.red = 0.0;
3751 im->gdes[im->gdes_c - 1].col.green = 0.0;
3752 im->gdes[im->gdes_c - 1].col.blue = 0.0;
3753 im->gdes[im->gdes_c - 1].col.alpha = 0.0;
3754 im->gdes[im->gdes_c - 1].legend[0] = '\0';
3755 im->gdes[im->gdes_c - 1].format[0] = '\0';
3756 im->gdes[im->gdes_c - 1].strftm = 0;
3757 im->gdes[im->gdes_c - 1].rrd[0] = '\0';
3758 im->gdes[im->gdes_c - 1].ds = -1;
3759 im->gdes[im->gdes_c - 1].cf_reduce = CF_AVERAGE;
3760 im->gdes[im->gdes_c - 1].cf = CF_AVERAGE;
3761 im->gdes[im->gdes_c - 1].yrule = DNAN;
3762 im->gdes[im->gdes_c - 1].xrule = 0;
3763 return 0;
3764 }
3766 /* copies input untill the first unescaped colon is found
3767 or until input ends. backslashes have to be escaped as well */
3768 int scan_for_col(
3769 const char *const input,
3770 int len,
3771 char *const output)
3772 {
3773 int inp, outp = 0;
3775 for (inp = 0; inp < len && input[inp] != ':' && input[inp] != '\0'; inp++) {
3776 if (input[inp] == '\\'
3777 && input[inp + 1] != '\0'
3778 && (input[inp + 1] == '\\' || input[inp + 1] == ':')) {
3779 output[outp++] = input[++inp];
3780 } else {
3781 output[outp++] = input[inp];
3782 }
3783 }
3784 output[outp] = '\0';
3785 return inp;
3786 }
3788 /* Now just a wrapper around rrd_graph_v */
3789 int rrd_graph(
3790 int argc,
3791 char **argv,
3792 char ***prdata,
3793 int *xsize,
3794 int *ysize,
3795 FILE * stream,
3796 double *ymin,
3797 double *ymax)
3798 {
3799 int prlines = 0;
3800 rrd_info_t *grinfo = NULL;
3801 rrd_info_t *walker;
3803 grinfo = rrd_graph_v(argc, argv);
3804 if (grinfo == NULL)
3805 return -1;
3806 walker = grinfo;
3807 (*prdata) = NULL;
3808 while (walker) {
3809 if (strcmp(walker->key, "image_info") == 0) {
3810 prlines++;
3811 if (((*prdata) =
3812 (char**)rrd_realloc((*prdata),
3813 (prlines + 1) * sizeof(char *))) == NULL) {
3814 rrd_set_error("realloc prdata");
3815 return 0;
3816 }
3817 /* imginfo goes to position 0 in the prdata array */
3818 (*prdata)[prlines - 1] = (char*)malloc((strlen(walker->value.u_str)
3819 + 2) * sizeof(char));
3820 strcpy((*prdata)[prlines - 1], walker->value.u_str);
3821 (*prdata)[prlines] = NULL;
3822 }
3823 /* skip anything else */
3824 walker = walker->next;
3825 }
3826 walker = grinfo;
3827 *xsize = 0;
3828 *ysize = 0;
3829 *ymin = 0;
3830 *ymax = 0;
3831 while (walker) {
3832 if (strcmp(walker->key, "image_width") == 0) {
3833 *xsize = walker->value.u_cnt;
3834 } else if (strcmp(walker->key, "image_height") == 0) {
3835 *ysize = walker->value.u_cnt;
3836 } else if (strcmp(walker->key, "value_min") == 0) {
3837 *ymin = walker->value.u_val;
3838 } else if (strcmp(walker->key, "value_max") == 0) {
3839 *ymax = walker->value.u_val;
3840 } else if (strncmp(walker->key, "print", 5) == 0) { /* keys are prdate[0..] */
3841 prlines++;
3842 if (((*prdata) =
3843 (char**)rrd_realloc((*prdata),
3844 (prlines + 1) * sizeof(char *))) == NULL) {
3845 rrd_set_error("realloc prdata");
3846 return 0;
3847 }
3848 (*prdata)[prlines - 1] = (char*)malloc((strlen(walker->value.u_str)
3849 + 2) * sizeof(char));
3850 (*prdata)[prlines] = NULL;
3851 strcpy((*prdata)[prlines - 1], walker->value.u_str);
3852 } else if (strcmp(walker->key, "image") == 0) {
3853 if ( fwrite(walker->value.u_blo.ptr, walker->value.u_blo.size, 1,
3854 (stream ? stream : stdout)) == 0 && ferror(stream ? stream : stdout)){
3855 rrd_set_error("writing image");
3856 return 0;
3857 }
3858 }
3859 /* skip anything else */
3860 walker = walker->next;
3861 }
3862 rrd_info_free(grinfo);
3863 return 0;
3864 }
3867 /* Some surgery done on this function, it became ridiculously big.
3868 ** Things moved:
3869 ** - initializing now in rrd_graph_init()
3870 ** - options parsing now in rrd_graph_options()
3871 ** - script parsing now in rrd_graph_script()
3872 */
3873 rrd_info_t *rrd_graph_v(
3874 int argc,
3875 char **argv)
3876 {
3877 image_desc_t im;
3878 rrd_info_t *grinfo;
3879 rrd_graph_init(&im);
3880 /* a dummy surface so that we can measure text sizes for placements */
3882 rrd_graph_options(argc, argv, &im);
3883 if (rrd_test_error()) {
3884 rrd_info_free(im.grinfo);
3885 im_free(&im);
3886 return NULL;
3887 }
3889 if (optind >= argc) {
3890 rrd_info_free(im.grinfo);
3891 im_free(&im);
3892 rrd_set_error("missing filename");
3893 return NULL;
3894 }
3896 if (strlen(argv[optind]) >= MAXPATH) {
3897 rrd_set_error("filename (including path) too long");
3898 rrd_info_free(im.grinfo);
3899 im_free(&im);
3900 return NULL;
3901 }
3903 strncpy(im.graphfile, argv[optind], MAXPATH - 1);
3904 im.graphfile[MAXPATH - 1] = '\0';
3906 if (strcmp(im.graphfile, "-") == 0) {
3907 im.graphfile[0] = '\0';
3908 }
3910 rrd_graph_script(argc, argv, &im, 1);
3911 if (rrd_test_error()) {
3912 rrd_info_free(im.grinfo);
3913 im_free(&im);
3914 return NULL;
3915 }
3917 /* Everything is now read and the actual work can start */
3919 if (graph_paint(&im) == -1) {
3920 rrd_info_free(im.grinfo);
3921 im_free(&im);
3922 return NULL;
3923 }
3926 /* The image is generated and needs to be output.
3927 ** Also, if needed, print a line with information about the image.
3928 */
3930 if (im.imginfo) {
3931 rrd_infoval_t info;
3932 char *path;
3933 char *filename;
3935 path = strdup(im.graphfile);
3936 filename = basename(path);
3937 info.u_str =
3938 sprintf_alloc(im.imginfo,
3939 filename,
3940 (long) (im.zoom *
3941 im.ximg), (long) (im.zoom * im.yimg));
3942 grinfo_push(&im, sprintf_alloc("image_info"), RD_I_STR, info);
3943 free(info.u_str);
3944 free(path);
3945 }
3946 if (im.rendered_image) {
3947 rrd_infoval_t img;
3949 img.u_blo.size = im.rendered_image_size;
3950 img.u_blo.ptr = im.rendered_image;
3951 grinfo_push(&im, sprintf_alloc("image"), RD_I_BLO, img);
3952 }
3953 grinfo = im.grinfo;
3954 im_free(&im);
3955 return grinfo;
3956 }
3958 static void
3959 rrd_set_font_desc (
3960 image_desc_t *im,int prop,char *font, double size ){
3961 if (font){
3962 strncpy(im->text_prop[prop].font, font, sizeof(text_prop[prop].font) - 1);
3963 im->text_prop[prop].font[sizeof(text_prop[prop].font) - 1] = '\0';
3964 im->text_prop[prop].font_desc = pango_font_description_from_string( font );
3965 };
3966 if (size > 0){
3967 im->text_prop[prop].size = size;
3968 };
3969 if (im->text_prop[prop].font_desc && im->text_prop[prop].size ){
3970 pango_font_description_set_size(im->text_prop[prop].font_desc, im->text_prop[prop].size * PANGO_SCALE);
3971 };
3972 }
3974 void rrd_graph_init(
3975 image_desc_t
3976 *im)
3977 {
3978 unsigned int i;
3979 char *deffont = getenv("RRD_DEFAULT_FONT");
3980 static PangoFontMap *fontmap = NULL;
3981 PangoContext *context;
3983 #ifdef HAVE_TZSET
3984 tzset();
3985 #endif
3986 #ifdef HAVE_SETLOCALE
3987 setlocale(LC_TIME, "");
3988 #ifdef HAVE_MBSTOWCS
3989 setlocale(LC_CTYPE, "");
3990 #endif
3991 #endif
3992 im->base = 1000;
3993 im->daemon_addr = NULL;
3994 im->draw_x_grid = 1;
3995 im->draw_y_grid = 1;
3996 im->draw_3d_border = 2;
3997 im->extra_flags = 0;
3998 im->font_options = cairo_font_options_create();
3999 im->forceleftspace = 0;
4000 im->gdes_c = 0;
4001 im->gdes = NULL;
4002 im->graph_antialias = CAIRO_ANTIALIAS_GRAY;
4003 im->grid_dash_off = 1;
4004 im->grid_dash_on = 1;
4005 im->gridfit = 1;
4006 im->grinfo = (rrd_info_t *) NULL;
4007 im->grinfo_current = (rrd_info_t *) NULL;
4008 im->imgformat = IF_PNG;
4009 im->imginfo = NULL;
4010 im->lazy = 0;
4011 im->legenddirection = TOP_DOWN;
4012 im->legendheight = 0;
4013 im->legendposition = SOUTH;
4014 im->legendwidth = 0;
4015 im->logarithmic = 0;
4016 im->maxval = DNAN;
4017 im->minval = 0;
4018 im->minval = DNAN;
4019 im->prt_c = 0;
4020 im->rigid = 0;
4021 im->rendered_image_size = 0;
4022 im->rendered_image = NULL;
4023 im->slopemode = 0;
4024 im->step = 0;
4025 im->symbol = ' ';
4026 im->tabwidth = 40.0;
4027 im->title[0] = '\0';
4028 im->unitsexponent = 9999;
4029 im->unitslength = 6;
4030 im->viewfactor = 1.0;
4031 im->watermark[0] = '\0';
4032 im->with_markup = 0;
4033 im->ximg = 0;
4034 im->xlab_user.minsec = -1;
4035 im->xorigin = 0;
4036 im->xOriginLegend = 0;
4037 im->xOriginLegendY = 0;
4038 im->xOriginLegendY2 = 0;
4039 im->xOriginTitle = 0;
4040 im->xsize = 400;
4041 im->ygridstep = DNAN;
4042 im->yimg = 0;
4043 im->ylegend[0] = '\0';
4044 im->second_axis_scale = 0; /* 0 disables it */
4045 im->second_axis_shift = 0; /* no shift by default */
4046 im->second_axis_legend[0] = '\0';
4047 im->second_axis_format[0] = '\0';
4048 im->yorigin = 0;
4049 im->yOriginLegend = 0;
4050 im->yOriginLegendY = 0;
4051 im->yOriginLegendY2 = 0;
4052 im->yOriginTitle = 0;
4053 im->ysize = 100;
4054 im->zoom = 1;
4056 im->surface = cairo_image_surface_create(CAIRO_FORMAT_ARGB32, 10, 10);
4057 im->cr = cairo_create(im->surface);
4059 for (i = 0; i < DIM(text_prop); i++) {
4060 im->text_prop[i].size = -1;
4061 rrd_set_font_desc(im,i, deffont ? deffont : text_prop[i].font,text_prop[i].size);
4062 }
4064 if (fontmap == NULL){
4065 fontmap = pango_cairo_font_map_get_default();
4066 }
4068 context = pango_cairo_font_map_create_context((PangoCairoFontMap*)fontmap);
4070 pango_cairo_context_set_resolution(context, 100);
4072 pango_cairo_update_context(im->cr,context);
4074 im->layout = pango_layout_new(context);
4076 // im->layout = pango_cairo_create_layout(im->cr);
4079 cairo_font_options_set_hint_style
4080 (im->font_options, CAIRO_HINT_STYLE_FULL);
4081 cairo_font_options_set_hint_metrics
4082 (im->font_options, CAIRO_HINT_METRICS_ON);
4083 cairo_font_options_set_antialias(im->font_options, CAIRO_ANTIALIAS_GRAY);
4087 for (i = 0; i < DIM(graph_col); i++)
4088 im->graph_col[i] = graph_col[i];
4091 }
4094 void rrd_graph_options(
4095 int argc,
4096 char *argv[],
4097 image_desc_t
4098 *im)
4099 {
4100 int stroff;
4101 char *parsetime_error = NULL;
4102 char scan_gtm[12], scan_mtm[12], scan_ltm[12], col_nam[12];
4103 time_t start_tmp = 0, end_tmp = 0;
4104 long long_tmp;
4105 rrd_time_value_t start_tv, end_tv;
4106 long unsigned int color;
4107 char *old_locale = "";
4109 /* defines for long options without a short equivalent. should be bytes,
4110 and may not collide with (the ASCII value of) short options */
4111 #define LONGOPT_UNITS_SI 255
4113 /* *INDENT-OFF* */
4114 struct option long_options[] = {
4115 { "alt-autoscale", no_argument, 0, 'A'},
4116 { "imgformat", required_argument, 0, 'a'},
4117 { "font-smoothing-threshold", required_argument, 0, 'B'},
4118 { "base", required_argument, 0, 'b'},
4119 { "color", required_argument, 0, 'c'},
4120 { "full-size-mode", no_argument, 0, 'D'},
4121 { "daemon", required_argument, 0, 'd'},
4122 { "slope-mode", no_argument, 0, 'E'},
4123 { "end", required_argument, 0, 'e'},
4124 { "force-rules-legend", no_argument, 0, 'F'},
4125 { "imginfo", required_argument, 0, 'f'},
4126 { "graph-render-mode", required_argument, 0, 'G'},
4127 { "no-legend", no_argument, 0, 'g'},
4128 { "height", required_argument, 0, 'h'},
4129 { "no-minor", no_argument, 0, 'I'},
4130 { "interlaced", no_argument, 0, 'i'},
4131 { "alt-autoscale-min", no_argument, 0, 'J'},
4132 { "only-graph", no_argument, 0, 'j'},
4133 { "units-length", required_argument, 0, 'L'},
4134 { "lower-limit", required_argument, 0, 'l'},
4135 { "alt-autoscale-max", no_argument, 0, 'M'},
4136 { "zoom", required_argument, 0, 'm'},
4137 { "no-gridfit", no_argument, 0, 'N'},
4138 { "font", required_argument, 0, 'n'},
4139 { "logarithmic", no_argument, 0, 'o'},
4140 { "pango-markup", no_argument, 0, 'P'},
4141 { "font-render-mode", required_argument, 0, 'R'},
4142 { "rigid", no_argument, 0, 'r'},
4143 { "step", required_argument, 0, 'S'},
4144 { "start", required_argument, 0, 's'},
4145 { "tabwidth", required_argument, 0, 'T'},
4146 { "title", required_argument, 0, 't'},
4147 { "upper-limit", required_argument, 0, 'u'},
4148 { "vertical-label", required_argument, 0, 'v'},
4149 { "watermark", required_argument, 0, 'W'},
4150 { "width", required_argument, 0, 'w'},
4151 { "units-exponent", required_argument, 0, 'X'},
4152 { "x-grid", required_argument, 0, 'x'},
4153 { "alt-y-grid", no_argument, 0, 'Y'},
4154 { "y-grid", required_argument, 0, 'y'},
4155 { "lazy", no_argument, 0, 'z'},
4156 { "units", required_argument, 0, LONGOPT_UNITS_SI},
4157 { "alt-y-mrtg", no_argument, 0, 1000}, /* this has no effect it is just here to save old apps from crashing when they use it */
4158 { "disable-rrdtool-tag",no_argument, 0, 1001},
4159 { "right-axis", required_argument, 0, 1002},
4160 { "right-axis-label", required_argument, 0, 1003},
4161 { "right-axis-format", required_argument, 0, 1004},
4162 { "legend-position", required_argument, 0, 1005},
4163 { "legend-direction", required_argument, 0, 1006},
4164 { "border", required_argument, 0, 1007},
4165 { 0, 0, 0, 0}
4166 };
4167 /* *INDENT-ON* */
4169 optind = 0;
4170 opterr = 0; /* initialize getopt */
4171 rrd_parsetime("end-24h", &start_tv);
4172 rrd_parsetime("now", &end_tv);
4173 while (1) {
4174 int option_index = 0;
4175 int opt;
4176 int col_start, col_end;
4178 opt = getopt_long(argc, argv,
4179 "Aa:B:b:c:Dd:Ee:Ff:G:gh:IiJjL:l:Nn:Bb:oPR:rS:s:T:t:u:v:W:w:X:x:Yy:z",
4180 long_options, &option_index);
4181 if (opt == EOF)
4182 break;
4183 switch (opt) {
4184 case 'I':
4185 im->extra_flags |= NOMINOR;
4186 break;
4187 case 'Y':
4188 im->extra_flags |= ALTYGRID;
4189 break;
4190 case 'A':
4191 im->extra_flags |= ALTAUTOSCALE;
4192 break;
4193 case 'J':
4194 im->extra_flags |= ALTAUTOSCALE_MIN;
4195 break;
4196 case 'M':
4197 im->extra_flags |= ALTAUTOSCALE_MAX;
4198 break;
4199 case 'j':
4200 im->extra_flags |= ONLY_GRAPH;
4201 break;
4202 case 'g':
4203 im->extra_flags |= NOLEGEND;
4204 break;
4205 case 1005:
4206 if (strcmp(optarg, "north") == 0) {
4207 im->legendposition = NORTH;
4208 } else if (strcmp(optarg, "west") == 0) {
4209 im->legendposition = WEST;
4210 } else if (strcmp(optarg, "south") == 0) {
4211 im->legendposition = SOUTH;
4212 } else if (strcmp(optarg, "east") == 0) {
4213 im->legendposition = EAST;
4214 } else {
4215 rrd_set_error("unknown legend-position '%s'", optarg);
4216 return;
4217 }
4218 break;
4219 case 1006:
4220 if (strcmp(optarg, "topdown") == 0) {
4221 im->legenddirection = TOP_DOWN;
4222 } else if (strcmp(optarg, "bottomup") == 0) {
4223 im->legenddirection = BOTTOM_UP;
4224 } else {
4225 rrd_set_error("unknown legend-position '%s'", optarg);
4226 return;
4227 }
4228 break;
4229 case 'F':
4230 im->extra_flags |= FORCE_RULES_LEGEND;
4231 break;
4232 case 1001:
4233 im->extra_flags |= NO_RRDTOOL_TAG;
4234 break;
4235 case LONGOPT_UNITS_SI:
4236 if (im->extra_flags & FORCE_UNITS) {
4237 rrd_set_error("--units can only be used once!");
4238 setlocale(LC_NUMERIC, old_locale);
4239 return;
4240 }
4241 if (strcmp(optarg, "si") == 0)
4242 im->extra_flags |= FORCE_UNITS_SI;
4243 else {
4244 rrd_set_error("invalid argument for --units: %s", optarg);
4245 return;
4246 }
4247 break;
4248 case 'X':
4249 im->unitsexponent = atoi(optarg);
4250 break;
4251 case 'L':
4252 im->unitslength = atoi(optarg);
4253 im->forceleftspace = 1;
4254 break;
4255 case 'T':
4256 old_locale = setlocale(LC_NUMERIC, "C");
4257 im->tabwidth = atof(optarg);
4258 setlocale(LC_NUMERIC, old_locale);
4259 break;
4260 case 'S':
4261 old_locale = setlocale(LC_NUMERIC, "C");
4262 im->step = atoi(optarg);
4263 setlocale(LC_NUMERIC, old_locale);
4264 break;
4265 case 'N':
4266 im->gridfit = 0;
4267 break;
4268 case 'P':
4269 im->with_markup = 1;
4270 break;
4271 case 's':
4272 if ((parsetime_error = rrd_parsetime(optarg, &start_tv))) {
4273 rrd_set_error("start time: %s", parsetime_error);
4274 return;
4275 }
4276 break;
4277 case 'e':
4278 if ((parsetime_error = rrd_parsetime(optarg, &end_tv))) {
4279 rrd_set_error("end time: %s", parsetime_error);
4280 return;
4281 }
4282 break;
4283 case 'x':
4284 if (strcmp(optarg, "none") == 0) {
4285 im->draw_x_grid = 0;
4286 break;
4287 };
4288 if (sscanf(optarg,
4289 "%10[A-Z]:%ld:%10[A-Z]:%ld:%10[A-Z]:%ld:%ld:%n",
4290 scan_gtm,
4291 &im->xlab_user.gridst,
4292 scan_mtm,
4293 &im->xlab_user.mgridst,
4294 scan_ltm,
4295 &im->xlab_user.labst,
4296 &im->xlab_user.precis, &stroff) == 7 && stroff != 0) {
4297 strncpy(im->xlab_form, optarg + stroff,
4298 sizeof(im->xlab_form) - 1);
4299 im->xlab_form[sizeof(im->xlab_form) - 1] = '\0';
4300 if ((int)
4301 (im->xlab_user.gridtm = tmt_conv(scan_gtm)) == -1) {
4302 rrd_set_error("unknown keyword %s", scan_gtm);
4303 return;
4304 } else if ((int)
4305 (im->xlab_user.mgridtm = tmt_conv(scan_mtm))
4306 == -1) {
4307 rrd_set_error("unknown keyword %s", scan_mtm);
4308 return;
4309 } else if ((int)
4310 (im->xlab_user.labtm = tmt_conv(scan_ltm)) == -1) {
4311 rrd_set_error("unknown keyword %s", scan_ltm);
4312 return;
4313 }
4314 im->xlab_user.minsec = 1;
4315 im->xlab_user.stst = im->xlab_form;
4316 } else {
4317 rrd_set_error("invalid x-grid format");
4318 return;
4319 }
4320 break;
4321 case 'y':
4323 if (strcmp(optarg, "none") == 0) {
4324 im->draw_y_grid = 0;
4325 break;
4326 };
4327 old_locale = setlocale(LC_NUMERIC, "C");
4328 if (sscanf(optarg, "%lf:%d", &im->ygridstep, &im->ylabfact) == 2) {
4329 setlocale(LC_NUMERIC, old_locale);
4330 if (im->ygridstep <= 0) {
4331 rrd_set_error("grid step must be > 0");
4332 return;
4333 } else if (im->ylabfact < 1) {
4334 rrd_set_error("label factor must be > 0");
4335 return;
4336 }
4337 } else {
4338 setlocale(LC_NUMERIC, old_locale);
4339 rrd_set_error("invalid y-grid format");
4340 return;
4341 }
4342 break;
4343 case 1007:
4344 im->draw_3d_border = atoi(optarg);
4345 break;
4346 case 1002: /* right y axis */
4348 if(sscanf(optarg,
4349 "%lf:%lf",
4350 &im->second_axis_scale,
4351 &im->second_axis_shift) == 2) {
4352 if(im->second_axis_scale==0){
4353 rrd_set_error("the second_axis_scale must not be 0");
4354 return;
4355 }
4356 } else {
4357 rrd_set_error("invalid right-axis format expected scale:shift");
4358 return;
4359 }
4360 break;
4361 case 1003:
4362 strncpy(im->second_axis_legend,optarg,150);
4363 im->second_axis_legend[150]='\0';
4364 break;
4365 case 1004:
4366 if (bad_format(optarg)){
4367 rrd_set_error("use either %le or %lf formats");
4368 return;
4369 }
4370 strncpy(im->second_axis_format,optarg,150);
4371 im->second_axis_format[150]='\0';
4372 break;
4373 case 'v':
4374 strncpy(im->ylegend, optarg, 150);
4375 im->ylegend[150] = '\0';
4376 break;
4377 case 'u':
4378 old_locale = setlocale(LC_NUMERIC, "C");
4379 im->maxval = atof(optarg);
4380 setlocale(LC_NUMERIC, old_locale);
4381 break;
4382 case 'l':
4383 old_locale = setlocale(LC_NUMERIC, "C");
4384 im->minval = atof(optarg);
4385 setlocale(LC_NUMERIC, old_locale);
4386 break;
4387 case 'b':
4388 im->base = atol(optarg);
4389 if (im->base != 1024 && im->base != 1000) {
4390 rrd_set_error
4391 ("the only sensible value for base apart from 1000 is 1024");
4392 return;
4393 }
4394 break;
4395 case 'w':
4396 long_tmp = atol(optarg);
4397 if (long_tmp < 10) {
4398 rrd_set_error("width below 10 pixels");
4399 return;
4400 }
4401 im->xsize = long_tmp;
4402 break;
4403 case 'h':
4404 long_tmp = atol(optarg);
4405 if (long_tmp < 10) {
4406 rrd_set_error("height below 10 pixels");
4407 return;
4408 }
4409 im->ysize = long_tmp;
4410 break;
4411 case 'D':
4412 im->extra_flags |= FULL_SIZE_MODE;
4413 break;
4414 case 'i':
4415 /* interlaced png not supported at the moment */
4416 break;
4417 case 'r':
4418 im->rigid = 1;
4419 break;
4420 case 'f':
4421 im->imginfo = optarg;
4422 break;
4423 case 'a':
4424 if ((int)
4425 (im->imgformat = if_conv(optarg)) == -1) {
4426 rrd_set_error("unsupported graphics format '%s'", optarg);
4427 return;
4428 }
4429 break;
4430 case 'z':
4431 im->lazy = 1;
4432 break;
4433 case 'E':
4434 im->slopemode = 1;
4435 break;
4436 case 'o':
4437 im->logarithmic = 1;
4438 break;
4439 case 'c':
4440 if (sscanf(optarg,
4441 "%10[A-Z]#%n%8lx%n",
4442 col_nam, &col_start, &color, &col_end) == 2) {
4443 int ci;
4444 int col_len = col_end - col_start;
4446 switch (col_len) {
4447 case 3:
4448 color =
4449 (((color & 0xF00) * 0x110000) | ((color & 0x0F0) *
4450 0x011000) |
4451 ((color & 0x00F)
4452 * 0x001100)
4453 | 0x000000FF);
4454 break;
4455 case 4:
4456 color =
4457 (((color & 0xF000) *
4458 0x11000) | ((color & 0x0F00) *
4459 0x01100) | ((color &
4460 0x00F0) *
4461 0x00110) |
4462 ((color & 0x000F) * 0x00011)
4463 );
4464 break;
4465 case 6:
4466 color = (color << 8) + 0xff /* shift left by 8 */ ;
4467 break;
4468 case 8:
4469 break;
4470 default:
4471 rrd_set_error("the color format is #RRGGBB[AA]");
4472 return;
4473 }
4474 if ((ci = grc_conv(col_nam)) != -1) {
4475 im->graph_col[ci] = gfx_hex_to_col(color);
4476 } else {
4477 rrd_set_error("invalid color name '%s'", col_nam);
4478 return;
4479 }
4480 } else {
4481 rrd_set_error("invalid color def format");
4482 return;
4483 }
4484 break;
4485 case 'n':{
4486 char prop[15];
4487 double size = 1;
4488 int end;
4490 old_locale = setlocale(LC_NUMERIC, "C");
4491 if (sscanf(optarg, "%10[A-Z]:%lf%n", prop, &size, &end) >= 2) {
4492 int sindex, propidx;
4494 setlocale(LC_NUMERIC, old_locale);
4495 if ((sindex = text_prop_conv(prop)) != -1) {
4496 for (propidx = sindex;
4497 propidx < TEXT_PROP_LAST; propidx++) {
4498 if (size > 0) {
4499 rrd_set_font_desc(im,propidx,NULL,size);
4500 }
4501 if ((int) strlen(optarg) > end+2) {
4502 if (optarg[end] == ':') {
4503 rrd_set_font_desc(im,propidx,optarg + end + 1,0);
4504 } else {
4505 rrd_set_error
4506 ("expected : after font size in '%s'",
4507 optarg);
4508 return;
4509 }
4510 }
4511 /* only run the for loop for DEFAULT (0) for
4512 all others, we break here. woodo programming */
4513 if (propidx == sindex && sindex != 0)
4514 break;
4515 }
4516 } else {
4517 rrd_set_error("invalid fonttag '%s'", prop);
4518 return;
4519 }
4520 } else {
4521 setlocale(LC_NUMERIC, old_locale);
4522 rrd_set_error("invalid text property format");
4523 return;
4524 }
4525 break;
4526 }
4527 case 'm':
4528 old_locale = setlocale(LC_NUMERIC, "C");
4529 im->zoom = atof(optarg);
4530 setlocale(LC_NUMERIC, old_locale);
4531 if (im->zoom <= 0.0) {
4532 rrd_set_error("zoom factor must be > 0");
4533 return;
4534 }
4535 break;
4536 case 't':
4537 strncpy(im->title, optarg, 150);
4538 im->title[150] = '\0';
4539 break;
4540 case 'R':
4541 if (strcmp(optarg, "normal") == 0) {
4542 cairo_font_options_set_antialias
4543 (im->font_options, CAIRO_ANTIALIAS_GRAY);
4544 cairo_font_options_set_hint_style
4545 (im->font_options, CAIRO_HINT_STYLE_FULL);
4546 } else if (strcmp(optarg, "light") == 0) {
4547 cairo_font_options_set_antialias
4548 (im->font_options, CAIRO_ANTIALIAS_GRAY);
4549 cairo_font_options_set_hint_style
4550 (im->font_options, CAIRO_HINT_STYLE_SLIGHT);
4551 } else if (strcmp(optarg, "mono") == 0) {
4552 cairo_font_options_set_antialias
4553 (im->font_options, CAIRO_ANTIALIAS_NONE);
4554 cairo_font_options_set_hint_style
4555 (im->font_options, CAIRO_HINT_STYLE_FULL);
4556 } else {
4557 rrd_set_error("unknown font-render-mode '%s'", optarg);
4558 return;
4559 }
4560 break;
4561 case 'G':
4562 if (strcmp(optarg, "normal") == 0)
4563 im->graph_antialias = CAIRO_ANTIALIAS_GRAY;
4564 else if (strcmp(optarg, "mono") == 0)
4565 im->graph_antialias = CAIRO_ANTIALIAS_NONE;
4566 else {
4567 rrd_set_error("unknown graph-render-mode '%s'", optarg);
4568 return;
4569 }
4570 break;
4571 case 'B':
4572 /* not supported curently */
4573 break;
4574 case 'W':
4575 strncpy(im->watermark, optarg, 100);
4576 im->watermark[99] = '\0';
4577 break;
4578 case 'd':
4579 {
4580 if (im->daemon_addr != NULL)
4581 {
4582 rrd_set_error ("You cannot specify --daemon "
4583 "more than once.");
4584 return;
4585 }
4587 im->daemon_addr = strdup(optarg);
4588 if (im->daemon_addr == NULL)
4589 {
4590 rrd_set_error("strdup failed");
4591 return;
4592 }
4594 break;
4595 }
4596 case '?':
4597 if (optopt != 0)
4598 rrd_set_error("unknown option '%c'", optopt);
4599 else
4600 rrd_set_error("unknown option '%s'", argv[optind - 1]);
4601 return;
4602 }
4603 } /* while (1) */
4605 { /* try to connect to rrdcached */
4606 int status = rrdc_connect(im->daemon_addr);
4607 if (status != 0) return;
4608 }
4610 pango_cairo_context_set_font_options(pango_layout_get_context(im->layout), im->font_options);
4611 pango_layout_context_changed(im->layout);
4615 if (im->logarithmic && im->minval <= 0) {
4616 rrd_set_error
4617 ("for a logarithmic yaxis you must specify a lower-limit > 0");
4618 return;
4619 }
4621 if (rrd_proc_start_end(&start_tv, &end_tv, &start_tmp, &end_tmp) == -1) {
4622 /* error string is set in rrd_parsetime.c */
4623 return;
4624 }
4626 if (start_tmp < 3600 * 24 * 365 * 10) {
4627 rrd_set_error
4628 ("the first entry to fetch should be after 1980 (%ld)",
4629 start_tmp);
4630 return;
4631 }
4633 if (end_tmp < start_tmp) {
4634 rrd_set_error
4635 ("start (%ld) should be less than end (%ld)", start_tmp, end_tmp);
4636 return;
4637 }
4639 im->start = start_tmp;
4640 im->end = end_tmp;
4641 im->step = max((long) im->step, (im->end - im->start) / im->xsize);
4642 }
4644 int rrd_graph_color(
4645 image_desc_t
4646 *im,
4647 char *var,
4648 char *err,
4649 int optional)
4650 {
4651 char *color;
4652 graph_desc_t *gdp = &im->gdes[im->gdes_c - 1];
4654 color = strstr(var, "#");
4655 if (color == NULL) {
4656 if (optional == 0) {
4657 rrd_set_error("Found no color in %s", err);
4658 return 0;
4659 }
4660 return 0;
4661 } else {
4662 int n = 0;
4663 char *rest;
4664 long unsigned int col;
4666 rest = strstr(color, ":");
4667 if (rest != NULL)
4668 n = rest - color;
4669 else
4670 n = strlen(color);
4671 switch (n) {
4672 case 7:
4673 sscanf(color, "#%6lx%n", &col, &n);
4674 col = (col << 8) + 0xff /* shift left by 8 */ ;
4675 if (n != 7)
4676 rrd_set_error("Color problem in %s", err);
4677 break;
4678 case 9:
4679 sscanf(color, "#%8lx%n", &col, &n);
4680 if (n == 9)
4681 break;
4682 default:
4683 rrd_set_error("Color problem in %s", err);
4684 }
4685 if (rrd_test_error())
4686 return 0;
4687 gdp->col = gfx_hex_to_col(col);
4688 return n;
4689 }
4690 }
4693 int bad_format(
4694 char *fmt)
4695 {
4696 char *ptr;
4697 int n = 0;
4699 ptr = fmt;
4700 while (*ptr != '\0')
4701 if (*ptr++ == '%') {
4703 /* line cannot end with percent char */
4704 if (*ptr == '\0')
4705 return 1;
4706 /* '%s', '%S' and '%%' are allowed */
4707 if (*ptr == 's' || *ptr == 'S' || *ptr == '%')
4708 ptr++;
4709 /* %c is allowed (but use only with vdef!) */
4710 else if (*ptr == 'c') {
4711 ptr++;
4712 n = 1;
4713 }
4715 /* or else '% 6.2lf' and such are allowed */
4716 else {
4717 /* optional padding character */
4718 if (*ptr == ' ' || *ptr == '+' || *ptr == '-')
4719 ptr++;
4720 /* This should take care of 'm.n' with all three optional */
4721 while (*ptr >= '0' && *ptr <= '9')
4722 ptr++;
4723 if (*ptr == '.')
4724 ptr++;
4725 while (*ptr >= '0' && *ptr <= '9')
4726 ptr++;
4727 /* Either 'le', 'lf' or 'lg' must follow here */
4728 if (*ptr++ != 'l')
4729 return 1;
4730 if (*ptr == 'e' || *ptr == 'f' || *ptr == 'g')
4731 ptr++;
4732 else
4733 return 1;
4734 n++;
4735 }
4736 }
4738 return (n != 1);
4739 }
4742 int vdef_parse(
4743 struct graph_desc_t
4744 *gdes,
4745 const char *const str)
4746 {
4747 /* A VDEF currently is either "func" or "param,func"
4748 * so the parsing is rather simple. Change if needed.
4749 */
4750 double param;
4751 char func[30];
4752 int n;
4753 char *old_locale;
4755 n = 0;
4756 old_locale = setlocale(LC_NUMERIC, "C");
4757 sscanf(str, "%le,%29[A-Z]%n", ¶m, func, &n);
4758 setlocale(LC_NUMERIC, old_locale);
4759 if (n == (int) strlen(str)) { /* matched */
4760 ;
4761 } else {
4762 n = 0;
4763 sscanf(str, "%29[A-Z]%n", func, &n);
4764 if (n == (int) strlen(str)) { /* matched */
4765 param = DNAN;
4766 } else {
4767 rrd_set_error
4768 ("Unknown function string '%s' in VDEF '%s'",
4769 str, gdes->vname);
4770 return -1;
4771 }
4772 }
4773 if (!strcmp("PERCENT", func))
4774 gdes->vf.op = VDEF_PERCENT;
4775 else if (!strcmp("PERCENTNAN", func))
4776 gdes->vf.op = VDEF_PERCENTNAN;
4777 else if (!strcmp("MAXIMUM", func))
4778 gdes->vf.op = VDEF_MAXIMUM;
4779 else if (!strcmp("AVERAGE", func))
4780 gdes->vf.op = VDEF_AVERAGE;
4781 else if (!strcmp("STDEV", func))
4782 gdes->vf.op = VDEF_STDEV;
4783 else if (!strcmp("MINIMUM", func))
4784 gdes->vf.op = VDEF_MINIMUM;
4785 else if (!strcmp("TOTAL", func))
4786 gdes->vf.op = VDEF_TOTAL;
4787 else if (!strcmp("FIRST", func))
4788 gdes->vf.op = VDEF_FIRST;
4789 else if (!strcmp("LAST", func))
4790 gdes->vf.op = VDEF_LAST;
4791 else if (!strcmp("LSLSLOPE", func))
4792 gdes->vf.op = VDEF_LSLSLOPE;
4793 else if (!strcmp("LSLINT", func))
4794 gdes->vf.op = VDEF_LSLINT;
4795 else if (!strcmp("LSLCORREL", func))
4796 gdes->vf.op = VDEF_LSLCORREL;
4797 else {
4798 rrd_set_error
4799 ("Unknown function '%s' in VDEF '%s'\n", func, gdes->vname);
4800 return -1;
4801 };
4802 switch (gdes->vf.op) {
4803 case VDEF_PERCENT:
4804 case VDEF_PERCENTNAN:
4805 if (isnan(param)) { /* no parameter given */
4806 rrd_set_error
4807 ("Function '%s' needs parameter in VDEF '%s'\n",
4808 func, gdes->vname);
4809 return -1;
4810 };
4811 if (param >= 0.0 && param <= 100.0) {
4812 gdes->vf.param = param;
4813 gdes->vf.val = DNAN; /* undefined */
4814 gdes->vf.when = 0; /* undefined */
4815 } else {
4816 rrd_set_error
4817 ("Parameter '%f' out of range in VDEF '%s'\n",
4818 param, gdes->vname);
4819 return -1;
4820 };
4821 break;
4822 case VDEF_MAXIMUM:
4823 case VDEF_AVERAGE:
4824 case VDEF_STDEV:
4825 case VDEF_MINIMUM:
4826 case VDEF_TOTAL:
4827 case VDEF_FIRST:
4828 case VDEF_LAST:
4829 case VDEF_LSLSLOPE:
4830 case VDEF_LSLINT:
4831 case VDEF_LSLCORREL:
4832 if (isnan(param)) {
4833 gdes->vf.param = DNAN;
4834 gdes->vf.val = DNAN;
4835 gdes->vf.when = 0;
4836 } else {
4837 rrd_set_error
4838 ("Function '%s' needs no parameter in VDEF '%s'\n",
4839 func, gdes->vname);
4840 return -1;
4841 };
4842 break;
4843 };
4844 return 0;
4845 }
4848 int vdef_calc(
4849 image_desc_t *im,
4850 int gdi)
4851 {
4852 graph_desc_t *src, *dst;
4853 rrd_value_t *data;
4854 long step, steps;
4856 dst = &im->gdes[gdi];
4857 src = &im->gdes[dst->vidx];
4858 data = src->data + src->ds;
4860 steps = (src->end - src->start) / src->step;
4861 #if 0
4862 printf
4863 ("DEBUG: start == %lu, end == %lu, %lu steps\n",
4864 src->start, src->end, steps);
4865 #endif
4866 switch (dst->vf.op) {
4867 case VDEF_PERCENT:{
4868 rrd_value_t *array;
4869 int field;
4870 if ((array = (rrd_value_t*)malloc(steps * sizeof(double))) == NULL) {
4871 rrd_set_error("malloc VDEV_PERCENT");
4872 return -1;
4873 }
4874 for (step = 0; step < steps; step++) {
4875 array[step] = data[step * src->ds_cnt];
4876 }
4877 qsort(array, step, sizeof(double), vdef_percent_compar);
4878 field = round((dst->vf.param * (double)(steps - 1)) / 100.0);
4879 dst->vf.val = array[field];
4880 dst->vf.when = 0; /* no time component */
4881 free(array);
4882 #if 0
4883 for (step = 0; step < steps; step++)
4884 printf("DEBUG: %3li:%10.2f %c\n",
4885 step, array[step], step == field ? '*' : ' ');
4886 #endif
4887 }
4888 break;
4889 case VDEF_PERCENTNAN:{
4890 rrd_value_t *array;
4891 int field;
4892 /* count number of "valid" values */
4893 int nancount=0;
4894 for (step = 0; step < steps; step++) {
4895 if (!isnan(data[step * src->ds_cnt])) { nancount++; }
4896 }
4897 /* and allocate it */
4898 if ((array = (rrd_value_t*)malloc(nancount * sizeof(double))) == NULL) {
4899 rrd_set_error("malloc VDEV_PERCENT");
4900 return -1;
4901 }
4902 /* and fill it in */
4903 field=0;
4904 for (step = 0; step < steps; step++) {
4905 if (!isnan(data[step * src->ds_cnt])) {
4906 array[field] = data[step * src->ds_cnt];
4907 field++;
4908 }
4909 }
4910 qsort(array, nancount, sizeof(double), vdef_percent_compar);
4911 field = round( dst->vf.param * (double)(nancount - 1) / 100.0);
4912 dst->vf.val = array[field];
4913 dst->vf.when = 0; /* no time component */
4914 free(array);
4915 }
4916 break;
4917 case VDEF_MAXIMUM:
4918 step = 0;
4919 while (step != steps && isnan(data[step * src->ds_cnt]))
4920 step++;
4921 if (step == steps) {
4922 dst->vf.val = DNAN;
4923 dst->vf.when = 0;
4924 } else {
4925 dst->vf.val = data[step * src->ds_cnt];
4926 dst->vf.when = src->start + (step + 1) * src->step;
4927 }
4928 while (step != steps) {
4929 if (finite(data[step * src->ds_cnt])) {
4930 if (data[step * src->ds_cnt] > dst->vf.val) {
4931 dst->vf.val = data[step * src->ds_cnt];
4932 dst->vf.when = src->start + (step + 1) * src->step;
4933 }
4934 }
4935 step++;
4936 }
4937 break;
4938 case VDEF_TOTAL:
4939 case VDEF_STDEV:
4940 case VDEF_AVERAGE:{
4941 int cnt = 0;
4942 double sum = 0.0;
4943 double average = 0.0;
4945 for (step = 0; step < steps; step++) {
4946 if (finite(data[step * src->ds_cnt])) {
4947 sum += data[step * src->ds_cnt];
4948 cnt++;
4949 };
4950 }
4951 if (cnt) {
4952 if (dst->vf.op == VDEF_TOTAL) {
4953 dst->vf.val = sum * src->step;
4954 dst->vf.when = 0; /* no time component */
4955 } else if (dst->vf.op == VDEF_AVERAGE) {
4956 dst->vf.val = sum / cnt;
4957 dst->vf.when = 0; /* no time component */
4958 } else {
4959 average = sum / cnt;
4960 sum = 0.0;
4961 for (step = 0; step < steps; step++) {
4962 if (finite(data[step * src->ds_cnt])) {
4963 sum += pow((data[step * src->ds_cnt] - average), 2.0);
4964 };
4965 }
4966 dst->vf.val = pow(sum / cnt, 0.5);
4967 dst->vf.when = 0; /* no time component */
4968 };
4969 } else {
4970 dst->vf.val = DNAN;
4971 dst->vf.when = 0;
4972 }
4973 }
4974 break;
4975 case VDEF_MINIMUM:
4976 step = 0;
4977 while (step != steps && isnan(data[step * src->ds_cnt]))
4978 step++;
4979 if (step == steps) {
4980 dst->vf.val = DNAN;
4981 dst->vf.when = 0;
4982 } else {
4983 dst->vf.val = data[step * src->ds_cnt];
4984 dst->vf.when = src->start + (step + 1) * src->step;
4985 }
4986 while (step != steps) {
4987 if (finite(data[step * src->ds_cnt])) {
4988 if (data[step * src->ds_cnt] < dst->vf.val) {
4989 dst->vf.val = data[step * src->ds_cnt];
4990 dst->vf.when = src->start + (step + 1) * src->step;
4991 }
4992 }
4993 step++;
4994 }
4995 break;
4996 case VDEF_FIRST:
4997 /* The time value returned here is one step before the
4998 * actual time value. This is the start of the first
4999 * non-NaN interval.
5000 */
5001 step = 0;
5002 while (step != steps && isnan(data[step * src->ds_cnt]))
5003 step++;
5004 if (step == steps) { /* all entries were NaN */
5005 dst->vf.val = DNAN;
5006 dst->vf.when = 0;
5007 } else {
5008 dst->vf.val = data[step * src->ds_cnt];
5009 dst->vf.when = src->start + step * src->step;
5010 }
5011 break;
5012 case VDEF_LAST:
5013 /* The time value returned here is the
5014 * actual time value. This is the end of the last
5015 * non-NaN interval.
5016 */
5017 step = steps - 1;
5018 while (step >= 0 && isnan(data[step * src->ds_cnt]))
5019 step--;
5020 if (step < 0) { /* all entries were NaN */
5021 dst->vf.val = DNAN;
5022 dst->vf.when = 0;
5023 } else {
5024 dst->vf.val = data[step * src->ds_cnt];
5025 dst->vf.when = src->start + (step + 1) * src->step;
5026 }
5027 break;
5028 case VDEF_LSLSLOPE:
5029 case VDEF_LSLINT:
5030 case VDEF_LSLCORREL:{
5031 /* Bestfit line by linear least squares method */
5033 int cnt = 0;
5034 double SUMx, SUMy, SUMxy, SUMxx, SUMyy, slope, y_intercept, correl;
5036 SUMx = 0;
5037 SUMy = 0;
5038 SUMxy = 0;
5039 SUMxx = 0;
5040 SUMyy = 0;
5041 for (step = 0; step < steps; step++) {
5042 if (finite(data[step * src->ds_cnt])) {
5043 cnt++;
5044 SUMx += step;
5045 SUMxx += step * step;
5046 SUMxy += step * data[step * src->ds_cnt];
5047 SUMy += data[step * src->ds_cnt];
5048 SUMyy += data[step * src->ds_cnt] * data[step * src->ds_cnt];
5049 };
5050 }
5052 slope = (SUMx * SUMy - cnt * SUMxy) / (SUMx * SUMx - cnt * SUMxx);
5053 y_intercept = (SUMy - slope * SUMx) / cnt;
5054 correl =
5055 (SUMxy -
5056 (SUMx * SUMy) / cnt) /
5057 sqrt((SUMxx -
5058 (SUMx * SUMx) / cnt) * (SUMyy - (SUMy * SUMy) / cnt));
5059 if (cnt) {
5060 if (dst->vf.op == VDEF_LSLSLOPE) {
5061 dst->vf.val = slope;
5062 dst->vf.when = 0;
5063 } else if (dst->vf.op == VDEF_LSLINT) {
5064 dst->vf.val = y_intercept;
5065 dst->vf.when = 0;
5066 } else if (dst->vf.op == VDEF_LSLCORREL) {
5067 dst->vf.val = correl;
5068 dst->vf.when = 0;
5069 };
5070 } else {
5071 dst->vf.val = DNAN;
5072 dst->vf.when = 0;
5073 }
5074 }
5075 break;
5076 }
5077 return 0;
5078 }
5080 /* NaN < -INF < finite_values < INF */
5081 int vdef_percent_compar(
5082 const void
5083 *a,
5084 const void
5085 *b)
5086 {
5087 /* Equality is not returned; this doesn't hurt except
5088 * (maybe) for a little performance.
5089 */
5091 /* First catch NaN values. They are smallest */
5092 if (isnan(*(double *) a))
5093 return -1;
5094 if (isnan(*(double *) b))
5095 return 1;
5096 /* NaN doesn't reach this part so INF and -INF are extremes.
5097 * The sign from isinf() is compatible with the sign we return
5098 */
5099 if (isinf(*(double *) a))
5100 return isinf(*(double *) a);
5101 if (isinf(*(double *) b))
5102 return isinf(*(double *) b);
5103 /* If we reach this, both values must be finite */
5104 if (*(double *) a < *(double *) b)
5105 return -1;
5106 else
5107 return 1;
5108 }
5110 void grinfo_push(
5111 image_desc_t *im,
5112 char *key,
5113 rrd_info_type_t type,
5114 rrd_infoval_t value)
5115 {
5116 im->grinfo_current = rrd_info_push(im->grinfo_current, key, type, value);
5117 if (im->grinfo == NULL) {
5118 im->grinfo = im->grinfo_current;
5119 }
5120 }