Home | History | Annotate | Download | only in MatrixFunctions

Lines Matching refs:blockIndex

393     for (Index blockIndex = 0; blockIndex < m_clusterSize.rows() - diagIndex; blockIndex++) {
394       // compute (blockIndex, blockIndex+diagIndex) block
395       DynMatrixType A = block(m_T, blockIndex, blockIndex);
396       DynMatrixType B = -block(m_T, blockIndex+diagIndex, blockIndex+diagIndex);
397       DynMatrixType C = block(m_fT, blockIndex, blockIndex) * block(m_T, blockIndex, blockIndex+diagIndex);
398       C -= block(m_T, blockIndex, blockIndex+diagIndex) * block(m_fT, blockIndex+diagIndex, blockIndex+diagIndex);
399       for (Index k = blockIndex + 1; k < blockIndex + diagIndex; k++) {
400 	C += block(m_fT, blockIndex, k) * block(m_T, k, blockIndex+diagIndex);
401 	C -= block(m_T, blockIndex, k) * block(m_fT, k, blockIndex+diagIndex);
403       block(m_fT, blockIndex, blockIndex+diagIndex) = solveTriangularSylvester(A, B, C);