DMARC-Filter:	OpenDMARC Filter v1.4.2 delorie.com 60J8sHmN4023012
Authentication-Results:	delorie.com; dmarc=pass (p=none dis=none) header.from=cygwin.com
Authentication-Results:	delorie.com; spf=pass smtp.mailfrom=cygwin.com
DKIM-Filter:	OpenDKIM Filter v2.11.0 delorie.com 60J8sHmN4023012
Authentication-Results:	delorie.com;
	dkim=pass (1024-bit key, unprotected) header.d=cygwin.com header.i=@cygwin.com header.a=rsa-sha256 header.s=default header.b=Q2LpaBbu
X-Recipient:	archive-cygwin AT delorie DOT com
DKIM-Filter:	OpenDKIM Filter v2.11.0 sourceware.org 2E2684BA9018
DKIM-Signature:	v=1; a=rsa-sha256; c=relaxed/relaxed; d=cygwin.com;
	s=default; t=1768812856;
	bh=zQZSrRfJxNCeoy0Pn3TVVVDaTaW8YEsVF1P5/zjzggs=;
	h=Date:To:In-Reply-To:Subject:References:List-Id:List-Unsubscribe:
	List-Archive:List-Post:List-Help:List-Subscribe:From:Reply-To:
	From;
	b=Q2LpaBbuyCmPAmRDu/1NrEURlJvvfX7hD0NkcnN5/CrdnBqT/4DXqOjhgNwmcA13P
	p8KiF6BlEM237dmr3zQzYGJlNUY6e0OYW9tayVq7DJNrpegK4rQlwftCpY6sVJviNj
	s5RggnA7OAEMaHlXVfKrLjrw8zftE9iKKG5gI3ZY=
X-Original-To:	cygwin AT cygwin DOT com
Delivered-To:	cygwin AT cygwin DOT com
DMARC-Filter:	OpenDMARC Filter v1.4.2 sourceware.org 0F2804BA2E22
ARC-Filter:	OpenARC Filter v1.0.0 sourceware.org 0F2804BA2E22
ARC-Seal:	i=1; a=rsa-sha256; d=sourceware.org; s=key; t=1768812806; cv=none;
	b=uBMuhY6PvWMJaF0OnI1UEmPLlmKmsql3uVV9KYKi5BkdFdjAjPhWEwvscRinh+oVxVpXHivYILoT7NTo1ZZTRdKQMDSKyKOq8Pl7OqGm3hVossXMAG6cOMGo2CJIXZth2P+cVyTN/WF6WQ3hngKgufTE3LEqYIhRuz5WCCgLxJ8=
ARC-Message-Signature:	i=1; a=rsa-sha256; d=sourceware.org; s=key;
	t=1768812806; c=relaxed/simple;
	bh=kRlW1IMkePVBdXrzvQImfbBsMqBVmXkW4MfsWuuuids=;
	h=DKIM-Signature:Date:Message-Id:From:To:Subject;
	b=rC+/iXHLwpMrUSsjnzyJsgKhR1//8Pga+Iz8QDvIzYssfkuI8fC3iabPVKSSWPSjdIOSKeiu2/T0JrWH4TI7ajmjJwmXL1yYPLh4gQyuy564o8LhZnRtnZrEHHLVNT9huT0f5Z5uIg4jVJkO2r+2lKfh3MeJyiG/M0jEU0TUmt4=
ARC-Authentication-Results:	i=1; server2.sourceware.org
DKIM-Filter:	OpenDKIM Filter v2.11.0 sourceware.org 0F2804BA2E22
X-CSE-ConnectionGUID:	NspuSPf4SZCMy/hqDRLPKQ==
X-CSE-MsgGUID:	6ZnK2ZhdSXqu/Bd6GwHbrw==
Received-SPF:	SoftFail (mail3-relais-sop.national.inria.fr:
	domain of Paul DOT Zimmermann AT inria DOT fr is inclined to not
	designate 152.81.2.78 as permitted sender) identity=mailfrom;
	client-ip=152.81.2.78;
	receiver=mail3-relais-sop.national.inria.fr;
	envelope-from="Paul DOT Zimmermann AT inria DOT fr";
	x-sender="Paul DOT Zimmermann AT inria DOT fr"; x-conformance=spf_only;
	x-record-type="v=spf1"; x-record-text="v=spf1
	include:mailout.safebrands.com a:basic-mail.safebrands.com
	a:basic-mail01.safebrands.com a:basic-mail02.safebrands.com
	ip4:128.93.142.0/24 ip4:192.134.164.0/24 ip4:128.93.162.160
	ip4:128.93.162.3 ip4:128.93.162.88 ip4:89.107.174.7 mx ~all"
