main/sc/source/core/tool/interpr5.cxx
// Same with transposed matrices A and Y.
Same with transposed matrices A and Y.
Same with transposed matrices A and Y .
void lcl_TApplyHouseholderTransformation(ScMatrixRef pMatA, SCSIZE nR,
                                          ScMatrixRef pMatY, SCSIZE nN)
{
    double fDenominator = lcl_TGetColumnSumProduct(pMatA, nR, pMatA, nR, nR, nN);
    double fNumerator = lcl_TGetColumnSumProduct(pMatA, nR, pMatY, 0, nR, nN);
    double fFactor = 2.0 * (fNumerator/fDenominator);
    for (SCSIZE col = nR; col < nN; col++)
        pMatY->PutDouble(
          pMatY->GetDouble(col) - fFactor * pMatA->GetDouble(col,nR), col);
}
void lcl_TApplyHouseholderTransformation ( ScMatrixRef pMatA , SCSIZE nR , ScMatrixRef pMatY , SCSIZE nN ) { double fDenominator = lcl_TGetColumnSumProduct ( pMatA , nR , pMatA , nR , nR , nN ) ; double fNumerator = lcl_TGetColumnSumProduct ( pMatA , nR , pMatY , 0 , nR , nN ) ; double fFactor = 2.0 * ( fNumerator / fDenominator ) ; for ( SCSIZE col = nR ; col < nN ; col ++ ) pMatY -> PutDouble ( pMatY -> GetDouble ( col ) - fFactor * pMatA -> GetDouble ( col , nR ) , col ) ; }
Same with transposed matrices A and Y.
Same with transposed matrices A and Y .
// ScMatrix matrices are zero based, index access (column,row)