musl - issues in double-extended functions

Message-Id: <p9u01pckstbt.fsf@thym>
Date: Thu, 30 Jul 2026 11:46:46 +0200
From: Paul Zimmermann <Paul.Zimmermann@...ia.fr>
To: musl@...ts.openwall.com
Subject: issues in double-extended functions

       Hi,

for x=-16384, expl returns 2^-16384 instead of 0, maybe the code was copied
from exp2l:

NEW exp 0 -1 -0x4p+12l [0] [2.31e+18] 2.30585e+18 2.305843009213694e+18
libm gives 0x1p-16384l
mpfr gives 0x0p+0l

and for x,y=-0x8p-16448,-0x1p+0, powl returns +Inf instead of -Inf:

NEW pow 0 -1 -0x8p-16448l,-0x1p+0l [inf] inf inf
libm gives infl
mpfr gives -infl

This was tested with 1.2.6, but maybe present in earlier versions.

Paul Zimmermann

Confused about mailing lists and their use? Read about mailing lists on Wikipedia and check out these guidelines on proper formatting of your messages.