From mboxrd@z Thu Jan 1 00:00:00 1970 X-Spam-Checker-Version: SpamAssassin 3.4.2 (2018-09-13) on inbox.vuxu.org X-Spam-Level: X-Spam-Status: No, score=-2.0 required=5.0 tests=MAILING_LIST_MULTI, RCVD_IN_DNSWL_MED,RCVD_IN_MSPIKE_H3,RCVD_IN_MSPIKE_WL,RDNS_NONE, SPF_PASS autolearn=ham autolearn_force=no version=3.4.2 Received: (qmail 18514 invoked from network); 22 Mar 2020 19:10:36 -0000 Received-SPF: pass (mother.openwall.net: domain of lists.openwall.com designates 195.42.179.200 as permitted sender) receiver=inbox.vuxu.org; client-ip=195.42.179.200 envelope-from= Received: from unknown (HELO mother.openwall.net) (195.42.179.200) by inbox.vuxu.org with ESMTP; 22 Mar 2020 19:10:36 -0000 Received: (qmail 5442 invoked by uid 550); 22 Mar 2020 19:10:32 -0000 Mailing-List: contact musl-help@lists.openwall.com; run by ezmlm Precedence: bulk List-Post: List-Help: List-Unsubscribe: List-Subscribe: List-ID: Reply-To: musl@lists.openwall.com Received: (qmail 5424 invoked from network); 22 Mar 2020 19:10:31 -0000 Date: Sun, 22 Mar 2020 15:10:19 -0400 From: Rich Felker To: musl@lists.openwall.com Message-ID: <20200322191019.GP11469@brightrain.aerifal.cx> References: <20200222195925.GK1663@brightrain.aerifal.cx> <20200223001942.GM1663@brightrain.aerifal.cx> <20200320181250.GC11469@brightrain.aerifal.cx> <20200322011958.GM11469@brightrain.aerifal.cx> <20200322175325.GO11469@brightrain.aerifal.cx> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In-Reply-To: User-Agent: Mutt/1.5.21 (2010-09-15) Subject: Re: [musl] Q: dealing with missing removal of excess precision On Sun, Mar 22, 2020 at 09:51:09PM +0300, Alexander Monakov wrote: > On Sun, 22 Mar 2020, Rich Felker wrote: > > > > Nack for sqrt 'unsigned short' fix, I recommend to consider > > > 'unsigned fpsr = 0' and "+a" constraint, the resulting assembly is > > > much better (GCC doesn't seem to do a good job optimizing the 'unsigned short' > > > variant at all). > > > > Sorry, I forgot to disassemble and check after making that change, and > > indeed it is very bad. > > > > How about just leaving it as-is? The value is masked such that upper > > bits don't matter, and "whatever the register happened to contain > > before" is a valid albeit ugly output from inline asm -- it doesn't > > admit any compiler transformations that would cause inconsistent value > > to be observed or other problems. > > Yes, I guess it's acceptable. > > > > Actually, may I ask why the initial commit did not mention that it relies > > > on this nontrivial property? > > > > The need for fixup of double was only realized later, in commit > > 809556e60a3359f646946879dd94c4760e5b8e84. It was discussed at the time > > that no action was needed for single, but it seems since there was no > > change there wasn't any mention of it in log. > > Are you sure? single-precision sqrtf received a change just two days > prior to that in commit e0a54e6725 ("correct rounding for i387 sqrtf function") > which is the same day as when Stephen Canon supplied his answer in > https://stackoverflow.com/questions/9678224/is-there-any-way-to-get-correct-rounding-with-the-i387-fsqrt-instruction > (and also the same day you asked the question) I just dug through the old gcc and glibc bz's it was mentioned on (52593 and 14032 resp.) and didn't find anything, but I'm pretty sure it was known in 2012 that it was a non-issue for single-precision. What I didn't understand then was that the callee was responsible for dropping the excess precision; I wrongly believed that "something" in the caller would make it not-matter/get collapsed down to nominal precision and rounded correctly. Commit e0a54e6725 and related commits were fixing that mistake, which I'd since recognized was wrong but hadn't yet recognized the urgency of fixing until I started seeing how bad it could break the compiler. Rich