[libsigrok.git] / src / analog.c

/*
 * This file is part of the libsigrok project.
 *
 * Copyright (C) 2014 Bert Vermeulen <bert@biot.com>
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */

#include <stdio.h>
#include <stdint.h>
#include <string.h>
#include <ctype.h>
#include <libsigrok/libsigrok.h>
#include "libsigrok-internal.h"

/** @cond PRIVATE */
#define LOG_PREFIX "analog"
/** @endcond */

/**
 * @file
 *
 * Handling and converting analog data.
 */

/**
 * @defgroup grp_analog Analog data handling
 *
 * Handling and converting analog data.
 *
 * @{
 */

struct unit_mq_string {
	uint64_t value;
	char *str;
};

static struct unit_mq_string unit_strings[] = {
	{ SR_UNIT_VOLT, "V" },
	{ SR_UNIT_AMPERE, "A" },
	{ SR_UNIT_OHM, "\xe2\x84\xa6" },
	{ SR_UNIT_FARAD, "F" },
	{ SR_UNIT_KELVIN, "K" },
	{ SR_UNIT_CELSIUS, "\xc2\xb0""C" },
	{ SR_UNIT_FAHRENHEIT, "\xc2\xb0""F" },
	{ SR_UNIT_HERTZ, "Hz" },
	{ SR_UNIT_PERCENTAGE, "%" },
	{ SR_UNIT_BOOLEAN, "" },
	{ SR_UNIT_SECOND, "s" },
	{ SR_UNIT_SIEMENS, "S" },
	{ SR_UNIT_DECIBEL_MW, "dBu" },
	{ SR_UNIT_DECIBEL_VOLT, "dBv" },
	{ SR_UNIT_UNITLESS, "" },
	{ SR_UNIT_DECIBEL_SPL, "dB" },
	{ SR_UNIT_CONCENTRATION, "ppm" },
	{ SR_UNIT_REVOLUTIONS_PER_MINUTE, "RPM" },
	{ SR_UNIT_VOLT_AMPERE, "VA" },
	{ SR_UNIT_WATT, "W" },
	{ SR_UNIT_WATT_HOUR, "Wh" },
	{ SR_UNIT_METER_SECOND, "m/s" },
	{ SR_UNIT_HECTOPASCAL, "hPa" },
	{ SR_UNIT_HUMIDITY_293K, "%rF" },
	{ SR_UNIT_DEGREE, "\xc2\xb0" },
	{ SR_UNIT_HENRY, "H" },
	{ SR_UNIT_GRAM, "g" },
	{ SR_UNIT_CARAT, "ct" },
	{ SR_UNIT_OUNCE, "oz" },
	{ SR_UNIT_TROY_OUNCE, "oz t" },
	{ SR_UNIT_POUND, "lb" },
	{ SR_UNIT_PENNYWEIGHT, "dwt" },
	{ SR_UNIT_GRAIN, "gr" },
	{ SR_UNIT_TAEL, "tael" },
	{ SR_UNIT_MOMME, "momme" },
	{ SR_UNIT_TOLA, "tola" },
	{ SR_UNIT_PIECE, "pcs" },
	ALL_ZERO
};

static struct unit_mq_string mq_strings[] = {
	{ SR_MQFLAG_AC, " AC" },
	{ SR_MQFLAG_DC, " DC" },
	{ SR_MQFLAG_RMS, " RMS" },
	{ SR_MQFLAG_DIODE, " DIODE" },
	{ SR_MQFLAG_HOLD, " HOLD" },
	{ SR_MQFLAG_MAX, " MAX" },
	{ SR_MQFLAG_MIN, " MIN" },
	{ SR_MQFLAG_AUTORANGE, " AUTO" },
	{ SR_MQFLAG_RELATIVE, " REL" },
	{ SR_MQFLAG_SPL_FREQ_WEIGHT_A, "(A)" },
	{ SR_MQFLAG_SPL_FREQ_WEIGHT_C, "(C)" },
	{ SR_MQFLAG_SPL_FREQ_WEIGHT_Z, "(Z)" },
	{ SR_MQFLAG_SPL_FREQ_WEIGHT_FLAT, "(SPL)" },
	{ SR_MQFLAG_SPL_TIME_WEIGHT_S, " S" },
	{ SR_MQFLAG_SPL_TIME_WEIGHT_F, " F" },
	{ SR_MQFLAG_SPL_LAT, " LAT" },
	/* Not a standard function for SLMs, so this is a made-up notation. */
	{ SR_MQFLAG_SPL_PCT_OVER_ALARM, "%oA" },
	{ SR_MQFLAG_DURATION, " DURATION" },
	{ SR_MQFLAG_AVG, " AVG" },
	{ SR_MQFLAG_REFERENCE, " REF" },
	{ SR_MQFLAG_UNSTABLE, " UNSTABLE" },
	ALL_ZERO
};

