324         m_alphas.coeffRef(i) = mS.diagonal().coeff(i);  in compute()
325         m_betas.coeffRef(i)  = mT.diagonal().coeff(i);  in compute()
329           v.coeffRef(i) = Scalar(1.0);  in compute()
331           if(abs(m_betas.coeffRef(i)) >= (std::numeric_limits<RealScalar>::min)())  in compute()
334             const Scalar alpha = real(m_alphas.coeffRef(i));  in compute()
335             const Scalar beta = m_betas.coeffRef(i);  in compute()
350 …coeffRef(j) = -v.segment(st,sz).transpose().cwiseProduct(beta*mS.block(j,st,1,sz) - alpha*mT.block…  in compute()
369         const RealScalar beta = m_betas.coeffRef(i) = m_betas.coeffRef(i+1) = a*b;  in compute()
377         m_alphas.coeffRef(i)   = conj(alpha);  in compute()
378         m_alphas.coeffRef(i+1) = alpha;  in compute()
383           cv.coeffRef(i+1) = Scalar(1.0);  in compute()
385 …        cv.coeffRef(i) = -(static_cast<Scalar>(beta*mS.coeffRef(i,i+1)) - alpha*mT.coeffRef(i,i+1))  in compute()
386                           / (static_cast<Scalar>(beta*mS.coeffRef(i,i))   - alpha*mT.coeffRef(i,i));  in compute()
399 …cv.coeffRef(j) =  cv.segment(st,sz).transpose().cwiseProduct(beta*mS.block(j,st,1,sz) - alpha*mT.b…  in compute()
400 …                           / (alpha*mT.coeffRef(j,j) - static_cast<Scalar>(beta*mS.coeffRef(j,j)));  in compute()