Jakub Jelinek 481ba4fb5f libquadmath: Use soft-fp for sqrtq finite positive arguments [PR114623] ...

sqrt should be 0.5ulp precise, but the current implementation is less
precise than that.
The following patch uses the soft-fp code (like e.g. glibc for x86) for it
if possible.  I didn't want to replicate the libgcc infrastructure for
choosing the right sfp-machine.h, so the patch just uses a single generic
implementation.  As the code is used solely for the finite positive arguments,
it shouldn't generate NaNs (so the exact form of canonical QNaN/SNaN is
irrelevant), and sqrt for these shouldn't produce underflows/overflows either,
for < 1.0 arguments it always returns larger values than the argument and for
> 1.0 smaller values than the argument.

2024-04-09  Jakub Jelinek  <jakub@redhat.com>

	PR libquadmath/114623
	* sfp-machine.h: New file.
	* math/sqrtq.c: Include from libgcc/soft-fp also soft-fp.h and quad.h
	if possible.
	(USE_SOFT_FP): Define in that case.
	(sqrtq): Use soft-fp based implementation for the finite positive
	arguments if possible.

2024-04-09 08:17:25 +02:00

..

acoshq.c

…

acosq.c

…

asinhq.c

…

asinq.c

…

atan2q.c

…

atanhq.c

…

atanq.c

…

cacoshq.c

…

cacosq.c

…

casinhq_kernel.c

…

casinhq.c

…

casinq.c

…

catanhq.c

…

catanq.c

…

cbrtq.c

…

ccoshq.c

…

ceilq.c

…

cexpq.c

…

cimagq.c

…

clog10q.c

…

clogq.c

…

complex.c

…

conjq.c

…

copysignq.c

…

coshq.c

…

cosq_kernel.c

…

cosq.c

…

cprojq.c

…

crealq.c

…

csinhq.c

…

csinq.c

…

csqrtq.c

…

ctanhq.c

…

ctanq.c

…

erfq.c

…

exp2q.c

…

expm1q.c

…

expq_table.h

…

expq.c

…

fabsq.c

…

fdimq.c

…

finiteq.c

…

floorq.c

…

fmaq.c

…

fmaxq.c

…

fminq.c

…

fmodq.c

…

frexpq.c

…

hypotq.c

…

ilogbq.c

…

isinfq.c

…

isnanq.c

…

issignalingq.c

…

j0q.c

…

j1q.c

…

jnq.c

…

ldexpq.c

…

lgammaq_neg.c

…

lgammaq_product.c

…

lgammaq.c

…

llrintq.c

…

llroundq.c

…

log1pq.c

…

log2q.c

…

log10q.c

…

logbq.c

…

logq.c

…

lrintq.c

…

lroundq.c

…

modfq.c

…

nanq.c

…

nearbyintq.c

…

nextafterq.c

…

powq.c

…

rem_pio2q.c

…

remainderq.c

…

remquoq.c

…

rintq.c

…

roundq.c

…

scalblnq.c

…

scalbnq.c

…

signbitq.c

…

sincos_table.c

…

sincosq_kernel.c

…

sincosq.c

…

sinhq.c

…

sinq_kernel.c

…

sinq.c

…

sqrtq.c

libquadmath: Use soft-fp for sqrtq finite positive arguments [PR114623]

2024-04-09 08:17:25 +02:00

tanhq.c

…

tanq_kernel.c

…

tanq.c

…

tgammaq_product.c

…

tgammaq.c

…

truncq.c

…

x2y2m1q.c

…