From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: by sourceware.org (Postfix, from userid 48) id A11613858C53; Sun, 19 Mar 2023 19:58:48 +0000 (GMT) DKIM-Filter: OpenDKIM Filter v2.11.0 sourceware.org A11613858C53 DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gcc.gnu.org; s=default; t=1679255928; bh=Hyvvg2iaavYB/ewR0w9NGaFwksgh2wbMFxSp2a7o6eM=; h=From:To:Subject:Date:In-Reply-To:References:From; b=jLRKxBwP8NIPDsUgXT9kBQyfrIew+KJn9WYElvJaSVaAJA5UjLTi8vEO6h3OE3Ee2 2RSP5Ebi2DATtucAQZAoAvLIonoCV1disgHbOiUlwEgurmmMtIK1KGBrne63aQuwK4 GfqlQw0P1FuDxagCf9udFULe3r7AYQAFYggrJcYc= From: "anlauf at gcc dot gnu.org" To: gcc-bugs@gcc.gnu.org Subject: [Bug fortran/109186] nearest(huge(x),-1.0_kind(x)) half of correct value Date: Sun, 19 Mar 2023 19:58:48 +0000 X-Bugzilla-Reason: CC X-Bugzilla-Type: changed X-Bugzilla-Watch-Reason: None X-Bugzilla-Product: gcc X-Bugzilla-Component: fortran X-Bugzilla-Version: 12.0 X-Bugzilla-Keywords: wrong-code X-Bugzilla-Severity: normal X-Bugzilla-Who: anlauf at gcc dot gnu.org X-Bugzilla-Status: NEW X-Bugzilla-Resolution: X-Bugzilla-Priority: P3 X-Bugzilla-Assigned-To: unassigned at gcc dot gnu.org X-Bugzilla-Target-Milestone: --- X-Bugzilla-Flags: X-Bugzilla-Changed-Fields: Message-ID: In-Reply-To: References: Content-Type: text/plain; charset="UTF-8" Content-Transfer-Encoding: quoted-printable X-Bugzilla-URL: http://gcc.gnu.org/bugzilla/ Auto-Submitted: auto-generated MIME-Version: 1.0 List-Id: https://gcc.gnu.org/bugzilla/show_bug.cgi?id=3D109186 --- Comment #4 from anlauf at gcc dot gnu.org --- The following patch fixes this: diff --git a/gcc/fortran/simplify.cc b/gcc/fortran/simplify.cc index 20ea38e0007..35ae637a483 100644 --- a/gcc/fortran/simplify.cc +++ b/gcc/fortran/simplify.cc @@ -6114,7 +6112,7 @@ gfc_simplify_nearest (gfc_expr *x, gfc_expr *s) kind =3D gfc_validate_kind (BT_REAL, x->ts.kind, 0); mpfr_set_emin ((mpfr_exp_t) gfc_real_kinds[kind].min_exponent - mpfr_get_prec(result->value.real) + 1); - mpfr_set_emax ((mpfr_exp_t) gfc_real_kinds[kind].max_exponent - 1); + mpfr_set_emax ((mpfr_exp_t) gfc_real_kinds[kind].max_exponent); mpfr_check_range (result->value.real, 0, MPFR_RNDU); if (mpfr_sgn (s->value.real) > 0) The above code was apparently introduced with r0-84566-gb6f63e898498e6=