math/Makefile

   1 # Copyright (C) 1996-2017 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 include ../Makeconfig
  23
  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 \
  32                    bits/long-double.h
  33
  34 # FPU support code.
  35 aux             := setfpucw fpu_control
  36
  37 # Build the -lm library.
  38
  39 extra-libs      := libm
  40 extra-libs-others = $(extra-libs)
  41
  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       \
  47                fesetmode
  48
  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
  58
  59 libm-calls =                                                              \
  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      \
  63         e_ilogbF                                                          \
  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                              \
  69         w_expF                                                            \
  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)
  79
  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
  87
  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))
  91
  92
  93 # Type specific routine support.
  94 #
  95 # The following three variables control what is included for each type:
  96 #
  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
 100 #
 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
 105
 106 # long double support
 107 type-ldouble-suffix := l
 108 type-ldouble-routines := t_sincosl k_sincosl s_iscanonicall
 109 type-ldouble-yes := ldouble
 110
 111 # double support
 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
 116
 117 # float support
 118 type-float-suffix := f
 119 type-float-routines := k_rem_pio2f
 120
 121
 122 # Apply suffix to each type in arg 1
 123 type-foreach = $(foreach t,$(types),$(subst F,$(type-$(t)-suffix),$(1)))
 124
 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)))  \
 129
 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)
 137 gen-calls = s_ldexpF
 138 generated += $(foreach s,.c .S,$(call type-foreach, $(calls:s_%=m_%$(s))))
 139 routines = $(call type-foreach, $(calls))
 140
 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 \
 147                         $(libm) \
 148                         $(libmvec) \
 149                         $(+force)
 150         (echo '/* GNU ld script'; echo '*/';\
 151          cat $<; \
 152          echo 'GROUP ( $(slibdir)/libm.so$(libm.so-version) ' \
 153               'AS_NEEDED ( $(libdir)/libmvec_nonshared.a $(slibdir)/libmvec.so$(libmvec.so-version) ) )' \
 154         ) > $@.new
 155         mv -f $@.new $@
 156
 157 $(inst_libdir)/libm-$(version).a: $(objpfx)libm.a \
 158                                   $(+force)
 159         $(do-install)
 160
 161 $(inst_libdir)/libm.a: $(common-objpfx)format.lds \
 162                        $(inst_libdir)/libm-$(version).a \
 163                        $(objpfx)../mathvec/libmvec.a \
 164                        $(+force)
 165         (echo '/* GNU ld script'; echo '*/';\
 166          cat $<; \
 167          echo 'GROUP ( $(libdir)/libm-$(version).a $(libdir)/libmvec.a )' \
 168         ) > $@.new
 169         mv -f $@.new $@
 170 endif
 171
 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
 193
 194 ifneq (,$(CXX))
 195 tests += test-math-isinff test-math-iszero
 196 endif
 197
 198 ifneq (no,$(PERL))
 199 libm-vec-tests = $(addprefix test-,$(libmvec-tests))
 200
 201 libm-tests = $(foreach t,$(types),test-$(t) test-$(t)-finite test-i$(t)) \
 202              $(libm-vec-tests)
 203
 204 libm-tests.o = $(addsuffix .o,$(libm-tests))
 205
 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)
 212
 213 ulps-file = $(firstword $(wildcard $(sysdirs:%=%/libm-test-ulps)))
 214
 215 $(objpfx)libm-test-ulps.h: $(ulps-file) gen-libm-test.pl
 216         $(make-target-directory)
 217         $(PERL) gen-libm-test.pl -u $< -H $@
 218
 219 $(objpfx)libm-test.c: libm-test.inc gen-libm-test.pl auto-libm-test-out
 220         $(make-target-directory)
 221         $(PERL) gen-libm-test.pl -c $< -a auto-libm-test-out -C $@
 222
 223 $(objpfx)libm-have-vector-test.h: libm-test.inc gen-libm-have-vector-test.sh
 224         $(make-target-directory)
 225         $(SHELL) gen-libm-have-vector-test.sh > $@
 226 endif
 227
 228 libm-test-fast-math-cflags = -fno-builtin -D__FAST_MATH__ -DTEST_FAST_MATH
 229 libm-test-vec-cflags = $(libm-test-fast-math-cflags) -fno-inline \
 230                        -ffloat-store -D_OPENMP=201307 -Wno-unknown-pragmas
 231
 232 CFLAGS-test-double-vlen2.c = $(libm-test-vec-cflags)
 233
 234 CFLAGS-test-double-vlen4.c = $(libm-test-vec-cflags)
 235 CFLAGS-test-double-vlen4-wrappers.c = $(double-vlen4-arch-ext-cflags)
 236
 237 CFLAGS-test-double-vlen8.c = $(libm-test-vec-cflags)
 238 CFLAGS-test-double-vlen8-wrappers.c = $(double-vlen8-arch-ext-cflags)
 239
 240 CFLAGS-test-float-vlen4.c = $(libm-test-vec-cflags)
 241
 242 CFLAGS-test-float-vlen8.c = $(libm-test-vec-cflags)
 243 CFLAGS-test-float-vlen8-wrappers.c = $(float-vlen8-arch-ext-cflags)
 244
 245 CFLAGS-test-float-vlen16.c = $(libm-test-vec-cflags)
 246 CFLAGS-test-float-vlen16-wrappers.c = $(float-vlen16-arch-ext-cflags)
 247
 248 # The no-inline tests use -fsignaling-nans, and thereby use the
 249 # versions of classification macros that support sNaNs.  The inline
 250 # function tests use the versions of classification macros that may
 251 # raise spurious exceptions for sNaNs, but also do not test for
 252 # exceptions.  Thus both versions of the classification macros are
 253 # validated.
 254 libm-test-no-inline-cflags = -fno-inline -ffloat-store -fno-builtin \
 255                              -fsignaling-nans
 256 libm-test-finite-cflags = $(libm-test-no-inline-cflags) \
 257                           -U__FINITE_MATH_ONLY__ -D__FINITE_MATH_ONLY__=1
 258 CFLAGS-test-float.c = $(libm-test-no-inline-cflags)
 259 CFLAGS-test-float-finite.c = $(libm-test-finite-cflags)
 260 CFLAGS-test-double.c = $(libm-test-no-inline-cflags)
 261 CFLAGS-test-double-finite.c = $(libm-test-finite-cflags)
 262 CFLAGS-test-ldouble.c = $(libm-test-no-inline-cflags)
 263 CFLAGS-test-ldouble-finite.c = $(libm-test-finite-cflags)
 264 CFLAGS-test-tgmath.c = -fno-builtin
 265 # The following testcase uses very long lines (>3 million), so it may take a
 266 # while to compile it. See: http://llvm.org/bugs/show_bug.cgi?id=14106 and
 267 # http://gcc.gnu.org/bugzilla/show_bug.cgi?id=54402
 268 CFLAGS-test-tgmath2.c = -fno-builtin
 269 CFLAGS-test-tgmath-ret.c = -fno-builtin
 270 CFLAGS-test-powl.c = -fno-builtin
 271
 272 CFLAGS-test-snan.c = -fsignaling-nans
 273 CPPFLAGS-test-ifloat.c = -U__LIBC_INTERNAL_MATH_INLINES \
 274                          $(libm-test-fast-math-cflags)
 275 CPPFLAGS-test-idouble.c = -U__LIBC_INTERNAL_MATH_INLINES \
 276                           $(libm-test-fast-math-cflags)
 277 CPPFLAGS-test-ildouble.c = -U__LIBC_INTERNAL_MATH_INLINES \
 278                           $(libm-test-fast-math-cflags)
 279
 280 CFLAGS-test-signgam-finite.c = -ffinite-math-only
 281 CFLAGS-test-signgam-finite-c99.c = -ffinite-math-only -std=c99
 282 CFLAGS-test-signgam-finite-c11.c = -ffinite-math-only -std=c11
 283 CFLAGS-test-signgam-uchar.c = -std=c99
 284 CFLAGS-test-signgam-uchar-init.c = -std=c99
 285 CFLAGS-test-signgam-uchar-static.c = -std=c99
 286 CFLAGS-test-signgam-uchar-init-static.c = -std=c99
 287 CFLAGS-test-signgam-uint.c = -std=c99
 288 CFLAGS-test-signgam-uint-init.c = -std=c99
 289 CFLAGS-test-signgam-uint-static.c = -std=c99
 290 CFLAGS-test-signgam-uint-init-static.c = -std=c99
 291 CFLAGS-test-signgam-ullong.c = -std=c99
 292 CFLAGS-test-signgam-ullong-init.c = -std=c99
 293 CFLAGS-test-signgam-ullong-static.c = -std=c99
 294 CFLAGS-test-signgam-ullong-init-static.c = -std=c99
 295
 296 CFLAGS-test-math-isinff.cc = -std=gnu++11
 297 CFLAGS-test-math-iszero.cc = -std=gnu++11
 298
 299 CFLAGS-test-iszero-excess-precision.c = -fexcess-precision=standard
 300 CFLAGS-test-iseqsig-excess-precision.c = -fexcess-precision=standard
 301 CFLAGS-test-flt-eval-method.c = -fexcess-precision=standard
 302
 303 CFLAGS-test-fe-snans-always-signal.c = -fsignaling-nans
 304
 305 # The -lieee module sets the _LIB_VERSION_ switch to IEEE mode
 306 # for error handling in the -lm functions.
 307 install-lib += libieee.a
 308 non-lib.a += libieee.a
 309 extra-objs += libieee.a ieee-math.o
 310
 311 include ../Rules
 312
 313 gen-all-calls = $(gen-libm-calls) $(gen-calls)
 314
 315 generated += $(addsuffix .c,$(call type-foreach,$(gen-all-calls))) \
 316              gen-libm-templates.stmp
 317
 318 # Create wrappers in the math build directory.
 319 $(objpfx)gen-libm-templates.stmp: Makefile
 320         $(make-target-directory)
 321         for gcall in $(gen-all-calls); do                                 \
 322           func=$${gcall%F*}$${gcall#*F};                                  \
 323           for type in $(foreach t,$(types),$(t)__$(type-$(t)-suffix)); do \
 324             suff=$${type#*__};                                            \
 325             type=$${type%__*};                                            \
 326             file=$(objpfx)$${gcall%F*}$${suff}$${gcall#*F}.c;             \
 327             (                                                             \
 328               echo "#include <math-type-macros-$${type}.h>";              \
 329               echo "#include <$${func}_template.c>";                      \
 330             ) > $${file};                                                 \
 331           done;                                                           \
 332         done;                                                             \
 333         echo > $(@)
 334
 335 # Add dependency to ensure the generator runs prior.
 336 $(foreach t, $(call type-foreach, $(gen-all-calls)), \
 337              $(objpfx)$(t).c): $(objpfx)gen-libm-templates.stmp
 338
 339 ifneq (no,$(PERL))
 340 # This must come after the inclusion of sysdeps Makefiles via Rules.
 341 $(addprefix $(objpfx), $(libm-tests.o)): $(addprefix $(objpfx),\
 342                                                      $(libm-tests-generated))
 343
 344 # Run the math programs to automatically generate ULPs files.
 345 .PHONY: regen-ulps
 346
 347 run-regen-ulps = $(test-wrapper-env) \
 348             $(run-program-env) \
 349             $($*-ENV) $(rtld-prefix) $${run}
 350
 351 regen-ulps: $(addprefix $(objpfx),$(libm-tests))
 352         rm -f $(objpfx)ULPs; rm -f $(objpfx)NewUlps; \
 353         cp $(ulps-file) $(objpfx)libm-test-ulps; \
 354         for run in $^; do \
 355           echo "Regenerating ULPs for $${run}"; \
 356           $(run-regen-ulps) -u -o $(objpfx); \
 357           cat $(objpfx)ULPs >> $(objpfx)libm-test-ulps; \
 358           rm $(objpfx)ULPs; \
 359         done; \
 360         $(PERL) gen-libm-test.pl -n $(objpfx)NewUlps \
 361                                  -u $(objpfx)libm-test-ulps; \
 362         echo "Automatic regeneration of ULPs complete."; \
 363         echo "Difference between the current baseline and the new baseline is:";\
 364         diff -urN $(ulps-file) $(objpfx)NewUlps; \
 365         echo "Copy $(objpfx)NewUlps to $(ulps-file) (relative to source)."
 366 else
 367 regen-ulps:
 368         @echo "Automatic regeneration of ULPs requires perl."; \
 369         exit 1;
 370 endif
 371
 372 # The generated sysd-rules file defines rules like this for sources
 373 # coming from sysdeps/ directories.  These rules find the generic sources.
 374 define o-iterator-doit
 375 $(objpfx)m_%$o: s_%.c $(before-compile); $$(compile-command.c)
 376 endef
 377 object-suffixes-left := $(all-object-suffixes)
 378 include $(o-iterator)
 379
 380 # Likewise, for those generated files shared with libc.
 381 define o-iterator-doit
 382 $(objpfx)m_%$o: $(objpfx)s_%.c $(before-compile); $$(compile-command.c)
 383 endef
 384 object-suffixes-left := $(all-object-suffixes)
 385 include $(o-iterator)
 386
 387 # This file defines the default _LIB_VERSION variable that controls
 388 # the error return conventions for the math functions.
 389 CPPFLAGS-s_lib_version.c := -D_POSIX_MODE
 390
 391 # We don't want the fdlibm code to use the inline math functions,
 392 # only the fdlibm code.
 393 math-CPPFLAGS += -D__NO_MATH_INLINES -D__LIBC_INTERNAL_MATH_INLINES
 394
 395 ifneq ($(long-double-fcts),yes)
 396 # The `double' and `long double' types are the same on this machine.
 397 # We won't compile the `long double' code at all.  Tell the `double' code
 398 # to define aliases for the `FUNCl' names.  To avoid type conflicts in
 399 # defining those aliases, tell <math.h> to declare the `FUNCl' names with
 400 # `double' instead of `long double'.
 401 math-CPPFLAGS += -DNO_LONG_DOUBLE -D_Mlong_double_=double
 402 endif
 403
 404 # These files quiet sNaNs in a way that is optimized away without
 405 # -fsignaling-nans.
 406 CFLAGS-s_modf.c += -fsignaling-nans
 407 CFLAGS-s_modff.c += -fsignaling-nans
 408 CFLAGS-s_modfl.c += -fsignaling-nans
 409
 410 # The -lieee library is actually an object file.
 411 # The module just defines the _LIB_VERSION_ variable.
 412 # It's not a library to make sure it is linked in instead of s_lib_version.o.
 413 $(objpfx)libieee.a: $(objpfx)ieee-math.o
 414         rm -f $@
 415         $(patsubst %/,cd % &&,$(objpfx)) \
 416         $(LN_S) $(<F) $(@F)
 417
 418 $(addprefix $(objpfx),\
 419             $(filter-out $(tests-static) $(libm-vec-tests),$(tests))): $(libm)
 420 $(addprefix $(objpfx),$(tests-static)): $(objpfx)libm.a
 421 $(addprefix $(objpfx), $(libm-vec-tests)): $(objpfx)%: $(objpfx)%-wrappers.o \
 422                                                        $(libm) $(libmvec)
 423
 424 gmp-objs = $(patsubst %,$(common-objpfx)stdlib/%.o,\
 425                       add_n sub_n cmp addmul_1 mul_1 mul_n divmod_1 \
 426                       lshift rshift mp_clz_tab udiv_qrnnd inlines \
 427                       $(gmp-sysdep_routines))
 428 $(objpfx)atest-exp: $(gmp-objs)
 429 $(objpfx)atest-sincos: $(gmp-objs)
 430 $(objpfx)atest-exp2: $(gmp-objs)
 431 $(objpfx)test-fenv-tls: $(shared-thread-library)