diff options
| author | Julian Weigt <juw@posteo.de> | 2025-12-28 11:59:45 +0000 |
|---|---|---|
| committer | Julian Weigt <juw@posteo.de> | 2026-02-04 15:55:48 +0100 |
| commit | 828654f79fa3e38e99319250de25350109786320 (patch) | |
| tree | 6f6c222af759531dc2e14e6a0619c3dfb1f25591 /double-error.c | |
| parent | 944f90c47ffcde862dfe5f258de0b1ebf229c20e (diff) | |
Handle error estimate for maximum properly.
Diffstat (limited to 'double-error.c')
| -rw-r--r-- | double-error.c | 15 |
1 files changed, 15 insertions, 0 deletions
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; |