SR_PRIV int sr_analog_init(struct sr_datafeed_analog2 *analog,
		struct sr_analog_encoding *encoding,
		struct sr_analog_meaning *meaning,
		struct sr_analog_spec *spec,
		int digits)
{
	memset(analog, 0, sizeof(*analog));
	memset(encoding, 0, sizeof(*encoding));
	memset(meaning, 0, sizeof(*meaning));
	memset(spec, 0, sizeof(*spec));

	analog->encoding = encoding;
	analog->meaning = meaning;
	analog->spec = spec;

	encoding->unitsize = sizeof(float);
	encoding->is_float = TRUE;
#ifdef WORDS_BIGENDIAN
	encoding->is_bigendian = TRUE;
#else
	encoding->is_bigendian = FALSE;
#endif
	encoding->digits = digits;
	encoding->is_digits_decimal = TRUE;
	encoding->scale.p = 1;
	encoding->scale.q = 1;
	encoding->offset.p = 0;
	encoding->offset.q = 1;

	spec->spec_digits = digits;

	return SR_OK;
}

SR_API int sr_analog_to_float(const struct sr_datafeed_analog2 *analog,
		float *outbuf)
{
	float offset;
	unsigned int b, i;
	gboolean bigendian;

#ifdef WORDS_BIGENDIAN
	bigendian = TRUE;
#else
	bigendian = FALSE;
#endif
	if (!analog->encoding->is_float) {
		/* TODO */
		sr_err("Only floating-point encoding supported so far.");
		return SR_ERR;
	}

	if (analog->encoding->unitsize == sizeof(float)
			&& analog->encoding->is_bigendian == bigendian
			&& (analog->encoding->scale.p == analog->encoding->scale.q)
			&& analog->encoding->offset.p / (float)analog->encoding->offset.q == 0) {
		/* The data is already in the right format. */
		memcpy(outbuf, analog->data, analog->num_samples * sizeof(float));
	} else {
		for (i = 0; i < analog->num_samples; i += analog->encoding->unitsize) {
			for (b = 0; b < analog->encoding->unitsize; b++) {
				if (analog->encoding->is_bigendian == bigendian)
					outbuf[i + b] = ((float *)analog->data)[i * analog->encoding->unitsize + b];
				else
					outbuf[i + (analog->encoding->unitsize - b)] = ((float *)analog->data)[i * analog->encoding->unitsize + b];
			}
			if (analog->encoding->scale.p != analog->encoding->scale.q)
				outbuf[i] = (outbuf[i] * analog->encoding->scale.p) / analog->encoding->scale.q;
			offset = ((float)analog->encoding->offset.p / (float)analog->encoding->offset.q);
			outbuf[i] += offset;
		}
	}

	return SR_OK;
}

/*
 * Convert a floating point value to a string, limited to the given
 * number of decimal digits.
 *
 * @param value The value to convert.
 * @param digits Number of digits after the decimal point to print.
 * @param result Pointer to store result.
 *
 * The string is allocated by the function and must be freed by the caller
 * after use by calling g_free().
 *
 * @retval SR_OK
 *
 * @since 0.4.0
 */
SR_API int sr_analog_float_to_string(float value, int digits, char **result)
{
	int cnt, i;

	/* This produces at least one too many digits */
	*result = g_strdup_printf("%.*f", digits, value);
	for (i = 0, cnt = 0; (*result)[i]; i++) {
		if (isdigit((*result)[i++]))
			cnt++;
		if (cnt == digits) {
			(*result)[i] = 0;
			break;
		}
	}

	return SR_OK;
}

/*
 * Convert the unit/MQ/MQ flags in the analog struct to a string.
 *
 * @param analog Struct containing the unit, MQ and MQ flags.
 * @param result Pointer to store result.
 *
 * The string is allocated by the function and must be freed by the caller
 * after use by calling g_free().
 *
 * @retval SR_OK
 *
 * @since 0.4.0
 */
