mirror of
https://github.com/justinethier/cyclone.git
synced 2025-05-23 20:15:05 +02:00
85 lines
2.7 KiB
Scheme
85 lines
2.7 KiB
Scheme
;;;; 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);")
|
|
))
|