• 제목/요약/키워드: Cholesky factorization

검색결과 32건 처리시간 0.028초

내부점방법을 위한 초마디 열촐레스키 분해의 실험적 고찰 (Experimental Study on Supernodal Column Choleksy Factorization in Interior-Point Methods)

  • 설동렬;정호원;박순달
    • 경영과학
    • /
    • 제15권1호
    • /
    • pp.87-95
    • /
    • 1998
  • The computational speed of interior point method depends on the speed of Cholesky factorization. The supernodal column Cholesky factorization is a fast method that performs Cholesky factorization of sparse matrices with exploiting computer's characteristics. Three steps are necessary to perform the supernodal column Cholesky factorization : symbolic factorization, creation of the elimination tree, ordering by a post-order of the elimination tree and creation of supernodes. We study performing sequences of these three steps and efficient implementation of them.

  • PDF

내부점 방법에서 촐레스키 분해의 수치적 안정성 (Numerical Stability of Cholesky Factorization in Interior Point Methods for Linear Programming)

  • 설동렬;성명기;안재근;박순달
    • 대한산업공학회지
    • /
    • 제25권3호
    • /
    • pp.290-297
    • /
    • 1999
  • In interior point methods for linear programming, we must solve a linear system with a symmetric positive definite matrix at every iteration, and Cholesky factorization is generally used to solve it. Therefore, if Cholesky factorization is not done successfully, many iterations are needed to find the optimal solution or we can not find it. We studied methods for improving the numerical stability of Cholesky factorization and the accuracy of the solution of the linear system.

  • PDF

내부점 방법에서 Augmented System의 촐레스키 분해 (Cholesky Factorization of the Augmented System in Interior Point Methods for Linear Programming)

  • 도승용;성명기;박순달
    • 한국경영과학회지
    • /
    • 제28권1호
    • /
    • pp.51-61
    • /
    • 2003
  • In the normal equations approach in which the ordering and factorization phases are separated, the factorization in the augmented system approach is computed dynamically. This means that in the augmented system the numerical factorization should be performed to obtain the non-zero structure of Cholesky factor L. This causes much time to set up the non-zero structure of Cholesky factor L. So, we present a method which can separate the ordering and numerical factorization in the augmented system. Experimental results show that the proposed method reduces the time for obtaining the non-zero structure of Cholesky factor L.

불완전분해법을 전처리로 하는 공액구배법의 안정화에 대한 연구 (Study on Robustness of Incomplete Cholesky Factorization using Preconditioning for Conjugate Gradient Method)

  • 고진환;이병채
    • 대한기계학회논문집A
    • /
    • 제27권2호
    • /
    • pp.276-284
    • /
    • 2003
  • The preconditioned conjugate gradient method is an efficient iterative solution scheme for large size finite element problems. As preconditioning method, we choose an incomplete Cholesky factorization which has efficiency and easiness in implementation in this paper. The incomplete Cholesky factorization mettled sometimes leads to breakdown of the computational procedure that means pivots in the matrix become minus during factorization. So, it is inevitable that a reduction process fur stabilizing and this process will guarantee robustness of the algorithm at the cost of a little computation. Recently incomplete factorization that enhances robustness through increasing diagonal dominancy instead of reduction process has been developed. This method has better efficiency for the problem that has rotational degree of freedom but is sensitive to parameters and the breakdown can be occurred occasionally. Therefore, this paper presents new method that guarantees robustness for this method. Numerical experiment shows that the present method guarantees robustness without further efficiency loss.

네트워크 문제에서 내부점 방법의 활용 (내부점 선형계획법에서 효율적인 공액경사법) (Interior Point Methods for Network Problems (An Efficient Conjugate Gradient Method for Interior Point Methods))

  • 설동렬
    • 한국국방경영분석학회지
    • /
    • 제24권1호
    • /
    • pp.146-156
    • /
    • 1998
  • Cholesky factorization is known to be inefficient to problems with dense column and network problems in interior point methods. We use the conjugate gradient method and preconditioners to improve the convergence rate of the conjugate gradient method. Several preconditioners were applied to LPABO 5.1 and the results were compared with those of CPLEX 3.0. The conjugate gradient method shows to be more efficient than Cholesky factorization to problems with dense columns and network problems. The incomplete Cholesky factorization preconditioner shows to be the most efficient among the preconditioners.

  • PDF

