;;;; Cyclone Scheme
;;;; https://github.com/justinethier/cyclone
;;;;
;;;; Copyright (c) 2014-2016, Justin Ethier
;;;; All rights reserved.
;;;;
;;;; This module contains the inexact library from r7rs.
;;;;
(define-library (scheme inexact)
  (export 
    acos 
    asin
    atan 
    cos
    exp 
    finite?
    infinite? 
    log
    nan? 
    sin
    sqrt 
    tan
  )
  (begin
    (define-c nan?
      "(void *data, int argc, closure _, object k, object z)"
      " Cyc_check_num(data, z);
        if (obj_is_int(z) || 
            type_of(z) == integer_tag ||
            type_of(z) == bignum_tag ||
            !isnan(((double_type *)z)->value)) 
        {
          return_closcall1(data, k, boolean_f);
        }
        return_closcall1(data, k, boolean_t);")
    (define-c infinite?
      "(void *data, int argc, closure _, object k, object z)"
      " Cyc_check_num(data, z);
        if (obj_is_int(z) || 
            type_of(z) == integer_tag ||
            type_of(z) == bignum_tag ||
            !isinf(((double_type *)z)->value)) 
        {
          return_closcall1(data, k, boolean_f);
        }
        return_closcall1(data, k, boolean_t);")
    (define (finite? z)
      (if (infinite? z) #f #t))
    (define-c acos
      "(void *data, int argc, closure _, object k, object z)"
      " return_inexact_double_op(data, k, acos, z);")
    (define-c asin
      "(void *data, int argc, closure _, object k, object z)"
      " return_inexact_double_op(data, k, asin, z);")
    (define-c atan
      "(void *data, int argc, closure _, object k, object z)"
      " return_inexact_double_op(data, k, atan, z);")
    (define-c cos
      "(void *data, int argc, closure _, object k, object z)"
      " return_inexact_double_op(data, k, cos, z);")
    (define-c exp
      "(void *data, int argc, closure _, object k, object z)"
      " return_inexact_double_op(data, k, exp, z);")
    (define (log z1 . z2)
      (if (null? z2)
          (c-log z1)
          (let ((z2* (car z2)))
            (/ (c-log z1) (c-log z2*)))))
    (define-c c-log
      "(void *data, int argc, closure _, object k, object z)"
      " return_inexact_double_op(data, k, log, z);"
      ;; TODO: this inline form is not stable yet
      ;"(void *data, object ptr, object z)"
      ;" return_inexact_double_op_no_cps(data, ptr, log, z);"
      )
    (define-c sin
      "(void *data, int argc, closure _, object k, object z)"
      " return_inexact_double_op(data, k, sin, z);")
    (define-c sqrt
      "(void *data, int argc, closure _, object k, object z)"
      " return_inexact_double_op(data, k, sqrt, z);")
    (define-c tan
      "(void *data, int argc, closure _, object k, object z)"
      " return_inexact_double_op(data, k, tan, z);")
))