perform the decomposition in place on h by zeroing out the sub sub diagonal elements in the first n - 2 columns of the matrix . concode_field_sep RectangularMatrix u concode_elem_sep RectangularMatrix h concode_field_sep RectangularMatrix getH concode_elem_sep RectangularMatrix getU
void function ( ) { for ( int loc0 = 0 ; loc0 < h . n ( ) - 2 ; loc0 ++ ) { Vector loc1 = h . getColumn ( loc0 ) ; Vector loc2 = loc1 . get ( loc0 + 1 , loc1 . size ( ) ) ; HouseholderReflection loc3 = new HouseholderReflection ( loc2 ) ; loc3 . applyLeft ( h , loc0 + 1 , h . m ( ) , loc0 , h . n ( ) ) ; loc3 . applyRight ( h , 0 , h . m ( ) , loc0 + 1 , h . n ( ) ) ; loc3 . applyRight ( u , 0 , u . m ( ) , loc0 + 1 , u . n ( ) ) ; } }