0.7.7.38:
[sbcl.git] / tests / irrat.pure.lisp
1 ;;;; tests of irrational floating point functions
2
3 ;;;; This software is part of the SBCL system. See the README file for
4 ;;;; more information.
5 ;;;;
6 ;;;; While most of SBCL is derived from the CMU CL system, the test
7 ;;;; files (like this one) were written from scratch after the fork
8 ;;;; from CMU CL.
9 ;;;; 
10 ;;;; This software is in the public domain and is provided with
11 ;;;; absolutely no warranty. See the COPYING and CREDITS files for
12 ;;;; more information.
13
14 (in-package :cl-user)
15 \f
16 ;;;; old bugs
17
18 ;;; This used to fail with
19 ;;;   The value -0.44579905382680446d0 is not of type (DOUBLE-FLOAT 0.0d0).
20 ;;; MNA's port of Raymond Toy work on CMU CL fixed this in sbcl-0.6.12.53.
21 (assert (equal (log #c(0.4 0.5)) #C(-0.44579905 0.8960554)))
22 \f
23 ;;;; other tests
24
25 ;;; expt
26 (assert (equal (expt #c(0 1) 2) -1))
27 (assert (equal (prin1-to-string (expt 2 #c(0 1))) "#C(0.7692389 0.63896126)"))
28
29 ;;; log
30 (assert (equal (prin1-to-string (log -3 10)) "#C(0.47712126 1.3643764)"))
31 (assert (= (log 3 0) 0))
32
33 ;;; sqrt, isqrt
34 (assert (= (sqrt 9) 3.0))
35 (assert (= (sqrt -9.0) #c(0.0 3.0)))
36 (assert (= (isqrt 9) 3))
37 (assert (= (isqrt 26) 5))
38
39
40 ;;; sin, sinh, asin, asinh
41 (assert (equal (prin1-to-string (sin (* 8 (/ pi 2)))) "-4.898425415289509d-16"))
42 (assert (equal (prin1-to-string (sin (expt 10 3))) "0.82687956"))
43 (assert (= (sinh 0) 0.0))
44 (assert (equal (prin1-to-string (sinh #c(5.0 -9.6)))
45                "#C(-73.06699 12.936809)"))
46 (assert (= (sin (* #c(0 1) 5)) (* #c(0 1) (sinh 5))))
47 (assert (= (sinh (* #c(0 1) 5)) (* #c(0 1) (sin 5))))
48 (assert (equal (prin1-to-string (asin -1)) "-1.5707964"))
49 (assert (= (asin 0) 0.0))
50 (assert (= (asin 2) #c(1.5707964 -1.3169578)))
51 (assert (equal (prin1-to-string (asinh 0.5)) "0.4812118"))
52 (assert (equal (prin1-to-string (asinh 3/7)) "0.41643077"))
53
54 ;;; cos, cosh, acos, acosh
55 (assert (=  (cos 0) 1.0))
56 (assert (equal (prin1-to-string (cos (/ pi 2))) "6.123031769111886d-17"))
57 (assert (= (cosh 0) 1.0))
58 (assert (equal (prin1-to-string (cosh 1)) "1.5430807"))
59 (assert (= (cos (* #c(0 1) 5)) (cosh 5)))
60 (assert (= (cosh (* #c(0 1) 5)) (cos 5)))
61 (assert (equal (prin1-to-string (acos 0)) "1.5707964"))
62 (assert (equal (prin1-to-string (acos -1)) "3.1415927"))
63 (assert (equal (prin1-to-string (acos 2)) "#C(0.0 1.3169578)"))
64 (assert (= (acos 1.00001) #c(0.0 0.0044751678)))
65 (assert (= (acosh 0) #c(0 1.5707964)))
66 (assert (= (acosh 1) 0))
67 (assert (= (acosh -1) #c(0 3.1415927)))
68
69 ;;; tan, tanh
70 (assert (equal (prin1-to-string (tan 1)) "1.5574077"))
71 (assert (equal (prin1-to-string (tan (/ pi 2))) "1.6331778728383844d+16"))
72 (assert (equal (prin1-to-string (tanh 0.00753)) "0.0075298576"))
73 (assert (= (tanh 50) 1.0))
74 (assert (= (tan (* #c(0 1) 5)) (* #c(0 1) (tanh 5))))
75 (assert (= (atan 1) 0.7853982))
76 (assert (equal (prin1-to-string (atanh 0.5) ) "0.54930615"))
77 (assert (equal (prin1-to-string (atanh 3/7)) "0.45814538"))