# HG changeset patch
# User nipkow
# Date 977147878 -3600
# Node ID 6077fd9335751978b26c7860012d2c33286f2670
# Parent  4ea37fba9c02d9de25afa6f4136f00ec18287b3c
added rational arithmetic

diff -r 4ea37fba9c02 -r 6077fd933575 src/Pure/library.ML
--- a/src/Pure/library.ML	Mon Dec 18 14:57:34 2000 +0100
+++ b/src/Pure/library.ML	Mon Dec 18 14:57:58 2000 +0100
@@ -117,6 +117,8 @@
   val take_suffix: ('a -> bool) -> 'a list -> 'a list * 'a list
 
   (*integers*)
+  val gcd: int * int -> int
+  val lcm: int * int -> int
   val inc: int ref -> int
   val dec: int ref -> int
   val upto: int * int -> int list
@@ -127,6 +129,16 @@
   val string_of_int: int -> string
   val string_of_indexname: string * int -> string
 
+  (*rational numbers*)
+  type rat
+  val rep_rat: rat -> int * int
+  val ratadd: rat * rat -> rat
+  val ratmul: rat * rat -> rat
+  val ratinv: rat -> rat
+  val int_ratdiv: int * int -> rat
+  val ratneg: rat -> rat
+  val rat_of_int: int -> rat
+
   (*strings*)
   val nth_elem_string: int * string -> string
   val foldl_string: ('a * string -> 'a) -> 'a * string -> 'a
@@ -635,6 +647,13 @@
 
 (** integers **)
 
+fun gcd(x,y) =
+  let fun gxd x y =
+    if y = 0 then x else gxd y (x mod y)
+  in if x < y then gxd y x else gxd x y end;
+
+fun lcm(x,y) = (x * y) div gcd(x,y);
+
 fun inc i = (i := ! i + 1; ! i);
 fun dec i = (i := ! i - 1; ! i);
 
@@ -679,7 +698,6 @@
   | string_of_indexname (a,i) = a ^ "_" ^ Int.toString i;
 
 
-
 (** strings **)
 
 (*functions tuned for strings, avoiding explode*)
@@ -1207,6 +1225,35 @@
 
 
 
+(** rational numbers **)
+
+datatype rat = Rat of bool * int * int
+
+fun rep_rat(Rat(a,p,q)) = (if a then p else ~p,q)
+
+fun ratnorm(a,p,q) = if p=0 then Rat(a,0,1) else
+  let val absp = abs p
+      val m = gcd(absp,q)
+  in Rat(a = (0 <= p), absp div m, q div m) end;
+
+fun ratadd(Rat(a,p,q),Rat(b,r,s)) =
+  let val den = lcm(q,s)
+      val p = p*(den div q) and r = r*(den div s)
+      val num = (if a then p else ~p) + (if b then r else ~r)
+  in ratnorm(true,num,den) end;
+
+fun ratmul(Rat(a,p,q),Rat(b,r,s)) = ratnorm(a=b,p*r,q*s)
+
+fun ratinv(Rat(a,p,q)) = if p=0 then error("ratinv") else Rat(a,q,p)
+
+fun int_ratdiv(p,q) =
+  if q=0 then error("int_ratdiv") else ratnorm(0<=q, p, abs q)
+
+fun ratneg(Rat(b,p,q)) = Rat(not b,p,q);
+
+fun rat_of_int i = if i < 0 then Rat(false,abs i,1) else Rat(true,i,1);
+
+
 (** misc **)
 
 fun overwrite_warn (args as (alist, (a, _))) msg =