MODIFLED INCOMPLETE CHOLESKY FACTORIZATION PRECONDITIONERS FOR A SYMMETRIC POSITIVE DEFINITE MATRIX

  • Yun, Jae-Heon;Han, Yu-Du
    • 대한수학회보
    • /
    • 제39권3호
    • /
    • pp.495-509
    • /
    • 2002
  • We propose variants of the modified incomplete Cho1esky factorization preconditioner for a symmetric positive definite (SPD) matrix. Spectral properties of these preconditioners are discussed, and then numerical results of the preconditioned CG (PCG) method using these preconditioners are provided to see the effectiveness of the preconditioners.

A MULTILEVEL BLOCK INCOMPLETE CHOLESKY PRECONDITIONER FOR SOLVING NORMAL EQUATIONS IN LINEAR LEAST SQUARES PROBLEMS

  • Jun, Zhang;Tong, Xiao
    • Journal of applied mathematics & informatics
    • /
    • 제11권1_2호
    • /
    • pp.59-80
    • /
    • 2003
  • An incomplete factorization method for preconditioning symmetric positive definite matrices is introduced to solve normal equations. The normal equations are form to solve linear least squares problems. The procedure is based on a block incomplete Cholesky factorization and a multilevel recursive strategy with an approximate Schur complement matrix formed implicitly. A diagonal perturbation strategy is implemented to enhance factorization robustness. The factors obtained are used as a preconditioner for the conjugate gradient method. Numerical experiments are used to show the robustness and efficiency of this preconditioning technique, and to compare it with two other preconditioners.

내부점 선형계획법의 밀집열 분할에 대하여 (On dence column splitting in interial point methods of linear programming)

  • 설동렬;박순달;정호원
    • 경영과학
    • /
    • 제14권2호
    • /
    • pp.69-79
    • /
    • 1997
  • The computational speed of interior point method of linear programming depends on the speed of Cholesky factorization. If the coefficient matrix A has dense columns then the matrix A.THETA. $A^{T}$ becomes a dense matrix. This causes Cholesky factorization to be slow. We study an efficient implementation method of the dense column splitting among dense column resolving technique and analyze the relation between dense column splitting and order methods to improve the sparsity of Cholesky factoror.

  • PDF

내부점 선형계획법의 쌍대문제 전환에 대하여 (On dual transformation in the interior point method of linear programming)

  • 설동렬;박순달;정호원
    • 한국경영과학회:학술대회논문집
    • /
    • 한국경영과학회 1996년도 추계학술대회발표논문집; 고려대학교, 서울; 26 Oct. 1996
    • /
    • pp.289-292
    • /
    • 1996
  • In Cholesky factorization of the interior point method, dense columns of A matrix make dense Cholesky factor L regardless of sparsity of A matrix. We introduce a method to transform a primal problem to a dual problem in order to preserve the sparsity.

  • PDF

FIR/IIR Lattice 필터의 설계를 위한 Circulant Matrix Factorization을 사용한 Spectral Factorization에 관한 연구 (Study of Spectral Factorization using Circulant Matrix Factorization to Design the FIR/IIR Lattice Filters)

  • 김상태;박종원
    • 한국정보통신학회논문지
    • /
    • 제7권3호
    • /
    • pp.437-447
    • /
    • 2003
  • Circulant Matrix Factorization (CMF)는 covariance 행렬의 spectral factorization된 결과를 얻을 수 있다. 우리는 얻어진 결과를 가지고 일반적으로 잘 알려진 방법인 Schur algorithm을 이용하여 finite impulse response(FIR)와 infinite impulse response (IIR) lattice 필터를 설계하는 방법을 제안하였다. CMF는 기존에 많이 사용되는 root finding을 사용하지 않고 covariance polynomial로부터 minimum phase 특성을 가지는 polynomial을 얻는데 유용한 방법이다. 그리고 Schur algorithm은 toeplitz matrix를 빠르게 Cholesky factorization하기 위한 방법으로 이 방법을 이용하면 FIR/IIR lattice 필터의 계수를 쉽게 찾아낼 수 있다. 본 논문에서는 이러한 방법들을 이용하여 FIR과 IIR lattice 필터의 설계의 계산적인 예제를 제시했으며, 제안된 방법과 다른 기존에 제시되었던 방법 (polynomial root finding과 cepstral deconvolution)들과 성능을 비교 평가하였다.