Changeset 9adb61d in mainline for uspace/lib/math/generic/mod.c

Timestamp:

2015-09-05T11:50:00Z (10 years ago)

Author:

Jiri Svoboda <jiri@…>

Branches:

lfn, master, serial, ticket/834-toolchain-update, topic/msim-upgrade, topic/simplify-dev-export

Children:

996dc042, ba8eecf

Parents:

e6f5766

Message:

Add single-precision variant for all functions. Allow generic implementations to call other functions while selecting the number of bits of precision, but not the implementation (generic or arch-specific).

File:

• uspace/lib/math/generic/mod.c (modified) (2 diffs)

uspace/lib/math/generic/mod.c

@@ -36 +36 @@
 #include <mod.h>
 
-/** Double precision modulo
+/** Remainder function (32-bit floating point)
+ *
+ * Calculate the modulo of dividend by divisor.
+ *
+ * This is a very basic implementation that uses
+ * division and multiplication (instead of exact
+ * arithmetics). Thus the result might be very
+ * imprecise (depending on the magnitude of the
+ * arguments).
+ *
+ * @param dividend Dividend.
+ * @param divisor  Divisor.
+ *
+ * @return Modulo.
+ *
+ */
+float32_t float32_mod(float32_t dividend, float32_t divisor)
+{
+        // FIXME: replace with exact arithmetics
+        
+        float32_t quotient = trunc_f32(dividend / divisor);
+        
+        return (dividend - quotient * divisor);
+}
+
+/** Remainder function (64-bit floating point)
  *
  * Calculate the modulo of dividend by divisor.
@@ -56 +81 @@
         // FIXME: replace with exact arithmetics
         
-        float64_t quotient = trunc(dividend / divisor);
+        float64_t quotient = trunc_f64(dividend / divisor);
         
         return (dividend - quotient * divisor);