cosmopolitan/libc/runtime/cosmo.S
Justine Tunney 5c6877b02b
Introduce support for trapping math
The feenableexcept() and fedisableexcept() APIs are now provided which
let you detect when NaNs appear the moment it happens from anywhere in
your program. Tests have also been added for the mission critical math
functions expf() and erff(), whose perfect operation has been assured.
See examples/trapping.c to see how to use this powerful functionality.
2024-04-30 13:38:43 -07:00

159 lines
5.2 KiB
ArmAsm

/*-*- mode:unix-assembly; indent-tabs-mode:t; tab-width:8; coding:utf-8 -*-│
vi: set noet ft=asm ts=8 sw=8 fenc=utf-8 :vi
Copyright 2020 Justine Alexandra Roberts Tunney
Permission to use, copy, modify, and/or distribute this software for
any purpose with or without fee is hereby granted, provided that the
above copyright notice and this permission notice appear in all copies.
THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL
WARRANTIES WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED
WARRANTIES OF MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE
AUTHOR BE LIABLE FOR ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL
DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR
PROFITS, WHETHER IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER
TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION WITH THE USE OR
PERFORMANCE OF THIS SOFTWARE.
*/
#include "libc/macros.internal.h"
#include "libc/sysv/consts/prot.h"
#include "libc/sysv/consts/map.h"
#include "libc/intrin/strace.internal.h"
#include "libc/sysv/consts/map.h"
#include "libc/dce.h"
#ifdef __x86_64__
.text.startup
// Cosmopolitan runtime.
//
// @param edi is argc
// @param rsi is argv
// @param rdx is environ
// @param rcx is auxv
// @noreturn
cosmo: push %rbp
mov %rsp,%rbp
mov %edi,%r12d
mov %rsi,%r13
mov %rdx,%r14
mov %rcx,%r15
#if SYSDEBUG
call __strace_init
mov %eax,%r12d
#endif /* SYSDEBUG */
#ifndef NOX87
// Windows always initializes FPU to douuble precision.
// WSL breaks Linux ABI by initializing FPU to double precision.
// This code makes long double long again.
//
// @see System V Application Binary Interface NexGen32e Architecture
// Processor Supplement, Version 1.0, December 5th, 2018
// Section 3.4.1: Initial Stack and Register State
fldcw 1f(%rip)
.rodata
.balign 2
// 8087 FPU Control Word
// IM: Invalid Operation
// DM: Denormal Operand
// ZM: Zero Divide
// OM: Overflow
// UM: Underflow
// PM: Precision
// PC: Precision Control
// {float,,double,long double}
// RC: Rounding Control
// {even, -, +, 0}
//
// drr
1: .short 0b00000000000000000001101111111
.previous
#endif
#ifdef __FAST_MATH__
push %rax
stmxcsr (%rsp)
//
// Enable hardware optimizations in violation of the IEEE standard.
//
// - 0x0040 enables "DAZ: Denormals Are Zeros" in MXCSR. This causes the
// processor to turn denormal inputs into zero, before computing them.
// See Intel Manual Vol. 1 §10.2.3.4
//
// - 0x8000 enables "FTZ: Flush To Zero" in MXCSR. This means a floating
// point operation that results in underflow will be set to zero, with
// the same sign, rather than producing a denormalized output. It will
// happen only if underflow trapping hasnt been enabled. See the Intel
// Manual Vol. 1 §10.2.3.3.
//
orl $0x8040,(%rsp)
ldmxcsr (%rsp)
pop %rax
#endif
// run assembly init
call _init
// call constructors
.weak __init_array_end
.weak __init_array_start
mov $__init_array_start,%eax
1: cmp $__init_array_end,%eax
je 2f
push %rax
push %rax
call .Largs
call *(%rax)
pop %rax
pop %rax
add $8,%eax
jmp 1b
// call main()
2: call .Largs
.weak main
call main
xchg %eax,%edi
call exit
.Largs: mov %r12d,%edi
mov %r13,%rsi
mov %r14,%rdx
mov %r15,%rcx
ret
.endfn cosmo,weak
// Enables Thread Local Storage.
.init.start 304,_init_tls
push %rdi
push %rsi
call __enable_tls
pop %rsi
pop %rdi
.init.end 304,_init_tls
#if IsAsan()
.init.start 305,_init_symbols
push %rdi
push %rsi
call __init_symbols
pop %rsi
pop %rdi
.init.end 305,_init_symbols
#endif
#ifdef FTRACE
.init.start 306,_init_ftrace
push %rdi
push %rsi
call ftrace_init
mov %eax,%r12d
pop %rsi
pop %rdi
.init.end 306,_init_ftrace
#endif
#endif /* __x86_64__ */