• Korean Management Science Review
• /
• v.15 no.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.

• Journal of Korean Institute of Industrial Engineers
• /
• v.25 no.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.

• Journal of the Korean Operations Research and Management Science Society
• /
• v.28 no.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.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.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.

• Journal of the military operations research society of Korea
• /
• v.24 no.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.

• Bulletin of the Korean Mathematical Society
• /
• v.39 no.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.

• Journal of applied mathematics & informatics
• /
• v.11 no.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.

• Korean Management Science Review
• /
• v.14 no.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.

• Proceedings of the Korean Operations and Management Science Society Conference
• /
• 1996.10a
• /
• 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.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.7 no.3
• /
• pp.437-447
• /
• 2003

We propose the methods to design the finite impulse response (FIR) and the infinite impulse response (IIR) lattice filters using Schur algorithm through the spectral factorization of the covariance matrix by circulant matrix factorization (CMF). Circulant matrix factorization is also very powerful tool used fur spectral factorization of the covariance polynomial in matrix domain to obtain the minimum phase polynomial without the polynomial root finding problem. Schur algorithm is the method for a fast Cholesky factorization of Toeplitz matrix, which easily determines the lattice filter parameters. Examples for the case of the FIR Inter and for the case of the IIR filter are included, and performance of our method check by comparing of our method and another methods (polynomial root finding and cepstral deconvolution).