• Title/Summary/Keyword: preconditioned conjugate gradient scheme

Search Result 12, Processing Time 0.019 seconds

A Deflation-Preconditioned Conjugate Gradient Method for Symmetric Eigenproblems

  • Jang, Ho-Jong
    • Journal of applied mathematics & informatics
    • /
    • v.9 no.1
    • /
    • pp.331-339
    • /
    • 2002
  • A preconditioned conjugate gradient(PCG) scheme with the aid of deflation for computing a few of the smallest eigenvalues arid their corresponding eigenvectors of the large generalized eigenproblems is considered. Topically there are two types of deflation techniques, the deflation with partial shifts and an arthogonal deflation. The efficient way of determining partial shifts is suggested and the deflation-PCG schemes with various partial shifts are investigated. Comparisons of theme schemes are made with orthogonal deflation-PCG, and their asymptotic behaviors with restart operation are also discussed.

AN ACCELERATED DEFLATION TECHNIQUE FOR LARGE SYMMETRIC GENERALIZED EIGENPROBLEMS

  • HYON, YUN-KYONG;JANG, HO-JONG
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • v.3 no.1
    • /
    • pp.99-106
    • /
    • 1999
  • An accelerated optimization technique combined with a stepwise deflation procedure is presented for the efficient evaluation of a few of the smallest eigenvalues and their corresponding eigenvectors of the generalized eigenproblems. The optimization is performed on the Rayleigh quotient of the deflated matrices by the aid of a preconditioned conjugate gradient scheme with the incomplete Cholesky factorization.

  • PDF

NUMERICAL STABILITY OF UPDATE METHOD FOR SYMMETRIC EIGENVALUE PROBLEM

  • Jang Ho-Jong;Lee Sung-Ho
    • Journal of applied mathematics & informatics
    • /
    • v.22 no.1_2
    • /
    • pp.467-474
    • /
    • 2006
  • We present and study the stability and convergence of a deflation-preconditioned conjugate gradient(PCG) scheme for the interior generalized eigenvalue problem $Ax = {\lambda}Bx$, where A and B are large sparse symmetric positive definite matrices. Numerical experiments are also presented to support our theoretical results.

PERTURBATION ANALYSIS OF DEFLATION TECHNIQUE FOR SYMMETRIC EIGENVALUE PROBLEM

  • JANG, HO-JONG
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • v.5 no.2
    • /
    • pp.17-23
    • /
    • 2001
  • The evaluation of a few of the smallest eigenpairs of large symmetric eigenvalue problem is of great interest in many physical and engineering applications. A deflation-preconditioned conjugate gradient(PCG) scheme for a such problem has been shown to be very efficient. In the present paper we provide the numerical stability of a deflation-PCG with partial shifts.

  • PDF

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

  • Ko, Jin-Hwan;Lee, Byung-Chai
    • 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.

An Efficient Adaptive Wavelet-Collocation Method Using Lifted Interpolating Wavelets (수정된 보간 웨이블렛응 이용한 적응 웨이블렛-콜로케이션 기법)

  • Kim, Yun-Yeong;Kim, Jae-Eun
    • Transactions of the Korean Society of Mechanical Engineers A
    • /
    • v.24 no.8 s.179
    • /
    • pp.2100-2107
    • /
    • 2000
  • The wavelet theory is relatively a new development and now acquires popularity and much interest in many areas including mathematics and engineering. This work presents an adaptive wavelet method for a numerical solution of partial differential equations in a collocation sense. Due to the multi-resolution nature of wavelets, an adaptive strategy can be easily realized it is easy to add or delete the wavelet coefficients as resolution levels progress. Typical wavelet-collocation methods use interpolating wavelets having no vanishing moment, but we propose a new wavelet-collocation method on modified interpolating wavelets having 2 vanishing moments. The use of the modified interpolating wavelets obtained by the lifting scheme requires a smaller number of wavelet coefficients as well as a smaller condition number of system matrices. The latter property makes a preconditioned conjugate gradient solver more useful for efficient analysis.

