From patchwork Thu Oct 6 15:01:37 2022 Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: 8bit X-Patchwork-Submitter: Letu Ren X-Patchwork-Id: 58476 Return-Path: X-Original-To: patchwork@sourceware.org Delivered-To: patchwork@sourceware.org Received: from server2.sourceware.org (localhost [IPv6:::1]) by sourceware.org (Postfix) with ESMTP id 67659383FB88 for ; Thu, 6 Oct 2022 15:03:13 +0000 (GMT) DKIM-Filter: OpenDKIM Filter v2.11.0 sourceware.org 67659383FB88 DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=sourceware.org; s=default; t=1665068593; bh=bDYxY1NT+IfDSf33XKfZvHiAJ8ywMEBBfE/tE5Ro2gI=; h=To:Subject:Date:List-Id:List-Unsubscribe:List-Archive:List-Post: List-Help:List-Subscribe:From:Reply-To:Cc:From; b=KIWzljUNv6z91qF9+Fe8wmMI4Ol+x7bAJeTZ7ovQCC6rHgpwsjh9kZF8M9rTj4+8T 23nRPnSO5KbX8ZlRp6UfUCpW2topDGRmn9f1JCR36dkSakivwe1hpTzpfPo8QJdzyz qbhx1irbyP5uZm0QmoVF8tRd669tabKDNxtqKcUk= X-Original-To: libc-alpha@sourceware.org Delivered-To: libc-alpha@sourceware.org Received: from mail-pj1-x102c.google.com (mail-pj1-x102c.google.com [IPv6:2607:f8b0:4864:20::102c]) by sourceware.org (Postfix) with ESMTPS id 4C7403898380 for ; Thu, 6 Oct 2022 15:01:54 +0000 (GMT) DMARC-Filter: OpenDMARC Filter v1.4.1 sourceware.org 4C7403898380 Received: by mail-pj1-x102c.google.com with SMTP id a5-20020a17090aa50500b002008eeb040eso5117690pjq.1 for ; Thu, 06 Oct 2022 08:01:54 -0700 (PDT) X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20210112; h=content-transfer-encoding:mime-version:message-id:date:subject:cc :to:from:x-gm-message-state:from:to:cc:subject:date; bh=bDYxY1NT+IfDSf33XKfZvHiAJ8ywMEBBfE/tE5Ro2gI=; b=u/JoRYnWa4xuRihaWS1jWtKFXPEqZGnGNAJXOJ+0YmTP58tH5FOVDUCsnP5Z6q5/5Z w1RoQq9gmXuGAacuBTi8iYk4cpw/l5uwB3YlQzOuXWaGl9eOJmMp+GCj+tUH5PnABecX oV+cNac+63IDwO9/wOeWBiDjhoobR7l7i1lO28ki1pTTeyVwU6Av+e3qMdpwlZeKDG8p QYqa2eCvC+cgSsF0HtTA+R6lMCkyzvxOfiX4JCAJXLWVp+CjpoiW8J8PrAch1TLZHHAl Qt6acBPfHOiLIXSHCIBb48GnY9NL14Z3FcRzwt8Lv91zhcFAXg3SlNwNJLuiCNcUYAwg Wkjg== X-Gm-Message-State: ACrzQf0rWXLevPsylOkL/3qlcoLosiUBhauDS0Ma6FJjZNCuIMouQyrM HvluX682uOttu6eZ9c9d7RRkUbYrfKrrAL9eA/A= X-Google-Smtp-Source: AMsMyM5cVqQ68TbyROOT2f4YtwJA4+0vCEXnJ/cVgWC8D2ovyMO37rLatb8EgdqJZklZmyppsLZptA== X-Received: by 2002:a17:90b:384a:b0:20a:dbe3:25d7 with SMTP id nl10-20020a17090b384a00b0020adbe325d7mr11241979pjb.85.1665068511452; Thu, 06 Oct 2022 08:01:51 -0700 (PDT) Received: from fanta-arch.tsinghua.edu.cn (ec2-18-162-45-232.ap-east-1.compute.amazonaws.com. [18.162.45.232]) by smtp.gmail.com with ESMTPSA id k15-20020a170902ce0f00b0017d7e5a9fa7sm10219944plg.92.2022.10.06.08.01.48 (version=TLS1_3 cipher=TLS_AES_256_GCM_SHA384 bits=256/256); Thu, 06 Oct 2022 08:01:50 -0700 (PDT) To: libc-alpha@sourceware.org Subject: [PATCH v5] stdlib/strfrom: Add copysign to fix NAN issue on riscv (BZ #29501) Date: Thu, 6 Oct 2022 23:01:37 +0800 Message-Id: <20221006150137.18691-1-fantasquex@gmail.com> X-Mailer: git-send-email 2.38.0 MIME-Version: 1.0 X-Spam-Status: No, score=-10.6 required=5.0 tests=BAYES_00, DKIM_SIGNED, DKIM_VALID, DKIM_VALID_AU, DKIM_VALID_EF, FREEMAIL_FROM, GIT_PATCH_0, RCVD_IN_DNSWL_NONE, SPF_HELO_NONE, SPF_PASS, TXREP autolearn=ham autolearn_force=no version=3.4.6 X-Spam-Checker-Version: SpamAssassin 3.4.6 (2021-04-09) on server2.sourceware.org X-BeenThere: libc-alpha@sourceware.org X-Mailman-Version: 2.1.29 Precedence: list List-Id: Libc-alpha mailing list List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-Patchwork-Original-From: Letu Ren via Libc-alpha From: Letu Ren Reply-To: Letu Ren Cc: fweimer@redhat.com, schwab@suse.de, Letu Ren , joseph@codesourcery.com Errors-To: libc-alpha-bounces+patchwork=sourceware.org@sourceware.org Sender: "Libc-alpha" According to the specification of ISO/IEC TS 18661-1:2014, The strfromd, strfromf, and strfroml functions are equivalent to snprintf(s, n, format, fp) (7.21.6.5), except the format string contains only the character %, an optional precision that does not contain an asterisk *, and one of the conversion specifiers a, A, e, E, f, F, g, or G, which applies to the type (double, float, or long double) indicated by the function suffix (rather than by a length modifier). Use of these functions with any other 20 format string results in undefined behavior. strfromf will convert the arguement with type float to double first. According to the latest version of IEEE754 which is published in 2019, Conversion of a quiet NaN from a narrower format to a wider format in the same radix, and then back to the same narrower format, should not change the quiet NaN payload in any way except to make it canonical. When either an input or result is a NaN, this standard does not interpret the sign of a NaN. However, operations on bit strings—copy, negate, abs, copySign—specify the sign bit of a NaN result, sometimes based upon the sign bit of a NaN operand. The logical predicates totalOrder and isSignMinus are also affected by the sign bit of a NaN operand. For all other operations, this standard does not specify the sign bit of a NaN result, even when there is only one input NaN, or when the NaN is produced from an invalid operation. converting NAN or -NAN with type float to double doesn't need to keep the signbit. As a result, this test case isn't mandatory. The problem is that according to RISC-V ISA manual in chapter 11.3 of riscv-isa-20191213, Except when otherwise stated, if the result of a floating-point operation is NaN, it is the canonical NaN. The canonical NaN has a positive sign and all significand bits clear except the MSB, a.k.a. the quiet bit. For single-precision floating-point, this corresponds to the pattern 0x7fc00000. which means that conversion -NAN from float to double won't keep the signbit. Since glibc ought to be consistent here between types and architectures, this patch adds copysign to fix this problem if the string is NAN. This patch adds two different functions under sysdeps directory to work around the issue. This patch has been tested on x86_64 and riscv64. Resolves: BZ #29501 v2: Change from macros to different inline functions. v3: Add unlikely check to isnan. v4: Fix wrong commit message header. v5: Fix style: add space before parentheses. Signed-off-by: Letu Ren --- I'm not sure about copyright of the new file I created and the name of the new function. Please give me some advice. Thanks. stdlib/strfrom-skeleton.c | 3 ++- .../generic/fix-float-double-convert-nan.h | 15 +++++++++++++ .../riscv/rvd/fix-float-double-convert-nan.h | 21 +++++++++++++++++++ 3 files changed, 38 insertions(+), 1 deletion(-) create mode 100644 sysdeps/generic/fix-float-double-convert-nan.h create mode 100644 sysdeps/riscv/rvd/fix-float-double-convert-nan.h diff --git a/stdlib/strfrom-skeleton.c b/stdlib/strfrom-skeleton.c index 1fba04bf6a..36e9adcad5 100644 --- a/stdlib/strfrom-skeleton.c +++ b/stdlib/strfrom-skeleton.c @@ -27,6 +27,7 @@ #include #include #include +#include #define UCHAR_T char #define L_(Str) Str @@ -61,7 +62,7 @@ STRFROM (char *dest, size_t size, const char *format, FLOAT f) because __printf_fp and __printf_fphex only accept double and long double as the floating-point argument. */ if (__builtin_types_compatible_p (FLOAT, float)) - fpnum.flt = f; + fpnum.flt = keep_sign_conversion (f); else fpnum.value = f; diff --git a/sysdeps/generic/fix-float-double-convert-nan.h b/sysdeps/generic/fix-float-double-convert-nan.h new file mode 100644 index 0000000000..0949dda102 --- /dev/null +++ b/sysdeps/generic/fix-float-double-convert-nan.h @@ -0,0 +1,15 @@ +/* Fix for conversion of float NAN to double. Generic version. */ + +#ifndef FIX_FLOAT_DOUBLE_CONVERT_NAN_H +#define FIX_FLOAT_DOUBLE_CONVERT_NAN_H + +/* This function aims to work around conversions of float -NAN + to double returning NAN instead of the correct -NAN in some + architectures. */ +static inline double __attribute__ ((always_inline)) +keep_sign_conversion (float flt) +{ + return flt; +} + +#endif diff --git a/sysdeps/riscv/rvd/fix-float-double-convert-nan.h b/sysdeps/riscv/rvd/fix-float-double-convert-nan.h new file mode 100644 index 0000000000..59806052a8 --- /dev/null +++ b/sysdeps/riscv/rvd/fix-float-double-convert-nan.h @@ -0,0 +1,21 @@ +/* Fix for conversion of float NAN to double. RISC-V version. */ + +#ifndef FIX_FLOAT_DOUBLE_CONVERT_NAN_H +#define FIX_FLOAT_DOUBLE_CONVERT_NAN_H + +#include + +/* RISC-V rvd instructions do not preserve the signbit of NAN + when converting from float to double. */ +static inline double +keep_sign_conversion (float flt) +{ + if (__glibc_unlikely (isnan (flt))) + { + float x = copysignf (1.f, flt); + return copysign ((double) flt, (double) x); + } + return flt; +} + +#endif