From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: from mail2-relais-roc.national.inria.fr (mail2-relais-roc.national.inria.fr [192.134.164.83]) by sourceware.org (Postfix) with ESMTPS id EA0CA385802F for ; Thu, 6 Jan 2022 15:33:12 +0000 (GMT) DMARC-Filter: OpenDMARC Filter v1.4.1 sourceware.org EA0CA385802F Authentication-Results: sourceware.org; dmarc=none (p=none dis=none) header.from=inria.fr Authentication-Results: sourceware.org; spf=pass smtp.mailfrom=inria.fr DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=inria.fr; s=dc; h=date:message-id:from:to:cc:in-reply-to:subject: references; bh=YWyEk+rew3qIyCjYn2JYcbbPrY1WKy30L8ybgI1RV2U=; b=B4Khawmbpnd2tgvuj8bDagB0uB5wkfa2EM8LtUbQ8vTiludjUL4osCEs Bv0rOYcoYbkR893pW3UEeXC5wBiBC8G2du34/39tzG66NJGZKaBBxtL8K EvORfZedZT5JgA/HC1tCEEUWe3BPan9owd429mg6Q+rIhX58FVKzWdW08 M=; IronPort-Data: =?us-ascii?q?A9a23=3Aq+wGyKIRjcqXBAvDFE+RW5clxSXFcZb7ZxGrkP8?= =?us-ascii?q?bfHC+12l20jdTm2saWmiBPv+KMWDxLdgkYIjj9UkB75/VzINqS1BcGVNFHysb8?= =?us-ascii?q?5KdbTi6Bh6tZH3KdpWroHqKXqzyU/GYRCwPZiKa9kfF3oTJ9yEmj/nRHuSkUoY?= =?us-ascii?q?oBwgqLeNaYHZ44f5cs75h6mJYqYDR7zKl4bsekeWHULOW82Ic3lYv1k62gEgHU?= =?us-ascii?q?MIeF98vlgdWifhj5DcynpSOZX4VDfnZw3DQGuG4EgMmLtsvwo1V/kuBl/ssIs+?= =?us-ascii?q?il7/nfyXmQJaLYFLI2iMQAvDk30EqSi8ai87XMNIkYFpTzQeImtV80tBEs5qYS?= =?us-ascii?q?AEzP6SKlv51vxxwSn0mZvcWqeSvzX+X9Jb7I1f9W2DtyfRiJFwoPIZe8e8fKUN?= =?us-ascii?q?D/PpdDTEGbh2HwcerxLa2UMFhgN9mKNjsOsURoBlIxjHDDO06RojDBaHD+MNVx?= =?us-ascii?q?h87htpPBurZf8pcbTNycBnaJRpVNT8/Dp8hmOairmL4bidTpU2cv6dx5HLcpCR?= =?us-ascii?q?y3bnxMML9edGLWcVN2E2fowru/WP0BjkAKNGa0nyC6H3EuwNltUsXQ6oOE7ix6?= =?us-ascii?q?vMshFqJx2VVBgd+aLdymtHh4mbWZj6VAxV8FvISkJUP?= IronPort-HdrOrdr: =?us-ascii?q?A9a23=3AVd8xiK3nDgLhGj6hDKDOyQqjBK8kLtp133Aq?= =?us-ascii?q?2lEZdPU1SL36qynApoV/6faZslgssRIb+exoRpPwI080nKQdieJ6AV7FZmjbUQ?= =?us-ascii?q?CTQL2Kg7GO/9SZIULDytI=3D?= X-IronPort-AV: E=Sophos;i="5.88,267,1635199200"; d="scan'208";a="14032893" Received: from tomate.loria.fr (HELO tomate) ([152.81.10.51]) by mail2-relais-roc.national.inria.fr with ESMTP/TLS/DHE-RSA-AES256-GCM-SHA384; 06 Jan 2022 16:33:12 +0100 Date: Thu, 06 Jan 2022 16:33:11 +0100 Message-Id: From: Paul Zimmermann To: joel@rtems.org Cc: keithp@keithp.com, newlib@sourceware.org, christoph.lauter@christoph-lauter.org, Jean-Michel.Muller@ens-lyon.fr, sibid@uvic.ca In-Reply-To: (message from Joel Sherrill on Thu, 6 Jan 2022 08:51:40 -0600) Subject: Re: correctly rounded mathematical functions References: <87mtkbbaiu.fsf@keithp.com> X-Spam-Status: No, score=-3.1 required=5.0 tests=BAYES_00, DKIM_SIGNED, DKIM_VALID, DKIM_VALID_AU, DKIM_VALID_EF, RCVD_IN_MSPIKE_H2, SPF_HELO_NONE, SPF_PASS, TXREP autolearn=ham autolearn_force=no version=3.4.4 X-Spam-Checker-Version: SpamAssassin 3.4.4 (2020-01-24) 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: Thu, 06 Jan 2022 15:33:14 -0000 Dear Joel, > > as far as I know, four libraries provide long double functions: > > GNU libc, Intel Math Library, OpenLibm and Musl. See Table 11 from > > https://members.loria.fr/PZimmermann/papers/accuracy.pdf. > > Is the FreeBSD long double support related to any of those libraries? > > https://github.com/freebsd/freebsd-src/tree/main/lib/msun after a quick look at the cosl code (https://github.com/freebsd/freebsd-src/blob/main/lib/msun/ld80/k_cosl.c) it seems to be the same as the OpenLibm implementation (openlibm-0.7.5/ld80/k_cosl.c). > Any comments on their implementation? the OpenLibm implementation has huge errors (up to 533 ulps for powl), but they have fixed some of the issues I have reported to them. See also https://github.com/JuliaMath/openlibm/issues/236. Best regards, Paul