math: Use accurate answers for cos and sincos.
[platform/upstream/glibc.git] / math / Makefile
1 # Copyright (C) 1996-2013 Free Software Foundation, Inc.
2 # This file is part of the GNU C Library.
3
4 # The GNU C Library is free software; you can redistribute it and/or
5 # modify it under the terms of the GNU Lesser General Public
6 # License as published by the Free Software Foundation; either
7 # version 2.1 of the License, or (at your option) any later version.
8
9 # The GNU C Library is distributed in the hope that it will be useful,
10 # but WITHOUT ANY WARRANTY; without even the implied warranty of
11 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
12 # Lesser General Public License for more details.
13
14 # You should have received a copy of the GNU Lesser General Public
15 # License along with the GNU C Library; if not, see
16 # <http://www.gnu.org/licenses/>.
17
18 # Makefile for the math library.
19
20 subdir          := math
21
22 # Installed header files.
23 headers         := math.h bits/mathcalls.h bits/mathinline.h bits/huge_val.h \
24                    bits/huge_valf.h bits/huge_vall.h bits/inf.h bits/nan.h \
25                    fpu_control.h complex.h bits/cmathcalls.h fenv.h \
26                    bits/fenv.h bits/fenvinline.h bits/mathdef.h tgmath.h \
27                    bits/math-finite.h
28
29 # FPU support code.
30 aux             := setfpucw fpu_control
31
32 # Build the -lm library.
33
34 extra-libs      := libm
35 extra-libs-others = $(extra-libs)
36
37 libm-support = k_standard s_lib_version s_matherr s_signgam             \
38                fclrexcpt fgetexcptflg fraiseexcpt fsetexcptflg          \
39                ftestexcept fegetround fesetround fegetenv feholdexcpt   \
40                fesetenv feupdateenv t_exp fedisblxcpt feenablxcpt       \
41                fegetexcept
42 libm-calls = e_acos e_acosh e_asin e_atan2 e_atanh e_cosh e_exp e_fmod  \
43              e_hypot e_j0 e_j1 e_jn e_lgamma_r e_log e_log10 e_pow      \
44              e_rem_pio2 e_remainder e_scalb e_sinh e_sqrt e_gamma_r     \
45              e_ilogb                                                    \
46              k_cos k_rem_pio2 k_sin k_tan s_asinh s_atan s_cbrt         \
47              s_ceil s_cos s_erf s_expm1 s_fabs                          \
48              s_floor s_log1p s_logb                                     \
49              s_nextafter s_nexttoward s_rint s_scalbln                  \
50              s_significand s_sin s_tan s_tanh w_acos w_acosh w_asin     \
51              w_atan2 w_atanh w_cosh w_drem w_exp w_exp2 w_exp10 w_fmod  \
52              w_tgamma w_hypot w_j0 w_j1 w_jn w_lgamma w_lgamma_r        \
53              w_log w_log10 w_pow w_remainder w_scalb w_sinh w_sqrt      \
54              w_ilogb                                                    \
55              s_fpclassify s_fmax s_fmin s_fdim s_nan s_trunc            \
56              s_remquo e_log2 e_exp2 s_round s_nearbyint s_sincos        \
57              conj cimag creal cabs carg s_cexp s_csinh s_ccosh s_clog   \
58              s_catan s_casin s_ccos s_csin s_ctan s_ctanh s_cacos       \
59              s_casinh s_cacosh s_catanh s_csqrt s_cpow s_cproj s_clog10 \
60              s_fma s_lrint s_llrint s_lround s_llround e_exp10 w_log2   \
61              s_isinf_ns s_issignaling $(calls:s_%=m_%) x2y2m1 k_casinh
62
63 include ../Makeconfig
64
65 dbl-only-routines := branred doasin dosincos halfulp mpa mpatan2        \
66                      mpatan mpexp mplog mpsqrt mptan sincos32 slowexp   \
67                      slowpow sincostab
68 libm-routines = $(strip $(libm-support) $(libm-calls) \
69                         $(patsubst %_rf,%f_r,$(libm-calls:=f))  \
70                         $(long-m-$(long-double-fcts))) \
71                 $(dbl-only-routines)
72 long-m-routines = $(patsubst %_rl,%l_r,$(libm-calls:=l))
73 long-m-support = t_sincosl k_sincosl
74 long-m-yes = $(long-m-routines) $(long-m-support)
75
76 # These functions are in libc instead of libm because __printf_fp
77 # calls them, so any program using printf will need them linked in,
78 # and we don't want to have to link every program with -lm.
79 # In libm-calls (above), list m_foo in place of s_foo for any
80 # routine that should be compiled separately for its libc and libm versions.
81 calls = s_isinf s_isnan s_finite s_copysign s_modf s_scalbn s_frexp s_ldexp \
82         s_signbit
83 generated += $(foreach s,.c .S l.c l.S f.c f.S,$(calls:s_%=m_%$s))
84 routines = $(calls) $(calls:=f) $(long-c-$(long-double-fcts))
85 long-c-yes = $(calls:=l)
86
87 # Rules for the test suite.
88 tests = test-matherr test-fenv atest-exp atest-sincos atest-exp2 basic-test \
89         test-misc test-fpucw tst-definitions test-tgmath test-tgmath-ret \
90         bug-nextafter bug-nexttoward bug-tgmath1 test-tgmath-int \
91         test-tgmath2 test-powl tst-CMPLX tst-CMPLX2 test-snan
92 # We do the `long double' tests only if this data type is available and
93 # distinct from `double'.
94 test-longdouble-yes = test-ldouble test-ildoubl
95
96 ifneq (no,$(PERL))
97 libm-tests = test-float test-double $(test-longdouble-$(long-double-fcts)) \
98         test-ifloat test-idouble
99 libm-tests.o = $(addsuffix .o,$(libm-tests))
100
101 tests += $(libm-tests)
102 libm-tests-generated = libm-test-ulps.h libm-test.c
103 generated += $(libm-tests-generated) libm-test.stmp
104
105 # This is needed for dependencies
106 before-compile += $(objpfx)libm-test.c
107
108 ulps-file = $(firstword $(wildcard $(sysdirs:%=%/libm-test-ulps)))
109
110 $(addprefix $(objpfx), $(libm-tests-generated)): $(objpfx)libm-test.stmp
111
112 $(objpfx)libm-test.stmp: $(ulps-file) libm-test.inc gen-libm-test.pl
113         $(make-target-directory)
114         $(PERL) gen-libm-test.pl -u $< -o "$(objpfx)"
115         @echo > $@
116
117 $(objpfx)test-float.o: $(objpfx)libm-test.stmp
118 $(objpfx)test-ifloat.o: $(objpfx)libm-test.stmp
119 $(objpfx)test-double.o: $(objpfx)libm-test.stmp
120 $(objpfx)test-idouble.o: $(objpfx)libm-test.stmp
121 $(objpfx)test-ldouble.o: $(objpfx)libm-test.stmp
122 $(objpfx)test-ildoubl.o: $(objpfx)libm-test.stmp
123 endif
124
125 CFLAGS-test-float.c = -fno-inline -ffloat-store -fno-builtin -frounding-math
126 CFLAGS-test-double.c = -fno-inline -ffloat-store -fno-builtin -frounding-math
127 CFLAGS-test-ldouble.c = -fno-inline -ffloat-store -fno-builtin -frounding-math
128 CFLAGS-test-tgmath.c = -fno-builtin
129 # The following testcase uses very long lines (>3 million), so it may take a
130 # while to compile it. See: http://llvm.org/bugs/show_bug.cgi?id=14106 and
131 # http://gcc.gnu.org/bugzilla/show_bug.cgi?id=54402
132 CFLAGS-test-tgmath2.c = -fno-builtin
133 CFLAGS-test-tgmath-ret.c = -fno-builtin
134 CFLAGS-test-powl.c = -fno-builtin
135 CFLAGS-test-test-fenv.c = -frounding-math
136 CFLAGS-test-misc.c = -frounding-math
137 CPPFLAGS-test-ifloat.c = -U__LIBC_INTERNAL_MATH_INLINES -D__FAST_MATH__ \
138                          -DTEST_FAST_MATH -fno-builtin -frounding-math
139 CPPFLAGS-test-idouble.c = -U__LIBC_INTERNAL_MATH_INLINES -D__FAST_MATH__ \
140                          -DTEST_FAST_MATH -fno-builtin -frounding-math
141 CPPFLAGS-test-ildoubl.c = -U__LIBC_INTERNAL_MATH_INLINES -D__FAST_MATH__ \
142                           -DTEST_FAST_MATH -fno-builtin -frounding-math
143
144
145 # The -lieee module sets the _LIB_VERSION_ switch to IEEE mode
146 # for error handling in the -lm functions.
147 install-lib += libieee.a
148 non-lib.a += libieee.a
149 extra-objs += libieee.a ieee-math.o
150
151 include ../Rules
152
153 ifneq (no,$(PERL))
154 # Run the math programs to automatically generate ULPs files.
155 .PHONY: regen-ulps
156
157 run-regen-ulps = $(test-wrapper-env) \
158             GCONV_PATH=$(common-objpfx)iconvdata LC_ALL=C \
159             $($*-ENV) $(rtld-prefix) $${run}
160
161 regen-ulps: $(addprefix $(objpfx),$(libm-tests))
162         rm -f $(objpfx)ULPs; rm -f $(objpfx)NewUlps; \
163         cp $(ulps-file) $(objpfx)libm-test-ulps; \
164         for run in $^; do \
165           echo "Regenerating ULPs for $${run}"; \
166           $(run-regen-ulps) -u -o $(objpfx); \
167           cat $(objpfx)ULPs >> $(objpfx)libm-test-ulps; \
168           rm $(objpfx)ULPs; \
169         done; \
170         $(PERL) gen-libm-test.pl -o $(objpfx) -n -u $(objpfx)libm-test-ulps; \
171         echo "Automatic regeneration of ULPs complete."; \
172         echo "Difference between the current baseline and the new baseline is:";\
173         diff -urN $(ulps-file) $(objpfx)NewUlps; \
174         echo "Copy $(objpfx)NewUlps to $(ulps-file) (relative to source)."
175 else
176 regen-ulps:
177         @echo "Automatic regeneration of ULPs requires perl."; \
178         exit 1;
179 endif
180
181 # The generated sysd-rules file defines rules like this for sources
182 # coming from sysdeps/ directories.  These rules find the generic sources.
183 define o-iterator-doit
184 $(objpfx)m_%$o: s_%.c $(before-compile); $$(compile-command.c)
185 endef
186 object-suffixes-left := $(all-object-suffixes)
187 include $(o-iterator)
188
189
190 # This file defines the default _LIB_VERSION variable that controls
191 # the error return conventions for the math functions.
192 CPPFLAGS-s_lib_version.c := -D_POSIX_MODE
193
194 # We don't want the fdlibm code to use the inline math functions,
195 # only the fdlibm code.
196 math-CPPFLAGS += -D__NO_MATH_INLINES -D__LIBC_INTERNAL_MATH_INLINES
197
198 ifneq ($(long-double-fcts),yes)
199 # The `double' and `long double' types are the same on this machine.
200 # We won't compile the `long double' code at all.  Tell the `double' code
201 # to define aliases for the `FUNCl' names.  To avoid type conflicts in
202 # defining those aliases, tell <math.h> to declare the `FUNCl' names with
203 # `double' instead of `long double'.
204 math-CPPFLAGS += -DNO_LONG_DOUBLE -D_Mlong_double_=double
205 endif
206
207 # The fdlibm code generates a lot of these warnings but is otherwise clean.
208 override CFLAGS += -Wno-uninitialized
209
210 # The -lieee library is actually an object file.
211 # The module just defines the _LIB_VERSION_ variable.
212 # It's not a library to make sure it is linked in instead of s_lib_version.o.
213 $(objpfx)libieee.a: $(objpfx)ieee-math.o
214         rm -f $@
215         $(patsubst %/,cd % &&,$(objpfx)) \
216         $(LN_S) $(<F) $(@F)
217
218 ifeq ($(build-shared),yes)
219 $(addprefix $(objpfx),$(tests)): $(objpfx)libm.so$(libm.so-version)
220 else
221 $(addprefix $(objpfx),$(tests)): $(objpfx)libm.a
222 endif
223
224 gmp-objs = $(patsubst %,$(common-objpfx)stdlib/%.o,\
225                       add_n sub_n cmp addmul_1 mul_1 mul_n divmod_1 \
226                       lshift rshift mp_clz_tab udiv_qrnnd inlines \
227                       $(gmp-sysdep_routines))
228 $(objpfx)atest-exp: $(gmp-objs)
229 $(objpfx)atest-sincos: $(gmp-objs)
230 $(objpfx)atest-exp2: $(gmp-objs)
231
232 # Depend on libc.so so a DT_NEEDED is generated in the shared objects.
233 # This ensures they will load libc.so for needed symbols if loaded by
234 # a statically-linked program that hasn't already loaded it.
235 $(objpfx)libm.so: $(common-objpfx)libc.so $(common-objpfx)libc_nonshared.a