SR_API int sr_analog_unit_to_string(const struct sr_datafeed_analog2 *analog,
		char **result)
{
	int i;
	GString *buf = g_string_new(NULL);

	for (i = 0; unit_strings[i].value; i++) {
		if (analog->meaning->unit == unit_strings[i].value) {
			g_string_assign(buf, unit_strings[i].str);
			break;
		}
	}

	/* More than one MQ flag may apply. */
	for (i = 0; mq_strings[i].value; i++)
		if (analog->meaning->mqflags & mq_strings[i].value)
			g_string_append(buf, mq_strings[i].str);

	*result = buf->str;
	g_string_free(buf, FALSE);

	return SR_OK;
}

/*
 * Set sr_rational r to the given value.
 *
 * @param p Numerator
 * @param q Denominator
 */
SR_API void sr_rational_set(struct sr_rational *r, uint64_t p, uint64_t q)
{
	r->p = p;
	r->q = q;
}

/** @} */
Commit	Line	Data
fb019a0e BV	1	/*
	2	* This file is part of the libsigrok project.
	3	*
	4	* Copyright (C) 2014 Bert Vermeulen <bert@biot.com>
	5	*
	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.
	10	*
	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.
	15	*
	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/>.
	18	*/
	19
c2a25ebb BV	20	#include <stdio.h>
c2a25ebb BV	21	#include <stdint.h>
fb019a0e	22	#include <string.h>
c2a25ebb	23	#include <ctype.h>
c1aae900	24	#include <libsigrok/libsigrok.h>
fb019a0e BV	25	#include "libsigrok-internal.h"
fb019a0e BV	26
e00b3f58	27	/** @cond PRIVATE */
fb019a0e	28	#define LOG_PREFIX "analog"
e00b3f58 UH	29	/** @endcond */
	30
	31	/**
	32	* @file
	33	*
	34	* Handling and converting analog data.
	35	*/
	36
	37	/**
	38	* @defgroup grp_analog Analog data handling
	39	*
	40	* Handling and converting analog data.
	41	*
	42	* @{
	43	*/
fb019a0e	44
a5892391 BV	45	struct unit_mq_string {
	46	uint64_t value;
	47	char *str;
	48	};
	49
	50	static struct unit_mq_string unit_strings[] = {
	51	{ SR_UNIT_VOLT, "V" },
	52	{ SR_UNIT_AMPERE, "A" },
	53	{ SR_UNIT_OHM, "\xe2\x84\xa6" },
	54	{ SR_UNIT_FARAD, "F" },
a5892391 BV	55	{ SR_UNIT_KELVIN, "K" },
	56	{ SR_UNIT_CELSIUS, "\xc2\xb0""C" },
	57	{ SR_UNIT_FAHRENHEIT, "\xc2\xb0""F" },
	58	{ SR_UNIT_HERTZ, "Hz" },
	59	{ SR_UNIT_PERCENTAGE, "%" },
f7bcc686	60	{ SR_UNIT_BOOLEAN, "" },
a5892391 BV	61	{ SR_UNIT_SECOND, "s" },
	62	{ SR_UNIT_SIEMENS, "S" },
	63	{ SR_UNIT_DECIBEL_MW, "dBu" },
	64	{ SR_UNIT_DECIBEL_VOLT, "dBv" },
f7bcc686	65	{ SR_UNIT_UNITLESS, "" },
a5892391 BV	66	{ SR_UNIT_DECIBEL_SPL, "dB" },
	67	{ SR_UNIT_CONCENTRATION, "ppm" },
	68	{ SR_UNIT_REVOLUTIONS_PER_MINUTE, "RPM" },
	69	{ SR_UNIT_VOLT_AMPERE, "VA" },
	70	{ SR_UNIT_WATT, "W" },
	71	{ SR_UNIT_WATT_HOUR, "Wh" },
	72	{ SR_UNIT_METER_SECOND, "m/s" },
	73	{ SR_UNIT_HECTOPASCAL, "hPa" },
	74	{ SR_UNIT_HUMIDITY_293K, "%rF" },
	75	{ SR_UNIT_DEGREE, "\xc2\xb0" },
f7bcc686 UH	76	{ SR_UNIT_HENRY, "H" },
	77	{ SR_UNIT_GRAM, "g" },
	78	{ SR_UNIT_CARAT, "ct" },
	79	{ SR_UNIT_OUNCE, "oz" },
	80	{ SR_UNIT_TROY_OUNCE, "oz t" },
	81	{ SR_UNIT_POUND, "lb" },
	82	{ SR_UNIT_PENNYWEIGHT, "dwt" },
	83	{ SR_UNIT_GRAIN, "gr" },
	84	{ SR_UNIT_TAEL, "tael" },
	85	{ SR_UNIT_MOMME, "momme" },
	86	{ SR_UNIT_TOLA, "tola" },
	87	{ SR_UNIT_PIECE, "pcs" },
a5892391 BV	88	ALL_ZERO
	89	};
	90
	91	static struct unit_mq_string mq_strings[] = {
a5892391 BV	92	{ SR_MQFLAG_AC, " AC" },
	93	{ SR_MQFLAG_DC, " DC" },
	94	{ SR_MQFLAG_RMS, " RMS" },
	95	{ SR_MQFLAG_DIODE, " DIODE" },
	96	{ SR_MQFLAG_HOLD, " HOLD" },
	97	{ SR_MQFLAG_MAX, " MAX" },
	98	{ SR_MQFLAG_MIN, " MIN" },
	99	{ SR_MQFLAG_AUTORANGE, " AUTO" },
	100	{ SR_MQFLAG_RELATIVE, " REL" },
f7bcc686 UH	101	{ SR_MQFLAG_SPL_FREQ_WEIGHT_A, "(A)" },
	102	{ SR_MQFLAG_SPL_FREQ_WEIGHT_C, "(C)" },
	103	{ SR_MQFLAG_SPL_FREQ_WEIGHT_Z, "(Z)" },
	104	{ SR_MQFLAG_SPL_FREQ_WEIGHT_FLAT, "(SPL)" },
	105	{ SR_MQFLAG_SPL_TIME_WEIGHT_S, " S" },
	106	{ SR_MQFLAG_SPL_TIME_WEIGHT_F, " F" },
	107	{ SR_MQFLAG_SPL_LAT, " LAT" },
	108	/* Not a standard function for SLMs, so this is a made-up notation. */
	109	{ SR_MQFLAG_SPL_PCT_OVER_ALARM, "%oA" },
	110	{ SR_MQFLAG_DURATION, " DURATION" },
a5892391 BV	111	{ SR_MQFLAG_AVG, " AVG" },
a5892391 BV	112	{ SR_MQFLAG_REFERENCE, " REF" },
f7bcc686	113	{ SR_MQFLAG_UNSTABLE, " UNSTABLE" },
a5892391 BV	114	ALL_ZERO
	115	};
	116
