From 828654f79fa3e38e99319250de25350109786320 Mon Sep 17 00:00:00 2001
From: Julian Weigt <juw@posteo.de>
Date: Sun, 28 Dec 2025 11:59:45 +0000
Subject: Handle error estimate for maximum properly.

---
 double-error.c | 15 +++++++++++++++
 1 file changed, 15 insertions(+)

(limited to 'double-error.c')

diff --git a/double-error.c b/double-error.c
index aaaa405..4cd1406 100644
--- a/double-error.c
+++ b/double-error.c
@@ -37,6 +37,21 @@ double_error difference(double_error de1, double_error de2){
 
 bool is_greater(double_error de1, double_error de2){ return (de1.d > de2.d); }
 
+double_error maximum(double_error de1, double_error de2){
+	double_error de;
+	double max1 = de1.d+de1.e;
+	double max2 = de2.d+de2.e;
+	if( max1 >= max2 ){
+		de.d = de1.d;
+		de.e = de1.e + max1*EPS;
+	} else {
+	de.d = de2.d;
+	de.e = de2.e + max2*EPS;
+	}
+
+	return de;
+}
+
 double_error product(double_error de1, double_error de2){
 	double_error de;
 	de.d = de1.d * de2.d;
-- 
cgit v1.2.3