exact abs

bksaiki · Oct 9, 2023 · ccff825 · ccff825
1 parent fe343fd
commit ccff825
Showing 1 changed file with 12 additions and 1 deletion.
diff --git a/src/real/ops.rs b/src/real/ops.rs
@@ -4,7 +4,7 @@ use std::{
 };
 
 use crate::{
-    ops::{RoundedAdd, RoundedMul, RoundedNeg, RoundedSub},
+    ops::{RoundedAdd, RoundedMul, RoundedNeg, RoundedSub, RoundedAbs},
     rfloat::RFloat,
     Real, RoundingContext,
 };
@@ -22,6 +22,17 @@ impl RoundedNeg for RealContext {
     }
 }
 
+impl RoundedAbs for RealContext {
+    fn abs<N: Real>(&self, src: &N) -> Self::Rounded {
+        let src = self.round(src); // convert (exactly) to RFloat
+        match src {
+            RFloat::Real(_, exp, c) => RFloat::Real(false, exp, c),
+            RFloat::Infinite(_) => RFloat::Infinite(false),
+            RFloat::Nan => RFloat::Nan,
+        }
+    }
+}
+
 impl RoundedAdd for RealContext {
     fn add<N1, N2>(&self, src1: &N1, src2: &N2) -> Self::Rounded
     where