Skip site navigation (1)Skip section navigation (2)
Date:      Sat, 02 May 2026 14:56:16 +0000
From:      Robert Clausecker <fuz@FreeBSD.org>
To:        src-committers@FreeBSD.org, dev-commits-src-all@FreeBSD.org, dev-commits-src-main@FreeBSD.org
Cc:        =?utf-8?Q?Jes=C3=BAs?= =?utf-8?Q?Bl=C3=A1zqu?=ez <jesuscblazquez@gmail.com>
Subject:   git: 7c20e15592a0 - main - lib/msun: Replaced pattern to force exception in _num families
Message-ID:  <69f61090.47462.d613e0d@gitrepo.freebsd.org>

index | next in thread | raw e-mail

The branch main has been updated by fuz:

URL: https://cgit.FreeBSD.org/src/commit/?id=7c20e15592a07ea457cacb0d6706948815c8420e

commit 7c20e15592a07ea457cacb0d6706948815c8420e
Author:     Jesús Blázquez <jesuscblazquez@gmail.com>
AuthorDate: 2026-04-27 16:43:41 +0000
Commit:     Robert Clausecker <fuz@FreeBSD.org>
CommitDate: 2026-05-02 14:50:49 +0000

    lib/msun: Replaced pattern to force exception in _num families
    
    Replaced the old pattern of using a ternary to force addition
    (raising exceptions for sNaN's) with a new one using a volatile
    variable. The _mag_num family was already implemented with this pattern
    
    PR:             294719
    Reviewed by:    fuz, kargl
    MFC after:      1 month
---
 lib/msun/src/s_fmaximum_num.c  | 14 ++++++++++----
 lib/msun/src/s_fmaximum_numf.c | 14 ++++++++++----
 lib/msun/src/s_fmaximum_numl.c | 14 ++++++++++----
 lib/msun/src/s_fminimum_num.c  | 14 ++++++++++----
 lib/msun/src/s_fminimum_numf.c | 14 ++++++++++----
 lib/msun/src/s_fminimum_numl.c | 14 ++++++++++----
 6 files changed, 60 insertions(+), 24 deletions(-)

diff --git a/lib/msun/src/s_fmaximum_num.c b/lib/msun/src/s_fmaximum_num.c
index c7cfcd26fd70..4d7c860d2ce9 100644
--- a/lib/msun/src/s_fmaximum_num.c
+++ b/lib/msun/src/s_fmaximum_num.c
@@ -53,12 +53,18 @@ fmaximum_num(double x, double y)
 	nan_y = isnan(y);
 
 	if (nan_x || nan_y) {
-		/* These ternary conditionals force (x+y), so that sNaN's raise exceptions */
+		/* If both are NaN, adding returns qNaN */
 		if (nan_x && nan_y)
-			return (x + y);
+		    return (x + y);
+
+		/* force_except makes sure sNaN's raise exceptions */
+		volatile double force_except = x + y;
+		force_except;
+
 		if (nan_x)
-			return ((x + y) != 0.0 ? y : y);
-		return ((x + y) != 0.0 ? x : x);
+			return (y);
+		else
+			return (x);
 	}
 
 	/* Handle comparisons of signed zeroes. */
diff --git a/lib/msun/src/s_fmaximum_numf.c b/lib/msun/src/s_fmaximum_numf.c
index 60a7ccae343e..b41f293a428d 100644
--- a/lib/msun/src/s_fmaximum_numf.c
+++ b/lib/msun/src/s_fmaximum_numf.c
@@ -52,12 +52,18 @@ fmaximum_numf(float x, float y)
 	nan_y = isnan(y);
 
 	if (nan_x || nan_y) {
-		/* These ternary conditionals force (x+y), so that sNaN's raise exceptions */
+		/* If both are NaN, adding returns qNaN */
 		if (nan_x && nan_y)
-			return (x + y);
+		    return (x + y);
+
+		/* force_except makes sure sNaN's raise exceptions */
+		volatile float force_except = x + y;
+		force_except;
+
 		if (nan_x)
-			return ((x + y) != 0.0 ? y : y);
-		return ((x + y) != 0.0 ? x : x);
+			return (y);
+		else
+			return (x);
 	}
 
 	/* Handle comparisons of signed zeroes. */
diff --git a/lib/msun/src/s_fmaximum_numl.c b/lib/msun/src/s_fmaximum_numl.c
index a2a5383ab85a..74a3d9ba2211 100644
--- a/lib/msun/src/s_fmaximum_numl.c
+++ b/lib/msun/src/s_fmaximum_numl.c
@@ -47,12 +47,18 @@ fmaximum_numl(long double x, long double y)
 	nan_y = isnan(y);
 
 	if (nan_x || nan_y) {
-		/* These ternary conditionals force (x+y), so that sNaN's raise exceptions */
+		/* If both are NaN, adding returns qNaN */
 		if (nan_x && nan_y)
-			return (x + y);
+		    return (x + y);
+
+		/* force_except makes sure sNaN's raise exceptions */
+		volatile long double force_except = x + y;
+		force_except;
+
 		if (nan_x)
-			return ((x + y) != 0.0 ? y : y);
-		return ((x + y) != 0.0 ? x : x);
+			return (y);
+		else
+			return (x);
 	}
 
 	/* Handle comparisons of signed zeroes. */
diff --git a/lib/msun/src/s_fminimum_num.c b/lib/msun/src/s_fminimum_num.c
index a5fd657955c5..67999c366a76 100644
--- a/lib/msun/src/s_fminimum_num.c
+++ b/lib/msun/src/s_fminimum_num.c
@@ -53,12 +53,18 @@ fminimum_num(double x, double y)
 	nan_y = isnan(y);
 
 	if (nan_x || nan_y) {
-		/* These ternary conditionals force (x+y), so that sNaN's raise exceptions */
+		/* If both are NaN, adding returns qNaN */
 		if (nan_x && nan_y)
-			return (x + y);
+		    return (x + y);
+
+		/* force_except makes sure sNaN's raise exceptions */
+		volatile double force_except = x + y;
+		force_except;
+
 		if (nan_x)
-			return ((x + y) != 0.0 ? y : y);
-		return ((x + y) != 0.0 ? x : x);
+			return (y);
+		else
+			return (x);
 	}
 
 	/* Handle comparisons of signed zeroes. */
diff --git a/lib/msun/src/s_fminimum_numf.c b/lib/msun/src/s_fminimum_numf.c
index d5e8f0812bdb..c4b686be158e 100644
--- a/lib/msun/src/s_fminimum_numf.c
+++ b/lib/msun/src/s_fminimum_numf.c
@@ -52,12 +52,18 @@ fminimum_numf(float x, float y)
 	nan_y = isnan(y);
 
 	if (nan_x || nan_y) {
-		/* These ternary conditionals force (x+y), so that sNaN's raise exceptions */
+		/* If both are NaN, adding returns qNaN */
 		if (nan_x && nan_y)
-			return (x + y);
+		    return (x + y);
+
+		/* force_except makes sure sNaN's raise exceptions */
+		volatile float force_except = x + y;
+		force_except;
+
 		if (nan_x)
-			return ((x + y) != 0.0 ? y : y);
-		return ((x + y) != 0.0 ? x : x);
+			return (y);
+		else
+			return (x);
 	}
 
 	/* Handle comparisons of signed zeroes. */
diff --git a/lib/msun/src/s_fminimum_numl.c b/lib/msun/src/s_fminimum_numl.c
index dc0884002499..b30fbc721b6b 100644
--- a/lib/msun/src/s_fminimum_numl.c
+++ b/lib/msun/src/s_fminimum_numl.c
@@ -47,12 +47,18 @@ fminimum_numl(long double x, long double y)
 	nan_y = isnan(y);
 
 	if (nan_x || nan_y) {
-		/* These ternary conditionals force (x+y), so that sNaN's raise exceptions */
+		/* If both are NaN, adding returns qNaN */
 		if (nan_x && nan_y)
-			return (x + y);
+		    return (x + y);
+
+		/* force_except makes sure sNaN's raise exceptions */
+		volatile long double force_except = x + y;
+		force_except;
+
 		if (nan_x)
-			return ((x + y) != 0.0 ? y : y);
-		return ((x + y) != 0.0 ? x : x);
+			return (y);
+		else
+			return (x);
 	}
 
 	/* Handle comparisons of signed zeroes. */


home | help

Want to link to this message? Use this
URL: <https://mail-archive.FreeBSD.org/cgi/mid.cgi?69f61090.47462.d613e0d>