Fcalva/cyclone
1
0
Fork 0
mirror of https://github.com/justinethier/cyclone.git synced 2025-05-19 05:39:17 +02:00
Code Issues Projects Releases Packages Wiki Activity Actions

Allow inline sqrt

Browse source
This commit is contained in:
Justin Ethier 2024-03-11 19:29:44 -07:00
parent fa6213b907
commit a2568d8589
1 changed files with 23 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

24
scheme/inexact.sld
View file

@ -120,6 +120,28 @@
} " } "
; "(void *data, object ptr, object z)" ; "(void *data, object ptr, object z)"
; " return_inexact_double_or_cplx_op_no_cps(data, ptr, sqrt, csqrt, z);" ; " return_inexact_double_or_cplx_op_no_cps(data, ptr, sqrt, csqrt, z);"
) "(void *data, object ptr, object z)"
" double complex unboxed;
Cyc_check_num(data, z);
if (obj_is_int(z)) {
unboxed = csqrt(obj_obj2int(z));
} else if (type_of(z) == integer_tag) {
unboxed = csqrt(((integer_type *)z)->value);
} else if (type_of(z) == bignum_tag) {
unboxed = csqrt(mp_get_double(&bignum_value(z)));
} else if (type_of(z) == complex_num_tag) {
unboxed = csqrt(complex_num_value(z));
assign_complex_num(ptr, unboxed);
return ptr;
} else {
unboxed = csqrt(((double_type *)z)->value);
}
if (cimag(unboxed) == 0.0) {
assign_double(ptr, creal(unboxed));
} else {
assign_double(ptr, unboxed);
}
return ptr; ")
)) ))