From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: from mx0a-00069f02.pphosted.com (mx0a-00069f02.pphosted.com [205.220.165.32]) by sourceware.org (Postfix) with ESMTPS id 30719383A831 for ; Thu, 12 Aug 2021 16:03:33 +0000 (GMT) DMARC-Filter: OpenDMARC Filter v1.4.1 sourceware.org 30719383A831 Received: from pps.filterd (m0246617.ppops.net [127.0.0.1]) by mx0b-00069f02.pphosted.com (8.16.0.43/8.16.0.43) with SMTP id 17CG1x0i021749; Thu, 12 Aug 2021 16:03:32 GMT Received: from aserp3030.oracle.com (aserp3030.oracle.com [141.146.126.71]) by mx0b-00069f02.pphosted.com with ESMTP id 3ad2ajgqux-1 (version=TLSv1.2 cipher=ECDHE-RSA-AES256-GCM-SHA384 bits=256 verify=OK); Thu, 12 Aug 2021 16:03:31 +0000 Received: from pps.filterd (aserp3030.oracle.com [127.0.0.1]) by aserp3030.oracle.com (8.16.0.42/8.16.0.42) with SMTP id 17CG0dGA081001; Thu, 12 Aug 2021 16:03:30 GMT Received: from nam10-mw2-obe.outbound.protection.outlook.com (mail-mw2nam10lp2109.outbound.protection.outlook.com [104.47.55.109]) by aserp3030.oracle.com with ESMTP id 3abx3y3m9f-1 (version=TLSv1.2 cipher=ECDHE-RSA-AES256-GCM-SHA384 bits=256 verify=OK); Thu, 12 Aug 2021 16:03:30 +0000 ARC-Seal: i=1; a=rsa-sha256; s=arcselector9901; d=microsoft.com; cv=none; b=PHGQAwaREoUlhIWh2/MK1ovFTrro0tS5VTbPfqs40Uo/j0nehXw8u1TXdEWrAyoAdVrcYNjd/dDESpatwPoBb1mM4MEnjFcpm+w69Sgbi6WvBM1HaibXBqLoHL2sHz8ocE6rAMNT6vwI68gvGTWhw5BuJI6CHBNqfBDTynOum+pn4w8ZzidTeXNIPro2bB2VgHP2GO98Ci4SYqSZ1tPiU757247gkmblOh2LJ1UffR4+MMhv9dYek5cl4831Qjv3sd5qDn1Qc+BBxlWvpDFWk+hZZi9Jp8LtJ6ftcMNvS8nHJTLCMxRf9oq4GH3Ci5+/gnruN9ywvvyK+bnTr38Ajw== ARC-Message-Signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=microsoft.com; s=arcselector9901; h=From:Date:Subject:Message-ID:Content-Type:MIME-Version:X-MS-Exchange-SenderADCheck; bh=xgvOl1N2Md9R/VevUtgl6PYFh4Q3E4mRCyOjlSZgu3o=; b=nByf69THIkl+HqsHHYCA27tMELOZWMKvQNSt04ZB92GpbGVK9aPyx8+F2bqhViXwdrhHjWOeepPgvAkvIANznu/3cMxqbYxxW+KsRwVqWrsnVKqKsqXLiVTdMJet0yQ3djEHGIwF4JYeOca638DhxUyPEqXbafhdIOgKRx14eLLWKk9KQ56BWGgTBFprxEMkTg4AbRDls4Pcz5/u+TMnBftqENoa1JiqhjYdvJuVYfrd3nSTs4PlefCxUxF21AizwKhJmkTdV2OowF6T05OlJejkMAP8CMGQ3t6Y0KObfC0UDzFKc4Pl7J903/dpBbTgEbikjqaDJjWU7DRW+TnoYg== ARC-Authentication-Results: i=1; mx.microsoft.com 1; spf=pass smtp.mailfrom=oracle.com; dmarc=pass action=none header.from=oracle.com; dkim=pass header.d=oracle.com; arc=none Received: from BYAPR10MB3208.namprd10.prod.outlook.com (2603:10b6:a03:159::10) by BY5PR10MB4065.namprd10.prod.outlook.com (2603:10b6:a03:1f6::15) with Microsoft SMTP Server (version=TLS1_2, cipher=TLS_ECDHE_RSA_WITH_AES_256_GCM_SHA384) id 15.20.4394.18; Thu, 12 Aug 2021 16:03:28 +0000 Received: from BYAPR10MB3208.namprd10.prod.outlook.com ([fe80::7c09:ef39:34f2:e79c]) by BYAPR10MB3208.namprd10.prod.outlook.com ([fe80::7c09:ef39:34f2:e79c%4]) with mapi id 15.20.4394.023; Thu, 12 Aug 2021 16:03:28 +0000 From: Patrick McGehearty To: gcc-patches@gcc.gnu.org Cc: segher@kernel.crashing.org Subject: [PATCH v3] Fix for powerpc64 long double complex divide failure Date: Thu, 12 Aug 2021 16:03:13 +0000 Message-Id: <1628784193-9675-1-git-send-email-patrick.mcgehearty@oracle.com> X-Mailer: git-send-email 1.8.3.1 Content-Type: text/plain X-ClientProxiedBy: SJ0PR03CA0087.namprd03.prod.outlook.com (2603:10b6:a03:331::32) To BYAPR10MB3208.namprd10.prod.outlook.com (2603:10b6:a03:159::10) MIME-Version: 1.0 X-MS-Exchange-MessageSentRepresentingType: 1 Received: from localhost.us.oracle.com (148.87.23.9) by SJ0PR03CA0087.namprd03.prod.outlook.com (2603:10b6:a03:331::32) with Microsoft SMTP Server (version=TLS1_2, cipher=TLS_ECDHE_RSA_WITH_AES_256_GCM_SHA384) id 15.20.4415.15 via Frontend Transport; Thu, 12 Aug 2021 16:03:27 +0000 X-MS-PublicTrafficType: Email X-MS-Office365-Filtering-Correlation-Id: af7d948b-bbc9-4d47-8c5e-08d95daab8fd X-MS-TrafficTypeDiagnostic: BY5PR10MB4065: X-Microsoft-Antispam-PRVS: X-MS-Oob-TLC-OOBClassifiers: OLM:9508; X-MS-Exchange-SenderADCheck: 1 X-MS-Exchange-AntiSpam-Relay: 0 X-Microsoft-Antispam: BCL:0; X-Microsoft-Antispam-Message-Info: 8rQFjQCUV5UwrwYJdHXLpQgt57V2gBUB3VD/uzwYlYZ1JaYJ0OxDJVOtytO9EgQPfyEa/Mz8AkiaMsvn+iobEE/rldGzJpfSsskGEbGJ0Fmoe8dedBURz35wiX9KJhVuGhN0M9qICu2exWR2IZ64YHvroMUCrY8BdwLHwernP4O49v5Dh3QJw3I1UpVbAov/s7wUkwgkxEXnf0jDZ4u4Q6+983JWuOPm6l9MhcXjaPp36M+gqfqLF25XxX0S0q5ec4DJBj8Of8pZgNXYAissTvA5hQXO7hmcMB39NsD1Pr0m8mZg8sX+wixC4Betp6BHwDJ4Tf53UGDczG4LrOB+WZDmdaiRG05LppYSWvvAF1N8kVwbP0rIrHRlc+uqrDDrgGQPlNfSIQtVZO8wQXo+G3LJwiD2CjZNI36ynJeFvTi/IHsb8gv4P6qpCQ9gCaO/W1sLGmFreP+Y+xMdz+H+ZTbaWaplSQ7d26k3IRk2I7S8XU2WK89+uVKAc2hHYKeFIau/tOJPETZ4Y8NsauMZyRxIqL/9eMdGg/7A+Rt4jl8eHEoWkMss2m4D/GafOMJxI6dD/QQ7JKRP0bkMMtL4trlPvon0oEm061HgvrxaytglzPrMDXZrJ0oLwRBH8OSnpvXbg73TBif1I0E2SvYFLH5jA58X4PF5KhCloms1H5RN6aDQRxARuU1M+HZjt7cAKjrUSKEGfUhDA+edKHl6y28KCckt9YXLYL8Yw4d1iry9PdzqQBkH9nVsckzIMzE6bfQp3x5p60NuUPWk86K9yOfE113Tj3hI4ge10J4boGc= X-Forefront-Antispam-Report: CIP:255.255.255.255; CTRY:; LANG:en; SCL:1; SRV:; IPV:NLI; SFV:NSPM; H:BYAPR10MB3208.namprd10.prod.outlook.com; PTR:; CAT:NONE; SFS:(136003)(376002)(366004)(39860400002)(346002)(396003)(36756003)(8936002)(966005)(66946007)(66476007)(38100700002)(7696005)(66556008)(38350700002)(478600001)(52116002)(8676002)(2906002)(2616005)(6666004)(956004)(316002)(83380400001)(86362001)(5660300002)(186003)(4326008)(6486002)(26005)(6916009); DIR:OUT; SFP:1101; X-MS-Exchange-AntiSpam-MessageData-ChunkCount: 1 X-MS-Exchange-AntiSpam-MessageData-0: =?us-ascii?Q?n7bq3ugIoHG4IGnq/K5SKxsEfaVz/vnltxowo0cHcu/FyAIquFT1riMl+SiM?= =?us-ascii?Q?i1wKGVB6DdJ6Nl39M4FDcw0EE2svHBTyl3W6leBgcuFBxj32hrTE23Za2r3P?= =?us-ascii?Q?Z0NMWMD7FL876dm0+PBdpLHKna1QGfxjyNLoZWHmFk0cgVOU5O2im5PycGYz?= =?us-ascii?Q?QhhkDowkTQlxBMVQWcD7WBj9EYKbGwlpi+xAf20OWBnBKMe0DOmcqHdN5PEr?= =?us-ascii?Q?VHbKo/H817eYniAzcT5DbieW3s4b671jZEGXvm/p5Oye1shNGep9IqvSzWLv?= =?us-ascii?Q?0EzjIQtyJcXKj/BStpdiwo4XR89MmbQ+qMipNIx3n8AS52mox2pnWQ4/gdji?= =?us-ascii?Q?7g4wmbUW1y+FPXFHTu4XTOgHp3AdEq7h34Nnx5L9J4IYVqKNhpRA8wa40MqS?= =?us-ascii?Q?NqCgKgs3nyV6EToYiiXk5wgnNWc4PrD64IICvASMB0u515+acV9lm82m4z/n?= =?us-ascii?Q?B0YfFDAoijQpzFWFeILsE+yXNUGZWwsLIXi9skpICPNhwa5anzLmlXL81bRm?= =?us-ascii?Q?SFAcWSC2kNds7hmkf5pGfoISSZWHySGWvsZUYJjXuKWhXwwHFNQxs8cuufNy?= =?us-ascii?Q?LRF2UKUyExCgH8hzNWlilwrmU3OBUTmym0CPCS8swTMOT11KqFUTGHvTugw/?= =?us-ascii?Q?bB+lCVsHkOTft4UENY4ETOrqwMXJ0uBxewvEFhS3MaJmKPVFaDz19q26ahuo?= =?us-ascii?Q?+0EQfUql0BsqY0OjnIjsHBMMoKok3yixJ5ASQ6Nydx3o4P3coDHgsjrkTQix?= =?us-ascii?Q?hIZifVU37RTZS+EQLyAsqAjxeLJrzRozywGyy63zB/+4Fldy9Q3+uCQrXMWe?= =?us-ascii?Q?YPwTXiKBmwV0BzEa7OdVaqG35P0OJxYYGTFHcDyAoWkrHPKC3N8+bDZUXSa4?= =?us-ascii?Q?lgYrf91w6LjdrXuwTGH1D7lIiQZbSdOJeK1a+y4TUNGPFVjTLHIsC7C/971p?= =?us-ascii?Q?zWYxuEjb6uuGpNlhmhL4efnu6SUn3nTBRWOkT8tVXP6TUtmOROvOuSHqXiQD?= =?us-ascii?Q?9L4N/1btUxDfmFgNGGnBQrYVmMJA1d5U/StXmnHdZP+qnRoZd5koGj9G+L6w?= =?us-ascii?Q?+dxe9gtyBAkeMgtgGSyk7xCIgoWI5keGAyR2qsB0nR8c/wmEFNgAFDbX0wPD?= =?us-ascii?Q?IBOFMOESOZI2RAPn8QRp/Y2zUwguoeq8c8h+LitEgpRw7ncQ10KaS4VcGcHp?= =?us-ascii?Q?D6NTS8MFQEYeGzG51sTqiM5CIB9UcBF+4pV2FDk++l/ERXQpL9tmBIXoCD4/?= =?us-ascii?Q?iDb/k5Ixv9plGZfhyWT/frvsPzCrC4Uoe7AnTzhJpXV1mDgQH6yKfCQEFN2N?= =?us-ascii?Q?WQliwS2QjR2o05JVq2vzVq7+?= X-OriginatorOrg: oracle.com X-MS-Exchange-CrossTenant-Network-Message-Id: af7d948b-bbc9-4d47-8c5e-08d95daab8fd X-MS-Exchange-CrossTenant-AuthSource: BYAPR10MB3208.namprd10.prod.outlook.com X-MS-Exchange-CrossTenant-AuthAs: Internal X-MS-Exchange-CrossTenant-OriginalArrivalTime: 12 Aug 2021 16:03:28.0852 (UTC) X-MS-Exchange-CrossTenant-FromEntityHeader: Hosted X-MS-Exchange-CrossTenant-Id: 4e2c6054-71cb-48f1-bd6c-3a9705aca71b X-MS-Exchange-CrossTenant-MailboxType: HOSTED X-MS-Exchange-CrossTenant-UserPrincipalName: 6KUjVbpbT4g2XxRfDCFoYm2wX5SkzZeyaYCDvOFICHUgVYDqzXJoTrA1OCEiE4paWAcX+gegzm7GoKrhQmFmPLIN5Kffk8h41xaSvZO+oxE= X-MS-Exchange-Transport-CrossTenantHeadersStamped: BY5PR10MB4065 X-Proofpoint-Virus-Version: vendor=nai engine=6200 definitions=10074 signatures=668682 X-Proofpoint-Spam-Details: rule=notspam policy=default score=0 adultscore=0 mlxscore=0 bulkscore=0 spamscore=0 phishscore=0 mlxlogscore=999 suspectscore=0 malwarescore=0 classifier=spam adjust=0 reason=mlx scancount=1 engine=8.12.0-2107140000 definitions=main-2108120104 X-Proofpoint-ORIG-GUID: iFEwu3LgfybJGGLixL4m5S8gHCGUaJOM X-Proofpoint-GUID: iFEwu3LgfybJGGLixL4m5S8gHCGUaJOM X-Spam-Status: No, score=-12.5 required=5.0 tests=BAYES_00, DKIM_SIGNED, DKIM_VALID, DKIM_VALID_AU, DKIM_VALID_EF, GIT_PATCH_0, KAM_SHORT, MSGID_FROM_MTA_HEADER, RCVD_IN_DNSWL_LOW, RCVD_IN_MSPIKE_BL, RCVD_IN_MSPIKE_L3, SPF_HELO_NONE, SPF_NONE, 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: gcc-patches@gcc.gnu.org X-Mailman-Version: 2.1.29 Precedence: list List-Id: Gcc-patches mailing list List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Thu, 12 Aug 2021 16:03:44 -0000 This patch resolves the failure of powerpc64 long double complex divide in native ibm long double format after the patch "Practical improvement to libgcc complex divide". The new code uses the following macros which are intended to be mapped to appropriate values according to the underlying hardware representation. See https://gcc.gnu.org/bugzilla/show_bug.cgi?id=101104 RBIG a value near the maximum representation RMIN a value near the minimum representation (but not in the subnormal range) RMIN2 a value moderately less than 1 RMINSCAL the inverse of RMIN2 RMAX2 RBIG * RMIN2 - a value to limit scaling to not overflow When "long double" values were not using the IEEE 128-bit format but the traditional IBM 128-bit, the previous code used the LDBL values which caused overflow for RMINSCAL. The new code uses the DBL values. RBIG LDBL_MAX = 0x1.fffffffffffff800p+1022 DBL_MAX = 0x1.fffffffffffff000p+1022 RMIN LDBL_MIN = 0x1.0000000000000000p-969 RMIN DBL_MIN = 0x1.0000000000000000p-1022 RMIN2 LDBL_EPSILON = 0x0.0000000000001000p-1022 = 0x1.0p-1074 RMIN2 DBL_EPSILON = 0x1.0000000000000000p-52 RMINSCAL 1/LDBL_EPSILON = inf (1.0p+1074 does not fit in IBM 128-bit). 1/DBL_EPSILON = 0x1.0000000000000000p+52 RMAX2 = RBIG * RMIN2 = 0x1.fffffffffffff800p-52 RBIG * RMIN2 = 0x1.fffffffffffff000p+970 The MAX and MIN values have only modest changes since the maximum and minimum values are about the same as for double precision. The EPSILON field is considerably different. Due to how very small values can be represented in the lower 64 bits of the IBM 128-bit floating point, EPSILON is extremely small, so far beyond the desired value that inversion of the value overflows and even without the overflow, the RMAX2 is so small as to eliminate most usage of the test. Instead of just replacing the use of KF_EPSILON with DF_EPSILON, we replace all uses of KF_* with DF_*. Since the exponent fields are essentially the same, we gain the positive benefits from the new formula while avoiding all under/overflow issues in the #defines. The change has been tested on gcc135.fsffrance.org and gains the expected improvements in accuracy for long double complex divide. libgcc/ PR target/101104 * config/rs6000/_divkc3.c (RBIG, RMIN, RMIN2, RMINSCAL, RMAX2): Use more correct values for native IBM 128-bit. --- libgcc/config/rs6000/_divkc3.c | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/libgcc/config/rs6000/_divkc3.c b/libgcc/config/rs6000/_divkc3.c index a1d29d2..2b229c8 100644 --- a/libgcc/config/rs6000/_divkc3.c +++ b/libgcc/config/rs6000/_divkc3.c @@ -38,10 +38,10 @@ see the files COPYING3 and COPYING.RUNTIME respectively. If not, see #endif #ifndef __LONG_DOUBLE_IEEE128__ -#define RBIG (__LIBGCC_KF_MAX__ / 2) -#define RMIN (__LIBGCC_KF_MIN__) -#define RMIN2 (__LIBGCC_KF_EPSILON__) -#define RMINSCAL (1 / __LIBGCC_KF_EPSILON__) +#define RBIG (__LIBGCC_DF_MAX__ / 2) +#define RMIN (__LIBGCC_DF_MIN__) +#define RMIN2 (__LIBGCC_DF_EPSILON__) +#define RMINSCAL (1 / __LIBGCC_DF_EPSILON__) #define RMAX2 (RBIG * RMIN2) #else #define RBIG (__LIBGCC_TF_MAX__ / 2) -- 1.8.3.1