From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: from elaine.keithp.com (home.keithp.com [63.227.221.253]) by sourceware.org (Postfix) with ESMTPS id AE092386103B for ; Wed, 16 Sep 2020 15:14:02 +0000 (GMT) DMARC-Filter: OpenDMARC Filter v1.3.2 sourceware.org AE092386103B Received: from localhost (localhost [127.0.0.1]) by elaine.keithp.com (Postfix) with ESMTP id E9EFB3F2D83A; Wed, 16 Sep 2020 08:14:00 -0700 (PDT) X-Virus-Scanned: Debian amavisd-new at keithp.com Received: from elaine.keithp.com ([127.0.0.1]) by localhost (elaine.keithp.com [127.0.0.1]) (amavisd-new, port 10024) with LMTP id jwlKD99bLi0x; Wed, 16 Sep 2020 08:13:59 -0700 (PDT) Received: from keithp.com (koto.keithp.com [10.0.0.2]) by elaine.keithp.com (Postfix) with ESMTPSA id 9634D3F2D2F1; Wed, 16 Sep 2020 08:13:59 -0700 (PDT) Received: by keithp.com (Postfix, from userid 1000) id 35EAE15821A6; Wed, 16 Sep 2020 08:13:59 -0700 (PDT) From: "Keith Packard" To: Paul Zimmermann Cc: newlib@sourceware.org Subject: Re: issue with tgammaf In-Reply-To: References: <87een3eyq4.fsf@keithp.com> Date: Wed, 16 Sep 2020 08:13:58 -0700 Message-ID: <87sgbheouh.fsf@keithp.com> MIME-Version: 1.0 Content-Type: multipart/signed; boundary="=-=-="; micalg=pgp-sha256; protocol="application/pgp-signature" X-Spam-Status: No, score=-4.8 required=5.0 tests=BAYES_00, DKIM_SIGNED, DKIM_VALID, DKIM_VALID_AU, DKIM_VALID_EF, SPF_HELO_NONE, SPF_PASS, TXREP autolearn=ham autolearn_force=no version=3.4.2 X-Spam-Checker-Version: SpamAssassin 3.4.2 (2018-09-13) on server2.sourceware.org X-BeenThere: newlib@sourceware.org X-Mailman-Version: 2.1.29 Precedence: list List-Id: Newlib mailing list List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 16 Sep 2020 15:14:05 -0000 --=-=-= Content-Type: text/plain Content-Transfer-Encoding: quoted-printable Paul Zimmermann writes: > thank you Keith. For your information I have updated my comparison with > Newlib 3.3.0: > > https://members.loria.fr/PZimmermann/papers/accuracy.pdf > > The differences with the previous (unknown) version of Newlib I used are > that j1 is less accurate (largest error goes from 2.66e6 ulps to 1.68e7 u= lps), > and y1 is more accurate (largest error goes from 4.65e8 ulps to 6.18e6 ul= ps). > (However, those huge errors simply tell us that j1f is completely off, in= the > first case for x=3D0x6.cfd78p+100 j1f returns 0 instead of -0x1p-52.) > > I am still unable to compile a program using lgammaf: You might have more luck with lgammaf_r as that avoids use of _impure_ptr in newlib. Alternatively, you can use the picolibc fork of newlib which uses native thread local storage instead of a giant struct full of library global data. That has built-in support for compiling on x86_64 for testing. =2D-=20 =2Dkeith --=-=-= Content-Type: application/pgp-signature; name="signature.asc" -----BEGIN PGP SIGNATURE----- iQIzBAEBCAAdFiEEw4O3eCVWE9/bQJ2R2yIaaQAAABEFAl9iK7YACgkQ2yIaaQAA ABFYlA/+LdGNFOq8ro4ccCz1l4CAXXZQvallqJLVzVnrfsuNylQPaWDOxRVtwmGw dC6NVhCMJC7MFJbE/ThgweZ9Dcu9QvrB9Li3LfcGBoMoutRq8RUOErJ65NuyZO00 TWs4JOSTJUsJAD8XmvSvbTaPvjGsj1xaVoe/DsCD0x75MqSiwi5lZ2OfXLzrKf2u udIblpoF2rMWAFPSAIPFqUODSRGMRhoSeHnyppyNGLp71wuAVAXwuynFXmDmWp+d YhCj8fkW0ckUTJX8n7RJirnlJxyucIaHcT3ap/Q0JVRM/C4PWYMRU/XU4j6YhrJu +LMT6JxXUmuPM7kyJgGZOtadvVoRn8XMyEtdyJpu963S1LdlT429ZI7RqFhfvODZ Bnsvr3vvHf3kQ13Dojbu8xS/9+6XEGM4F8W1MrmwSyvHjecbefJ5ulGxsHYniBQW ug1IGvuNhKjgiWYJsc7aWG7OwyQyqrg9WGZTjJ1sbfocdRmV4E9MH/Y9lyMVbds+ E3a3bP54zhP9ZxZIzi9YlBiFq8kfr0cmmPSeIH9PsYLGvbfHUg1P89z+5MOMqag0 VGYRZCjRLgDGNXeYHWylC3wTFwyXeZ0pg1hUKXpN7e6m+beyZyTfYUVp6DnYObzc pPslqTGfdHv7glbAhe0NIpZJEAzQNRdZAEA+yEw3NWymDsi72Rs= =VfWQ -----END PGP SIGNATURE----- --=-=-=--