Effective condition number for finite difference method

被引:49
|
作者
Li, Zi-Cai [1 ]
Chien, Cheng-Sheng
Huang, Hung-Tsai
机构
[1] Natl Sun Yat Sen Univ, Dept Appl Math, Kaohsiung 804, Taiwan
[2] Natl Sun Yat Sen Univ, Dept Comp Sci & Engn, Kaohsiung 804, Taiwan
[3] Natl Ctr Theoret Sci, Hsinchu, Taiwan
[4] Natl Chung Hsing Univ, Dept Appl Math, Taichung 402, Taiwan
关键词
condition number; effective condition number; finite difference method; Poisson's equation;
D O I
10.1016/j.cam.2005.11.037
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
For solving the linear algebraic equations Ax = b with the symmetric and positive definite matrix A, from elliptic equations, the traditional condition number in the 2-norm is defined by Cond. = lambda(1)/lambda(n) where lambda(1) and lambda(n) are the maximal and minimal eigenvalues of the matrix A, respectively. The condition number is used to provide the bounds of the relative errors from the perturbation of both A and b. Such a Cond. can only be reached by the worst situation of all rounding errors and all b. For the given b, the true relative errors may be smaller, or even much smaller than the Cond., which is called the effective condition number in Chan and Foulser [Effectively well-conditioned linear systems, SIAM J. Sci. Statist. Comput. 9 (1988) 963-969] and Christiansen and Hansen [The effective condition number applied to error analysis of certain boundary collocation methods, J. Comput. Appl. Math. 54(1) (1994) 15-36]. In this paper, we propose the new computational formulas for effective condition number Cond_eff, and define the new simplified effective condition number Cond-E. For the latter, we only need the eigenvector corresponding to the minimal eigenvalue of A, which can be easily obtained by the inverse power method. In this paper, we also apply the effective condition number for the finite difference method for Poisson's equation. The difference grids are not supposed to be quasiuniform. Under a non-orthogonality assumption, the effective condition number is proven to be O(1) for the homogeneous boundary conditions. Such a result is extraordinary, compared with the traditional Cond. = O(h(min)(-2)), where h(min) is the minimal meshspacing of the difference grids used. For the non-homogeneous Neumann and Dirichlet boundary conditions, the effective condition number is proven to be O(h(-1/2)) and O(h(-1/2)h(min)(-1)), respectively, where It is the maximal meshspacing of the difference grids. Numerical experiments are carried out to verify the analysis made. (c) 2006 Elsevier B.V. All rights reserved.
引用
收藏
页码:208 / 235
页数:28
相关论文
共 50 条
  • [1] Effective Condition Number of Finite Difference Method for Poisson's Equation Involving Boundary Singularities
    Li, Zi-Cai
    Huang, Hung-Tsai
    Huang, Chien-Sen
    Lu, Tzon-Tzer
    Fang, Qing
    NUMERICAL FUNCTIONAL ANALYSIS AND OPTIMIZATION, 2011, 32 (06) : 659 - 681
  • [2] Effective condition number for the finite element method using local mesh refinements
    Li, Zi-Cai
    Huang, Hung-Tsai
    APPLIED NUMERICAL MATHEMATICS, 2009, 59 (08) : 1779 - 1795
  • [3] Effective Condition Number for Boundary Knot Method
    Wang, F. Z.
    Ling, L.
    Chen, W.
    CMC-COMPUTERS MATERIALS & CONTINUA, 2009, 12 (01): : 57 - 70
  • [4] Condition number analysis and preconditioning of the finite cell method
    de Prenter, F.
    Verhoosel, C. V.
    van Zwieten, G. J.
    van Brummelen, E. H.
    COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2017, 316 : 297 - 327
  • [5] Effective condition number of the Hermite finite element methods for biharmonic equations
    Li, Zi Cai
    Huang, Hung-Tsai
    Huang, Jin
    APPLIED NUMERICAL MATHEMATICS, 2008, 58 (09) : 1291 - 1308
  • [6] Finite Difference Method for Hyperbolic Equations with the Nonlocal Integral Condition
    Ashyralyev, Allaberen
    Aggez, Necmettin
    DISCRETE DYNAMICS IN NATURE AND SOCIETY, 2011, 2011
  • [7] Finite Difference Method for The Reverse Parabolic Problem with Neumann Condition
    Ashyralyyev, Charyyar
    Dural, Ayfer
    Sozen, Yasar
    FIRST INTERNATIONAL CONFERENCE ON ANALYSIS AND APPLIED MATHEMATICS (ICAAM 2012), 2012, 1470 : 102 - 105
  • [8] A FINITE DIFFERENCE METHOD FOR ELLIPTIC PROBLEMS WITH IMPLICIT JUMP CONDITION
    Cao, Fujun
    Yuan, Dongfang
    Sheng, Zhiqiang
    Yuan, Guangwei
    He, Limin
    INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, 2022, 19 (04) : 439 - 457
  • [9] A Rigorous Condition Number Estimate of an Immersed Finite Element Method
    Saihua Wang
    Feng Wang
    Xuejun Xu
    Journal of Scientific Computing, 2020, 83
  • [10] On the condition number of the finite element method for the Laplace–Beltrami operator
    Marcial Nguemfouo
    Michaël Ndjinga
    Journal of Elliptic and Parabolic Equations, 2024, 10 : 59 - 86