acos.lsp 2.78 KB
Newer Older
pfdietz's avatar
pfdietz committed
1
2
3
4
5
6
7
8
9
;-*- Mode:     Lisp -*-
;;;; Author:   Paul Dietz
;;;; Created:  Tue Feb 10 05:39:24 2004
;;;; Contains: Tess of ACOS

(in-package :cl-test)

(deftest acos.1
  (loop for i from -1000 to 1000
sds's avatar
sds committed
10
11
12
13
        for rlist = (multiple-value-list (acos i))
        for y = (car rlist)
        always (and (null (cdr rlist))
                    (numberp y)))
pfdietz's avatar
pfdietz committed
14
15
16
17
  t)

(deftest acos.2
  (loop for type in '(short-float single-float double-float long-float)
sds's avatar
sds committed
18
19
20
21
22
23
24
25
26
        collect
        (let ((a (coerce 2000 type))
              (b (coerce -1000 type)))
          (loop for x = (- (random a) b)
                for rlist = (multiple-value-list (acos x))
                for y = (car rlist)
                repeat 1000
                always (and (null (cdr rlist))
                            (numberp y)))))
pfdietz's avatar
pfdietz committed
27
28
29
30
  (t t t t))

(deftest acos.3
  (loop for type in '(integer short-float single-float double-float long-float)
sds's avatar
sds committed
31
32
33
34
35
36
37
38
39
        collect
        (let ((a (coerce 2000 type))
              (b (coerce -1000 type)))
          (loop for x = (- (random a) b)
                for rlist = (multiple-value-list (acos (complex 0 x)))
                for y = (car rlist)
                repeat 1000
                always (and (null (cdr rlist))
                            (numberp y)))))
pfdietz's avatar
pfdietz committed
40
41
42
43
  (t t t t t))

(deftest acos.4
  (loop for type in '(integer short-float single-float double-float long-float)
sds's avatar
sds committed
44
45
46
47
48
49
50
51
52
53
        collect
        (let ((a (coerce 2000 type))
              (b (coerce -1000 type)))
          (loop for x1 = (- (random a) b)
                for x2 = (- (random a) b)
                for rlist = (multiple-value-list (acos (complex x1 x2)))
                for y = (car rlist)
                repeat 1000
                always (and (null (cdr rlist))
                            (numberp y)))))
pfdietz's avatar
pfdietz committed
54
55
56
57
58
59
60
61
  (t t t t t))

(deftest acos.5
  (approx= (acos 0) (coerce (/ pi 2) 'single-float))
  t)

(deftest acos.6
  (loop for type in '(single-float short-float double-float long-float)
sds's avatar
sds committed
62
63
64
        unless (approx= (acos (coerce 0 type))
                        (coerce (/ pi 2) type))
        collect type)
pfdietz's avatar
pfdietz committed
65
66
67
68
  nil)

(deftest acos.7
  (loop for type in '(single-float short-float double-float long-float)
sds's avatar
sds committed
69
70
71
        unless (approx= (acos (coerce 1 type))
                        (coerce 0 type))
        collect type)
pfdietz's avatar
pfdietz committed
72
73
74
75
  nil)

(deftest acos.8
  (loop for type in '(single-float short-float double-float long-float)
sds's avatar
sds committed
76
77
78
        unless (approx= (acos (coerce -1 type))
                        (coerce pi type))
        collect type)
pfdietz's avatar
pfdietz committed
79
80
  nil)

81
82
83
84
(deftest acos.9
  (macrolet ((%m (z) z)) (not (not (> (acos (expand-in-current-env (%m 0))) 0))))
  t)

pfdietz's avatar
pfdietz committed
85
86
87
88
89
90
91
92
93
94
95
96
97
98
;;; FIXME
;;; Add accuracy tests

;;; Error tests

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

(deftest acos.error.2
  (signals-error (acos 0.0 0.0) program-error)
  t)

(deftest acos.error.3
99
  (check-type-error #'acos #'numberp)
pfdietz's avatar
pfdietz committed
100
  nil)