Changes in Makefile

fa3f69fa4cabf50d667631b4b8083673bccc874f - libmsun: Fix incorrect MLINK for sincosl(3)

Thu, 04 Jun 2026 18:31:26 +0200

libmsun: Fix incorrect MLINK for sincosl(3)PR: 295704MFC after: 1 week List of files: /freebsd/lib/msun/Makefile

3085fc9d97bd83785ba3ba43e0378d7d67987d1f - [libm] implementation of rsqrt, rsqrtf, and rsqrtl

Fri, 08 May 2026 16:06:08 +0200

[libm] implementation of rsqrt, rsqrtf, and rsqrtlFrom the PR:The attached diff implements the inverse square root function, i.e,rsqrt(x) = 1 / sqrt(x). Exhaustive testing of the float versionsuggests that it is correctly rounded in round-to-nearest for alltest values in the range [0x1p-127,0x1p126].Exhaustive testing of rsqrt and rsqrtl cannot be done, but 1100Mvalues of x for rsqrt and 400M values for rsqrtl were tested. Alltested values were correctly rounded.I do not have access to LD128 (i.e., IEEE 128-bit floating point)hardware, so the implementation of rsqrtl() is untested.The following is a summary of changes to source code.* lib/msun/Makefile: . Add s_rsqrt.c and s_rsqrtf.c to COMMON_SRCS. . For non-53-bit long double targets, add s_rsqrtl.c to COMMON_SRCS. . Add MLINKS for rsqrt.3, rsqrtf.3, and rsqrtl.3 to sqrt.3.* lib/msun/Symbol.map: . Add rsqrt, rsqrtf, and rsqrtl to the Symbol map for shared libm.so.* lib/msun/man/sqrt.3: . Update the sqrt.3 manual page to include information for rsqrt[fl]. . Note, these function come from ISO C23 (and IEEE-754 2008).* lib/msun/src/math.h: . Add prototypes for new functions.* lib/msun/src/math_private.h: . Add _SPLIT, _FAST2SUM, _SLOW2SUM, _XADD, _MUL, and _XMUL macros to perform type-type arthimetic (i.e., float-float).* src/s_rsqrt.c: . New file with the implementation of 'double rsqrt(double)'. . For 53-bit long double targets, add a weak reference for rsqrtl.* src/s_rsqrtf.c: . New file with the implementation of 'float rsqrt(float)'.* src/s_rsqrtl.c . New file with the implementation of 'long double rsqrt(long double)'. Note, the LD80 version uses bit twiddling and LD128 version is a straight C language implementation. The LD128 is untested due to lack of hardware.PR: 295089MFC after: 1 week List of files: /freebsd/lib/msun/Makefile

f62d826a6f5b9022b0cedfe22a698998ad9cb7f4 - lib/msun: fmaximum_mag_num family. Tests and man page

Tue, 21 Apr 2026 19:17:00 +0200

lib/msun: fmaximum_mag_num family. Tests and man pageAdded the fmaximum_mag_num{,f,l} and fminimum_mag_num{,f,l} functions.PR: 294719Reviewed by: fuz, karglMFC after: 1 month List of files: /freebsd/lib/msun/Makefile

4e30c12973c73539bef01be3870a90b6346141e8 - lib/msun: Added fmaximum_mag and fmaximum_num families

Sat, 11 Apr 2026 10:11:44 +0200

lib/msun: Added fmaximum_mag and fmaximum_num familiesAdded support for the f{maximum,minimum}_{mag,num} families, the newC23 standard functions for maximum magnitude and number-preferringmaximum. This includes modifying fmax.3, on top of D56230, torecommend the use of fmaximum_num and fminimum_num.Reviewed by: fuz, karglMFC after: 1 monthDifferential Revision: https://reviews.freebsd.org/D56236 List of files: /freebsd/lib/msun/Makefile

c913dce86e7b6ff9c0c3265631529586aaf8520a - lib/msun: centralize addition of the arch-specific symbol map files

Mon, 06 Apr 2026 19:46:45 +0200

lib/msun: centralize addition of the arch-specific symbol map filesThis also adds the aarch64 symbols, exporting them.Reported and tested by: fluffyReviewed by: emasteFixes: 3a01e1e1a50cb9a9594aac2148dc920a6b295428Sponsored by: The FreeBSD FoundationMFC after: 1 weekDifferential revision: https://reviews.freebsd.org/D56283 List of files: /freebsd/lib/msun/Makefile

e55db843ef45a8788f69e110d97210fb3968b92f - lib/msun: Added fmaximum and fminimum family. Tests and man page

Tue, 31 Mar 2026 18:31:15 +0200