Received-SPF:	None (mail3-relais-sop.national.inria.fr: no sender
	authenticity information available from domain of
	postmaster AT araignee) identity=helo; client-ip=152.81.2.78;
	receiver=mail3-relais-sop.national.inria.fr;
	envelope-from="Paul DOT Zimmermann AT inria DOT fr";
	x-sender="postmaster AT araignee"; x-conformance=spf_only
X-IronPort-AV:	E=Sophos;i="6.21,237,1763420400"; d="scan'208";a="136359244"
Date:	Mon, 19 Jan 2026 09:53:24 +0100
Message-Id:	<p9u0ikcyugor.fsf@araignee.loria.fr>
To:	cygwin AT cygwin DOT com
In-Reply-To:	<3c56778a-870c-4bae-a855-6b96392870c8@SystematicSW.ab.ca>
	(message from Brian Inglis on Sun, 18 Jan 2026 15:43:21 -0700)
Subject:	Re: Log2 seems inefficient
References:	<226987210 DOT 591128 DOT 1768766440723 AT connect DOT xfinity DOT com>
	<3c56778a-870c-4bae-a855-6b96392870c8 AT SystematicSW DOT ab DOT ca>
X-BeenThere:	cygwin AT cygwin DOT com
X-Mailman-Version:	2.1.30
List-Id:	General Cygwin discussions and problem reports <cygwin.cygwin.com>
List-Unsubscribe:	<https://cygwin.com/mailman/options/cygwin>,
	<mailto:cygwin-request AT cygwin DOT com?subject=unsubscribe>
List-Archive:	<https://cygwin.com/pipermail/cygwin/>
List-Post:	<mailto:cygwin AT cygwin DOT com>
List-Help:	<mailto:cygwin-request AT cygwin DOT com?subject=help>
List-Subscribe:	<https://cygwin.com/mailman/listinfo/cygwin>,
	<mailto:cygwin-request AT cygwin DOT com?subject=subscribe>
From:	Paul Zimmermann via Cygwin <cygwin AT cygwin DOT com>
Reply-To:	Paul Zimmermann <Paul DOT Zimmermann AT inria DOT fr>
MIME-Version:	1.0
Errors-To:	cygwin-bounces~archive-cygwin=delorie DOT com AT cygwin DOT com
Sender:	"Cygwin" <cygwin-bounces~archive-cygwin=delorie DOT com AT cygwin DOT com>

       Hi,

> With default settings the worst case error in nearest rounding mode is 0.547 ULP 
> with inlined fma and fma contraction.
> It uses a 1 KB lookup table,
> ...
> Note that the math.h header defines log2(x) to be log(x)/Ln2, this is
> not changed, so the new code is only used if that macro is suppressed."
>
> Paul Zimmerman regularly posts library math function comparisons on the newlib 
> and comparable library mailing lists, available at:
> 
> 	https://members.loria.fr/PZimmermann/papers/accuracy.pdf
> 
> that should give you an idea of the current state of the functions, if he is 
> even aware of the issue.

indeed, the known worst case error for Newlib log2 is 2.06 ulp:

log2 0 -1 0x1.68d778f076021p+0 [2] [2.06] 2.05526 2.055255219183752
libm gives 0x1.fb1b88680e7adp-2
mpfr gives 0x1.fb1b88680e7abp-2

The 0.547 value above could match the GNU libc worst error:

log2 0 -1 0x1.1406d79e1b574p+0 [1] [0.548] 0.547538 0.5475371438132311
libm gives 0x1.bd16a11d3404ep-4
mpfr gives 0x1.bd16a11d3404fp-4

If your goal is accuracy, you can try the CORE-MATH implementation:

https://gitlab.inria.fr/core-math/core-math/-/blob/master/src/binary64/log2/log2.c

It takes 19 cycles on average on an Intel(R) Core(TM) i7-8700, and is
correctly rounded (for any rounding mode).

Paul



-- 
Problem reports:      https://cygwin.com/problems.html
FAQ:                  https://cygwin.com/faq/
Documentation:        https://cygwin.com/docs.html
Unsubscribe info:     https://cygwin.com/ml/#unsubscribe-simple

- Raw text -