tests/trig.lsp

   1 ; test trigonometric functions
   2 (setf eps 1.e-6 deg2rad (/ pi 180) rad2deg (/ 180 pi))
   3
   4 (check #'< (abs (- (* 2 pi rad2deg) 360)) eps)
   5 (check #'< (abs (- (* (atan (/ 1 1)) rad2deg) 45)) eps)
   6 (check #'< (abs (- (atan (/ 5 5)) (/ pi 4))) eps)
   7 (check #'<
   8        (abs (- (sqrt (+ (^ 2 2) (^ 3 2)))  (/ 3 (sin (atan (/ 3 2)))))) eps)
   9 (check #'<
  10        (abs (- (cos (* pi .37)) (/ (sin (* pi .37)) (tan (* pi .37))))) eps)
  11 (check #'< (abs (- (cos (/ pi 4)) (sin (/ pi 4)))) eps)
  12 (check #'< (abs (- (acos -1) pi)) eps)
  13 (check #'< (abs (- (asin 1) (/ pi 2))) eps)
  14
  15 ; fft test
  16 (let* ((n 1000)
  17        (x (normal-rand n)))
  18   (check #'< (max (abs (- x (/ (fft (fft x) t) n)))) eps))