lib/msun: Added fmaximum and fminimum family. Tests and man pageStarting from the existing fmax{,f,l} functions I've added the fmaximumfamily, which handles NaN according to the newest standard (propagatingit).This commit is a PoC for GSoC 2026.Reviewed by: fuz, karglMFC after: 1 monthDifferential Revision: https://reviews.freebsd.org/D55834 List of files: /freebsd/lib/msun/Makefile

e38f2308273c8a51ec45f013d22c963590917cca - lib/msun: Fix tgammal(3) on IEEE 128-bit platforms

Fri, 01 Mar 2024 16:53:58 +0100

lib/msun: Fix tgammal(3) on IEEE 128-bit platformsUndo the 80-bit "stub" implementation of the 128-bit long doubletgammal(3) function. The latest (as of Feb 2024) version of thesrc/contrib/arm-optimised-routines library includes a standalone,full 128-bit replacement. This needs a small bit of wrapping tofit it in, but is otherwise a drop-in replacement.Testing this is hard, as most maths packages blow up as soon astheir 80-bit floating-point capability is exceeded. With 128-bittgammal(), this is easy to do, and this is the range that needs tobe checked the most carefully. Using my copy of Maple, I was ableto check that the output was within a few ULP of the correct answer,right up to the point of 128-bit over- and underflow. Additionally,the results are no worse, and indeed better than the 80-bit version.Steve Kargl sent me his libm testing code, which I used to verifythat the excpetions for certain key values were correct. Tested inthis case were +-Inf, +-NaN, +-1 and +-0.Differential Revision: https://reviews.freebsd.org/D44168Reviewed by: theraven, andrew, imp List of files: /freebsd/lib/msun/Makefile

d04e03c19a8d7660ea36f62720eef9bc96b69889 - msun: remove fabs from Symbol.map, and adjust comment

Sat, 27 Jan 2024 22:24:38 +0100

msun: remove fabs from Symbol.map, and adjust commentWe have s_fabs.c, but fabs(3) is already provided by libc due tohistorical reasons, so it is not compiled into libm. When the linkerdoes not use --undefined-version, this leads to a complaint about thesymbol being nonexistent, so remove it from Symbol.map.While here, adjust the comment about some functions being supplied bylibc: while it is true that all these are indeed in libc, libm stillincludes its own versions of frexp(3), isnan(3), isnanf(3), andisnanl(3).Reported by: Steve Kargl MFC after: 3 days List of files: /freebsd/lib/msun/Makefile

dc36d6f9bb1753f3808552f3afd30eda9a7b206a - lib: Remove ancient SCCS tags.

Thu, 23 Nov 2023 18:21:37 +0100

lib: Remove ancient SCCS tags.Remove ancient SCCS tags from the tree, automated scripting, with twominor fixup to keep things compiling. All the common forms in the treewere removed with a perl script.Sponsored by: Netflix List of files: /freebsd/lib/msun/Makefile

b19d8afe4dd205f7d4bfa3282f58b73d504926db - msun: LIBCSRCDIR is too fragile, use ${SRCTOP}/lib/libc instead

Tue, 05 Sep 2023 22:58:06 +0200

msun: LIBCSRCDIR is too fragile, use ${SRCTOP}/lib/libc insteadLIBCSRCDIR is defined in bsd.libnames.mk, which is read in later in theMakefile than the line:.if exists(${LIBCSRCDIR}/${MACHINE_ARCH})so we test to see if /${MARCHIN_ARCH} exists which it usually doesn't(but did for me since I mounted 13.2R SD image there). Move to definingour own LIBC_SRCTOP in terms of SRCTOP to treat these uniformily.Sponsored by: NetflixReviewed by: sjgDifferential Revision: https://reviews.freebsd.org/D41661 List of files: /freebsd/lib/msun/Makefile

d0b2dbfa0ecf2bbc9709efc5e20baf8e4b44bbbf - Remove $FreeBSD$: one-line sh pattern

Wed, 16 Aug 2023 19:55:03 +0200

Remove $FreeBSD$: one-line sh patternRemove /^\s*#[#!]?\s*\$FreeBSD\$.*$\n/ List of files: /freebsd/lib/msun/Makefile

c44ca4c227f7e49bbaf42c014b8bd307b1a94094 - msun/Makefile: Fix conditional and unbreak the build.

Fri, 09 Sep 2022 09:48:48 +0200

msun/Makefile: Fix conditional and unbreak the build.Approved by: impFixes: 20adba8bc1beb125d5e5ed8f12e747ae79ca6a10 List of files: /freebsd/lib/msun/Makefile

20adba8bc1beb125d5e5ed8f12e747ae79ca6a10 - src.conf: Introduce WITHOUT_MACHDEP knob.

