2 * This file is part of the libsigrok project.
4 * Copyright (C) 2014 Bert Vermeulen <bert@biot.com>
6 * This program is free software: you can redistribute it and/or modify
7 * it under the terms of the GNU General Public License as published by
8 * the Free Software Foundation, either version 3 of the License, or
9 * (at your option) any later version.
11 * This program is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 * GNU General Public License for more details.
16 * You should have received a copy of the GNU General Public License
17 * along with this program. If not, see <http://www.gnu.org/licenses/>.
26 #include <libsigrok/libsigrok.h>
27 #include "libsigrok-internal.h"
30 #define LOG_PREFIX "analog"
36 * Handling and converting analog data.
40 * @defgroup grp_analog Analog data handling
42 * Handling and converting analog data.
47 struct unit_mq_string {
52 /* Please use the same order as in enum sr_unit (libsigrok.h). */
53 static struct unit_mq_string unit_strings[] = {
54 { SR_UNIT_VOLT, "V" },
55 { SR_UNIT_AMPERE, "A" },
56 { SR_UNIT_OHM, "\xe2\x84\xa6" },
57 { SR_UNIT_FARAD, "F" },
58 { SR_UNIT_KELVIN, "K" },
59 { SR_UNIT_CELSIUS, "\xc2\xb0""C" },
60 { SR_UNIT_FAHRENHEIT, "\xc2\xb0""F" },
61 { SR_UNIT_HERTZ, "Hz" },
62 { SR_UNIT_PERCENTAGE, "%" },
63 { SR_UNIT_BOOLEAN, "" },
64 { SR_UNIT_SECOND, "s" },
65 { SR_UNIT_SIEMENS, "S" },
66 { SR_UNIT_DECIBEL_MW, "dBu" },
67 { SR_UNIT_DECIBEL_VOLT, "dBv" },
68 { SR_UNIT_UNITLESS, "" },
69 { SR_UNIT_DECIBEL_SPL, "dB" },
70 { SR_UNIT_CONCENTRATION, "ppm" },
71 { SR_UNIT_REVOLUTIONS_PER_MINUTE, "RPM" },
72 { SR_UNIT_VOLT_AMPERE, "VA" },
73 { SR_UNIT_WATT, "W" },
74 { SR_UNIT_WATT_HOUR, "Wh" },
75 { SR_UNIT_METER_SECOND, "m/s" },
76 { SR_UNIT_HECTOPASCAL, "hPa" },
77 { SR_UNIT_HUMIDITY_293K, "%rF" },
78 { SR_UNIT_DEGREE, "\xc2\xb0" },
79 { SR_UNIT_HENRY, "H" },
80 { SR_UNIT_GRAM, "g" },
81 { SR_UNIT_CARAT, "ct" },
82 { SR_UNIT_OUNCE, "oz" },
83 { SR_UNIT_TROY_OUNCE, "oz t" },
84 { SR_UNIT_POUND, "lb" },
85 { SR_UNIT_PENNYWEIGHT, "dwt" },
86 { SR_UNIT_GRAIN, "gr" },
87 { SR_UNIT_TAEL, "tael" },
88 { SR_UNIT_MOMME, "momme" },
89 { SR_UNIT_TOLA, "tola" },
90 { SR_UNIT_PIECE, "pcs" },
94 /* Please use the same order as in enum sr_mqflag (libsigrok.h). */
95 static struct unit_mq_string mq_strings[] = {
96 { SR_MQFLAG_AC, " AC" },
97 { SR_MQFLAG_DC, " DC" },
98 { SR_MQFLAG_RMS, " RMS" },
99 { SR_MQFLAG_DIODE, " DIODE" },
100 { SR_MQFLAG_HOLD, " HOLD" },
101 { SR_MQFLAG_MAX, " MAX" },
102 { SR_MQFLAG_MIN, " MIN" },
103 { SR_MQFLAG_AUTORANGE, " AUTO" },
104 { SR_MQFLAG_RELATIVE, " REL" },
105 { SR_MQFLAG_SPL_FREQ_WEIGHT_A, "(A)" },
106 { SR_MQFLAG_SPL_FREQ_WEIGHT_C, "(C)" },
107 { SR_MQFLAG_SPL_FREQ_WEIGHT_Z, "(Z)" },
108 { SR_MQFLAG_SPL_FREQ_WEIGHT_FLAT, "(SPL)" },
109 { SR_MQFLAG_SPL_TIME_WEIGHT_S, " S" },
110 { SR_MQFLAG_SPL_TIME_WEIGHT_F, " F" },
111 { SR_MQFLAG_SPL_LAT, " LAT" },
112 /* Not a standard function for SLMs, so this is a made-up notation. */
113 { SR_MQFLAG_SPL_PCT_OVER_ALARM, "%oA" },
114 { SR_MQFLAG_DURATION, " DURATION" },
115 { SR_MQFLAG_AVG, " AVG" },
116 { SR_MQFLAG_REFERENCE, " REF" },
117 { SR_MQFLAG_UNSTABLE, " UNSTABLE" },
118 { SR_MQFLAG_FOUR_WIRE, " 4-WIRE" },
122 SR_PRIV int sr_analog_init(struct sr_datafeed_analog *analog,
123 struct sr_analog_encoding *encoding,
124 struct sr_analog_meaning *meaning,
125 struct sr_analog_spec *spec,
128 memset(analog, 0, sizeof(*analog));
129 memset(encoding, 0, sizeof(*encoding));
130 memset(meaning, 0, sizeof(*meaning));
131 memset(spec, 0, sizeof(*spec));
133 analog->encoding = encoding;
134 analog->meaning = meaning;
137 encoding->unitsize = sizeof(float);
138 encoding->is_float = TRUE;
139 #ifdef WORDS_BIGENDIAN
140 encoding->is_bigendian = TRUE;
142 encoding->is_bigendian = FALSE;
144 encoding->digits = digits;
145 encoding->is_digits_decimal = TRUE;
146 encoding->scale.p = 1;
147 encoding->scale.q = 1;
148 encoding->offset.p = 0;
149 encoding->offset.q = 1;
151 spec->spec_digits = digits;
157 * Convert an analog datafeed payload to an array of floats.
159 * @param[in] analog The analog payload to convert. Must not be NULL.
160 * analog->data, analog->meaning, and analog->encoding
162 * @param[out] outbuf Memory where to store the result. Must not be NULL.
164 * Sufficient memory for outbuf must have been pre-allocated by the caller,
165 * who is also responsible for freeing it when no longer needed.
167 * @retval SR_OK Success.
168 * @retval SR_ERR Unsupported encoding.
169 * @retval SR_ERR_ARG Invalid argument.
173 SR_API int sr_analog_to_float(const struct sr_datafeed_analog *analog,
177 unsigned int b, i, count;
180 if (!analog || !(analog->data) || !(analog->meaning)
181 || !(analog->encoding) || !outbuf)
184 count = analog->num_samples * g_slist_length(analog->meaning->channels);
186 #ifdef WORDS_BIGENDIAN
191 if (!analog->encoding->is_float) {
192 float offset = analog->encoding->offset.p / (float)analog->encoding->offset.q;
193 float scale = analog->encoding->scale.p / (float)analog->encoding->scale.q;
194 gboolean is_signed = analog->encoding->is_signed;
195 gboolean is_bigendian = analog->encoding->is_bigendian;
196 int8_t *data8 = (int8_t *)(analog->data);
197 int16_t *data16 = (int16_t *)(analog->data);
198 int32_t *data32 = (int32_t *)(analog->data);
200 switch (analog->encoding->unitsize) {
203 for (unsigned int i = 0; i < count; i++) {
204 outbuf[i] = scale * data8[i];
208 for (unsigned int i = 0; i < count; i++) {
209 outbuf[i] = scale * R8(data8 + i);
215 if (is_signed && is_bigendian) {
216 for (unsigned int i = 0; i < count; i++) {
217 outbuf[i] = scale * RB16S(&data16[i]);
220 } else if (is_bigendian) {
221 for (unsigned int i = 0; i < count; i++) {
222 outbuf[i] = scale * RB16(&data16[i]);
225 } else if (is_signed) {
226 for (unsigned int i = 0; i < count; i++) {
227 outbuf[i] = scale * RL16S(&data16[i]);
231 for (unsigned int i = 0; i < count; i++) {
232 outbuf[i] = scale * RL16(&data16[i]);
238 if (is_signed && is_bigendian) {
239 for (unsigned int i = 0; i < count; i++) {
240 outbuf[i] = scale * RB32S(&data32[i]);
243 } else if (is_bigendian) {
244 for (unsigned int i = 0; i < count; i++) {
245 outbuf[i] = scale * RB32(&data32[i]);
248 } else if (is_signed) {
249 for (unsigned int i = 0; i < count; i++) {
250 outbuf[i] = scale * RL32S(&data32[i]);
254 for (unsigned int i = 0; i < count; i++) {
255 outbuf[i] = scale * RL32(&data32[i]);
261 sr_err("Unsupported unit size '%d' for analog-to-float conversion.",
262 analog->encoding->unitsize);
268 if (analog->encoding->unitsize == sizeof(float)
269 && analog->encoding->is_bigendian == bigendian
270 && analog->encoding->scale.p == 1
271 && analog->encoding->scale.q == 1
272 && analog->encoding->offset.p / (float)analog->encoding->offset.q == 0) {
273 /* The data is already in the right format. */
274 memcpy(outbuf, analog->data, count * sizeof(float));
276 for (i = 0; i < count; i += analog->encoding->unitsize) {
277 for (b = 0; b < analog->encoding->unitsize; b++) {
278 if (analog->encoding->is_bigendian == bigendian)
279 ((uint8_t *)outbuf)[i + b] =
280 ((uint8_t *)analog->data)[i * analog->encoding->unitsize + b];
282 ((uint8_t *)outbuf)[i + (analog->encoding->unitsize - b)] =
283 ((uint8_t *)analog->data)[i * analog->encoding->unitsize + b];
285 if (analog->encoding->scale.p != 1
286 || analog->encoding->scale.q != 1)
287 outbuf[i] = (outbuf[i] * analog->encoding->scale.p) / analog->encoding->scale.q;
288 offset = ((float)analog->encoding->offset.p / (float)analog->encoding->offset.q);
297 * Scale a float value to the appropriate SI prefix.
299 * @param[in,out] value The float value to convert to appropriate SI prefix.
300 * @param[in,out] digits The number of significant decimal digits in value.
302 * @return The SI prefix to which value was scaled, as a printable string.
306 SR_API const char *sr_analog_si_prefix(float *value, int *digits)
308 #define NEG_PREFIX_COUNT 5 /* number of prefixes below unity */
309 #define POS_PREFIX_COUNT (int)(ARRAY_SIZE(prefixes) - NEG_PREFIX_COUNT - 1)
310 static const char *prefixes[] = { "f","p","n","µ","m","","k","M","G","T" };
312 if (value == NULL || digits == NULL || isnan(*value))
313 return prefixes[NEG_PREFIX_COUNT];
315 float logval = log10f(fabsf(*value));
316 int prefix = (logval / 3) - (logval < 1);
318 if (prefix < -NEG_PREFIX_COUNT) prefix = -NEG_PREFIX_COUNT;
319 if (3 * prefix < -*digits) prefix = (-*digits + 2 * (*digits < 0)) / 3;
320 if (prefix > POS_PREFIX_COUNT) prefix = POS_PREFIX_COUNT;
322 *value *= powf(10, -3 * prefix);
323 *digits += 3 * prefix;
324 return prefixes[prefix + NEG_PREFIX_COUNT];
328 * Convert the unit/MQ/MQ flags in the analog struct to a string.
330 * @param[in] analog Struct containing the unit, MQ and MQ flags.
331 * Must not be NULL. analog->meaning must not be NULL.
332 * @param[out] result Pointer to store result. Must not be NULL.
334 * The string is allocated by the function and must be freed by the caller
335 * after use by calling g_free().
337 * @retval SR_OK Success.
338 * @retval SR_ERR_ARG Invalid argument.
342 SR_API int sr_analog_unit_to_string(const struct sr_datafeed_analog *analog,
348 if (!analog || !(analog->meaning) || !result)
351 buf = g_string_new(NULL);
353 for (i = 0; unit_strings[i].value; i++) {
354 if (analog->meaning->unit == unit_strings[i].value) {
355 g_string_assign(buf, unit_strings[i].str);
360 /* More than one MQ flag may apply. */
361 for (i = 0; mq_strings[i].value; i++)
362 if (analog->meaning->mqflags & mq_strings[i].value)
363 g_string_append(buf, mq_strings[i].str);
366 g_string_free(buf, FALSE);
372 * Set sr_rational r to the given value.
374 * @param[out] r Rational number struct to set. Must not be NULL.
375 * @param[in] p Numerator.
376 * @param[in] q Denominator.
380 SR_API void sr_rational_set(struct sr_rational *r, int64_t p, uint64_t q)
389 #ifndef HAVE___INT128_T
395 struct sr_uint128_t {
400 static void mult_int64(struct sr_int128_t *res, const int64_t a,
403 uint64_t t1, t2, t3, t4;
405 t1 = (UINT32_MAX & a) * (UINT32_MAX & b);
406 t2 = (UINT32_MAX & a) * (b >> 32);
407 t3 = (a >> 32) * (UINT32_MAX & b);
408 t4 = (a >> 32) * (b >> 32);
410 res->low = t1 + (t2 << 32) + (t3 << 32);
411 res->high = (t1 >> 32) + (uint64_t)((uint32_t)(t2)) + (uint64_t)((uint32_t)(t3));
413 res->high += ((int64_t)t2 >> 32) + ((int64_t)t3 >> 32) + t4;
416 static void mult_uint64(struct sr_uint128_t *res, const uint64_t a,
419 uint64_t t1, t2, t3, t4;
421 // (x1 + x2) * (y1 + y2) = x1*y1 + x1*y2 + x2*y1 + x2*y2
422 t1 = (UINT32_MAX & a) * (UINT32_MAX & b);
423 t2 = (UINT32_MAX & a) * (b >> 32);
424 t3 = (a >> 32) * (UINT32_MAX & b);
425 t4 = (a >> 32) * (b >> 32);
427 res->low = t1 + (t2 << 32) + (t3 << 32);
428 res->high = (t1 >> 32) + (uint64_t)((uint32_t)(t2)) + (uint64_t)((uint32_t)(t3));
430 res->high += ((int64_t)t2 >> 32) + ((int64_t)t3 >> 32) + t4;
435 * Compare two sr_rational for equality
437 * @param[in] a First value
438 * @param[in] b Second value
440 * The values are compared for numerical equality, i.e. 2/10 == 1/5
442 * @retval 1 if both values are equal
443 * @retval 0 otherwise
447 SR_API int sr_rational_eq(const struct sr_rational *a, const struct sr_rational *b)
449 #ifdef HAVE___INT128_T
452 /* p1/q1 = p2/q2 <=> p1*q2 = p2*q1 */
453 m1 = ((__int128_t)(b->p)) * ((__uint128_t)a->q);
454 m2 = ((__int128_t)(a->p)) * ((__uint128_t)b->q);
459 struct sr_int128_t m1, m2;
461 mult_int64(&m1, a->q, b->p);
462 mult_int64(&m2, a->p, b->q);
464 return (m1.high == m2.high) && (m1.low == m2.low);
469 * Multiply two sr_rational
471 * @param[in] a First value
472 * @param[in] b Second value
473 * @param[out] res Result
475 * The resulting nominator/denominator are reduced if the result would not fit
476 * otherwise. If the resulting nominator/denominator are relatively prime,
477 * this may not be possible.
479 * It is save to use the same variable for result and input values
481 * @retval SR_OK Success.
482 * @retval SR_ERR_ARG Resulting value to large
486 SR_API int sr_rational_mult(struct sr_rational *res, const struct sr_rational *a,
487 const struct sr_rational *b)
489 #ifdef HAVE___INT128_T
493 p = (__int128_t)(a->p) * (__int128_t)(b->p);
494 q = (__uint128_t)(a->q) * (__uint128_t)(b->q);
496 if ((p > INT64_MAX) || (p < INT64_MIN) || (q > UINT64_MAX)) {
497 while (!((p & 1) || (q & 1))) {
503 if ((p > INT64_MAX) || (p < INT64_MIN) || (q > UINT64_MAX)) {
504 // TODO: determine gcd to do further reduction
508 res->p = (int64_t)(p);
509 res->q = (uint64_t)(q);
514 struct sr_int128_t p;
515 struct sr_uint128_t q;
517 mult_int64(&p, a->p, b->p);
518 mult_uint64(&q, a->q, b->q);
520 while (!(p.low & 1) && !(q.low & 1)) {
522 if (p.high & 1) p.low |= (1ll << 63);
525 if (q.high & 1) q.low |= (1ll << 63);
531 if ((p.high >= 0) && (p.low > INT64_MAX))
536 res->p = (int64_t)p.low;
544 * Divide rational a by rational b
546 * @param[in] num numerator
547 * @param[in] div divisor
548 * @param[out] res Result
550 * The resulting nominator/denominator are reduced if the result would not fit
551 * otherwise. If the resulting nominator/denominator are relatively prime,
552 * this may not be possible.
554 * It is save to use the same variable for result and input values
556 * @retval SR_OK Success.
557 * @retval SR_ERR_ARG Division by zero
558 * @retval SR_ERR_ARG Denominator of divisor to large
559 * @retval SR_ERR_ARG Resulting value to large
563 SR_API int sr_rational_div(struct sr_rational *res, const struct sr_rational *num,
564 const struct sr_rational *div)
566 struct sr_rational t;
568 if (div->q > INT64_MAX)
581 return sr_rational_mult(res, num, &t);