Add support for xsd:double and xsd:float

These can be used to serialise a float or double in the shortest normalised
form that can be read back in to the exact same floating point value.

2 files changed, 83 insertions, 0 deletions

diff --git a/src/node.c b/src/node.c
index 98c8e0ac..3d63361f 100644
--- a/src/node.c
+++ b/src/node.c
@@ -34,6 +34,9 @@
 #include <string.h>
 
 // Define C11 numeric constants if the compiler hasn't already
+#ifndef FLT_DECIMAL_DIG
+#    define FLT_DECIMAL_DIG 9
+#endif
 #ifndef DBL_DECIMAL_DIG
 #    define DBL_DECIMAL_DIG 17
 #endif
@@ -696,6 +699,84 @@ serd_new_decimal(const double    d,
 	return node;
 }
 
+static SerdNode*
+serd_new_scientific(const double    d,
+                    const unsigned  precision,
+                    const SerdNode* datatype)
+{
+	const size_t type_len = serd_node_total_size(datatype);
+	const int    fpclass  = fpclassify(d);
+
+	if (fpclass == FP_INFINITE && d < 0) {
+		return serd_new_typed_literal("-INF", datatype);
+	} else if (fpclass == FP_INFINITE && d > 0) {
+		return serd_new_typed_literal("INF", datatype);
+	} else if (fpclass == FP_NAN) {
+		return serd_new_typed_literal("NaN", datatype);
+	} else if (fpclass == FP_ZERO) {
+		return signbit(d) ? serd_new_typed_literal("-0.0E0", datatype)
+		                  : serd_new_typed_literal("0.0E0", datatype);
+	}
+
+	// Get decimal digits
+	const double           abs_d      = fabs(d);
+	char                   digits[32] = {0};
+	const SerdDecimalCount count      = serd_decimals(abs_d, digits, precision);
+	assert(count.count == 1 || digits[count.count - 1] != '0');
+
+	// Calculate string length and allocate node
+	const size_t   len            = count.count;
+	const int      expt           = count.expt;
+	unsigned       abs_expt       = (unsigned)abs(expt);
+	const unsigned abs_exp_digits = (unsigned)serd_count_digits(abs_expt);
+
+	SerdNode* node = serd_node_malloc(type_len + len + abs_exp_digits + 4,
+	                                  SERD_HAS_DATATYPE,
+	                                  SERD_LITERAL);
+
+	char* const buf = serd_node_buffer(node);
+	if (d < 0.0) {
+		buf[node->n_bytes++] = '-';
+	}
+
+	// Write mantissa, with decimal point after the first (normal form)
+	buf[node->n_bytes++] = digits[0];
+	buf[node->n_bytes++] = '.';
+	if (len > 1) {
+		node->n_bytes += copy_digits(buf + node->n_bytes, digits + 1, len - 1);
+	} else {
+		buf[node->n_bytes++] = '0';
+	}
+
+	// Write exponent
+	buf[node->n_bytes++] = 'E';
+	if (expt < 0) {
+		buf[node->n_bytes++] = '-';
+	}
+	char* s = buf + node->n_bytes + abs_exp_digits - 1;
+	do {
+		*s-- = (char)('0' + (abs_expt % 10));
+	} while ((abs_expt /= 10) > 0);
+	node->n_bytes += abs_exp_digits;
+
+	memcpy(serd_node_meta(node), datatype, type_len);
+	serd_node_check_padding(node);
+	return node;
+}
+
+SerdNode*
+serd_new_double(const double d)
+{
+	return serd_new_scientific(d, DBL_DECIMAL_DIG, &serd_xsd_double.node);
+}
+
+SerdNode*
+serd_new_float(const float f)
+{
+	return serd_new_scientific(
+	        (double)f, FLT_DECIMAL_DIG, &serd_xsd_float.node);
+}
+
 SerdNode*
 serd_new_integer(int64_t i, const SerdNode* datatype)
 {
diff --git a/src/static_nodes.h b/src/static_nodes.h
index f4794d72..b7770432 100644
--- a/src/static_nodes.h
+++ b/src/static_nodes.h
@@ -32,6 +32,8 @@ typedef struct StaticNode {
 DEFINE_XSD_NODE(base64Binary)
 DEFINE_XSD_NODE(boolean)
 DEFINE_XSD_NODE(decimal)
+DEFINE_XSD_NODE(double)
+DEFINE_XSD_NODE(float)
 DEFINE_XSD_NODE(integer)
 
 #endif // SERD_STATIC_NODES_H