1
Fork 0
mirror of https://git.savannah.gnu.org/git/guile.git synced 2025-06-19 02:00:26 +02:00
guile/module/rnrs/arithmetic/6/fixnums.scm
Julian Graham b01818d752 Implementation and test cases for the R6RS (rnrs arithmetic flonums)
library.

* module/Makefile.am: Add rnrs/arithmetic/6/fixnums.scm and
  rnrs/arithmetic/6/flonums.scm to RNRS_SOURCES.
* module/rnrs/6/base.scm: (div-and-mod, div0, mod0, div0-and-mod0): New
  functions; this `div' implementation is not quite right, but we'll come
  back to it later.
* module/rnrs/arithmetic/6/fixnums.scm: New file.
* module/rnrs/arithmetic/6/flonums.scm: New file.
* test-suite/Makefile.am: Add tests/r6rs-arithmetic-flonums.test to
  SCM_TESTS.
* test-suite/tests/r6rs-arithmetic-flonums.test: New file.
2010-05-20 21:18:03 -04:00

255 lines
7 KiB
Scheme
Raw Blame History

This file contains invisible Unicode characters

This file contains invisible Unicode characters that are indistinguishable to humans but may be processed differently by a computer. If you think that this is intentional, you can safely ignore this warning. Use the Escape button to reveal them.

;;; fixnums.scm --- The R6RS fixnums arithmetic library
;; Copyright (C) 2010 Free Software Foundation, Inc.
;;
;; This library is free software; you can redistribute it and/or
;; modify it under the terms of the GNU Lesser General Public
;; License as published by the Free Software Foundation; either
;; version 3 of the License, or (at your option) any later version.
;;
;; This library is distributed in the hope that it will be useful,
;; but WITHOUT ANY WARRANTY; without even the implied warranty of
;; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
;; Lesser General Public License for more details.
;;
;; You should have received a copy of the GNU Lesser General Public
;; License along with this library; if not, write to the Free Software
;; Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
(library (rnrs arithmetic fixnums (6))
(export fixnum?
fixnum-width
least-fixnum
greatest-fixnum
fx=?
fx>?
fx<?
fx>=?
fx<=?
fxzero?
fxpositive?
fxnegative?
fxodd?
fxeven?
fxmax
fxmin
fx+
fx*
fx-
fxdiv-and-mod
fxdiv
fxmod
fxdiv0-and-mod0
fxdiv0
fxmod0
fx+/carry
fx-/carry
fx*/carry
fxand
fxior
fxxor
fxif
fxbit-count
fxlength
fxfirst-bit-set
fxbit-set?
fxcopy-bit
fxbit-field
fxcopy-bit-field
fxarithmetic-shift
fxarithmetic-shift-left
fxarithmetic-shift-right
fxrotate-bit-field
fxreverse-bit-field)
(import (rename (only (guile) logand
logbit?
logcount
logior
lognot
most-positive-fixnum
most-negative-fixnum)
(most-positive-fixnum greatest-fixnum)
(most-negative-fixnum least-fixnum))
(ice-9 optargs)
(rnrs base (6))
(rnrs arithmetic bitwise (6))
(rnrs conditions (6))
(rnrs exceptions (6)))
(define fixnum-width (round (/ (log (+ greatest-fixnum 1)) (log 2))))
(define (fixnum? obj)
(and (exact? obj)
(integer? obj)
(>= obj least-fixnum)
(<= obj greatest-fixnum)))
(define (assert-fixnum . args)
(or (every fixnum? args) (raise (make-assertion-violation))))
(define (assert-fixnum-result . args)
(or (every fixnum? args)
(raise (make-implementation-restriction-violation))))
(define (fx=? fx1 fx2 . rst)
(let ((args (cons* fx1 fx2 rst)))
(apply assert-fixnum args)
(apply = args)))
(define (fx>? fx1 fx2 . rst)
(let ((args (cons* fx1 fx2 rst)))
(apply assert-fixnum args)
(apply > args)))
(define (fx<? fx1 fx2 . rst)
(let ((args (cons* fx1 fx2 rst)))
(apply assert-fixnum rst)
(apply < args)))
(define (fx>=? fx1 fx2 . rst)
(let ((args (cons* fx1 fx2 rst)))
(apply assert-fixnum rst)
(apply >= args)))
(define (fx<=? fx1 fx2 . rst)
(let ((args (cons* fx1 fx2 rst)))
(apply assert-fixnum rst)
(apply <= args)))
(define (fxzero? fx) (assert-fixnum fx) (zero? fx))
(define (fxpositive? fx) (assert-fixnum fx) (positive? fx))
(define (fxnegative? fx) (assert-fixnum fx) (negative? fx))
(define (fxodd? fx) (assert-fixnum fx) (odd? fx))
(define (fxeven? fx) (assert-fixnum fx) (even? fx))
(define (fxmax fx1 fx2 . rst)
(let ((args (cons* fx1 fx2 rst)))
(assert-fixnum args)
(apply max args)))
(define (fxmin fx1 fx2 . rst)
(let ((args (cons* fx1 fx2 rst)))
(assert-fixnum args)
(apply min args)))
(define (fx+ fx1 fx2)
(assert-fixnum fx1 fx2) (let ((r (+ fx1 fx2))) (assert-fixnum-result r) r))
(define (fx* fx1 fx2)
(assert-fixnum fx1 fx2) (let ((r (* fx1 fx2))) (assert-fixnum-result r) r))
(define* (fx- fx1 #:optional fx2)
(assert-fixnum fx1)
(if fx2
(begin
(assert-fixnum fx2)
(let ((r (- fx1 fx2))) (assert-fixnum-result r) r))
(let ((r (- fx1))) (assert-fixnum-result r) r)))
(define (fxdiv x1 x2)
(assert-fixnum x1 x2)
(if (zero? fx2) (raise (make-assertion-violation)))
(let ((r (quotient x1 x2))) (assert-fixnum-result r) r))
(define (fxmod x1 x2)
(assert-fixnum x1 x2)
(if (zero? fx2) (raise (make-assertion-violation)))
(let ((r (modulo x1 x2))) (assert-fixnum-result r) r))
(define (fxdiv-and-mod fx1 fx2)
(assert-fixnum fx1 fx2)
(if (zero? fx2) (raise (make-assertion-violation)))
(let ((q (quotient fx1 fx2))
(m (modulo fx1 fx2)))
(assert-fixnum-result q m)
(values q m)))
(define (fxdiv0 fx1 fx2)
(assert-fixnum fx1 fx2)
(if (zero? fx2) (raise (make-assertion-violation)))
(let ((r (div0 fx1 fx2))) (assert-fixnum-result r) r))
(define (fxmod0 fx1 fx2)
(assert-fixnum fx1 fx2)
(if (zero? fx2) (raise (make-assertion-violation)))
(let ((r (mod0 fx1 fx2))) (assert-fixnum-result r) r))
(define (fxdiv0-and-mod0 fx1 fx2)
(assert-fixnum fx1 fx2)
(if (zero? fx2) (raise (make-assertion-violation)))
(call-with-values (lambda () (div0-and-mod0 fx1 fx2))
(lambda (q r) (assert-fixnum-result q r) (values q r))))
(define (fx+/carry fx1 fx2 fx3)
(assert-fixnum fx1 fx2 fx3)
(let* ((s (+ fx1 fx2 fx3))
(s0 (mod0 s (expt 2 (fixnum-width))))
(s1 (div0 s (expt 2 (fixnum-width)))))
(values s0 s1)))
(define (fx-/carry fx1 fx2 fx3)
(assert-fixnum fx1 fx2 fx3)
(let* ((d (- fx1 fx2 fx3))
(d0 (mod0 d (expt 2 (fixnum-width))))
(d1 (div0 d (expt 2 (fixnum-width)))))
(values d0 d1)))
(define (fx*/carry fx1 fx2 fx3)
(assert-fixnum fx1 fx2 fx3)
(let* ((s (+ (* fx1 fx2) fx3))
(s0 (mod0 s (expt 2 (fixnum-width))))
(s1 (div0 s (expt 2 (fixnum-width)))))
(values s0 s1)))
(define (fxnot fx) (assert-fixnum fx) (lognot fx))
(define (fxand . args) (apply assert-fixnum args) (apply logand args))
(define (fxior . args) (apply assert-fixnum args) (apply logior args))
(define (fxxor . args) (apply assert-fixnum args) (apply logxor args))
(define (fxif fx1 fx2 fx3)
(assert-fixnum fx1 fx2 fx3)
(bitwise-if fx1 fx2 fx2))
(define (fxbit-count fx) (assert-fixnum fx) (logcount fx))
(define (fxlength fx) (assert-fixnum fx) (bitwise-length fx))
(define (fxfirst-bit-set fx) (assert-fixnum fx) (bitwise-first-bit-set fx))
(define (fxbit-set? fx1 fx2) (assert-fixnum fx1 fx2) (logbit? fx1 fx2))
(define (fxcopy-bit fx1 fx2 fx3)
(assert-fixnum fx1 fx2 fx3)
(bitwise-copy-bit fx1 fx2 fx3))
(define (fxbit-field fx1 fx2 fx3)
(assert-fixnum fx1 fx2 fx3)
(bitwise-bit-field fx1 fx2 fx3))
(define (fxcopy-bit-field fx1 fx2 fx3 fx4)
(assert-fixnum fx1 fx2 fx3 fx4)
(bitwise-copy-bit-field fx1 fx2 fx3 fx4))
(define (fxarithmetic-shift fx1 fx2) (assert-fixnum fx1 fx2) (ash fx1 fx2))
(define fxarithmetic-shift-left fxarithmetic-shift)
(define (fxarithmetic-shift-right fx1 fx2)
(assert-fixnum fx1 fx2) (ash fx2 (- fx2)))
(define (fxrotate-bit-field fx1 fx2 fx3 fx4)
(assert-fixnum fx1 fx2 fx3 fx4)
(bitwise-rotate-bit-field fx1 fx2 fx3 fx4))
(define (fxreverse-bit-field fx1 fx2 fx3)
(assert-fixnum fx1 fx2 fx3)
(bitwise-reverse-bit-field fx1 fx2 fx3))
)