1 /*
2 * Copyright (C) 2015 The Android Open Source Project
3 *
4 * Licensed under the Apache License, Version 2.0 (the "License");
5 * you may not use this file except in compliance with the License.
6 * You may obtain a copy of the License at
7 *
8 * http://www.apache.org/licenses/LICENSE-2.0
9 *
10 * Unless required by applicable law or agreed to in writing, software
11 * distributed under the License is distributed on an "AS IS" BASIS,
12 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 * See the License for the specific language governing permissions and
14 * limitations under the License.
15 */
16
17 #include <math.h>
18
19 #include "fpmath.h"
20
fabs(double x)21 double fabs(double x) {
22 #if __arm__
23 // Both Clang and GCC insist on moving r0/r1 into a double register
24 // and using fabs where bit-twiddling would be a better choice.
25 // They get fabsf right, but we need to be careful in fabsl too.
26 IEEEd2bits u;
27 u.d = x;
28 u.bits.sign = 0;
29 return u.d;
30 #else
31 return __builtin_fabs(x);
32 #endif
33 }
34
fabsf(float x)35 float fabsf(float x) {
36 return __builtin_fabsf(x);
37 }
38
39 #if defined(__LP64__)
fabsl(long double x)40 long double fabsl(long double x) { return __builtin_fabsl(x); }
41 #else
fabsl(long double x)42 long double fabsl(long double x) {
43 // Don't use __builtin_fabs here because of ARM. (See fabs above.)
44 return fabs(x);
45 }
46 #endif
47
48 #if defined(__aarch64__)
ceilf(float x)49 float ceilf(float x) { return __builtin_ceilf(x); }
ceil(double x)50 double ceil(double x) { return __builtin_ceil(x); }
51
floorf(float x)52 float floorf(float x) { return __builtin_floorf(x); }
floor(double x)53 double floor(double x) { return __builtin_floor(x); }
54
fmaf(float x,float y,float z)55 float fmaf(float x, float y, float z) { return __builtin_fmaf(x, y, z); }
fma(double x,double y,double z)56 double fma(double x, double y, double z) { return __builtin_fma(x, y, z); }
57
fmaxf(float x,float y)58 float fmaxf(float x, float y) { return __builtin_fmaxf(x, y); }
fmax(double x,double y)59 double fmax(double x, double y) { return __builtin_fmax(x, y); }
60
fminf(float x,float y)61 float fminf(float x, float y) { return __builtin_fminf(x, y); }
fmin(double x,double y)62 double fmin(double x, double y) { return __builtin_fmin(x, y); }
63
rintf(float x)64 float rintf(float x) { return __builtin_rintf(x); }
rint(double x)65 double rint(double x) { return __builtin_rint(x); }
66
roundf(float x)67 float roundf(float x) { return __builtin_roundf(x); }
round(double x)68 double round(double x) { return __builtin_round(x); }
69
truncf(float x)70 float truncf(float x) { return __builtin_truncf(x); }
trunc(double x)71 double trunc(double x) { return __builtin_trunc(x); }
72 #endif
73