From 044403657d09aed9c56a9bca6decc864ed987f69 Mon Sep 17 00:00:00 2001
From: Robert Norton
Date: Thu, 22 Jun 2017 16:22:22 +0100
Subject: fix three different copies of the hardware_quot function to do proper
 trucation towards zero. Previous version was incorrect if result was exact
 and a<0 and b>0.

---
 src/lem_interp/interp_lib.lem | 15 +++++++++------
 1 file changed, 9 insertions(+), 6 deletions(-)

(limited to 'src/lem_interp')

diff --git a/src/lem_interp/interp_lib.lem b/src/lem_interp/interp_lib.lem
index fe87c9e1..dd39373d 100644
--- a/src/lem_interp/interp_lib.lem
+++ b/src/lem_interp/interp_lib.lem
@@ -66,13 +66,16 @@ let hardware_mod (a: integer) (b:integer) : integer =
  then (a mod b) - b
  else a mod b
 
+(* There are different possible answers for integer divide regarding
+rounding behaviour on negative operands. Positive operands always
+round down so derive the one we want (trucation towards zero) from
+that *)
 let hardware_quot (a:integer) (b:integer) : integer = 
-  if a < 0 && b < 0
-  then (abs a) / (abs b)
-  else if (a < 0 && b > 0) 
-  then (a/b) + 1
-  else a/b
-
+  let q = (abs a) / (abs b) in
+  if ((a<0) = (b<0)) then
+    q  (* same sign -- result positive *)
+  else
+    ~q (* different sign -- result negative *)
 
 let (max_64u : integer) = integerFromNat ((natPow 2 64) - 1)
 let (max_64  : integer) = integerFromNat ((natPow 2 63) - 1)
-- 
cgit v1.2.3