41caa319 AJ	117	SR_PRIV int sr_analog_init(struct sr_datafeed_analog2 *analog,
	118	struct sr_analog_encoding *encoding,
	119	struct sr_analog_meaning *meaning,
	120	struct sr_analog_spec *spec,
	121	int digits)
	122	{
	123	memset(analog, 0, sizeof(*analog));
	124	memset(encoding, 0, sizeof(*encoding));
	125	memset(meaning, 0, sizeof(*meaning));
	126	memset(spec, 0, sizeof(*spec));
	127
	128	analog->encoding = encoding;
	129	analog->meaning = meaning;
	130	analog->spec = spec;
	131
	132	encoding->unitsize = sizeof(float);
	133	encoding->is_float = TRUE;
	134	#ifdef WORDS_BIGENDIAN
	135	encoding->is_bigendian = TRUE;
	136	#else
	137	encoding->is_bigendian = FALSE;
	138	#endif
	139	encoding->digits = digits;
	140	encoding->is_digits_decimal = TRUE;
	141	encoding->scale.p = 1;
	142	encoding->scale.q = 1;
	143	encoding->offset.p = 0;
	144	encoding->offset.q = 1;
	145
	146	spec->spec_digits = digits;
	147
	148	return SR_OK;
	149	}
	150
fb019a0e	151	SR_API int sr_analog_to_float(const struct sr_datafeed_analog2 *analog,
4b4fdeea	152	float *outbuf)
fb019a0e BV	153	{
	154	float offset;
	155	unsigned int b, i;
	156	gboolean bigendian;
	157
	158	#ifdef WORDS_BIGENDIAN
	159	bigendian = TRUE;
	160	#else
	161	bigendian = FALSE;
	162	#endif
	163	if (!analog->encoding->is_float) {
	164	/* TODO */
4b4fdeea	165	sr_err("Only floating-point encoding supported so far.");
fb019a0e BV	166	return SR_ERR;
	167	}
	168
	169	if (analog->encoding->unitsize == sizeof(float)
	170	&& analog->encoding->is_bigendian == bigendian
	171	&& (analog->encoding->scale.p == analog->encoding->scale.q)
4b4fdeea	172	&& analog->encoding->offset.p / (float)analog->encoding->offset.q == 0) {
fb019a0e	173	/* The data is already in the right format. */
4b4fdeea	174	memcpy(outbuf, analog->data, analog->num_samples * sizeof(float));
fb019a0e BV	175	} else {
	176	for (i = 0; i < analog->num_samples; i += analog->encoding->unitsize) {
	177	for (b = 0; b < analog->encoding->unitsize; b++) {
	178	if (analog->encoding->is_bigendian == bigendian)
4b4fdeea	179	outbuf[i + b] = ((float )analog->data)[i analog->encoding->unitsize + b];
fb019a0e	180	else
4b4fdeea	181	outbuf[i + (analog->encoding->unitsize - b)] = ((float )analog->data)[i analog->encoding->unitsize + b];
fb019a0e BV	182	}
fb019a0e BV	183	if (analog->encoding->scale.p != analog->encoding->scale.q)
4b4fdeea BV	184	outbuf[i] = (outbuf[i] * analog->encoding->scale.p) / analog->encoding->scale.q;
	185	offset = ((float)analog->encoding->offset.p / (float)analog->encoding->offset.q);
	186	outbuf[i] += offset;
fb019a0e BV	187	}
	188	}
	189
	190	return SR_OK;
	191	}
