1- // ===-- Collection of utils for sinf16/cosf16 ------------------*- C++ -*-===//
1+ // ===-- Collection of utils for sinf16/cosf16 ------------------- *- C++ -*-===//
22//
33// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
44// See https://llvm.org/LICENSE.txt for license information.
55// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
66//
7- // ===----------------------------------------------------------------------===/
7+ // ===----------------------------------------------------------------------===//
88#ifndef LLVM_LIBC_SRC_MATH_GENERIC_SINCOSF16_UTILS_H
99#define LLVM_LIBC_SRC_MATH_GENERIC_SINCOSF16_UTILS_H
1010
1414#include " src/__support/common.h"
1515#include " src/__support/macros/config.h"
1616
17-
1817namespace LIBC_NAMESPACE_DECL {
1918
2019// Lookup table for sin(k * pi / 32) with k = 0, ..., 63.
@@ -46,10 +45,11 @@ LIBC_INLINE int32_t range_reduction(float x, float &y) {
4645 return static_cast <int32_t >(kf);
4746}
4847
49- LIBC_INLINE void sincospif16_eval (float xf, float &sin_k, float &cos_k, float &sin_y, float &cosm1_y) {
48+ LIBC_INLINE void sincospif16_eval (float xf, float &sin_k, float &cos_k,
49+ float &sin_y, float &cosm1_y) {
5050 float y;
5151 int32_t k = range_reduction (xf, y);
52-
52+
5353 sin_k = SIN_K_PI_OVER_32[k & 63 ];
5454 cos_k = SIN_K_PI_OVER_32[(k + 16 ) & 63 ];
5555
@@ -65,7 +65,7 @@ LIBC_INLINE void sincospif16_eval(float xf, float &sin_k, float &cos_k, float &s
6565 // Degree-6 minimax even polynomial for sin(y*pi/32)/y generated by Sollya
6666 // with: > Q = fpminimax(sin(y*pi/32)/y, [|0, 2, 4, 6|], [|SG...|], [0, 0.5]);
6767 sin_y = y * fputil::polyeval (ysq, 0x1 .921fb6p-4f , -0x1 .4aeabcp-13f ,
68- 0x1 .a03354p -21f , -0x1 .ad02d2p -20f );
68+ 0x1 .a03354p -21f , -0x1 .ad02d2p -20f );
6969
7070 // Note that cosm1_y = cos(y*pi/32) - 1 = cos_y - 1
7171 // Derivation:
@@ -75,7 +75,7 @@ LIBC_INLINE void sincospif16_eval(float xf, float &sin_k, float &cos_k, float &s
7575 // Degree-6 minimax even polynomial for cos(y*pi/32) generated by Sollya with:
7676 // > P = fpminimax(cos(y*pi/32), [|0, 2, 4, 6|],[|1, SG...|], [0, 0.5]);
7777 cosm1_y = ysq * fputil::polyeval (ysq, -0x1 .3bd3ccp-8f , 0x1 .03a61ap-18f ,
78- 0x1 .a6f7a2p -29f );
78+ 0x1 .a6f7a2p -29f );
7979}
8080
8181} // namespace LIBC_NAMESPACE_DECL
0 commit comments