269     DenseMatrix refS = DenseMatrix::Zero(rows, rows);  in sparse_product()  local
281     refS = refUp + refLo;  in sparse_product()
282     refS.diagonal() *= 0.5;  in sparse_product()
290     VERIFY_IS_APPROX(refS.adjoint(), refS);  in sparse_product()
292     VERIFY_IS_APPROX(mS, refS);  in sparse_product()
293     VERIFY_IS_APPROX(x=mS*b, refX=refS*b);  in sparse_product()
296     VERIFY_IS_APPROX(x=mUp.template selfadjointView<Upper>()*b, refX=refS*b);  in sparse_product()
297     VERIFY_IS_APPROX(x=mLo.template selfadjointView<Lower>()*b, refX=refS*b);  in sparse_product()
298     VERIFY_IS_APPROX(x=mS.template selfadjointView<Upper|Lower>()*b, refX=refS*b);  in sparse_product()
300     VERIFY_IS_APPROX(x=b * mUp.template selfadjointView<Upper>(),       refX=b*refS);  in sparse_product()
301     VERIFY_IS_APPROX(x=b * mLo.template selfadjointView<Lower>(),       refX=b*refS);  in sparse_product()
302     VERIFY_IS_APPROX(x=b * mS.template selfadjointView<Upper|Lower>(),  refX=b*refS);  in sparse_product()
304     VERIFY_IS_APPROX(x.noalias()+=mUp.template selfadjointView<Upper>()*b, refX+=refS*b);  in sparse_product()
305     VERIFY_IS_APPROX(x.noalias()-=mLo.template selfadjointView<Lower>()*b, refX-=refS*b);  in sparse_product()
306     VERIFY_IS_APPROX(x.noalias()+=mS.template selfadjointView<Upper|Lower>()*b, refX+=refS*b);  in sparse_product()
311                      refX = refLo.template selfadjointView<Lower>()*refS);  in sparse_product()
313                      refX = refS * refLo.template selfadjointView<Lower>());  in sparse_product()
322 …res = mA.template triangularView<Lower>()*mS,   refX = refA.template triangularView<Lower>()*refS);  in sparse_product()
323 …VERIFY_IS_APPROX(mSres = mS * mA.template triangularView<Lower>(), refX = refS * refA.template tri…  in sparse_product()
324 …res = mA.template triangularView<Upper>()*mS,   refX = refA.template triangularView<Upper>()*refS);  in sparse_product()
325 …VERIFY_IS_APPROX(mSres = mS * mA.template triangularView<Upper>(), refX = refS * refA.template tri…  in sparse_product()