Changes in Makefile

b45fd24a3ca666fbb1c754d9de94a9b33fa7949e - lib/msun: move the calculation of LIBC_SRCTOP earlier

Mon, 07 Sep 2026 01:23:45 +0200

lib/msun: move the calculation of LIBC_SRCTOP earlierwhich makes the variable available for machine/Makefile.inc usage.Reviewed by: fuz, kfvSponsored by: The FreeBSD FoundationMFC after: 1 weekDifferential revision: https://reviews.freebsd.org/D59462 List of files: /freebsd/lib/msun/Makefile

d8b6461f8c56237a00ddc32678b3928bf76f5440 - libm: Implement femode_t, fegetmode(), and fesetmode() as per C23

Thu, 03 Sep 2026 01:35:22 +0200

libm: Implement femode_t, fegetmode(), and fesetmode() as per C23Reviewed by: kargl, kibApproved by: fuz (mentor)MFC after: 1 monthDifferential Revision: https://reviews.freebsd.org/D59288 List of files: /freebsd/lib/msun/Makefile

ae417b3194e76ce26065dc20281493ee83619879 - msun: add asinpi, acospi, and atanpi

Wed, 29 Jul 2026 13:45:42 +0200

msun: add asinpi, acospi, and atanpiThis commit implements the inverse half-cycletrigonometric functions: asinpi(x) = asin(x) / pi Eq. (1) acospi(x) = acos(x) / pi atanpi(x) = atan(x) / piImplemention details are contained in src/s_asinpi.c andsrc/a_atanpi.c, where the details for acospi(x) appear inthe former.*************CAVEAT EMPTOR: The ld128 code has been only compiled. It hasnot been tested for correctness due to lack of hardware.*************Code compiled on AMD Ryzen 7 7700X system run FreeBSD 16.0-CURRENTmain-n284956-de9fe28ab847.Exhaustive testing of acospif(x), asinpif(x), and atanpif(x)on the indicated intervals yields% ./tlibm acospi -fPE -x 0x1p-120 -X 1Interval tested for acospif: [7.52316e-37,1] ulp <= 0.5: 99.627% 1002878299 | 99.627% 10028782990.5 < ulp <= 0.6: 0.277% 2789599 | 99.904% 10056678980.6 < ulp <= 0.7: 0.096% 965062 | 100.000% 1006632960Max ulp: 0.63661975 at 5.96046412e-08 0x1.fffffep-25% ./tlibm asinpi -fPED -x 0x1p-120f -X 1.fInterval tested for asinpif: [7.52316e-37,1] ulp <= 0.5: 99.851% 1005129353 | 99.851% 10051293530.5 < ulp <= 0.6: 0.149% 1501097 | 100.000% 10066304500.6 < ulp <= 0.7: 0.000% 2510 | 100.000% 1006632960Max ulp: 0.68957579 at 5.04878759e-01 0x1.027f78p-1% ./tlibm atanpi -fPE -x 0x1p-120 -X max > zatanpif.txt &Interval tested for atanpif: [7.52316e-37,3.40282e+38] ulp <= 0.5: 99.865% 2077574602 | 99.865% 20775746020.5 < ulp <= 0.6: 0.131% 2735011 | 99.997% 20803096130.6 < ulp <= 0.7: 0.003% 65170 | 100.000% 2080374783Max ulp: 0.68433094 at 5.01186252e-01 0x1.009b7cp-1Testing the double and long double version cannot be donein an exhaustive manner. For 300 M values, uniformilydistributed in the indicated interals, one finds the max ULP:Interval tested for acospi: [9.31323e-10,0.25] xm = 2.4423788416892520e-01, /* 0x3fcf432f, 0xde79920f */libm = 4.2146222480005391e-01, /* 0x3fdaf93c, 0xb201001c */mpfr = 4.2146222480005396e-01, /* 0x3fdaf93c, 0xb201001d */ ULP = 0.50499351466286857Interval tested for acospi: [0.25,0.5] xm = 4.9689430915631438e-01, /* 0x3fdfcd1d, 0xc9d945c6 */libm = 3.3447366122373884e-01, /* 0x3fd56804, 0x371513ef */mpfr = 3.3447366122373889e-01, /* 0x3fd56804, 0x371513f0 */ ULP = 0.57195275455053829Interval tested for acospi: [0.5,0.75] xm = 5.0238623667462079e-01, /* 0x3fe0138c, 0x4d0f4be0 */libm = 3.3245556599062825e-01, /* 0x3fd546f3, 0xb5d36303 */mpfr = 3.3245556599062820e-01, /* 0x3fd546f3, 0xb5d36302 */ ULP = 0.63427929243758807Interval tested for acospi: [0.75,1] xm = 7.5853651919512177e-01, /* 0x3fe845ee, 0x60d8789f */libm = 2.2591472240382732e-01, /* 0x3fcceac6, 0x0c3465ce */mpfr = 2.2591472240382729e-01, /* 0x3fcceac6, 0x0c3465cd */ ULP = 0.56915750216472161Interval tested for asinpi: [9.31323e-10,0.25] xm = 1.9502362835488171e-01, /* 0x3fc8f688, 0xc4dda0fb */libm = 6.2478354989018887e-02, /* 0x3faffd29, 0xb6c57c61 */mpfr = 6.2478354989018881e-02, /* 0x3faffd29, 0xb6c57c60 */ ULP = 0.52347765415885006Interval tested for asinpi: [0.25,0.5] xm = 4.9937103583123676e-01, /* 0x3fdff5b1, 0xeeddbf62 */libm = 1.6643553767987129e-01, /* 0x3fc54dc2, 0x7b9d15a4 */mpfr = 1.6643553767987126e-01, /* 0x3fc54dc2, 0x7b9d15a3 */ ULP = 0.66214688371031072Interval tested for asinpi: [0.5,0.75] xm = 5.0228515250761718e-01, /* 0x3fe012b8, 0x4fe92bbb */libm = 1.6750722213679006e-01, /* 0x3fc570e0, 0x6c75edd5 */mpfr = 1.6750722213679009e-01, /* 0x3fc570e0, 0x6c75edd6 */ ULP = 0.78223048105528226Interval tested for asinpi: [0.75,1] xm = 7.5425933001419776e-01, /* 0x3fe822e4, 0x7663a4aa */libm = 2.7200385380185182e-01, /* 0x3fd16882, 0xda1dc13b */mpfr = 2.7200385380185188e-01, /* 0x3fd16882, 0xda1dc13c */ ULP = 0.53747973176773822Interval tested for atanpi: [9.31323e-10,0.25] xm = 1.9666113418757322e-01, /* 0x3fc92c31, 0x29dd6d2f */libm = 6.1810387818117797e-02, /* 0x3fafa59c, 0x7476baa5 */mpfr = 6.1810387818117804e-02, /* 0x3fafa59c, 0x7476baa6 */ ULP = 0.54674297446584263Interval tested for atanpi: [0.25,0.5] xm = 4.1312119637707068e-01, /* 0x3fda7093, 0xe2ee5494 */libm = 1.2470309560460152e-01, /* 0x3fbfec8a, 0xc554ebec */mpfr = 1.2470309560460154e-01, /* 0x3fbfec8a, 0xc554ebed */ ULP = 0.73116638175113347Interval tested for atanpi: [0.5,0.75] xm = 5.0018949583396499e-01, /* 0x3fe0018d, 0x66cd1b82 */libm = 1.4763186871058706e-01, /* 0x3fc2e599, 0xdffacb8f */mpfr = 1.4763186871058709e-01, /* 0x3fc2e599, 0xdffacb90 */ ULP = 0.69192753950764663Interval tested for atanpi: [0.75,1] xm = 7.5007880583359599e-01, /* 0x3fe800a5, 0x448f4c03 */libm = 2.0484881828445453e-01, /* 0x3fca387c, 0x6f93f71f */mpfr = 2.0484881828445450e-01, /* 0x3fca387c, 0x6f93f71e */ ULP = 0.65765471872064396Interval tested for atanpi: [1,2] xm = 1.0103228000344093e+00, /* 0x3ff02a48, 0x3d88d0a2 */libm = 2.5163447403817019e-01, /* 0x3fd01ac7, 0x7b229108 */mpfr = 2.5163447403817013e-01, /* 0x3fd01ac7, 0x7b229107 */ ULP = 0.67409519689166042Interval tested for atanpi: [2,4] xm = 2.0231383267437946e+00, /* 0x40002f63, 0x25a530a9 */libm = 3.5387589538123299e-01, /* 0x3fd6a5e7, 0x156053c6 */mpfr = 3.5387589538123293e-01, /* 0x3fd6a5e7, 0x156053c5 */ ULP = 0.69695587476021503Interval tested for atanpi: [4,1.79769e+308] xm = 4.0000000000000000e+00, /* 0x40100000, 0x00000000 */libm = 4.2202086962263069e-01, /* 0x3fdb0263, 0xd2508e31 */mpfr = 4.2202086962263069e-01, /* 0x3fdb0263, 0xd2508e31 */ ULP = 0.27709400511686716PR: 295884MFC after: 1 monthReviewed by: fuz List of files: /freebsd/lib/msun/Makefile

c9991e01149b5cd5a6710d70921b203ef7f16fc2 - powerpc64le: switch from calling wrappers to weak references

Thu, 25 Jun 2026 19:33:42 +0200

powerpc64le: switch from calling wrappers to weak referencesFixes: 255538cd906045095d0c2113ae6c4731ce36c0cfDifferential Revision: https://reviews.freebsd.org/D57850Reviewed by: adrian List of files: /freebsd/lib/msun/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