;-*- Mode:     Lisp -*-
;;;; Author:   Paul Dietz
;;;; Created:  Wed Feb 11 19:53:33 2004
;;;; Contains: Tests of LOG



(deftest log.1
  (let ((result (log 1)))
    (or (eqlt result 0)
        (eqlt result 0.0)))
  t)

(deftest log.2
  (mapcar #'log '(1.0s0 1.0f0 1.0d0 1.0l0))
  (0.0s0 0.0f0 0.0d0 0.0l0))

(deftest log.3
  (loop for type in '(short-float single-float double-float long-float)
        nconc
        (loop
         for x = (+ (random (coerce 1 type)) (/ 1 1000))
         for rlist = (multiple-value-list (log x))
         for y = (car rlist)
         repeat 1000
         unless (and (null (cdr rlist))
                     (typep y type))
         collect (list x rlist)))
  nil)

(deftest log.4
  (loop for type in '(short-float single-float double-float long-float)
        nconc
        (loop
         for x = (1+ (random (coerce 1000000 type)))
         for rlist = (multiple-value-list (log x))
         for y = (car rlist)
         repeat 1000
         unless (and (null (cdr rlist))
                     (typep y type))
         collect (list x rlist)))
  nil)

(deftest log.5
  (loop for type in '(short-float single-float double-float long-float)
        for zero = (coerce 0 type)
        nconc
        (loop
         for x = (- (random (coerce 1 type)))
         for rlist = (and (/= x zero) (multiple-value-list (log x)))
         for y = (car rlist)
         repeat 1000
         unless (or (= x zero)
                    (and (null (cdr rlist))
                         (typep y `(complex ,type))))
         collect (list x rlist)))
  nil)

(deftest log.6
  (loop for type in '(short-float single-float double-float long-float)
        for zero = (coerce 0 type)
        nconc
        (loop
         for x = (- (random (coerce 1000000 type)))
         for rlist = (and (/= x zero) (multiple-value-list (log x)))
         for y = (car rlist)
         repeat 1000
         unless (or (= x zero)
                    (and (null (cdr rlist))
                         (typep y `(complex ,type))))
         collect (list x rlist)))
  nil)

(deftest log.7
  (loop for type in '(short-float single-float double-float long-float)
        for zero = (coerce 0 type)
        nconc
        (loop
         for x1 = (- (random (coerce 2000 type)) 1000)
         for x2 = (1+ (random (coerce 1000 type)))
         for rlist = (and (/= x1 zero)
                          (multiple-value-list (log (complex x1 x2))))
         for y = (car rlist)
         repeat 1000
         unless (or (= x1 zero)
                    (and (null (cdr rlist))
                         (typep y `(complex ,type))))
         collect (list x1 x2 rlist)))
  nil)

(deftest log.8
  (loop for type in '(short-float single-float double-float long-float)
        for zero = (coerce 0 type)
        nconc
        (loop
         for x1 = (- (random (coerce 2000 type)) 1000)
         for x2 = (- -1 (random (coerce 1000 type)))
         for rlist = (and (/= x1 zero)
                          (multiple-value-list (log (complex x1 x2))))
         for y = (car rlist)
         repeat 1000
         unless (or (= x1 zero)
                    (and (null (cdr rlist))
                         (typep y `(complex ,type))))
         collect (list x1 x2 rlist)))
  nil)

;;; FIXME
;;; Add tests for two-arg calls

;;; FIXME
;;; More accuracy tests here

;;; Error tests

(deftest log.error.1
  (signals-error (log) program-error)
  t)

(deftest log.error.2
  (signals-error (log 1.0 2.0 3.0) program-error)
  t)