From patchwork Sat Sep 9 20:41:45 2023 Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: 8bit X-Patchwork-Submitter: Richard Henderson X-Patchwork-Id: 721100 Delivered-To: patch@linaro.org Received: by 2002:adf:f64d:0:b0:31d:da82:a3b4 with SMTP id x13csp412745wrp; Sat, 9 Sep 2023 13:42:44 -0700 (PDT) X-Google-Smtp-Source: AGHT+IEefvjulOyV9ideNoQKfppnj86vhE11fuT8ZpW6cYdbWEY/ETTMSiJr8C989aepqza4Mw6E X-Received: by 2002:a05:620a:bc7:b0:76f:16e9:bea2 with SMTP id s7-20020a05620a0bc700b0076f16e9bea2mr6814183qki.30.1694292164650; Sat, 09 Sep 2023 13:42:44 -0700 (PDT) ARC-Seal: i=1; a=rsa-sha256; t=1694292164; cv=none; d=google.com; s=arc-20160816; b=gRaRmCd2EX+VPGznr9VQuICnQHcYGNW8S2fnTqjLh3oWI8ohkx/zp4srvIlxLl8FPX bAYx21AzVTplIugSZvzMc+a9t8kJ/iZGn/VkH+4X/vRMs8ksZhb03YEm88fr2PangP5J pZdCk+GJZUlz6UW6la+dZpf8/20VECBMKl4kkXXUxBra+Vc/1OvrDeWRn2gNcH32Ei4r nZLJ43uzuRAL/bpkpR77y21kZf2YiTl61nezb3dpOJZp7AEvksIvAmW1x2gupopXYuLm sW6AdYDqXIEvgPRDchEhSEv/AgUtLVClVo1szbNkLQ+8z6oP3eE0eeyexJ6hWYWR2Yd7 2PEg== ARC-Message-Signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=sender:errors-to:list-subscribe:list-help:list-post:list-archive :list-unsubscribe:list-id:precedence:content-transfer-encoding :mime-version:message-id:date:subject:cc:to:from:dkim-signature; bh=TjSfPst91wYzmd/Lt7HYTR7Bfmwse9WbD2zkkLX5J/U=; fh=gmHlaZIMMboJ2MqxTgcMNAggxV68YaGaSfcBwyCSkVg=; b=ZiUHCCjbiqtSK1CxAESLlbcGRl4UuGHZP00nHw7cdO1iGzaSenBDBgBS8/WKNU6fVU d3qlo5M3JMnKUavWvN6Y/ODc37UVfns23dLGhqOGaW0MUo5f6yBpJ502qR8spWA/4vsz qqbniSE3QFegnEUFMBCxl36R18i4RizGGabFB10gYQRtrN6QLbwTWqfo6ruGAXNW6dKe /8Lzbt69ZKC0+p7s4p+bZ9C1r2dJ1q1NZiZsSo7NjD80MCUMleSFO48q9/D1LDyq8TF+ xQLZ46wR+W1vGv8+Gu2S4cKnKVqe1bF9Qv3lfrBiij9VW9N3ORvndAL3tl8Ylx3KNsUm i0HA== ARC-Authentication-Results: i=1; mx.google.com; dkim=pass header.i=@linaro.org header.s=google header.b=EINq9nxq; spf=pass (google.com: domain of qemu-devel-bounces+patch=linaro.org@nongnu.org designates 209.51.188.17 as permitted sender) smtp.mailfrom="qemu-devel-bounces+patch=linaro.org@nongnu.org"; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=linaro.org Return-Path: Received: from lists.gnu.org (lists.gnu.org. [209.51.188.17]) by mx.google.com with ESMTPS id qx6-20020a05620a8b8600b0076cf352b3f9si2630751qkn.761.2023.09.09.13.42.44 for (version=TLS1_2 cipher=ECDHE-ECDSA-CHACHA20-POLY1305 bits=256/256); Sat, 09 Sep 2023 13:42:44 -0700 (PDT) Received-SPF: pass (google.com: domain of qemu-devel-bounces+patch=linaro.org@nongnu.org designates 209.51.188.17 as permitted sender) client-ip=209.51.188.17; Authentication-Results: mx.google.com; dkim=pass header.i=@linaro.org header.s=google header.b=EINq9nxq; spf=pass (google.com: domain of qemu-devel-bounces+patch=linaro.org@nongnu.org designates 209.51.188.17 as permitted sender) smtp.mailfrom="qemu-devel-bounces+patch=linaro.org@nongnu.org"; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=linaro.org Received: from localhost ([::1] helo=lists1p.gnu.org) by lists.gnu.org with esmtp (Exim 4.90_1) (envelope-from ) id 1qf4m9-0005Nj-Vh; Sat, 09 Sep 2023 16:41:54 -0400 Received: from eggs.gnu.org ([2001:470:142:3::10]) by lists.gnu.org with esmtps (TLS1.2:ECDHE_RSA_AES_256_GCM_SHA384:256) (Exim 4.90_1) (envelope-from ) id 1qf4m8-0005CT-3x for qemu-devel@nongnu.org; Sat, 09 Sep 2023 16:41:52 -0400 Received: from mail-oi1-x22c.google.com ([2607:f8b0:4864:20::22c]) by eggs.gnu.org with esmtps (TLS1.2:ECDHE_RSA_AES_128_GCM_SHA256:128) (Exim 4.90_1) (envelope-from ) id 1qf4m5-0005rU-M6 for qemu-devel@nongnu.org; Sat, 09 Sep 2023 16:41:51 -0400 Received: by mail-oi1-x22c.google.com with SMTP id 5614622812f47-3ab29af398aso2296849b6e.1 for ; Sat, 09 Sep 2023 13:41:48 -0700 (PDT) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=linaro.org; s=google; t=1694292107; x=1694896907; darn=nongnu.org; h=content-transfer-encoding:mime-version:message-id:date:subject:cc :to:from:from:to:cc:subject:date:message-id:reply-to; bh=TjSfPst91wYzmd/Lt7HYTR7Bfmwse9WbD2zkkLX5J/U=; b=EINq9nxqpDjWnrSD9aMpz5WJXJsN9garayHBnKtjxEb2y+0qF1zRIOQOEZzNu+qIKu fY2KsDXIZ0wA/Asm6xM4IYm2L9SpCYFO6MSTWWjZYqG83vY7cPVv4vcIAcC06lFKtq7Q /eaht8JFkDLEqJt62qx5Zp4PNSwFpVzAo/Gwe93YELMj1DkMPQRJ/aPtcSbbRi/81xeu wGbO9KvjoQMMUUN/GpJSoTur08ai+3iUBQ9PVys0/2W0orp5JI+6ptk53OCfoMgm82FA 8l/8xeaifgqkZ/WHRmcW8mAWi19Toa1owR691+Phg09xAZ/m/c6LqKa+qZfPvDpt/7p2 USeA== X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20230601; t=1694292107; x=1694896907; h=content-transfer-encoding:mime-version:message-id:date:subject:cc :to:from:x-gm-message-state:from:to:cc:subject:date:message-id :reply-to; bh=TjSfPst91wYzmd/Lt7HYTR7Bfmwse9WbD2zkkLX5J/U=; b=gY0s0P/Doe/PirHJkYHLZLBQJS+qTcQ7E9tagcatPjvyPWLkSYX4Gkd/9XDXXlVU9M 1iNSQ4UoLzYgO0K4KKlzPpdFtqdX6OHs7NeNC6JA72z6zBZszobIwJ6n35pLQ51DQhKn MY56IMTxeTHFo+rOZa9mfdPQzKXcaIc6k9HuPhYYkWIHLudKOl+6LUuZbUOV833eMPj/ j9vgZmn0lIJ+fW/FnjmvHugba2ZvK8x3GJJ4R2X1sNoXeDCRQDVKpAMjmfMaVZXAB0ei aiIkb36dmaoZ5WzsYl/BnvwLGgMzgDW2MlJ0FB/JNRC/mL0/tL/KlDEit9fx5O0hrwoj 6rlA== X-Gm-Message-State: AOJu0YyFlv2nswTnaLfZG5omtin9nXXHRylL5Jn6//8Fk3Z5SyscyITf gOkIOS8BP222uIkBdXgvu2jEzlAK5FFGigN8y/k= X-Received: by 2002:a54:4015:0:b0:3a3:e6d4:69d6 with SMTP id x21-20020a544015000000b003a3e6d469d6mr6300618oie.7.1694292107396; Sat, 09 Sep 2023 13:41:47 -0700 (PDT) Received: from stoup.. ([71.212.131.115]) by smtp.gmail.com with ESMTPSA id y17-20020a637d11000000b00569fd817b73sm3049155pgc.48.2023.09.09.13.41.46 (version=TLS1_3 cipher=TLS_AES_256_GCM_SHA384 bits=256/256); Sat, 09 Sep 2023 13:41:46 -0700 (PDT) From: Richard Henderson To: qemu-devel@nongnu.org Cc: Keith Packard , =?utf-8?q?Philippe_Mathieu-Daud=C3=A9?= Subject: [PATCH v2] softfloat: Handle m68k extended precision denormals properly Date: Sat, 9 Sep 2023 13:41:45 -0700 Message-Id: <20230909204145.37548-1-richard.henderson@linaro.org> X-Mailer: git-send-email 2.34.1 MIME-Version: 1.0 Received-SPF: pass client-ip=2607:f8b0:4864:20::22c; envelope-from=richard.henderson@linaro.org; helo=mail-oi1-x22c.google.com X-Spam_score_int: -20 X-Spam_score: -2.1 X-Spam_bar: -- X-Spam_report: (-2.1 / 5.0 requ) BAYES_00=-1.9, DKIM_SIGNED=0.1, DKIM_VALID=-0.1, DKIM_VALID_AU=-0.1, DKIM_VALID_EF=-0.1, RCVD_IN_DNSWL_NONE=-0.0001, SPF_HELO_NONE=0.001, SPF_PASS=-0.001 autolearn=ham autolearn_force=no X-Spam_action: no action X-BeenThere: qemu-devel@nongnu.org X-Mailman-Version: 2.1.29 Precedence: list List-Id: List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , Errors-To: qemu-devel-bounces+patch=linaro.org@nongnu.org Sender: qemu-devel-bounces+patch=linaro.org@nongnu.org Motorola treats denormals with explicit integer bit set as having unbiased exponent 0, unlike Intel which treats it as having unbiased exponent 1 (more like all other IEEE formats that have no explicit integer bit). Add a flag on FloatFmt to differentiate the behaviour. Reported-by: Keith Packard Reviewed-by: Philippe Mathieu-Daudé Signed-off-by: Richard Henderson --- v2: Incorporate Keith's identity test and fix. --- fpu/softfloat.c | 9 +++++- tests/tcg/m68k/denormal.c | 53 ++++++++++++++++++++++++++++++++++ fpu/softfloat-parts.c.inc | 7 +++-- tests/tcg/m68k/Makefile.target | 2 +- 4 files changed, 66 insertions(+), 5 deletions(-) create mode 100644 tests/tcg/m68k/denormal.c diff --git a/fpu/softfloat.c b/fpu/softfloat.c index 0cc130ae9b..3adfa8cee0 100644 --- a/fpu/softfloat.c +++ b/fpu/softfloat.c @@ -517,6 +517,7 @@ typedef struct { * round_mask: bits below lsb which must be rounded * The following optional modifiers are available: * arm_althp: handle ARM Alternative Half Precision + * m68k_denormal: explicit integer bit for extended precision may be 1 */ typedef struct { int exp_size; @@ -526,6 +527,7 @@ typedef struct { int frac_size; int frac_shift; bool arm_althp; + bool m68k_denormal; uint64_t round_mask; } FloatFmt; @@ -576,7 +578,12 @@ static const FloatFmt float128_params = { static const FloatFmt floatx80_params[3] = { [floatx80_precision_s] = { FLOATX80_PARAMS(23) }, [floatx80_precision_d] = { FLOATX80_PARAMS(52) }, - [floatx80_precision_x] = { FLOATX80_PARAMS(64) }, + [floatx80_precision_x] = { + FLOATX80_PARAMS(64), +#ifdef TARGET_M68K + .m68k_denormal = true, +#endif + }, }; /* Unpack a float to parts, but do not canonicalize. */ diff --git a/tests/tcg/m68k/denormal.c b/tests/tcg/m68k/denormal.c new file mode 100644 index 0000000000..20bd8c7332 --- /dev/null +++ b/tests/tcg/m68k/denormal.c @@ -0,0 +1,53 @@ +/* + * Test m68k extended double denormals. + */ + +#include +#include + +#define TEST(X, Y) { X, Y, X * Y } + +static volatile long double test[][3] = { + TEST(0x1p+16383l, 0x1p-16446l), + TEST(0x1.1p-8223l, 0x1.1p-8224l), + TEST(1.0l, 0x1p-16383l), +}; + +#undef TEST + +static void dump_ld(const char *label, long double ld) +{ + union { + long double d; + struct { + uint32_t exp:16; + uint32_t space:16; + uint32_t h; + uint32_t l; + }; + } u; + + u.d = ld; + printf("%12s: % -27La 0x%04x 0x%08x 0x%08x\n", label, u.d, u.exp, u.h, u.l); +} + +int main(void) +{ + int i, n = sizeof(test) / sizeof(test[0]), err = 0; + + for (i = 0; i < n; ++i) { + long double x = test[i][0]; + long double y = test[i][1]; + long double build_mul = test[i][2]; + long double runtime_mul = x * y; + + if (runtime_mul != build_mul) { + dump_ld("x", x); + dump_ld("y", y); + dump_ld("build_mul", build_mul); + dump_ld("runtime_mul", runtime_mul); + err = 1; + } + } + return err; +} diff --git a/fpu/softfloat-parts.c.inc b/fpu/softfloat-parts.c.inc index 527e15e6ab..a44649f4f4 100644 --- a/fpu/softfloat-parts.c.inc +++ b/fpu/softfloat-parts.c.inc @@ -118,7 +118,8 @@ static void partsN(canonicalize)(FloatPartsN *p, float_status *status, } else { int shift = frac_normalize(p); p->cls = float_class_normal; - p->exp = fmt->frac_shift - fmt->exp_bias - shift + 1; + p->exp = fmt->frac_shift - fmt->exp_bias + - shift + !fmt->m68k_denormal; } } else if (likely(p->exp < fmt->exp_max) || fmt->arm_althp) { p->cls = float_class_normal; @@ -256,7 +257,7 @@ static void partsN(uncanon_normal)(FloatPartsN *p, float_status *s, is_tiny = !frac_addi(&discard, p, inc); } - frac_shrjam(p, 1 - exp); + frac_shrjam(p, !fmt->m68k_denormal - exp); if (p->frac_lo & round_mask) { /* Need to recompute round-to-even/round-to-odd. */ @@ -287,7 +288,7 @@ static void partsN(uncanon_normal)(FloatPartsN *p, float_status *s, p->frac_lo &= ~round_mask; } - exp = (p->frac_hi & DECOMPOSED_IMPLICIT_BIT) != 0; + exp = (p->frac_hi & DECOMPOSED_IMPLICIT_BIT) && !fmt->m68k_denormal; frac_shr(p, frac_shift); if (is_tiny && (flags & float_flag_inexact)) { diff --git a/tests/tcg/m68k/Makefile.target b/tests/tcg/m68k/Makefile.target index 1163c7ef03..6ff214e60a 100644 --- a/tests/tcg/m68k/Makefile.target +++ b/tests/tcg/m68k/Makefile.target @@ -4,7 +4,7 @@ # VPATH += $(SRC_PATH)/tests/tcg/m68k -TESTS += trap +TESTS += trap denormal # On m68k Linux supports 4k and 8k pages (but 8k is currently broken) EXTRA_RUNS+=run-test-mmap-4096 # run-test-mmap-8192