From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: by sourceware.org (Postfix, from userid 48) id 64D1F3870874; Tue, 16 Feb 2021 11:47:43 +0000 (GMT) DKIM-Filter: OpenDKIM Filter v2.11.0 sourceware.org 64D1F3870874 From: "redi at gcc dot gnu.org" To: gcc-bugs@gcc.gnu.org Subject: [Bug libstdc++/99117] [9/10/11 Regression] cannot accumulate std::valarray Date: Tue, 16 Feb 2021 11:47:43 +0000 X-Bugzilla-Reason: CC X-Bugzilla-Type: changed X-Bugzilla-Watch-Reason: None X-Bugzilla-Product: gcc X-Bugzilla-Component: libstdc++ X-Bugzilla-Version: 10.1.0 X-Bugzilla-Keywords: needs-reduction, wrong-code X-Bugzilla-Severity: normal X-Bugzilla-Who: redi at gcc dot gnu.org X-Bugzilla-Status: NEW X-Bugzilla-Resolution: X-Bugzilla-Priority: P2 X-Bugzilla-Assigned-To: unassigned at gcc dot gnu.org X-Bugzilla-Target-Milestone: 9.4 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 X-BeenThere: gcc-bugs@gcc.gnu.org X-Mailman-Version: 2.1.29 Precedence: list List-Id: Gcc-bugs mailing list List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 16 Feb 2021 11:47:43 -0000 https://gcc.gnu.org/bugzilla/show_bug.cgi?id=3D99117 --- Comment #7 from Jonathan Wakely --- (In reply to Jakub Jelinek from comment #5) > int* p =3D sum._M_data; > int* e1 =3D sum._M_data; >=20 > If p and e1 aren't __restrict__ too, shouldn't that be fine? p (called __p below) doesn't use __restrict__: template void __valarray_copy(const _Expr<_Dom, _Tp>& __e, size_t __n, _Array<_Tp> __= a) { _Tp* __p (__a._M_data); for (size_t __i =3D 0; __i < __n; ++__i, ++__p) *__p =3D __e[__i]; } Here __e is the expression template which has two const valarray& memb= ers, so maybe more accurately it's: struct valarray { int* __restrict__ _M_data; size_t _M_size; }; valarray sum{ new int[2]{1,1}, 2 }; valarray rhs{ new int[2]{2,2}, 2 }; int* p =3D sum._M_data; const valarray& e1 =3D sum; const valarray& e2 =3D rhs; for (size_t i =3D 0; i < sum._M_size; ++i, ++p) *p =3D e1._M_data[i] + e2._M_data[i]; So e1._M_data is marked with __restrict__.=