From mboxrd@z Thu Jan  1 00:00:00 1970
Return-Path: <gcc-bugzilla@gcc.gnu.org>
Received: by sourceware.org (Postfix, from userid 48)
	id 43FF63858C5E; Fri,  3 Mar 2023 12:54:55 +0000 (GMT)
DKIM-Filter: OpenDKIM Filter v2.11.0 sourceware.org 43FF63858C5E
DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gcc.gnu.org;
	s=default; t=1677848095;
	bh=UJ729SV3LmVxmz78ROHRbFuObsEKR3gXYraAdVeBKlg=;
	h=From:To:Subject:Date:In-Reply-To:References:From;
	b=ERPRYvrw5DYWEQyCQVT8uHtVW948eGHynoX1Sx+UpK9sxkrVadINEL4yH3H3un5xO
	 Q890753iFMgyjUPrXjV3sIJzVPrAXtTkZTwOJAT2T5Tdk0mx/BlvWBQex20ToUI2NF
	 yo8HGxZTo22c62gtyLDzEChhZEhnL+eCDn3aKq+c=
From: "jakub at gcc dot gnu.org" <gcc-bugzilla@gcc.gnu.org>
To: gcc-bugs@gcc.gnu.org
Subject: [Bug tree-optimization/109008] [13 Regression] Wrong code in scipy
 package since r13-3926-gd4c2f1d376da6f
Date: Fri, 03 Mar 2023 12:54:55 +0000
X-Bugzilla-Reason: CC
X-Bugzilla-Type: changed
X-Bugzilla-Watch-Reason: None
X-Bugzilla-Product: gcc
X-Bugzilla-Component: tree-optimization
X-Bugzilla-Version: 13.0
X-Bugzilla-Keywords: wrong-code
X-Bugzilla-Severity: normal
X-Bugzilla-Who: jakub at gcc dot gnu.org
X-Bugzilla-Status: NEW
X-Bugzilla-Resolution: 
X-Bugzilla-Priority: P1
X-Bugzilla-Assigned-To: unassigned at gcc dot gnu.org
X-Bugzilla-Target-Milestone: 13.0
X-Bugzilla-Flags: 
X-Bugzilla-Changed-Fields: 
Message-ID: <bug-109008-4-YeLJlufFK6@http.gcc.gnu.org/bugzilla/>
In-Reply-To: <bug-109008-4@http.gcc.gnu.org/bugzilla/>
References: <bug-109008-4@http.gcc.gnu.org/bugzilla/>
Content-Type: text/plain; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable
X-Bugzilla-URL: http://gcc.gnu.org/bugzilla/
Auto-Submitted: auto-generated
MIME-Version: 1.0
List-Id: <gcc-bugs.sourceware.org>

https://gcc.gnu.org/bugzilla/show_bug.cgi?id=3D109008

--- Comment #6 from Jakub Jelinek <jakub at gcc dot gnu.org> ---
Note, the ulps frange_arithmetic are ulps of the result, result is in this =
case
0.0,
so 1ulp is the smallest subnormal number.
That is something completely different from what we need here though for the
reverse
operations.
[1.0, 1.0] + x =3D=3D [1.0, 1.0] is true with round to nearest for x in
[-0x1.0p-54, 0x1.0p-53], so that is ulp of 1.0, not the result.
Now, is this problem just around the 0 result or other values too?
Say [1.0, 1.0] + x =3D=3D [2.0, 2.0]?=