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 #ifdef __declare_extern_inline
43
44 #ifdef _WANT_MATH_ERRNO
45 #include <errno.h>
46 #endif
47
48 __declare_extern_inline(double)
sqrt(double x)49 sqrt (double x)
50 {
51 double result;
52 #ifdef _WANT_MATH_ERRNO
53 if (isless(x, 0.0))
54 errno = EDOM;
55 #endif
56 __asm__ __volatile__ ("fsqrt\t%d0, %d1" : "=w" (result) : "w" (x));
57 return result;
58 }
59
60 __declare_extern_inline(float)
sqrtf(float x)61 sqrtf (float x)
62 {
63 float result;
64 #ifdef _WANT_MATH_ERRNO
65 if (isless(x, 0.0f))
66 errno = EDOM;
67 #endif
68 __asm__ __volatile__ ("fsqrt\t%s0, %s1" : "=w" (result) : "w" (x));
69 return result;
70 }
71
72 __declare_extern_inline(double)
fma(double x,double y,double z)73 fma (double x, double y, double z)
74 {
75 double result;
76 __asm__ __volatile__ ("fmadd\t%d0, %d1, %d2, %d3" : "=w" (result) : "w" (x), "w" (y), "w" (z));
77 return result;
78 }
79
80 __declare_extern_inline(float)
fmaf(float x,float y,float z)81 fmaf (float x, float y, float z)
82 {
83 float result;
84 __asm__ __volatile__ ("fmadd\t%s0, %s1, %s2, %s3" : "=w" (result) : "w" (x), "w" (y), "w" (z));
85 return result;
86 }
87
88 #endif
89
90 #endif /* _MACHINE_MATH_H_ */
91
