examples/coma/product.coma

   1 use int.Int
   2 use int.ComputerDivision
   3 use coma.Std
   4
   5 let product (a b: int) {b >= 0} (return (c: int) {c = a * b}) =
   6   loop
   7   [ loop {p * q + r = a * b /\ q >= 0}
   8     = if {q > 0} next (-> return {r})
   9       [ next ->
  10           if {mod q 2 = 1} (-> [&r <- r + p] step) step
  11           [ step -> [&p <- p + p | &q <- div q 2] loop ] ] ]
  12   [ &p: int = a | &q: int = b | &r: int = 0 ]
  13
  14 let product_pure (a b: int) {b >= 0} (return (c: int) {c = a * b}) =
  15   loop {a} {b} {0}
  16   [ loop (p q r: int) {p * q + r = a * b /\ q >= 0}
  17     = if {q > 0} (-> if {mod q 2 = 1} (-> next {r+p}) (-> next {r})
  18                      [ next (s: int) -> loop {p+p} {div q 2} {s} ])
  19                  (-> return {r})
  20   ]
  21
  22 let product_inline (a b: int) {b >= 0} (return (c: int) {c = a * b}) =
  23   loop {a} {b} {0}
  24   [ loop (p q r: int) {p * q + r = a * b /\ q >= 0}
  25     = if {q > 0} (-> if {mod q 2 = 1} (-> loop {p+p} {div q 2} {r+p})
  26                                       (-> loop {p+p} {div q 2} {r}))
  27                  (-> return {r})
  28   ]