Mon, 08 Aug 2022 17:16:00 +0200

src.conf: Introduce WITHOUT_MACHDEP knob.Summary:This knob can be used to make buildsystem prefer generic C implentations ofvarious functions, instead of machine-specific assembler ones.Test Plan: `make buildworld` on amd64Reviewed by: imp, emasteDifferential Revision: https://reviews.freebsd.org/D36076MFC after: 3 days List of files: /freebsd/lib/msun/Makefile

03a88e3de9c68182d21df94b1c8c7ced930dbd1f - * lib/msun/Makefile b/lib/msun/Makefile:

Tue, 14 Dec 2021 10:08:57 +0100

* lib/msun/Makefile b/lib/msun/Makefile: . Disconnect imprecise.c from the build. This file can be deleted. . Add b_tgammal.c to the build for ld80 and ld128 targets. The ld128 is a 'git mv' of imprecise.c to ld128/b_tgammal.c.* lib/msun/ld80/b_expl.c: . New file. Implement __exp__D for ld80 targets. This is based on bsdsrc/b_exp.c.* lib/msun/ld80/b_logl.c: . New file. Implement __log__D for ld80 targets. This is based on bsdsrc/b_log.c.* lib/msun/ld80/b_tgammal.c b/lib/msun/ld80/b_tgammal.c . New file. Implement tgammal(x) for ld80 targets.Submitted by: Steve KarglDifferential Revision: https://reviews.freebsd.org/D33444Reviewed by: pfg List of files: /freebsd/lib/msun/Makefile

455b2ccda3df35f31a167f8a35f11ec31fac89bc - * lib/msun/Makefile:

Tue, 14 Dec 2021 10:02:45 +0100

* lib/msun/Makefile: . Disconnect b_exp.c and b_log.c from the build.* lib/msun/bsdsrc/b_exp.c: . Replace scalb() usage with C99's ldexp(). . Replace finite(x) usage with C99's isfinite(). . Whitespace changes towards style(9). . Remove include of "mathimpl.h". It is no longer needed. . Remove #if 0 ... #endif code, which has been present since svn r93211 (2002-03-26). . New minimax polynomial coefficients. . Add comments to explain origins of some constants. . Use ansi-C prototype. Remove K&R prototype. Add static to prototype.* lib/msun/bsdsrc/b_log.c: . Remove include of "mathimpl.h". It is no longer needed. . Fix comments to actually describe the code. . Reduce minimax polynomial from degree 4 to degree 3. This uses newly computed coefficients. . Use ansi-C prototype. Remove K&R prototype. Add static to prototype. . Remove volatile in declaration of u1. . Alphabetize decalaration list. . Whitespace changes towards style(9). . In argument reduction of x to g and m, replace use of logb() and ldexp() with a single call to frexp(). Add code to get 1 <= g < 2. . Remove #if 0 ... #endif code, which has been present since svn r93211 (2002-03-26). . The special case m == -1022, replace logb() with ilogb().* lib/msun/bsdsrc/b_tgamma.c: . Update comments. Fix comments where needed. . Add float.h to get LDBL_MANT_DIG for weak reference of tgammal to tgamma. . Remove include of "mathimpl.h". It is no longer needed. . Use "math.h" instead of . . Add '#include math_private.h" . Add struct Double from mathimpl.h and include b_log.c and b_exp.c. . Remove forward declarations of neg_gam(), small_gam(), smaller_gam, large_gam() and ratfun_gam() by re-arranging the code to move these function above their first reference. . New minimax coefficients for polynomial in large_gam(). . New splitting of a0 into a0hi nd a0lo, which include additional bits of precision. . Use ansi-C prototype. Remove K&R prototype. . Replace the TRUNC() macro with a simple cast of a double entities to float before assignment (functional changes). . Replace sin(M_PI*z) with sinpi(z) and cos(M_PI*(0.5-z)) with cospi(0.5-z).Submitted by: Steve KarglDifferential Revision: https://reviews.freebsd.org/D33444Reviewed by: pfg List of files: /freebsd/lib/msun/Makefile

b2e843161dc3b79777e873183447c92ed9c3703a - Use a builtin where possible in msun

Tue, 02 Nov 2021 12:31:17 +0100

Use a builtin where possible in msunSome of the functions in msun can be implemented using a compilerbuiltin function to generate a small number of instructions. Implementthis support in fma, fmax, fmin, and sqrt on arm64.Care must be taken as the builtin can be implemented as a functioncall on some architectures that lack direct support. In these caseswe need to use the original code path.As we don't set errno on failure build with -fno-math-errno so thetoolchain doesn't convert a builtin into a function call when itdetects a failure, e.g. gcc will add a call to sqrt when the inputis negative leading to an infinite loop.Sponsored by: The FreeBSD FoundationDifferential Revision: https://reviews.freebsd.org/D32801 List of files: /freebsd/lib/msun/Makefile