c2a25ebb BV	192
	193	/*
	194	* Convert a floating point value to a string, limited to the given
	195	* number of decimal digits.
	196	*
	197	* @param value The value to convert.
	198	* @param digits Number of digits after the decimal point to print.
a24da9a8 ML	199	* @param result Pointer to store result.
	200	*
	201	* The string is allocated by the function and must be freed by the caller
	202	* after use by calling g_free().
c2a25ebb BV	203	*
	204	* @retval SR_OK
	205	*
	206	* @since 0.4.0
	207	*/
a24da9a8	208	SR_API int sr_analog_float_to_string(float value, int digits, char **result)
c2a25ebb BV	209	{
	210	int cnt, i;
	211
	212	/* This produces at least one too many digits */
a24da9a8	213	result = g_strdup_printf("%.f", digits, value);
a9b2283f BV	214	for (i = 0, cnt = 0; (*result)[i]; i++) {
a9b2283f BV	215	if (isdigit((*result)[i++]))
c2a25ebb BV	216	cnt++;
c2a25ebb BV	217	if (cnt == digits) {
a9b2283f	218	(*result)[i] = 0;
c2a25ebb BV	219	break;
	220	}
	221	}
	222
	223	return SR_OK;
	224	}
	225
a5892391 BV	226	/*
	227	* Convert the unit/MQ/MQ flags in the analog struct to a string.
	228	*
	229	* @param analog Struct containing the unit, MQ and MQ flags.
a24da9a8 ML	230	* @param result Pointer to store result.
	231	*
	232	* The string is allocated by the function and must be freed by the caller
	233	* after use by calling g_free().
a5892391 BV	234	*
	235	* @retval SR_OK
	236	*
	237	* @since 0.4.0
	238	*/
	239	SR_API int sr_analog_unit_to_string(const struct sr_datafeed_analog2 *analog,
a24da9a8	240	char **result)
a5892391	241	{
a24da9a8 ML	242	int i;
a24da9a8 ML	243	GString *buf = g_string_new(NULL);
a5892391	244
a5892391 BV	245	for (i = 0; unit_strings[i].value; i++) {
a5892391 BV	246	if (analog->meaning->unit == unit_strings[i].value) {
a24da9a8	247	g_string_assign(buf, unit_strings[i].str);
a5892391 BV	248	break;
	249	}
	250	}
	251
	252	/* More than one MQ flag may apply. */
a24da9a8 ML	253	for (i = 0; mq_strings[i].value; i++)
	254	if (analog->meaning->mqflags & mq_strings[i].value)
	255	g_string_append(buf, mq_strings[i].str);
	256
	257	*result = buf->str;
	258	g_string_free(buf, FALSE);
a5892391 BV	259
	260	return SR_OK;
	261	}
	262
90cefe0c BV	263	/*
	264	* Set sr_rational r to the given value.
	265	*
	266	* @param p Numerator
	267	* @param q Denominator
	268	*/
2617c81a	269	SR_API void sr_rational_set(struct sr_rational *r, uint64_t p, uint64_t q)
90cefe0c BV	270	{
	271	r->p = p;
	272	r->q = q;
	273	}
	274
e00b3f58	275	/** @} */