;;;; mathh-test
;;;; Kon Lovett, May '17

;;;; Issues
;;;;

(require-extension test)

;;;

(require-extension mathh)

;;

(test-group "ISO C Functions"

	(test 1.0 (bessel-j0 0.0))
	(test 0.0 (bessel-j1 0.0))
	(test 1.0 (bessel-jn 0 0.0))
	(test 0.0 (bessel-jn 1 0.0))

	(test 0.088256964215677 (bessel-y0 1.0))
	(test -0.781212821300289 (bessel-y1 1.0))
	(test 0.088256964215677 (bessel-yn 0 1.0))
	(test -0.781212821300289 (bessel-yn 1 1.0))

	(test 1.0 (cosh 0.0))
	(test 0.0 (sinh 0.0))
	(test 0.0 (tanh 0.0))

	(test 5.0 (hypot -5.0 0))

	(test 1.0 (tgamma 1.0))
	(test 0.0 (lgamma 1.0))

	(test 2.0 (log10 100.0))
	(test 3.0 (log2 8.0))
	(test 0.0 (log1p 0.0))

	(test 0.0 (fpmod 0.0 1.0))

	(test (values 5.0 0.5) (modf 5.5))

	(test 20.0 (ldexp 5.0 2))

	(test 20.0 (scalbn 5.0 2))

	(test (values 0.536870912 -30) (frexp 5.0e-10))
)

;;

(test-group "BSD Functions"

	(test-assert (signbit -1.0))
	(test-assert (not (signbit 1.0)))
	(test-assert (signbit -0.0))

	(test -1.0 (copysign 1.0 -1.0))
	(test 1.0 (copysign -1.0 1.0))

	(test 1.0 (nextafter 1.0 -1.0))
	(test -1.0 (nextafter -1.0 1.0))

	(test 2.4662 (cbrt 15.0))
)

;;

(test-group "Function fpclass"

	(test 'negative-infinite (fpclass -inf.0))
	(test 'signaling-nan (fpclass -nan.0))
	(test 'negative-zero (fpclass -0.0))
	(test 'positive-normal (fpclass 0.741573033707865))
	(test 'positive-normal (fpclass (fp/ 33.0 44.5)))
)

;;

(test-group "Function fpclassify"

	(test 'infinite (fpclassify -inf.0))
	(test 'nan (fpclassify -nan.0))
	(test 'zero (fpclassify -0.0))
	(test 'normal (fpclassify 0.741573033707865))
	(test 'normal (fpclassify (fp/ 33.0 44.5)))
)

;;;

(require-extension fp-utils)

(define-constant 5eps (fp/ 9.0 1e06))
(define-constant 4eps (fp/ 9.0 1e05))

(test-group "FP Utils"

  (test-assert (fpzero? 0.0))
  (test-assert (not (fpzero? 1.0)))
  (test-assert (not (fpzero? maximum-flonum)))
  (test-assert (not (fpzero? minimum-flonum)))

  (test-assert (not (fppositive? 0.0)))
  (test-assert (not (fppositive? (fpneg minimum-flonum))))
  (test-assert (fppositive? maximum-flonum))

  (test-assert (fpcardinal? 0.0))
  (test-assert (not (fpcardinal? (fpneg minimum-flonum))))
  (test-assert (fpcardinal? maximum-flonum))

  (test-assert (not (fpnegative? 0.0)))
  (test-assert (fpnegative? (fpneg minimum-flonum)))
  (test-assert (not (fpnegative? maximum-flonum)))

	(test-assert (not (fpeven? 7.0)))
	(test-assert (fpeven? 6.0))
	(test-assert (not (fpodd? 6.0)))
	(test-assert (fpodd? 7.0))

	(test-assert (flonum? (fprandom)))
	(test-assert (flonum? (fprandom 2456)))

	(test 4.0 (fpadd1 3.0))
	(test 2.0 (fpsub1 3.0))

	(test 27.0 (fpcub 3.0))

	(test 1.0 (fpmodulo 5.0 2.0))
	(test 0.0 (fpmodulo 0.0 1.0))

	(test 2.0 (fpquotient 5.0 2.0))
	(test 1.0 (fpremainder 5.0 2.0))

	(test-assert (fp~= 0.123456 0.123457 5eps))
	(test-assert (fp~<= 0.123456 0.123457 5eps))
	(test-assert (fp~>= 0.123456 0.123457 5eps))
	(test-assert (fp~<= 0.123456 0.12346 5eps))
	(test-assert (fp~>= 0.123456 0.12344 5eps))

  (parameterize ((current-test-epsilon 4eps))
	  (test 5.6568 (fpdistance 1.0 1.0 5.0 5.0)) )

	(receive (mx mn) (fpmax-and-min 1.0 -1.0 -16.0 13.0 2.0 16.0 7.0 -8.0)
	  (test "fpmax-and-min max" 16.0 mx)
	  (test "fpmax-and-min min" -16.0 mn) )
)

;;;

(require-extension fx-utils)

(test-group "FX Utils"

  (test-assert (fxzero? 0))
  (test-assert (not (fxzero? 1)))
  (test-assert (not (fxzero? most-positive-fixnum)))
  (test-assert (not (fxzero? most-negative-fixnum)))

  (test-assert (not (fxpositive? 0)))
  (test-assert (not (fxpositive? most-negative-fixnum)))
  (test-assert (fxpositive? most-positive-fixnum))

  (test-assert (fxcardinal? 0))
  (test-assert (not (fxcardinal? most-negative-fixnum)))
  (test-assert (fxcardinal? most-positive-fixnum))

  (test-assert (not (fxnegative? 0)))
  (test-assert (fxnegative? most-negative-fixnum))
  (test-assert (not (fxnegative? most-positive-fixnum)))

	(test-assert (fixnum? (fxrandom)))
	(test-assert (fixnum? (fxrandom 2456)))

	(test 4 (fxadd1 3))
	(test 2 (fxsub1 3))

	(test 27 (fxcub 3))

	(test 8 (fxpow2log2 3))

	(test 16 (fxdistance 1 1 5 5))

	(receive (mx mn) (fxmax-and-min 1 -1 -16 13 2 16 7 -8)
	  (test "fxmax-and-min max" 16 mx)
	  (test "fxmax-and-min min" -16 mn) )
)

;;;

;(import (prefix mathh-consts C:))
;(require-library mathh-consts)
;=> C:sqrt2 C:degree C:ln2 C:log2e C:e
(require-extension mathh-consts)

(test-group "Math Constants"

	; Well, some
	(test sqrt2 (sqrt 2.0))
	(test degree (/ pi 180.0))
	(test ln2 (log 2.0))
	(test log2e (log2 e))
)

;;;

(test-exit)