046e2d5db1e8afd2d09ea28e5d2a7550535d4b77 - Implementations of cexpl()

Fri, 05 Nov 2021 03:04:01 +0100

Implementations of cexpl()The change implements cexpl() for both ld80 and ld128 architectures.Testing was done on x86_64 and aarch64 systems.Along the way sincos[fl]() use an optimization that reduces the argumentto being done one rather than twice. This optimization actually pointedto a bug in the ld128 version of sincosl(), which is now fixed. Inaddition, the minmax polynomial coefficients for sincosl() have beenupdated.A concise log of the file-by-file changes follows.* include/complex.h: . Add a prototype for cexpl().* lib/msun/Makefile: . Add s_cexpl.c to the build. . Setup a link for cexpl.3 to cexp.3.* lib/msun/Symbol.map: . Expose cexpl symbol in libm shared library.* lib/msun/ld128/s_cexpl.c: * Implementation of cexpl() for 128-bit long double architectures. Tested on an aarch64 system.* lib/msun/ld80/s_cexpl.c: * Implementation of cexpl() for Intel 80-bit long double.* lib/msun/man/cexp.3: . Document cexpl().* lib/msun/man/complex.3: . Add a BUGS section about cpow[fl].* lib/msun/src/s_cexp.c: . Include float.h for weak references on 53-bit long double targets. . Use sincos() to reduce argument reduction cost.* lib/msun/src/s_cexpf.c: . Use sincosf() to reduce argument reduction cost.* lib/msun/src/k_sincosl.h: . Catch up with the new minmax polynomial coefficients for the kernel for the 128-bit cosl() implementation. . BUG FIX: *cs was used where *sn should have been. This means that sinl() was no computed correctly when iy != 0.* lib/msun/src/s_cosl.c: . Include fpmath.h to get access to IEEEl2bits. . Replace M_PI_4 with pio4, a 64-bit or 113-bit approximation for pi / 4.PR: 216862MFC after: 1 week List of files: /freebsd/lib/msun/Makefile

ca3d8cb087cd5b40369478b1693f3e4038b5fa23 - lib/msun: Move the files to appropriate locations in the Makefile

Tue, 26 Oct 2021 22:53:51 +0200

lib/msun: Move the files to appropriate locations in the MakefileFixes: dce5f3abed7181cc533ca5edPR: 218514MFC after: 1 week List of files: /freebsd/lib/msun/Makefile

dce5f3abed7181cc533ca5ed3de44517775e78dd - [LIBM] implementations of sinpi[fl], cospi[fl], and tanpi[fl]

Mon, 25 Oct 2021 15:13:52 +0200

[LIBM] implementations of sinpi[fl], cospi[fl], and tanpi[fl]Both IEEE-754 2008 and ISO/IEC TS 18661-4 define the half-cycletrignometric functions cospi, sinpi, and tanpi. The attachedpatch implements cospi[fl], sinpi[fl], and tanpi[fl]. Limitedtesting on the cospi and sinpi reveal a max ULP less than 0.89;while tanpi is more problematic with a max ULP less than 2.01in the interval [0,0.5]. The algorithms used in these functionsare documented in {ks}_cospi.c, {ks}_sinpi.c, and s_tanpi.c.Note. I no longer have access to a system with ld128 andadequate support to compile and test the ld128 implementationsof these functions. Given the almost complete lack of input fromothers on improvements to libm, I doubt that anyone cares. Ifsomeone does care, the ld128 files contain a number of FIXME comments,and in particular, while the polynomial coefficients are givenI did not update the polynomial algorithms to properly use thecoefficients.PR: 218514MFC after: 2 weeks List of files: /freebsd/lib/msun/Makefile

168234fa67c38f898b784b3265dd77ace0b0a2f9 - lib/msun: Exclude ignored-pragmas from -Werror

Thu, 15 Apr 2021 17:27:52 +0200

lib/msun: Exclude ignored-pragmas from -WerrorThis avoids build failures due to the clang 12 warning: '#pragma FENV_ACCESS' is not supported on this target - ignoredClang 12 currently emits this warning for all non-x86 architectures.While this can result in incorrect code generation (e.g. on AArch64 someexceptions are not raised as expected), this is a pre-existing issue andwe should not fail the build due to this warning.Reviewed By: dim, emasteMFC after: 1 weekDifferential Revision: https://reviews.freebsd.org/D29743 List of files: /freebsd/lib/msun/Makefile