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 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 35 float fabsf(float x) { 36 return __builtin_fabsf(x); 37 } 38 39 #if defined(__LP64__) 40 long double fabsl(long double x) { return __builtin_fabsl(x); } 41 #else 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__) 49 float ceilf(float x) { return __builtin_ceilf(x); } 50 double ceil(double x) { return __builtin_ceil(x); } 51 52 float floorf(float x) { return __builtin_floorf(x); } 53 double floor(double x) { return __builtin_floor(x); } 54 55 float fmaf(float x, float y, float z) { return __builtin_fmaf(x, y, z); } 56 double fma(double x, double y, double z) { return __builtin_fma(x, y, z); } 57 58 float fmaxf(float x, float y) { return __builtin_fmaxf(x, y); } 59 double fmax(double x, double y) { return __builtin_fmax(x, y); } 60 61 float fminf(float x, float y) { return __builtin_fminf(x, y); } 62 double fmin(double x, double y) { return __builtin_fmin(x, y); } 63 64 float rintf(float x) { return __builtin_rintf(x); } 65 double rint(double x) { return __builtin_rint(x); } 66 67 float roundf(float x) { return __builtin_roundf(x); } 68 double round(double x) { return __builtin_round(x); } 69 70 float truncf(float x) { return __builtin_truncf(x); } 71 double trunc(double x) { return __builtin_trunc(x); } 72 #endif 73
