1 # Copyright (C) 1996-2017 Free Software Foundation, Inc.
2 # This file is part of the GNU C Library.
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.
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.
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/>.
18 # Makefile for the math library.
24 # Installed header files.
25 headers
:= math.h bits
/mathcalls.h bits
/mathinline.h bits
/huge_val.h \
26 bits
/huge_valf.h bits
/huge_vall.h bits
/inf.h bits
/nan.h \
27 fpu_control.h complex.h bits
/cmathcalls.h fenv.h \
28 bits
/fenv.h bits
/fenvinline.h bits
/mathdef.h tgmath.h \
29 bits
/math-finite.h bits
/math-vector.h \
30 bits
/libm-simd-decl-stubs.h bits
/iscanonical.h \
31 bits
/flt-eval-method.h bits
/fp-fast.h bits
/fp-logb.h \
35 aux
:= setfpucw fpu_control
37 # Build the -lm library.
40 extra-libs-others
= $(extra-libs
)
42 libm-support
= s_lib_version s_matherr s_signgam \
43 fclrexcpt fgetexcptflg fraiseexcpt fsetexcptflg \
44 ftestexcept fegetround fesetround fegetenv feholdexcpt \
45 fesetenv feupdateenv t_exp fedisblxcpt feenablxcpt \
46 fegetexcept fesetexcept fetestexceptflag fegetmode \
49 # Wrappers for these functions generated per type using a file named
50 # <func>_template.c and the appropriate math-type-macros-<TYPE>.h.
51 gen-libm-calls
= cargF conjF cimagF crealF cabsF s_cacosF \
52 s_cacoshF s_ccosF s_ccoshF s_casinF s_csinF s_casinhF \
53 k_casinhF s_csinhF k_casinhF s_csinhF s_catanhF s_catanF \
54 s_ctanF s_ctanhF s_cexpF s_clogF s_cprojF s_csqrtF \
55 s_cpowF s_clog10F s_fdimF s_nextdownF s_fmaxF s_fminF \
56 s_nanF s_iseqsigF s_canonicalizeF w_ilogbF w_llogbF \
57 w_log1pF w_scalblnF s_fmaxmagF s_fminmagF
60 e_acosF e_acoshF e_asinF e_atan2F e_atanhF e_coshF e_expF e_fmodF \
61 e_hypotF e_j0F e_j1F e_jnF e_lgammaF_r e_logF e_log10F e_powF \
62 e_rem_pio2F e_remainderF e_scalbF e_sinhF e_sqrtF e_gammaF_r \
64 k_cosF k_sinF k_tanF s_asinhF s_atanF s_cbrtF \
65 s_ceilF s_cosF s_erfF s_expm1F s_fabsF \
66 s_floorF s_log1pF s_logbF \
67 s_nextafterF s_nexttowardF s_rintF s_scalblnF \
68 s_significandF s_sinF s_tanF s_tanhF \
70 w_lgammaF w_lgammaF_r \
71 s_fpclassifyF s_truncF \
72 s_remquoF e_log2F e_exp2F s_roundF s_nearbyintF s_sincosF \
73 s_fmaF s_lrintF s_llrintF s_lroundF s_llroundF e_exp10F \
74 s_issignalingF
$(calls
:s_
%=m_
%) x2y2m1F \
75 gamma_productF lgamma_negF lgamma_productF \
76 s_nextupF s_totalorderF s_totalordermagF s_getpayloadF \
77 s_setpayloadF s_setpayloadsigF s_roundevenF s_fromfpF s_ufromfpF \
78 s_fromfpxF s_ufromfpxF
$(gen-libm-calls
)
80 libm-compat-calls-auto
= \
81 w_acosF_compat w_acoshF_compat w_asinF_compat w_atan2F_compat \
82 w_atanhF_compat w_coshF_compat w_exp2F_compat w_exp10F_compat \
83 w_fmodF_compat w_hypotF_compat w_j0F_compat w_j1F_compat \
84 w_jnF_compat w_log2F_compat w_log10F_compat w_logF_compat \
85 w_powF_compat w_remainderF_compat w_scalbF_compat \
86 w_sinhF_compat w_sqrtF_compat w_tgammaF_compat
88 libm-compat-calls-ldouble-yes
= w_lgamma_compatl k_standardl
89 libm-compat-calls
= w_lgamma_compatf w_lgamma_compat k_standard k_standardf \
90 $(libm-compat-calls-ldouble-
$(long-double-fcts
))
93 # Type specific routine support.
95 # The following three variables control what is included for each type:
97 # type-floatN-suffix = The suffix of the type
98 # type-floatN-routines = Type specific support objects
99 # type-floatN-yes = If the type is supported, evaluates to floatN
101 # Finally, note that types is an intentionally recursive variable.
102 # We only know the full set of supported types for the target machine
103 # after the Rules makefile has been parsed.
104 types
= $(type-ldouble-
$(long-double-fcts
)) double float
106 # long double support
107 type-ldouble-suffix
:= l
108 type-ldouble-routines
:= t_sincosl k_sincosl s_iscanonicall
109 type-ldouble-yes
:= ldouble
112 type-double-suffix
:=
113 type-double-routines
:= branred doasin dosincos halfulp mpa mpatan2 \
114 mpatan mpexp mplog mpsqrt mptan sincos32 slowexp \
115 slowpow sincostab k_rem_pio2
118 type-float-suffix
:= f
119 type-float-routines
:= k_rem_pio2f
122 # Apply suffix to each type in arg 1
123 type-foreach
= $(foreach t
,$(types
),$(subst F
,$(type-
$(t
)-suffix),$(1)))
125 libm-routines
= $(strip $(libm-support
) $(libm-compat-calls
) \
126 $(call type-foreach
, $(libm-compat-calls-auto
)) \
127 $(call type-foreach
, $(libm-calls
)) \
128 $(foreach t
, $(types
), $(type-
$(t
)-routines
))) \
130 # These functions are in libc instead of libm because __printf_fp
131 # calls them, so any program using printf will need them linked in,
132 # and we don't want to have to link every program with -lm.
133 # In libm-calls (above), list m_foo in place of s_foo for any
134 # routine that should be compiled separately for its libc and libm versions.
135 calls
= s_isinfF s_isnanF s_finiteF s_copysignF s_modfF s_scalbnF s_frexpF \
136 s_signbitF
$(gen-calls
)
138 generated
+= $(foreach s
,.c .S
,$(call type-foreach
, $(calls
:s_
%=m_
%$(s
))))
139 routines
= $(call type-foreach
, $(calls
))
141 ifeq ($(build-mathvec
),yes
)
142 # We need to install libm.so and libm.a as linker scripts
143 # for transparent use of vector math library.
144 install-lib-ldscripts
:= libm.so libm.a
145 install-others
= $(inst_libdir
)/libm.so
$(inst_libdir
)/libm.a
146 $(inst_libdir
)/libm.so
: $(common-objpfx
)format.lds \
150 (echo
'/* GNU ld script'; echo
'*/';\
152 echo
'GROUP ( $(slibdir)/libm.so$(libm.so-version) ' \
153 'AS_NEEDED ( $(libdir)/libmvec_nonshared.a $(slibdir)/libmvec.so$(libmvec.so-version) ) )' \
157 $(inst_libdir
)/libm-
$(version
).a
: $(objpfx
)libm.a \
161 $(inst_libdir
)/libm.a
: $(common-objpfx
)format.lds \
162 $(inst_libdir
)/libm-
$(version
).a \
163 $(objpfx
)..
/mathvec
/libmvec.a \
165 (echo
'/* GNU ld script'; echo
'*/';\
167 echo
'GROUP ( $(libdir)/libm-$(version).a $(libdir)/libmvec.a )' \
172 # Rules for the test suite.
173 tests
= test-matherr test-fenv atest-exp atest-sincos atest-exp2 basic-test \
174 test-misc test-fpucw test-fpucw-ieee tst-definitions test-tgmath \
175 test-tgmath-ret bug-nextafter bug-nexttoward bug-tgmath1 \
176 test-tgmath-int test-tgmath2 test-powl tst-CMPLX tst-CMPLX2 test-snan \
177 test-fenv-tls test-fenv-preserve test-fenv-return test-fenvinline \
178 test-nearbyint-except test-fenv-clear test-signgam-finite \
179 test-signgam-finite-c99 test-signgam-finite-c11 \
180 test-nearbyint-except-2 test-signgam-uchar test-signgam-uchar-init \
181 test-signgam-uint test-signgam-uint-init test-signgam-ullong \
182 test-signgam-ullong-init test-nan-overflow test-nan-payload \
183 test-fexcept test-fexcept-traps test-fesetexcept \
184 test-fesetexcept-traps test-fetestexceptflag test-femode \
185 test-femode-traps test-iszero-excess-precision \
186 test-iseqsig-excess-precision test-flt-eval-method \
187 test-fp-ilogb-constants test-fp-llogb-constants \
188 test-fe-snans-always-signal
$(tests-static
)
189 tests-static
= test-fpucw-static test-fpucw-ieee-static \
190 test-signgam-uchar-static test-signgam-uchar-init-static \
191 test-signgam-uint-static test-signgam-uint-init-static \
192 test-signgam-ullong-static test-signgam-ullong-init-static
195 tests
+= test-math-isinff test-math-iszero
199 libm-vec-tests
= $(addprefix test-
,$(libmvec-tests
))
201 libm-tests
= $(foreach t
,$(types
),test-
$(t
) test-
$(t
)-finite test-i
$(t
)) \
204 libm-tests.o
= $(addsuffix .o
,$(libm-tests
))
206 tests
+= $(libm-tests
)
207 libm-vec-test-wrappers
= $(addsuffix -wrappers
, $(libm-vec-tests
))
208 test-extras
+= $(libm-vec-test-wrappers
)
209 extra-test-objs
+= $(addsuffix .o
, $(libm-vec-test-wrappers
))
210 libm-tests-generated
= libm-test-ulps.h libm-have-vector-test.h libm-test.c
211 generated
+= $(libm-tests-generated
) auto-libm-test-out
213 ulps-file
= $(firstword $(wildcard $(sysdirs
:%=%/libm-test-ulps
)))
215 $(objpfx
)libm-test-ulps.h
: $(ulps-file
) gen-libm-test.pl
216 $(make-target-directory
)
217 $(PERL
) gen-libm-test.pl
-u
$< -H
$@
219 libm-test-funcs-auto
= acos acosh asin asinh atan atan2 atanh cabs carg cbrt \
220 ccos ccosh cexp clog clog10 cos cosh cpow csin csinh \
221 csqrt ctan ctanh erf erfc exp exp10 exp2 expm1 fma \
222 hypot j0 j1 jn lgamma log log10 log1p log2 pow sin \
223 sincos sinh sqrt tan tanh tgamma y0 y1 yn
224 auto-libm-test-out-files
= $(addprefix auto-libm-test-out-
,\
225 $(libm-test-funcs-auto
))
227 $(objpfx
)libm-test.c
: libm-test.inc gen-libm-test.pl
$(auto-libm-test-out-files
)
228 $(make-target-directory
)
229 cat
$(auto-libm-test-out-files
) > $(objpfx
)auto-libm-test-out
230 $(PERL
) gen-libm-test.pl
-c
$< -a
$(objpfx
)auto-libm-test-out
-C
$@
232 $(objpfx
)libm-have-vector-test.h
: libm-test.inc gen-libm-have-vector-test.sh
233 $(make-target-directory
)
234 $(SHELL
) gen-libm-have-vector-test.sh
> $@
237 libm-test-fast-math-cflags
= -fno-builtin
-D__FAST_MATH__
-DTEST_FAST_MATH
238 libm-test-vec-cflags
= $(libm-test-fast-math-cflags
) -fno-inline \
239 -ffloat-store
-D_OPENMP
=201307 -Wno-unknown-pragmas
241 CFLAGS-test-double-vlen2.c
= $(libm-test-vec-cflags
)
243 CFLAGS-test-double-vlen4.c
= $(libm-test-vec-cflags
)
244 CFLAGS-test-double-vlen4-wrappers.c
= $(double-vlen4-arch-ext-cflags
)
246 CFLAGS-test-double-vlen8.c
= $(libm-test-vec-cflags
)
247 CFLAGS-test-double-vlen8-wrappers.c
= $(double-vlen8-arch-ext-cflags
)
249 CFLAGS-test-float-vlen4.c
= $(libm-test-vec-cflags
)
251 CFLAGS-test-float-vlen8.c
= $(libm-test-vec-cflags
)
252 CFLAGS-test-float-vlen8-wrappers.c
= $(float-vlen8-arch-ext-cflags
)
254 CFLAGS-test-float-vlen16.c
= $(libm-test-vec-cflags
)
255 CFLAGS-test-float-vlen16-wrappers.c
= $(float-vlen16-arch-ext-cflags
)
257 # The no-inline tests use -fsignaling-nans, and thereby use the
258 # versions of classification macros that support sNaNs. The inline
259 # function tests use the versions of classification macros that may
260 # raise spurious exceptions for sNaNs, but also do not test for
261 # exceptions. Thus both versions of the classification macros are
263 libm-test-no-inline-cflags
= -fno-inline
-ffloat-store
-fno-builtin \
265 libm-test-finite-cflags
= $(libm-test-no-inline-cflags
) \
266 -U__FINITE_MATH_ONLY__
-D__FINITE_MATH_ONLY__
=1
267 CFLAGS-test-float.c
= $(libm-test-no-inline-cflags
)
268 CFLAGS-test-float-finite.c
= $(libm-test-finite-cflags
)
269 CFLAGS-test-double.c
= $(libm-test-no-inline-cflags
)
270 CFLAGS-test-double-finite.c
= $(libm-test-finite-cflags
)
271 CFLAGS-test-ldouble.c
= $(libm-test-no-inline-cflags
)
272 CFLAGS-test-ldouble-finite.c
= $(libm-test-finite-cflags
)
273 CFLAGS-test-tgmath.c
= -fno-builtin
274 # The following testcase uses very long lines (>3 million), so it may take a
275 # while to compile it. See: http://llvm.org/bugs/show_bug.cgi?id=14106 and
276 # http://gcc.gnu.org/bugzilla/show_bug.cgi?id=54402
277 CFLAGS-test-tgmath2.c
= -fno-builtin
278 CFLAGS-test-tgmath-ret.c
= -fno-builtin
279 CFLAGS-test-powl.c
= -fno-builtin
281 CFLAGS-test-snan.c
= -fsignaling-nans
282 CPPFLAGS-test-ifloat.c
= -U__LIBC_INTERNAL_MATH_INLINES \
283 $(libm-test-fast-math-cflags
)
284 CPPFLAGS-test-idouble.c
= -U__LIBC_INTERNAL_MATH_INLINES \
285 $(libm-test-fast-math-cflags
)
286 CPPFLAGS-test-ildouble.c
= -U__LIBC_INTERNAL_MATH_INLINES \
287 $(libm-test-fast-math-cflags
)
289 CFLAGS-test-signgam-finite.c
= -ffinite-math-only
290 CFLAGS-test-signgam-finite-c99.c
= -ffinite-math-only
-std
=c99
291 CFLAGS-test-signgam-finite-c11.c
= -ffinite-math-only
-std
=c11
292 CFLAGS-test-signgam-uchar.c
= -std
=c99
293 CFLAGS-test-signgam-uchar-init.c
= -std
=c99
294 CFLAGS-test-signgam-uchar-static.c
= -std
=c99
295 CFLAGS-test-signgam-uchar-init-static.c
= -std
=c99
296 CFLAGS-test-signgam-uint.c
= -std
=c99
297 CFLAGS-test-signgam-uint-init.c
= -std
=c99
298 CFLAGS-test-signgam-uint-static.c
= -std
=c99
299 CFLAGS-test-signgam-uint-init-static.c
= -std
=c99
300 CFLAGS-test-signgam-ullong.c
= -std
=c99
301 CFLAGS-test-signgam-ullong-init.c
= -std
=c99
302 CFLAGS-test-signgam-ullong-static.c
= -std
=c99
303 CFLAGS-test-signgam-ullong-init-static.c
= -std
=c99
305 CFLAGS-test-math-isinff.
cc = -std
=gnu
++11
306 CFLAGS-test-math-iszero.
cc = -std
=gnu
++11
308 CFLAGS-test-iszero-excess-precision.c
= -fexcess-precision
=standard
309 CFLAGS-test-iseqsig-excess-precision.c
= -fexcess-precision
=standard
310 CFLAGS-test-flt-eval-method.c
= -fexcess-precision
=standard
312 CFLAGS-test-fe-snans-always-signal.c
= -fsignaling-nans
314 # The -lieee module sets the _LIB_VERSION_ switch to IEEE mode
315 # for error handling in the -lm functions.
316 install-lib
+= libieee.a
317 non-lib.a
+= libieee.a
318 extra-objs
+= libieee.a ieee-math.o
322 gen-all-calls
= $(gen-libm-calls
) $(gen-calls
)
324 generated
+= $(addsuffix .c
,$(call type-foreach
,$(gen-all-calls
))) \
325 gen-libm-templates.stmp
327 # Create wrappers in the math build directory.
328 $(objpfx
)gen-libm-templates.stmp
: Makefile
329 $(make-target-directory
)
330 for gcall in
$(gen-all-calls
); do \
331 func
=$${gcall
%F
*}$${gcall
#*F}; \
332 for type in
$(foreach t
,$(types
),$(t
)__
$(type-
$(t
)-suffix)); do \
335 file
=$(objpfx
)$${gcall
%F
*}$${suff}$${gcall
#*F}.c; \
337 echo
"#include <math-type-macros-$${type}.h>"; \
338 echo
"#include <$${func}_template.c>"; \
344 # Add dependency to ensure the generator runs prior.
345 $(foreach t
, $(call type-foreach
, $(gen-all-calls
)), \
346 $(objpfx
)$(t
).c
): $(objpfx
)gen-libm-templates.stmp
349 # This must come after the inclusion of sysdeps Makefiles via Rules.
350 $(addprefix $(objpfx
), $(libm-tests.o
)): $(addprefix $(objpfx
),\
351 $(libm-tests-generated
))
353 # Run the math programs to automatically generate ULPs files.
356 run-regen-ulps
= $(test-wrapper-env
) \
358 $($*-ENV
) $(rtld-prefix
) $${run}
360 regen-ulps
: $(addprefix $(objpfx
),$(libm-tests
))
361 rm -f
$(objpfx
)ULPs
; rm -f
$(objpfx
)NewUlps
; \
362 cp
$(ulps-file
) $(objpfx
)libm-test-ulps
; \
364 echo
"Regenerating ULPs for $${run}"; \
365 $(run-regen-ulps
) -u
-o
$(objpfx
); \
366 cat
$(objpfx
)ULPs
>> $(objpfx
)libm-test-ulps
; \
369 $(PERL
) gen-libm-test.pl
-n
$(objpfx
)NewUlps \
370 -u
$(objpfx
)libm-test-ulps
; \
371 echo
"Automatic regeneration of ULPs complete."; \
372 echo
"Difference between the current baseline and the new baseline is:";\
373 diff
-urN
$(ulps-file
) $(objpfx
)NewUlps
; \
374 echo
"Copy $(objpfx)NewUlps to $(ulps-file) (relative to source)."
377 @echo
"Automatic regeneration of ULPs requires perl."; \
381 # The generated sysd-rules file defines rules like this for sources
382 # coming from sysdeps/ directories. These rules find the generic sources.
383 define o-iterator-doit
384 $(objpfx
)m_
%$o: s_
%.c
$(before-compile
); $$(compile-command.c
)
386 object-suffixes-left
:= $(all-object-suffixes
)
387 include $(o-iterator
)
389 # Likewise, for those generated files shared with libc.
390 define o-iterator-doit
391 $(objpfx
)m_
%$o: $(objpfx
)s_
%.c
$(before-compile
); $$(compile-command.c
)
393 object-suffixes-left
:= $(all-object-suffixes
)
394 include $(o-iterator
)
396 # This file defines the default _LIB_VERSION variable that controls
397 # the error return conventions for the math functions.
398 CPPFLAGS-s_lib_version.c
:= -D_POSIX_MODE
400 # We don't want the fdlibm code to use the inline math functions,
401 # only the fdlibm code.
402 math-CPPFLAGS
+= -D__NO_MATH_INLINES
-D__LIBC_INTERNAL_MATH_INLINES
404 ifneq ($(long-double-fcts
),yes
)
405 # The `double' and `long double' types are the same on this machine.
406 # We won't compile the `long double' code at all. Tell the `double' code
407 # to define aliases for the `FUNCl' names. To avoid type conflicts in
408 # defining those aliases, tell <math.h> to declare the `FUNCl' names with
409 # `double' instead of `long double'.
410 math-CPPFLAGS
+= -DNO_LONG_DOUBLE
-D_Mlong_double_
=double
413 # These files quiet sNaNs in a way that is optimized away without
415 CFLAGS-s_modf.c
+= -fsignaling-nans
416 CFLAGS-s_modff.c
+= -fsignaling-nans
417 CFLAGS-s_modfl.c
+= -fsignaling-nans
419 # The -lieee library is actually an object file.
420 # The module just defines the _LIB_VERSION_ variable.
421 # It's not a library to make sure it is linked in instead of s_lib_version.o.
422 $(objpfx
)libieee.a
: $(objpfx
)ieee-math.o
424 $(patsubst %/,cd
% &&,$(objpfx
)) \
427 $(addprefix $(objpfx
),\
428 $(filter-out $(tests-static
) $(libm-vec-tests
),$(tests
))): $(libm
)
429 $(addprefix $(objpfx
),$(tests-static
)): $(objpfx
)libm.a
430 $(addprefix $(objpfx
), $(libm-vec-tests
)): $(objpfx
)%: $(objpfx
)%-wrappers.o \
433 gmp-objs
= $(patsubst %,$(common-objpfx
)stdlib
/%.o
,\
434 add_n sub_n cmp addmul_1 mul_1 mul_n divmod_1 \
435 lshift rshift mp_clz_tab udiv_qrnnd inlines \
436 $(gmp-sysdep_routines
))
437 $(objpfx
)atest-exp
: $(gmp-objs
)
438 $(objpfx
)atest-sincos
: $(gmp-objs
)
439 $(objpfx
)atest-exp2
: $(gmp-objs
)
440 $(objpfx
)test-fenv-tls
: $(shared-thread-library
)