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

검색결과 19건 처리시간 0.027초

BLOCK INCOMPLETE FACTORIZATION PRECONDITIONERS FOR A SYMMETRIC H-MATRIX

  • Yun, Jae-Heon;Kim, Sang-Wook
    • 대한수학회보
    • /
    • 제37권3호
    • /
    • pp.551-568
    • /
    • 2000
  • We propose new parallelizable block incomplete factorization preconditioners for a symmetric block-tridiagonal H-matrix. Theoretical properties of these block preconditioners are compared with those of block incomplete factorization preconditioners for the corresponding comparison matrix. Numerical results of the preconditioned CG(PCG) method using these block preconditioners are compared with those of PCG method using a standard incomplete factorization preconditioner to see the effectiveness of the block incomplete factorization preconditioners.

  • PDF

A VARIANT OF BLOCK INCOMPLETE FACTORIZATION PRECONDITIONERS FOR A SYMMETRIC H-MATRIX

  • Yun, Jae-Heon;Kim, Sang-Wook
    • Journal of applied mathematics & informatics
    • /
    • 제8권3호
    • /
    • pp.705-720
    • /
    • 2001
  • We propose a variant of parallel block incomplete factorization preconditioners for a symmetric block-tridiagonal H-matrix. Theoretical properties of these block preconditioners are compared with those of block incomplete factoriztion preconditioners for the corresponding somparison matrix. Numerical results of the preconditioned CG(PCG) method using these block preconditioners are compared with those of PCG using other types of block incomplete factorization preconditioners. Lastly, parallel computations of the block incomplete factorization preconditioners are carried out on the Cray C90.

불완전분해법을 전처리로 하는 공액구배법의 안정화에 대한 연구 (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.

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.

PARALLEL BLOCK ILU PRECONDITIONERS FOR A BLOCK-TRIDIAGONAL M-MATRIX

  • Yun, Jae-Heon;Kim, Sang-Wook
    • 대한수학회지
    • /
    • 제36권1호
    • /
    • pp.209-227
    • /
    • 1999
  • We propose new parallel block ILU (Incomplete LU) factorization preconditioners for a nonsymmetric block-tridiagonal M-matrix. Theoretial properties of these block preconditioners are studied to see the convergence rate of the preconditioned iterative methods, Lastly, numerical results of the right preconditioned GMRES and BiCGSTAB methods using the block ILU preconditioners are compared with those of these two iterative methods using a standard ILU preconditioner to see the effectiveness of the block ILU preconditioners.

  • PDF

네트워크 문제에서 내부점 방법의 활용 (내부점 선형계획법에서 효율적인 공액경사법) (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

A CLASS OF MULTILEVEL RECURSIVE INCOMPLETE LU PRECONDITIONING TECHNIQUES

  • Zhang, Jun
    • Journal of applied mathematics & informatics
    • /
    • 제8권2호
    • /
    • pp.305-326
    • /
    • 2001
  • We introduce a class of multilevel recursive incomplete LU preconditioning techniques (RILUM) for solving general sparse matrices. This techniques is based on a recursive two by two block incomplete LU factorization on the coefficient martix. The coarse level system is constructed as an (approximate) Schur complement. A dynamic preconditioner is obtained by solving the Schur complement matrix approximately. The novelty of the proposed techniques is to solve the Schur complement matrix by a preconditioned Krylov subspace method. Such a reduction process is repeated to yield a multilevel recursive preconditioner.

HYBRID REORDERING STRATEGIES FOR ILU PRECONDITIONING OF INDEFINITE SPARSE MATRICES

  • Lee Eun-Joo;Zgang Jun
    • Journal of applied mathematics & informatics
    • /
    • 제22권1_2호
    • /
    • pp.307-316
    • /
    • 2006
  • Incomplete LU factorization preconditioning techniques often have difficulty on indefinite sparse matrices. We present hybrid reordering strategies to deal with such matrices, which include new diagonal reorderings that are in conjunction with a symmetric nondecreasing degree algorithm. We first use the diagonal reorderings to efficiently search for entries of single element rows and columns and/or the maximum absolute value to be placed on the diagonal for computing a nonsymmetric permutation. To augment the effectiveness of the diagonal reorderings, a nondecreasing degree algorithm is applied to reduce the amount of fill-in during the ILU factorization. With the reordered matrices, we achieve a noticeable improvement in enhancing the stability of incomplete LU factorizations. Consequently, we reduce the convergence cost of the preconditioned Krylov subspace methods on solving the reordered indefinite matrices.

RECURSIVE TWO-LEVEL ILU PRECONDITIONER FOR NONSYMMETRIC M-MATRICES

  • Guessous, N.;Souhar, O.
    • Journal of applied mathematics & informatics
    • /
    • 제16권1_2호
    • /
    • pp.19-35
    • /
    • 2004
  • We develop in this paper some preconditioners for sparse non-symmetric M-matrices, which combine a recursive two-level block I LU factorization with multigrid method, we compare these preconditioners on matrices arising from discretized convection-diffusion equations using up-wind finite difference schemes and multigrid orderings, some comparison theorems and experiment results are demonstrated.