Lines Matching refs:LU
87 inline static void LU_factor(gene_matrix & LU, Pivot_Vector & pivot, int N)
90 ASSERT( LU.rows()==LU.cols() ) ;
100 if( abs( LU( i, j ) )>=big ) big = abs( LU( i, j ) ) ;
111 theSum = LU( i, j ) ;
112 theSum -= matrix_matrix_product_sliced(LU, i, 0, i-1, LU, 0, j) ;
113 // theSum -= sum( LU( i, Range( fromStart, i-1 ) )*LU( Range( fromStart, i-1 ), j ) ) ;
114 LU( i, j ) = theSum ;
120 theSum = LU( i, j ) ;
121 theSum -= matrix_matrix_product_sliced(LU, i, 0, j-1, LU, 0, j) ;
122 // theSum -= sum( LU( i, Range( fromStart, j-1 ) )*LU( Range( fromStart, j-1 ), j ) ) ;
123 LU( i, j ) = theSum ;
133 dum = LU( index_max, k ) ;
134 LU( index_max, k ) = LU( j, k ) ;
135 LU( j, k ) = dum ;
140 if ( LU( j, j )==0. ) LU( j, j ) = 1.e-20 ;
143 dum = 1./LU( j, j ) ;
144 for( int i=j+1; i<N; i++ ) LU( i, j ) *= dum ;
150 inline static void LU_solve(const gene_matrix & LU, const Pivot_Vector pivot, gene_vector &B, gene_vector X, int N)
155 ASSERT( LU.rows()==LU.cols() ) ;
167 theSum -= matrix_vector_product_sliced(LU, X, i, ii-1, i-1) ;
168 // theSum -= sum( LU( i, Range( ii-1, i-1 ) )*X( Range( ii-1, i-1 ) ) ) ;
169 // theSum -= sum( LU( i, Range( ii-1, i-1 ) )*B( Range( ii-1, i-1 ) ) ) ;
180 theSum -= matrix_vector_product_sliced(LU, X, i, i+1, N) ;
181 // theSum -= sum( LU( i, Range( i+1, toEnd ) )*X( Range( i+1, toEnd ) ) ) ;
182 // theSum -= sum( LU( i, Range( i+1, toEnd ) )*B( Range( i+1, toEnd ) ) ) ;
184 X( i ) = theSum/LU( i, i ) ;
185 // B( i ) = theSum/LU( i, i ) ;