AN EFFICIENT ALGORITHM FOR INCOMPRESSIBLE FREE SURFACE FLOW ON CARTESIAN MESHES (직교격자상에서 효율적인 비압축성 자유표면유동 해법)

  • Go, G.S.;Ahn, H.T.
    • Journal of computational fluids engineering
    • /
    • v.19 no.4
    • /
    • pp.20-28
    • /
    • 2014
  • An efficient solution algorithm for simulating free surface problem is presented. Navier-Stokes equations for variable density incompressible flow are employed as the governing equation on Cartesian meshes. In order to describe the free surface motion efficiently, VOF(Volume Of Fluid) method utilizing THINC(Tangent of Hyperbola for Interface Capturing) scheme is employed. The most time-consuming part of the current free surface flow simulations is the solution step of the linear system, derived by the pressure Poisson equation. To solve a pressure Poisson equation efficiently, the PCG(Preconditioned Conjugate Gradient) method is utilized. This study showed that the proper application of the preconditioner is the key for the efficient solution of the free surface flow when its pressure Poisson equation is solved by the CG method. To demonstrate the efficiency of the current approach, we compared the convergence histories of different algorithms for solving the pressure Poisson equation.

Numerical Analysis of Shallow Water Equation with Fully Implicit Method (음해법을 이용한 천수방정식의 수치해석)

  • Kang, Ju Whan;Park, Sang Hyun;Lee, Kil Seong
    • KSCE Journal of Civil and Environmental Engineering Research
    • /
    • v.13 no.3
    • /
    • pp.119-127
    • /
    • 1993
  • Recently, ADI scheme has been a most common tool for solving shallow water equation numerically. But ADI models of tidal flow is likely to cause so called ADI effect in such a region of the Yellow Sea which shows complex topography and has submarine canyons especially. To overcome this, a finite difference algorithm is developed which adopts fully implicit method and preconditioned conjugate gradient squared method. Applying the algorithm including simulation of intertidal zone to Sae-Man-Keum. velocity fields and flooding/drying phenomena are simulated well in spite of complex topography.

  • PDF

Convergence of Nonlocal Integral Operator in Peridynamics (비국부 적분 연산기로 표현되는 페리다이나믹 방정식의 수렴성)

  • Jo, Gwanghyun;Ha, Youn Doh
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.34 no.3
    • /
    • pp.151-157
    • /
    • 2021
  • This paper is devoted to a convergence study of the nonlocal integral operator in peridynamics. The implicit formulation can be an efficient approach to obtain the static/quasi-static solution of crack propagation problems. Implicit methods require constly large-matrix operations. Therefore, convergence is important for improving computational efficiency. When the radial influence function is utilized in the nonlocal integral equation, the fractional Laplacian integral equation is obtained. It has been mathematically proved that the condition number of the system matrix is affected by the order of the radial influence function and nonlocal horizon size. We formulate the static crack problem with peridynamics and utilize Newton-Raphson methods with a preconditioned conjugate gradient scheme to solve this nonlinear stationary system. The convergence behavior and the computational time for solving the implicit algebraic system have been studied with respect to the order of the radial influence function and nonlocal horizon size.

Diffusion synthetic acceleration with the fine mesh rebalance of the subcell balance method with tetrahedral meshes for SN transport calculations

  • Muhammad, Habib;Hong, Ser Gi
    • Nuclear Engineering and Technology
    • /
    • v.52 no.3
    • /
    • pp.485-498
    • /
    • 2020
  • A diffusion synthetic acceleration (DSA) technique for the SN transport equation discretized with the linear discontinuous expansion method with subcell balance (LDEM-SCB) on unstructured tetrahedral meshes is presented. The LDEM-SCB scheme solves the transport equation with the discrete ordinates method by using the subcell balances and linear discontinuous expansion of the flux. Discretized DSA equations are derived by consistently discretizing the continuous diffusion equation with the LDEM-SCB method, however, the discretized diffusion equations are not fully consistent with the discretized transport equations. In addition, a fine mesh rebalance (FMR) method is devised to accelerate the discretized diffusion equation coupled with the preconditioned conjugate gradient (CG) method. The DSA method is applied to various test problems to show its effectiveness in speeding up the iterative convergence of the transport equation. The results show that the DSA method gives small spectral radii for the tetrahedral meshes having various minimum aspect ratios even in highly scattering dominant mediums for the homogeneous test problems. The numerical tests for the homogeneous and heterogeneous problems show that DSA with FMR (with preconditioned CG) gives significantly higher speedups and robustness than the one with the Gauss-Seidel-like iteration.