about summary refs log tree commit diff
path: root/sysdeps/x86_64/fpu/svml_d_cos_data.S
blob: c9bfd63840acfd7ac9d5ae5511c76463bbc2d6c8 (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
/* Data for vectorized cos.
   Copyright (C) 2014-2015 Free Software Foundation, Inc.
   This file is part of the GNU C Library.

   The GNU C Library is free software; you can redistribute it and/or
   modify it under the terms of the GNU Lesser General Public
   License as published by the Free Software Foundation; either
   version 2.1 of the License, or (at your option) any later version.

   The GNU C Library is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
   Lesser General Public License for more details.

   You should have received a copy of the GNU Lesser General Public
   License along with the GNU C Library; if not, see
   <http://www.gnu.org/licenses/>.  */

#include "svml_d_cos_data.h"

.macro double_vector offset value
.if .-__svml_dcos_data != \offset
.err
.endif
.rept 8
.quad \value
.endr
.endm

	.section .rodata, "a"
	.align 64

/* Data table for vector implementations of function cos.
   The table may contain polynomial, reduction, lookup
   coefficients and other constants obtained through different
   methods of research and experimental work.
 */
	.globl __svml_dcos_data
__svml_dcos_data:

/* General purpose constants:
   absolute value mask
 */
double_vector __dAbsMask 0x7fffffffffffffff

/* working range threshold */
double_vector __dRangeVal 0x4160000000000000

/* PI/2 */
double_vector __dHalfPI 0x3ff921fb54442d18

/* 1/PI */
double_vector __dInvPI 0x3fd45f306dc9c883

/* right-shifter constant */
double_vector __dRShifter 0x4338000000000000

/* 0.5 */
double_vector __dOneHalf 0x3fe0000000000000

/* Range reduction PI-based constants:
   PI high part
 */
double_vector __dPI1 0x400921fb40000000

/* PI mid  part 1 */
double_vector __dPI2 0x3e84442d00000000

/* PI mid  part 2 */
double_vector __dPI3 0x3d08469880000000

/* PI low  part */
double_vector __dPI4 0x3b88cc51701b839a

/* Range reduction PI-based constants if FMA available:
   PI high part (FMA available)
 */
double_vector __dPI1_FMA 0x400921fb54442d18

/* PI mid part  (FMA available) */
double_vector __dPI2_FMA 0x3ca1a62633145c06

/* PI low part  (FMA available) */
double_vector __dPI3_FMA 0x395c1cd129024e09

/* Polynomial coefficients (relative error 2^(-52.115)): */
double_vector __dC1 0xbfc55555555554a7
double_vector __dC2 0x3f8111111110a4a8
double_vector __dC3 0xbf2a01a019a5b86d
double_vector __dC4 0x3ec71de38030fea0
double_vector __dC5 0xbe5ae63546002231
double_vector __dC6 0x3de60e6857a2f220
double_vector __dC7 0xbd69f0d60811aac8

/*
   Additional constants:
   absolute value mask
 */
double_vector __dAbsMask_la 0x7fffffffffffffff

/* 1/PI */
double_vector __dInvPI_la 0x3fd45f306dc9c883

/* right-shifer for low accuracy version */
double_vector __dRShifter_la 0x4330000000000000

/* right-shifer-1.0 for low accuracy version */
double_vector __dRShifterm5_la 0x432fffffffffffff

/* right-shifer with low mask for low accuracy version */
double_vector __dRXmax_la 0x43300000007ffffe

	.type	__svml_dcos_data,@object
	.size	__svml_dcos_data,.-__svml_dcos_data