diff options
author | Daniel Baumann <daniel.baumann@progress-linux.org> | 2024-05-05 12:08:03 +0000 |
---|---|---|
committer | Daniel Baumann <daniel.baumann@progress-linux.org> | 2024-05-05 12:08:18 +0000 |
commit | 5da14042f70711ea5cf66e034699730335462f66 (patch) | |
tree | 0f6354ccac934ed87a2d555f45be4c831cf92f4a /libnetdata/storage_number/storage_number.c | |
parent | Releasing debian version 1.44.3-2. (diff) | |
download | netdata-5da14042f70711ea5cf66e034699730335462f66.tar.xz netdata-5da14042f70711ea5cf66e034699730335462f66.zip |
Merging upstream version 1.45.3+dfsg.
Signed-off-by: Daniel Baumann <daniel.baumann@progress-linux.org>
Diffstat (limited to 'libnetdata/storage_number/storage_number.c')
-rw-r--r-- | libnetdata/storage_number/storage_number.c | 233 |
1 files changed, 0 insertions, 233 deletions
diff --git a/libnetdata/storage_number/storage_number.c b/libnetdata/storage_number/storage_number.c deleted file mode 100644 index 6468951bd..000000000 --- a/libnetdata/storage_number/storage_number.c +++ /dev/null @@ -1,233 +0,0 @@ -// SPDX-License-Identifier: GPL-3.0-or-later - -#include "../libnetdata.h" - -bool is_system_ieee754_double(void) { - static bool logged = false; - - struct { - NETDATA_DOUBLE original; - - union { - uint64_t i; - NETDATA_DOUBLE d; - }; - } tests[] = { - { .original = 1.25, .i = 0x3FF4000000000000 }, - { .original = 1.0, .i = 0x3FF0000000000000 }, - { .original = 2.0, .i = 0x4000000000000000 }, - { .original = 4.0, .i = 0x4010000000000000 }, - { .original = 8.8, .i = 0x402199999999999A }, - { .original = 16.16, .i = 0x403028F5C28F5C29 }, - { .original = 32.32, .i = 0x404028F5C28F5C29 }, - { .original = 64.64, .i = 0x405028F5C28F5C29 }, - { .original = 128.128, .i = 0x406004189374BC6A }, - { .original = 32768.32768, .i = 0x40E0000A7C5AC472 }, - { .original = 65536.65536, .i = 0x40F0000A7C5AC472 }, - { .original = -65536.65536, .i = 0xC0F0000A7C5AC472 }, - { .original = 65535.65535, .i = 0x40EFFFF4F8A0902E }, - { .original = -65535.65535, .i = 0xC0EFFFF4F8A0902E }, - { .original = 4.503599627e15, .i = 0x432FFFFFFFF4B180 }, - { .original = -4.503599627e15, .i = 0xC32FFFFFFFF4B180 }, - { .original = 1.25e25, .i = 0x4524ADF4B7320335 }, - { .original = 1.25e307, .i = 0x7FB1CCF385EBC8A0 }, - { .original = 1.25e-25, .i = 0x3AC357C299A88EA7 }, - { .original = 1.25e-100, .i = 0x2B317F7D4ED8C33E }, - { .original = NAN, .i = 0x7FF8000000000000 }, - { .original = -INFINITY, .i = 0xFFF0000000000000 }, - { .original = INFINITY, .i = 0x7FF0000000000000 }, - { .original = 1.25e-132, .i = 0x248C6463225AB7EC }, - { .original = 0.0, .i = 0x0000000000000000 }, - { .original = -0.0, .i = 0x8000000000000000 }, - { .original = DBL_MIN, .i = 0x0010000000000000 }, - { .original = DBL_MAX, .i = 0x7FEFFFFFFFFFFFFF }, - { .original = -DBL_MIN, .i = 0x8010000000000000 }, - { .original = -DBL_MAX, .i = 0xFFEFFFFFFFFFFFFF }, - }; - - size_t errors = 0; - size_t elements = sizeof(tests) / sizeof(tests[0]); - for(size_t i = 0; i < elements ; i++) { - uint64_t *ptr = (uint64_t *)&tests[i].original; - - if(*ptr != tests[i].i && (tests[i].original == tests[i].d || (isnan(tests[i].original) && isnan(tests[i].d)))) { - if(!logged) - netdata_log_info("IEEE754: test #%zu, value " NETDATA_DOUBLE_FORMAT_G " is represented in this system as %lX, but it was expected as %lX", - i+1, tests[i].original, *ptr, tests[i].i); - errors++; - } - } - - if(!errors && sizeof(NETDATA_DOUBLE) == sizeof(uint64_t)) { - if(!logged) - netdata_log_info("IEEE754: system is using IEEE754 DOUBLE PRECISION values"); - - logged = true; - return true; - } - else { - if(!logged) - netdata_log_info("IEEE754: system is NOT compatible with IEEE754 DOUBLE PRECISION values"); - - logged = true; - return false; - } -} - -storage_number pack_storage_number(NETDATA_DOUBLE value, SN_FLAGS flags) { - // bit 32 = sign 0:positive, 1:negative - // bit 31 = 0:divide, 1:multiply - // bit 30, 29, 28 = (multiplier or divider) 0-7 (8 total) - // bit 27 SN_EXISTS_100 - // bit 26 SN_EXISTS_RESET - // bit 25 SN_ANOMALY_BIT = 0: anomalous, 1: not anomalous - // bit 24 to bit 1 = the value - - if(unlikely(fpclassify(value) == FP_NAN || fpclassify(value) == FP_INFINITE)) - return SN_EMPTY_SLOT; - - storage_number r = flags & SN_USER_FLAGS; - - if(unlikely(fpclassify(value) == FP_ZERO || fpclassify(value) == FP_SUBNORMAL)) - return r; - - int m = 0; - NETDATA_DOUBLE n = value, factor = 10; - - // if the value is negative - // add the sign bit and make it positive - if(n < 0) { - r += SN_FLAG_NEGATIVE; // the sign bit 32 - n = -n; - } - - if(n / 10000000.0 > 0x00ffffff) { - factor = 100; - r |= SN_FLAG_NOT_EXISTS_MUL100; - } - - // make its integer part fit in 0x00ffffff - // by dividing it by 10 up to 7 times - // and increasing the multiplier - while(m < 7 && n > (NETDATA_DOUBLE)0x00ffffff) { - n /= factor; - m++; - } - - if(m) { - // the value was too big, and we divided it - // so, we add a multiplier to unpack it - r += SN_FLAG_MULTIPLY + (m << 27); // the multiplier m - - if(n > (NETDATA_DOUBLE)0x00ffffff) { - #ifdef NETDATA_INTERNAL_CHECKS - netdata_log_error("Number " NETDATA_DOUBLE_FORMAT " is too big.", value); - #endif - r += 0x00ffffff; - return r; - } - } - else { - // 0x0019999e is the number that can be multiplied - // by 10 to give 0x00ffffff - // while the value is below 0x0019999e we can - // multiply it by 10, up to 7 times, increasing - // the multiplier - while(m < 7 && n < (NETDATA_DOUBLE)0x0019999e) { - n *= 10; - m++; - } - - if (unlikely(n > (NETDATA_DOUBLE)0x00ffffff)) { - n /= 10; - m--; - } - // the value was small enough, and we multiplied it - // so, we add a divider to unpack it - r += (m << 27); // the divider m - } - -#ifdef STORAGE_WITH_MATH - // without this there are rounding problems - // example: 0.9 becomes 0.89 - r += lrint((double) n); -#else - r += (storage_number)n; -#endif - - return r; -} - -// Lookup table to make storage number unpacking efficient. -NETDATA_DOUBLE unpack_storage_number_lut10x[4 * 8]; - -__attribute__((constructor)) void initialize_lut(void) { - // The lookup table is partitioned in 4 subtables based on the - // values of the factor and exp bits. - for (int i = 0; i < 8; i++) { - // factor = 0 - unpack_storage_number_lut10x[0 * 8 + i] = 1 / pow(10, i); // exp = 0 - unpack_storage_number_lut10x[1 * 8 + i] = pow(10, i); // exp = 1 - - // factor = 1 - unpack_storage_number_lut10x[2 * 8 + i] = 1 / pow(100, i); // exp = 0 - unpack_storage_number_lut10x[3 * 8 + i] = pow(100, i); // exp = 1 - } -} - -/* -int print_netdata_double(char *str, NETDATA_DOUBLE value) -{ - char *wstr = str; - - int sign = (value < 0) ? 1 : 0; - if(sign) value = -value; - -#ifdef STORAGE_WITH_MATH - // without llrintl() there are rounding problems - // for example 0.9 becomes 0.89 - unsigned long long uvalue = (unsigned long long int) llrintl(value * (NETDATA_DOUBLE)100000); -#else - unsigned long long uvalue = value * (NETDATA_DOUBLE)100000; -#endif - - wstr = print_number_llu_r_smart(str, uvalue); - - // make sure we have 6 bytes at least - while((wstr - str) < 6) *wstr++ = '0'; - - // put the sign back - if(sign) *wstr++ = '-'; - - // reverse it - char *begin = str, *end = --wstr, aux; - while (end > begin) aux = *end, *end-- = *begin, *begin++ = aux; - // wstr--; - // strreverse(str, wstr); - - // remove trailing zeros - int decimal = 5; - while(decimal > 0 && *wstr == '0') { - *wstr-- = '\0'; - decimal--; - } - - // terminate it, one position to the right - // to let space for a dot - wstr[2] = '\0'; - - // make space for the dot - int i; - for(i = 0; i < decimal ;i++) { - wstr[1] = wstr[0]; - wstr--; - } - - // put the dot - if(wstr[2] == '\0') { wstr[1] = '\0'; decimal--; } - else wstr[1] = '.'; - - // return the buffer length - return (int) ((wstr - str) + 2 + decimal ); -} -*/ |