median.c 2.2 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879
  1. // SPDX-License-Identifier: GPL-3.0-or-later
  2. #include "median.h"
  3. // ----------------------------------------------------------------------------
  4. // median
  5. struct grouping_median {
  6. size_t series_size;
  7. size_t next_pos;
  8. LONG_DOUBLE series[];
  9. };
  10. void *grouping_create_median(RRDR *r) {
  11. long entries = r->group;
  12. if(entries < 0) entries = 0;
  13. struct grouping_median *g = (struct grouping_median *)callocz(1, sizeof(struct grouping_median) + entries * sizeof(LONG_DOUBLE));
  14. g->series_size = (size_t)entries;
  15. return g;
  16. }
  17. // resets when switches dimensions
  18. // so, clear everything to restart
  19. void grouping_reset_median(RRDR *r) {
  20. struct grouping_median *g = (struct grouping_median *)r->internal.grouping_data;
  21. g->next_pos = 0;
  22. }
  23. void grouping_free_median(RRDR *r) {
  24. freez(r->internal.grouping_data);
  25. r->internal.grouping_data = NULL;
  26. }
  27. void grouping_add_median(RRDR *r, calculated_number value) {
  28. struct grouping_median *g = (struct grouping_median *)r->internal.grouping_data;
  29. if(unlikely(g->next_pos >= g->series_size)) {
  30. error("INTERNAL ERROR: median buffer overflow on chart '%s' - next_pos = %zu, series_size = %zu, r->group = %ld.", r->st->name, g->next_pos, g->series_size, r->group);
  31. }
  32. else {
  33. if(calculated_number_isnumber(value))
  34. g->series[g->next_pos++] = (LONG_DOUBLE)value;
  35. }
  36. }
  37. calculated_number grouping_flush_median(RRDR *r, RRDR_VALUE_FLAGS *rrdr_value_options_ptr) {
  38. struct grouping_median *g = (struct grouping_median *)r->internal.grouping_data;
  39. calculated_number value;
  40. if(unlikely(!g->next_pos)) {
  41. value = 0.0;
  42. *rrdr_value_options_ptr |= RRDR_VALUE_EMPTY;
  43. }
  44. else {
  45. if(g->next_pos > 1) {
  46. sort_series(g->series, g->next_pos);
  47. value = (calculated_number)median_on_sorted_series(g->series, g->next_pos);
  48. }
  49. else
  50. value = (calculated_number)g->series[0];
  51. if(!calculated_number_isnumber(value)) {
  52. value = 0.0;
  53. *rrdr_value_options_ptr |= RRDR_VALUE_EMPTY;
  54. }
  55. //log_series_to_stderr(g->series, g->next_pos, value, "median");
  56. }
  57. g->next_pos = 0;
  58. return value;
  59. }