87         m_S(size, size),  in m_S()  function
105         m_S(A.rows(),A.cols()),
141         return m_S;  in matrixS()
191       MatrixType m_S, m_T, m_Q, m_Z;
220       const Index dim = m_S.cols();  in hessenbergTriangular()
228       m_S.applyOnTheLeft(m_Q.adjoint());  in hessenbergTriangular()
237           if(m_S.coeff(i,j) != 0)  in hessenbergTriangular()
239             G.makeGivens(m_S.coeff(i-1,j), m_S.coeff(i,j), &m_S.coeffRef(i-1, j));  in hessenbergTriangular()
240             m_S.coeffRef(i,j) = Scalar(0.0);  in hessenbergTriangular()
241             m_S.rightCols(dim-j-1).applyOnTheLeft(i-1,i,G.adjoint());  in hessenbergTriangular()
252             m_S.applyOnTheRight(i,i-1,G);  in hessenbergTriangular()
266       const Index size = m_S.cols();  in computeNorms()
271         m_normOfS += m_S.col(j).segment(0, (std::min)(size,j+2)).cwiseAbs().sum();  in computeNorms()
285         Scalar s = abs(m_S.coeff(res-1,res-1)) + abs(m_S.coeff(res,res));  in findSmallSubdiagEntry()
288         if (abs(m_S.coeff(res,res-1)) < NumTraits<Scalar>::epsilon() * s)  in findSmallSubdiagEntry()
315       const Index dim=m_S.cols();  in splitOffTwoRows()
316       if (abs(m_S.coeff(i+1,i))==Scalar(0))  in splitOffTwoRows()
323           template solve<OnTheRight>(m_S.template block<2,2>(i,i));  in splitOffTwoRows()
336           m_S.rightCols(dim-i).applyOnTheLeft(i,i+1,G.adjoint());  in splitOffTwoRows()
343           m_S.topRows(i+2).applyOnTheRight(i+1,i,G);  in splitOffTwoRows()
349           m_S.coeffRef(i+1,i) = Scalar(0.0);  in splitOffTwoRows()
364       const Index dim = m_S.cols();  in pushDownZero()
370         m_S.rightCols(dim-firstColS).applyOnTheLeft(zz,zz+1,G.adjoint());  in pushDownZero()
379           G.makeGivens(m_S.coeff(zz+1, zz), m_S.coeff(zz+1,zz-1));  in pushDownZero()
380           m_S.topRows(zz+2).applyOnTheRight(zz, zz-1,G);  in pushDownZero()
382           m_S.coeffRef(zz+1,zz-1) = Scalar(0.0);  in pushDownZero()
389       G.makeGivens(m_S.coeff(l,l), m_S.coeff(l,l-1));  in pushDownZero()
390       m_S.applyOnTheRight(l,l-1,G);  in pushDownZero()
392       m_S.coeffRef(l,l-1)=Scalar(0.0);  in pushDownZero()
403       const Index dim = m_S.cols();  in step()
411           a11=m_S.coeff(f+0,f+0), a12=m_S.coeff(f+0,f+1),  in step()
412           a21=m_S.coeff(f+1,f+0), a22=m_S.coeff(f+1,f+1), a32=m_S.coeff(f+2,f+1),  in step()
416           a87=m_S.coeff(l-1,l-2),  in step()
417           a98=m_S.coeff(l-0,l-1),  in step()
432 …x = m_S.coeff(f,f)/m_T.coeff(f,f)-m_S.coeff(l,l)/m_T.coeff(l,l) + m_S.coeff(l,l-1)*m_T.coeff(l-1,l…  in step()
434         y = m_S.coeff(f+1,f)/m_T.coeff(f,f);  in step()
453           a11 = m_S.coeff(f,f),     a12 = m_S.coeff(f,f+1),  in step()
454           a21 = m_S.coeff(f+1,f),   a22 = m_S.coeff(f+1,f+1),  in step()
455                                     a32 = m_S.coeff(f+2,f+1),  in step()
457           a88 = m_S.coeff(l-1,l-1), a89 = m_S.coeff(l-1,l),  in step()
458           a98 = m_S.coeff(l,l-1),   a99 = m_S.coeff(l,l),  in step()
486 …m_S.template middleRows<3>(k).rightCols(dim-fc).applyHouseholderOnTheLeft(essential2, tau, m_works…  in step()
491           m_S.coeffRef(k+2,k-1) = m_S.coeffRef(k+1,k-1) = Scalar(0.0);  in step()
501           tmp = m_S.template middleCols<2>(k).topRows(lr) * essential2;  in step()
502           tmp += m_S.col(k+2).head(lr);  in step()
503           m_S.col(k+2).head(lr) -= tau*tmp;  in step()
504           m_S.template middleCols<2>(k).topRows(lr) -= (tau*tmp) * essential2.adjoint();  in step()
524         m_S.applyOnTheRight(k+1,k,G);  in step()
532         x = m_S.coeff(k+1,k);  in step()
533         y = m_S.coeff(k+2,k);  in step()
535           z = m_S.coeff(k+3,k);  in step()
540       m_S.applyOnTheLeft(l-1,l,G.adjoint());  in step()
544       m_S.coeffRef(l,l-2) = Scalar(0.0);  in step()
548       m_S.applyOnTheRight(l,l-1,G);  in step()
568       m_S = A_in; m_T = B_in;  in compute()
585         if (f>0) m_S.coeffRef(f,f-1) = Scalar(0.0);  in compute()