core: Fixes for int-max, int-min.
[mes.git] / src / math.c
index c0b6190b434375fb23d4985f83509414afe00c51..2a90f93326f7a71ca4d25544221e04b9814b7df1 100644 (file)
@@ -1,6 +1,6 @@
 /* -*-comment-start: "//";comment-end:""-*-
  * Mes --- Maxwell Equations of Software
- * Copyright © 2016,2017 Jan Nieuwenhuizen <janneke@gnu.org>
+ * Copyright © 2016,2017,2018 Jan (janneke) Nieuwenhuizen <janneke@gnu.org>
  *
  * This file is part of Mes.
  *
 
 #include <limits.h>
 
+void
+assert_number (char const* name, SCM x)
+{
+  if (TYPE (x) != TNUMBER)
+    {
+      eputs (name);
+      error (cell_symbol_not_a_number, x);
+    }
+}
+
 SCM
 greater_p (SCM x) ///((name . ">") (arity . n))
 {
-  int n = INT_MAX;
+  if (x == cell_nil)
+    return cell_t;
+  assert_number ("greater_p", CAR (x));
+  int n = VALUE (CAR (x));
+  x = CDR (x);
   while (x != cell_nil)
     {
-      assert (TYPE (car (x)) == TNUMBER);
-      if (VALUE (car (x)) >= n) return cell_f;
+      assert_number ("greater_p", CAR (x));
+      if (VALUE (car (x)) >= n)
+        return cell_f;
       n = VALUE (car (x));
       x = cdr (x);
     }
@@ -37,11 +52,16 @@ greater_p (SCM x) ///((name . ">") (arity . n))
 SCM
 less_p (SCM x) ///((name . "<") (arity . n))
 {
-  int n = INT_MIN;
+  if (x == cell_nil)
+    return cell_t;
+  assert_number ("less_p", CAR (x));
+  int n = VALUE (CAR (x));
+  x = CDR (x);
   while (x != cell_nil)
     {
-      assert (TYPE (car (x)) == TNUMBER);
-      if (VALUE (car (x)) <= n) return cell_f;
+      assert_number ("less_p", CAR (x));
+      if (VALUE (car (x)) <= n)
+        return cell_f;
       n = VALUE (car (x));
       x = cdr (x);
     }
@@ -51,13 +71,15 @@ less_p (SCM x) ///((name . "<") (arity . n))
 SCM
 is_p (SCM x) ///((name . "=") (arity . n))
 {
-  if (x == cell_nil) return cell_t;
-  assert (TYPE (car (x)) == TNUMBER);
-  int n = VALUE (car (x));
+  if (x == cell_nil)
+    return cell_t;
+  assert_number ("is_p", CAR (x));
+  int n = VALUE (CAR (x));
   x = cdr (x);
   while (x != cell_nil)
     {
-      if (VALUE (car (x)) != n) return cell_f;
+      if (VALUE (car (x)) != n)
+        return cell_f;
       x = cdr (x);
     }
   return cell_t;
@@ -66,15 +88,14 @@ is_p (SCM x) ///((name . "=") (arity . n))
 SCM
 minus (SCM x) ///((name . "-") (arity . n))
 {
-  SCM a = car (x);
-  assert (TYPE (a) == TNUMBER);
-  int n = VALUE (a);
+  assert_number ("minus", CAR (x));
+  int n = VALUE (CAR (x));
   x = cdr (x);
   if (x == cell_nil)
     n = -n;
   while (x != cell_nil)
     {
-      assert (TYPE (car (x)) == TNUMBER);
+      assert_number ("minus", CAR (x));
       n -= VALUE (car (x));
       x = cdr (x);
     }
@@ -87,7 +108,7 @@ plus (SCM x) ///((name . "+") (arity . n))
   int n = 0;
   while (x != cell_nil)
     {
-      assert (TYPE (car (x)) == TNUMBER);
+      assert_number ("plus", CAR (x));
       n += VALUE (car (x));
       x = cdr (x);
     }
@@ -98,14 +119,15 @@ SCM
 divide (SCM x) ///((name . "/") (arity . n))
 {
   int n = 1;
-  if (x != cell_nil) {
-    assert (TYPE (car (x)) == TNUMBER);
-    n = VALUE (car (x));
-    x = cdr (x);
-  }
+  if (x != cell_nil)
+    {
+      assert_number ("divide", CAR (x));
+      n = VALUE (car (x));
+      x = cdr (x);
+    }
   while (x != cell_nil)
     {
-      assert (TYPE (car (x)) == TNUMBER);
+      assert_number ("divide", CAR (x));
       n /= VALUE (car (x));
       x = cdr (x);
     }
@@ -115,8 +137,8 @@ divide (SCM x) ///((name . "/") (arity . n))
 SCM
 modulo (SCM a, SCM b)
 {
-  assert (TYPE (a) == TNUMBER);
-  assert (TYPE (b) == TNUMBER);
+  assert_number ("modulo", a);
+  assert_number ("modulo", b);
   int x = VALUE (a);
   while (x < 0) x += VALUE (b);
   return MAKE_NUMBER (x % VALUE (b));
@@ -128,31 +150,52 @@ multiply (SCM x) ///((name . "*") (arity . n))
   int n = 1;
   while (x != cell_nil)
     {
-      assert (TYPE (car (x)) == TNUMBER);
+      assert_number ("multiply", CAR (x));
       n *= VALUE (car (x));
       x = cdr (x);
     }
   return MAKE_NUMBER (n);
 }
 
+SCM
+logand (SCM x) ///((arity . n))
+{
+  int n = 0;
+  while (x != cell_nil)
+    {
+      assert_number ("multiply", CAR (x));
+      n &= VALUE (car (x));
+      x = cdr (x);
+    }
+  return MAKE_NUMBER (n);
+}
+
 SCM
 logior (SCM x) ///((arity . n))
 {
   int n = 0;
   while (x != cell_nil)
     {
-      assert (TYPE (car (x)) == TNUMBER);
+      assert_number ("logior", CAR (x));
       n |= VALUE (car (x));
       x = cdr (x);
     }
   return MAKE_NUMBER (n);
 }
 
+SCM
+lognot (SCM x)
+{
+  assert_number ("lognot", x);
+  int n = ~VALUE (x);
+  return MAKE_NUMBER (n);
+}
+
 SCM
 ash (SCM n, SCM count)
 {
-  assert (TYPE (n) == TNUMBER);
-  assert (TYPE (count) == TNUMBER);
+  assert_number ("ash", n);
+  assert_number ("ash", count);
   int cn = VALUE (n);
   int ccount = VALUE (count);
   return MAKE_NUMBER ((ccount < 0) ? cn >> -ccount : cn << ccount);