1 /* Software floating-point emulation.
2 Helper routine for _Q_* routines.
3 Simulate exceptions using double arithmetics.
4 Copyright (C) 1999-2018 Free Software Foundation, Inc.
5 This file is part of the GNU C Library.
6 Contributed by Jakub Jelinek (jj@ultra.linux.cz).
8 The GNU C Library is free software; you can redistribute it and/or
9 modify it under the terms of the GNU Lesser General Public
10 License as published by the Free Software Foundation; either
11 version 2.1 of the License, or (at your option) any later version.
13 The GNU C Library is distributed in the hope that it will be useful,
14 but WITHOUT ANY WARRANTY; without even the implied warranty of
15 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 Lesser General Public License for more details.
18 You should have received a copy of the GNU Lesser General Public
19 License along with the GNU C Library; if not, see
20 <http://www.gnu.org/licenses/>. */
27 unsigned long long ___Q_zero
= 0x0000000000000000ULL
;
29 void ___Q_simulate_exceptions(int exceptions
)
31 if (exceptions
& FP_EX_INVALID
)
34 __asm__
__volatile__ ("fdivs %0, %0, %0" : "+f" (f
));
36 if (exceptions
& FP_EX_DIVZERO
)
38 float f
= 1.0, g
= 0.0;
39 __asm__
__volatile__ ("fdivs %0, %1, %0"
43 if (exceptions
& FP_EX_OVERFLOW
)
46 __asm__
__volatile__("fmuls %0, %0, %0" : "+f" (f
));
47 exceptions
&= ~FP_EX_INEXACT
;
49 if (exceptions
& FP_EX_UNDERFLOW
)
52 __asm__
__volatile__("fmuls %0, %0, %0" : "+f" (f
));
53 exceptions
&= ~FP_EX_INEXACT
;
55 if (exceptions
& FP_EX_INEXACT
)
57 double d
= 1.0, e
= M_PI
;
58 __asm__
__volatile__ ("fdivd %0, %1, %0"