Defer set-to-zero in triangular = product so that no aliasing issue occur in the common:
  A.triangularView() = B*A.sefladjointView()*B.adjoint()
case that used to work in 3.2.
(grafted from df279134beee3771e025160762e976680de90503)

--HG--
branch : 3.3
extra : source : df279134beee3771e025160762e976680de90503
3 files changed