Files
l2/stdlib/float.sl
2026-02-18 13:58:08 +01:00

111 lines
1.9 KiB
Plaintext

# L2 Floating Point Library (Double Precision)
# Arithmetic
#f+ [*, x1 | x2] -> [* | x3]
:asm f+ {
movq xmm0, [r12]
add r12, 8
movq xmm1, [r12]
addsd xmm1, xmm0
movq [r12], xmm1
} ;
#f- [*, x1 | x2] -> [* | x3]
:asm f- {
movq xmm0, [r12]
add r12, 8
movq xmm1, [r12]
subsd xmm1, xmm0
movq [r12], xmm1
} ;
#f* [*, x1 | x2] -> [* | x3]
:asm f* {
movq xmm0, [r12]
add r12, 8
movq xmm1, [r12]
mulsd xmm1, xmm0
movq [r12], xmm1
} ;
#f/ [*, x1 | x2] -> [* | x3]
:asm f/ {
movq xmm0, [r12]
add r12, 8
movq xmm1, [r12]
divsd xmm1, xmm0
movq [r12], xmm1
} ;
#fneg [* | x] -> [* | -x]
:asm fneg {
movq xmm0, [r12]
mov rax, 0x8000000000000000
movq xmm1, rax
xorpd xmm0, xmm1
movq [r12], xmm0
} ;
# Comparison
#f== [*, x1 | x2] -> [* | flag]
:asm f== {
movq xmm0, [r12]
add r12, 8
movq xmm1, [r12]
ucomisd xmm0, xmm1
mov rax, 0
setz al
mov [r12], rax
} ;
#f< [*, x1 | x2] -> [* | flag]
:asm f< {
movq xmm0, [r12] ; a
add r12, 8
movq xmm1, [r12] ; b
ucomisd xmm0, xmm1
mov rax, 0
seta al ; Above (a > b) -> b < a
mov [r12], rax
} ;
#f> [*, x1 | x2] -> [* | flag]
:asm f> {
movq xmm0, [r12] ; a
add r12, 8
movq xmm1, [r12] ; b
ucomisd xmm1, xmm0
mov rax, 0
seta al ; b > a
mov [r12], rax
} ;
# Conversion
#int>float [* | x] -> [* | xf]
:asm int>float {
cvtsi2sd xmm0, [r12]
movq [r12], xmm0
} ;
#float>int [* | xf] -> [* | x]
:asm float>int {
cvttsd2si rax, [r12]
mov [r12], rax
} ;
# Output
# extern declarations are required for runtime linking
extern int printf(char* fmt, double x)
extern int fflush(void* stream)
#fput [* | xf] -> [*]
word fput
"%f" drop swap printf drop
0 fflush drop
end
#fputln [* | xf] -> [*]
word fputln
"%f\n" drop swap printf drop
0 fflush drop
end