140 likes | 257 Views
„Discretization and Rounding Errors: ODE´s and Beyond“. Alan Krautstengl Institute of Finance and Administration Czech Republic Institute of Computer Science of the Czech Academy of Sciences Czech Republic. Errors the Mathematical Modeling and Scientific Computing. Real World Problem.
E N D
„Discretization and Rounding Errors: ODE´s and Beyond“ Alan Krautstengl Institute of Finance and Administration Czech Republic Institute of Computer Science of the Czech Academy of Sciences Czech Republic
Real World Problem • Mathematical model (Integro-differencial equations) • Infinite – dimensional problem • Existence and uniqueness of solution • Numerical model (Galerkin, Ritz, FEM, FD) • Finite – dimensional approximation • Convergence to analytic solution • Linear algebraic model (System of lin. equations) • Forming of the linear system • Matrix computation
Error Analysis • Mathematical model (Errors of the model) • Numerical model (Discretization errors) 3. Linear algebraic model (Computational errors)
Total Error TOTAL ERROR = MODEL ERROR + DISCRETIZATION ERROR + COMPUTATION ERROR • Not independent • Importance of understanding how they may influence each other • Good solution process carries fundamental properties (e.g. positive definiteness) throughout the embedding scheme
ODE´s • Mutual relationship between discretization and computational errors is well-understood • interesting influence
Ingredients • Algebra • Analysis Similar holds for general one-step and multi-step method • Higher dimesions? • General result non-existent • Marginal treatment in leading texts (e. g. Jim Demmel – „Applied Numerical Linear Algebra“)
Algebraic and Analytical Treatment Method of lines Finite differences Finite elements • Analysis get complicated • Intuition and special cases suggest similarities
Other Problems • Other problems • Choice of norm • Algorithm dependence of rounding errors Acknowledgement Author wishes to acknowledge Zdeněk Strakoš for fruitful discussions and mathematical insight contributing to this presentation