Finite Difference Method Error Estimation
Contents |
tour help Tour Start here for a quick overview of the site Help Center Detailed answers to any questions you might have Meta Discuss the workings and policies of this site About Us Learn more about Stack Overflow the company Business Learn more about hiring developers or posting l2 error matlab ads with us Computational Science beta Questions Tags Users Badges Unanswered Ask Question _ Computational Science l2 error norm Stack Exchange is a question and answer site for scientists using computers to solve scientific problems. Join them; it only takes a minute: Sign error norms up Here's how it works: Anybody can ask a question Anybody can answer The best answers are voted up and rise to the top how to measure the error of a finite difference method up vote 0 down vote favorite 2
L2 Norm
Suppose I am solving a pde with a solution known with a finite-difference method. I can represent it as $A_hu_h=f$ for some approximating matrix $A_h$. And I define the discrete norm in which I will analyze the convergence to be $$||e||^2=h\sum e_i^2$$. I have the estimate for the error $$||e||\leq C ||\tau||$$ for $\tau$ local error if first order in space and $C$ a stability bound for $A^{-1}$. Then, I want to see which error gives my program to confirm the theory. I pick a point on the grid and watch what happens as a double the mesh. Assume I see the error gets decreased by 2, so I suspect it is a linear convergence at that point. Thus, I can see only rate of convergence at one point. How can I guarantee that the convergence is uniform "everywhere" on the grid? to confirm the theory results, should not I measure the error at each point and see I have a linear convergence in order to claim that the error goes to zero linearly in a discrete norm defined above? Because I could state results in a similar norms such as discrete $L^{\infty}$ or discrete $L^1$, however, what I measure by the computer is the same: difference between numerical solution and the function value at the point. What confuses me is the fact that before I implement the method, I can state error estimates in a variety of norms, however, how do I relate these theoretical approximations to the error that I can actually measure, since this is just a value at the point on the grid and is independent from the way I do my theoretical analysis? Edit: I should rephrase the question: There are a number of papers, where the convergence is measured for a particular point on the grid and the results are stated to confirm the theoretical estimates that are done in discrete $L^{\infty}$ or $L^{2}$ for the error vector. I
has multiple issues. Please help improve it or discuss these issues on the talk page. (Learn how and when to remove these template messages) This article's lead section may not adequately summarize key points of its contents. Please consider expanding the lead to provide an accessible overview of all important aspects of the article. Please discuss this issue on the article's talk page. (April 2015) This article may be http://scicomp.stackexchange.com/questions/2775/how-to-measure-the-error-of-a-finite-difference-method too technical for most readers to understand. Please help improve this article to make it understandable to non-experts, without removing the technical details. The talk page may contain suggestions. (April 2015) (Learn how and when to remove this template message) (Learn how and when to remove this template message) Differential equations https://en.wikipedia.org/wiki/Finite_difference_method Navier–Stokes differential equations used to simulate airflow around an obstruction. Scope Natural sciences Engineering Astronomy Physics Chemistry Biology Geology Applied mathematics Continuum mechanics Chaos theory Dynamical systems Social sciences Economics Population dynamics Classification Types Ordinary Partial Differential-algebraic Integro-differential Fractional Linear Non-linear By variable type Dependent and independent variables Autonomous Complex Coupled/ Decoupled Exact Homogeneous/ Nonhomogeneous Features Order Operator Notation Relation to processes Difference (discrete analogue) Stochastic Delay Solution General topics Picard–Lindelöf theorem Wronskian Phase portrait Phase space Lyapunov/ Asymptotic/ Exponential stability Rate of convergence Series/ Integral solutions Numerical integration Dirac delta function Solution methods Inspection Method of characteristics Euler Finite difference(Crank–Nicolson) Finite element Finite volume Galerkin Integrating factor Integral transforms Perturbation theory Runge–Kutta Separation of variables Undetermined coefficients People Isaac Newton Leonhard Euler Émile Picard Józef Maria Hoene-Wroński Ernst Lindelöf Rudolf Lipschitz Augustin-Louis Cauchy John Crank Phyllis Nicolson Carl David Tolmé Runge Martin Wilhelm Kutta v t e In mathematics,
of error estimates for finite-difference methodsAuthorsAuthors and affiliationsM. N. SpijkerArticleReceived: 01 January 1971DOI: 10.1007/BF01398460Cite this article as: Spijker, M.N. Numer. Math. (1971) 18: 73. doi:10.1007/BF01398460AbstractIn this paper we study in http://link.springer.com/article/10.1007/BF01398460 an abstract setting the structure of estimates for the global (accumulated) error in semilinear finite-difference methods. We derive error estimates, which are the most refined ones (in a sense specified precisely in this paper) that are possible for the difference methods considered. Applications and (numerical) examples are presented in the following fields: 1. Numerical solution of ordinary as well as partial l2 error differential equations with prescribed initial or boundary values. 2. Accumulation of local round-off error as well as of local discretization error. 3. The problem of fixing which methods out of a given class of finite-difference methods are “most stable”. 4. The construction of finite-difference methods which are convergent but not consistent with respect to a given differential equation.References1.Ceschino, F., Kuntzmann, finite difference method J.: Problèmes différentiels de conditions initiales. Paris: Dunod 1963.Google Scholar2.Forsythe, G. E., Wasow, W. R.: Finite-difference methods for partial differential equations. New York: J. Wiley & Sons 1960.Google Scholar3.Godunov, S. K., Ryabenki, V. S.: Theory of difference schemes. Amsterdam: North-Holland Publishing Company 1964.Google Scholar4.Gragg, W. B., Stetter, H. J.: Generalized multistep predictor-corrector methods. J. Assoc. Comput. Mach.11, 188–209 (1964).Google Scholar5.Henrici, P.: Discrete variable methods in ordinary differential equations. New York: J. Wiley & Sons 1962.Google Scholar6.Hull, T. E., Luxemburg, W. A. J.: Numerical methods and existence theorems for ordinary differential equations. Numer. Math.2, 30–41 (1960).Google Scholar7.Isaacson, E., Keller, H. B.: Analysis of numerical methods. New York: J. Wiley & Sons 1966.Google Scholar8.Lees, M.: Discrete methods for nonlinear two-point boundary value problems. In: Numerical solution of partial differential equations, e. d. J. H. Bramble. New York: Academic Press 1966.Google Scholar9.Metté, A.: Essai de résolution du problème de Goursat par la methode de Runge-Kutta pour une equation aux dérivées partielles du type hyperbolique. Rev. Francaise Informat. Recherche Opérationnelle1, 67–90 (1967).Google Scholar10.Spijker, M. N.: Convergence and stability of step-by-step methods for the