1 //===-- sinf/cosf data tables ---------------------------------------------===// 2 // 3 // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions. 4 // See https://llvm.org/LICENSE.txt for license information. 5 // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception 6 // 7 //===----------------------------------------------------------------------===// 8 9 #include "math_utils.h" 10 #include "sincosf_utils.h" 11 12 #include <stdint.h> 13 14 namespace __llvm_libc { 15 16 // The constants and polynomials for sine and cosine. The 2nd entry 17 // computes -cos (x) rather than cos (x) to get negation for free. 18 const sincos_t __sincosf_table[2] = { 19 {{1.0, -1.0, -1.0, 1.0}, 20 as_double(0x41645f306dc9c883), 21 as_double(0x3ff921fb54442d18), 22 as_double(0x3ff0000000000000), 23 as_double(0xbfdffffffd0c621c), 24 as_double(0x3fa55553e1068f19), 25 as_double(0xbf56c087e89a359d), 26 as_double(0x3ef99343027bf8c3), 27 as_double(0xbfc555545995a603), 28 as_double(0x3f81107605230bc4), 29 as_double(0xbf2994eb3774cf24)}, 30 {{1.0, -1.0, -1.0, 1.0}, 31 as_double(0x41645f306dc9c883), 32 as_double(0x3ff921fb54442d18), 33 as_double(0xbff0000000000000), 34 as_double(0x3fdffffffd0c621c), 35 as_double(0xbfa55553e1068f19), 36 as_double(0x3f56c087e89a359d), 37 as_double(0xbef99343027bf8c3), 38 as_double(0xbfc555545995a603), 39 as_double(0x3f81107605230bc4), 40 as_double(0xbf2994eb3774cf24)}, 41 }; 42 43 // Table with 4/PI to 192 bit precision. To avoid unaligned accesses 44 // only 8 new bits are added per entry, making the table 4 times larger. 45 const uint32_t __inv_pio4[24] = { 46 0xa2, 0xa2f9, 0xa2f983, 0xa2f9836e, 0xf9836e4e, 0x836e4e44, 47 0x6e4e4415, 0x4e441529, 0x441529fc, 0x1529fc27, 0x29fc2757, 0xfc2757d1, 48 0x2757d1f5, 0x57d1f534, 0xd1f534dd, 0xf534ddc0, 0x34ddc0db, 0xddc0db62, 49 0xc0db6295, 0xdb629599, 0x6295993c, 0x95993c43, 0x993c4390, 0x3c439041}; 50 51 } // namespace __llvm_libc 52