git: 25cca51ed294 - main - libm: remainder: make sure x is zero

From: Ahmad Khalifa <vexeduxr_at_FreeBSD.org>
Date: Sat, 11 Oct 2025 10:17:31 UTC
The branch main has been updated by vexeduxr:

URL: https://cgit.FreeBSD.org/src/commit/?id=25cca51ed294890d20a3c0290814cd26875db686

commit 25cca51ed294890d20a3c0290814cd26875db686
Author:     Ahmad Khalifa <vexeduxr@FreeBSD.org>
AuthorDate: 2025-10-10 09:30:52 +0000
Commit:     Ahmad Khalifa <vexeduxr@FreeBSD.org>
CommitDate: 2025-10-11 10:13:19 +0000

    libm: remainder: make sure x is zero
    
    Make sure the entirety of x is zero before flipping the sign bit.
    Otherwise the sign would be wrong for small values of x when x is
    negative and |n*y| > |x|
    
    Reported by:    alfredo
    PR:             251091
    Reviewed by:    kargl
    MFC after:      3 days
    Differential Revision: https://reviews.freebsd.org/D53023
---
 lib/msun/src/e_remainder.c | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/lib/msun/src/e_remainder.c b/lib/msun/src/e_remainder.c
index a5fb7141d01a..cc6cd320073e 100644
--- a/lib/msun/src/e_remainder.c
+++ b/lib/msun/src/e_remainder.c
@@ -64,8 +64,8 @@ remainder(double x, double p)
 		if(x>=p_half) x -= p;
 	    }
 	}
-	GET_HIGH_WORD(hx,x);
-	if ((hx&0x7fffffff)==0) hx = 0;
+	EXTRACT_WORDS(hx, lx, x);
+	if (((hx&0x7fffffff)|lx) == 0) hx = 0;
 	SET_HIGH_WORD(x,hx^sx);
 	return x;
 }