From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: X-Spam-Checker-Version: SpamAssassin 3.4.4 (2020-01-24) on inbox.vuxu.org X-Spam-Level: X-Spam-Status: No, score=-2.8 required=5.0 tests=DKIM_ADSP_CUSTOM_MED, DKIM_INVALID,DKIM_SIGNED,FREEMAIL_FORGED_FROMDOMAIN,FREEMAIL_FROM, HEADER_FROM_DIFFERENT_DOMAINS,MAILING_LIST_MULTI,RCVD_IN_DNSWL_MED, RCVD_IN_MSPIKE_H4,RCVD_IN_MSPIKE_WL autolearn=ham autolearn_force=no version=3.4.4 Received: from second.openwall.net (second.openwall.net [193.110.157.125]) by inbox.vuxu.org (Postfix) with SMTP id 57FAE22285 for ; Wed, 29 May 2024 02:38:38 +0200 (CEST) Received: (qmail 18179 invoked by uid 550); 29 May 2024 00:38: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 18147 invoked from network); 29 May 2024 00:38:31 -0000 DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20230601; t=1716943104; x=1717547904; darn=lists.openwall.com; h=content-transfer-encoding:in-reply-to:from:content-language :references:to:subject:user-agent:mime-version:date:message-id:from :to:cc:subject:date:message-id:reply-to; bh=32r0+6GcorC58CivE+0coT/0anBBneoeDiJE+EQIqnw=; b=lTaHRhSE1mvGZjDIs4XZ+b3Mau6u23FLIKo1JKGsJ2jepHxKyM0t1rKJhGvG5BoUXL d5wVGDoaaTwonaj2V0vpzawH8eRcWZdE1zxF3MixzfOcQKRVNDbj5oElZs5DgxBfAPzF KObtZURkB5sRbKEeYdt+V96Hzx5D5M+UBJ80YyQpxNkTLDoioF2X05FTJL2wlGLYvfyj +M+ryu04xYleejdgOzTDz9rLIlY9Jf0AmdKjEkc2ICuN0ETYAb4dhnFLP2FeA5Ui4ezC UA+K0MErgbo9Hn24nSTicTKHbQzYp1BUd81mhzPPr+4l0hG6q5KazhDEYiuPbS9i2CjA aJhg== X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20230601; t=1716943104; x=1717547904; h=content-transfer-encoding:in-reply-to:from:content-language :references:to:subject:user-agent:mime-version:date:message-id :x-gm-message-state:from:to:cc:subject:date:message-id:reply-to; bh=32r0+6GcorC58CivE+0coT/0anBBneoeDiJE+EQIqnw=; b=j6PPeUVVWHb2yrYyNJPBwnd46zvHY3+sA8yKSFYlu1D+CimoEUHUU3IRYm0mhfrZ4H 48zcdYVTzeSYw2PpYveSmplCRnPKOSwCVRMAtd8hlqch0k796bx7+s7IygCTAMWIsEZo Vfm81c/bGj+WiiVlnYZZlqQZl9ATh7GzBS1wCVRi6pIZ5uVXcxAildIY4mNk35Y+12yI 3XB73hEthBbij4/QN/JGgoFUOi48AKukneTWcjCnGIoFeHE8f4PJ5Oy7y5vaN5vLuwkP o3zrddbux79aXM3zVIx8BUG839/VyTX9wXd2Br5YE7noky1iq2SuWgYObLIUp2n8lz+T 8FzQ== X-Gm-Message-State: AOJu0YzP8Hi48kgKcAjMLttLxgcmkpnuxMMOeIW1vIVIgXf7AMWY38Ba T6UetNHGl/W50UItvEAswqRt4bG/1J/xMBPBKQ24WdG4c+L4DHSXkNdX65IA X-Google-Smtp-Source: AGHT+IFe8KiDOfFTMe3T4zPUR4nlOZx22dAi9FeTZxVV9tbhDJNsj9y98Wa0kgBH9JQujM4WXR5hdg== X-Received: by 2002:a17:906:36d5:b0:a62:832a:e2c9 with SMTP id a640c23a62f3a-a62832ae3bcmr1149847466b.67.1716943103668; Tue, 28 May 2024 17:38:23 -0700 (PDT) Message-ID: <724090dc-83e7-4f40-9ff2-1a60196b9966@gmail.com> Date: Wed, 29 May 2024 02:38:22 +0200 MIME-Version: 1.0 User-Agent: Mozilla Thunderbird To: musl@lists.openwall.com, Nikolaos Chatzikonstantinou References: Content-Language: en-US From: Gabriel Ravier In-Reply-To: Content-Type: text/plain; charset=UTF-8; format=flowed Content-Transfer-Encoding: 8bit Subject: Re: [musl] Re: Implementing csqrtl() On 5/28/24 10:48 PM, Nikolaos Chatzikonstantinou wrote: > On Mon, Jul 4, 2022 at 5:35 AM Nikolaos Chatzikonstantinou > wrote: >> Hello list, >> >> I wanted to implement some function from >> >> which is an open issue in the wiki. >> >> One of the missing complex functions is csqrtl(), the long double >> version of complex square root. I was able to find a 1987 article from >> W. Kahan, titled "Branch cuts for complex elementary functions." that >> contained an implementation for complex square root for arbitrary >> floating-point numbers. In this e-mail you'll find an attached git >> patch with the implementation. >> >> As a very basic test, I wrote a program that produces random >> complex numbers in the square [0, N] x [0, N] for N=1,100 and >> calculates csqrt{,f,l}() with my implementation, glibc and the >> arbitrary precision mpc_sqrt() from MPC, >> . >> >> Glibc stays almost within 1 ulp in float and double, but my >> implementation wasn't so good with float. The double >> implementation seems to get the exact same results as glibc does. >> >> I was not able to even test the long double version with this >> method, because I did not write the code that produces random >> long double complex numbers yet. >> >> There's a few things that I don't quite understand here. One is, I'm >> not sure why Kahan's implementation is accurate. For another, I don't >> know how to do any sort of speed tests; I've read online that >> microbenchmarking is not reproducible for math functions, and so >> google's benchmark does not seem >> to help. Of course I'm aware that the C implementation would not be >> the one used in most systems, as the assembly implementations are >> usually better. Finally, I don't know what the right way to test an >> implementation for accuracy is: whether by using automation or writing >> proofs. It seems the state of the art has evolved quite a bit from >> 1987, and yet I don't know where to look for information on this >> topic, as it seems very specific to chips & the C std lib. >> >> Feel free to provide any sort of criticism. I'm e-mailing this >> implementation for the purpose of starting a discussion, but I'm >> hoping to be able to contribute something in the near future. > Hello, I want to revive the discussion on this message. Last time I > attached the patch I was ignored, apart from some nitpicks on the C > source in the patch, which was a draft anyway. Missing complex math > functions are still an open issue. I can't help but think that last > time I was ignored apart from some nitpick on some signed/unsigned > type, as if getting schooled on C is the purpose of the mailing list > -- bizarre! the patch was a draft to start a discussion on resolving > the complex math issue. What I concluded from the responses I received > is that the musl project does not have a maintainer who is capable of > understanding the underlying issues with the complex math algorithms > in libc. Is that still the case and is there any interest in resolving > it perhaps by finding a maintainer who can work on this issue? I cannot read people's minds to determine exactly why it is that they did not thoroughly review the implementation but my guess would be that, unable to understand the mathematical details of how the code worked, they only reviewed the parts of it they were more familiar with (e.g. the "nitpick on some signed/unsigned type") - I guess it makes sense it would seem odd as the only thing to get reviewed at all, but normally those parts would also get reviewed regardless even if people were better able to understand the mathematics of the patch. > > Regards, > Nikolaos Chatzikonstantinou