1 /*
2 * SPDX-License-Identifier: BSD-3-Clause
3 *
4 * Copyright © 2022 Keith Packard
5 *
6 * Redistribution and use in source and binary forms, with or without
7 * modification, are permitted provided that the following conditions
8 * are met:
9 *
10 * 1. Redistributions of source code must retain the above copyright
11 * notice, this list of conditions and the following disclaimer.
12 *
13 * 2. Redistributions in binary form must reproduce the above
14 * copyright notice, this list of conditions and the following
15 * disclaimer in the documentation and/or other materials provided
16 * with the distribution.
17 *
18 * 3. Neither the name of the copyright holder nor the names of its
19 * contributors may be used to endorse or promote products derived
20 * from this software without specific prior written permission.
21 *
22 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
23 * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
24 * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
25 * FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
26 * COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT,
27 * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
28 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
29 * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
30 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
31 * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
32 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED
33 * OF THE POSSIBILITY OF SUCH DAMAGE.
34 */
35
36 #ifndef _MACHINE_MATH_H_
37 #define _MACHINE_MATH_H_
38
39 #define _HAVE_FAST_FMA 1
40 #define _HAVE_FAST_FMAF 1
41
42 #if defined(_HAVE_ATTRIBUTE_ALWAYS_INLINE) && defined(_HAVE_ATTRIBUTE_GNU_INLINE)
43 #define __declare_aarch64_macro(type) extern __inline type __attribute((gnu_inline, always_inline))
44
45 #ifdef _WANT_MATH_ERRNO
46 #include <errno.h>
47 #endif
48
49 __declare_aarch64_macro(double)
sqrt(double x)50 sqrt (double x)
51 {
52 double result;
53 #ifdef _WANT_MATH_ERRNO
54 if (isless(x, 0.0))
55 errno = EDOM;
56 #endif
57 __asm__ __volatile__ ("fsqrt\t%d0, %d1" : "=w" (result) : "w" (x));
58 return result;
59 }
60
61 __declare_aarch64_macro(float)
sqrtf(float x)62 sqrtf (float x)
63 {
64 float result;
65 #ifdef _WANT_MATH_ERRNO
66 if (isless(x, 0.0f))
67 errno = EDOM;
68 #endif
69 __asm__ __volatile__ ("fsqrt\t%s0, %s1" : "=w" (result) : "w" (x));
70 return result;
71 }
72
73 __declare_aarch64_macro(double)
fma(double x,double y,double z)74 fma (double x, double y, double z)
75 {
76 double result;
77 __asm__ __volatile__ ("fmadd\t%d0, %d1, %d2, %d3" : "=w" (result) : "w" (x), "w" (y), "w" (z));
78 return result;
79 }
80
81 __declare_aarch64_macro(float)
fmaf(float x,float y,float z)82 fmaf (float x, float y, float z)
83 {
84 float result;
85 __asm__ __volatile__ ("fmadd\t%s0, %s1, %s2, %s3" : "=w" (result) : "w" (x), "w" (y), "w" (z));
86 return result;
87 }
88
89 #endif
90
91 #endif /* _MACHINE_MATH_H_ */
92
