93     real theSum = 0. ;  in LU_factor()  local
111 	theSum = LU( i, j ) ;  in LU_factor()
112 	theSum -= matrix_matrix_product_sliced(LU, i, 0, i-1, LU, 0, j) ;  in LU_factor()
114 	LU( i, j ) = theSum ;  in LU_factor()
120 	theSum = LU( i, j ) ;  in LU_factor()
121 	theSum -= matrix_matrix_product_sliced(LU, i, 0, j-1, LU, 0, j) ;  in LU_factor()
123 	LU( i, j ) = theSum ;  in LU_factor()
124 	if( (ImplicitScaling( i )*abs( theSum ))>=big ) {  in LU_factor()
125 	  dum = ImplicitScaling( i )*abs( theSum ) ;  in LU_factor()
159     real theSum = 0. ;  in LU_solve()  local
162       theSum = X( ip ) ;  in LU_solve()
167 	theSum -= matrix_vector_product_sliced(LU, X, i, ii-1, i-1) ;  in LU_solve()
170       } else if( theSum ) {  in LU_solve()
173       X( i ) = theSum ;  in LU_solve()
178       theSum = X( i ) ;  in LU_solve()
180       theSum -= matrix_vector_product_sliced(LU, X, i, i+1, N) ;  in LU_solve()
184       X( i ) = theSum/LU( i, i ) ;  in LU_solve()