FiniteDifferenceMethods幻灯片.pptVIP

Copyright ? The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Copyright ? The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Parallel Programmingin C with MPI and OpenMP Michael J. Quinn Chapter 13 Finite Difference Methods Outline Ordinary and partial differential equations Finite difference methods Vibrating string problem Steady state heat distribution problem Ordinary and Partial Differential Equations Ordinary differential equation: equation containing derivatives of a function of one variable Partial differential equation: equation containing derivatives of a function of two or more variables Examples of Phenomena Modeled by PDEs Air flow over an aircraft wing Blood circulation in human body Water circulation in an ocean Bridge deformations as its carries traffic Evolution of a thunderstorm Oscillations of a skyscraper hit by earthquake Strength of a toy Model of Sea Surface Temperaturein Atlantic Ocean Courtesy MICOM group at the Rosenstiel School of Marine and Atmospheric Science, University of Miami Solving PDEs Finite element method Finite difference method (our focus) Converts PDE into matrix equation Result is usually a sparse matrix Matrix-based algorithms represent matrices explicitly Matrix-free algorithms represent matrix values implicitly (our focus) Linear Second-order PDEs Linear second-order PDEs are of the form where A - H are functions of x and y only Elliptic PDEs: B2 - AC 0 Parabolic PDEs: B2 - AC = 0 Hyperbolic PDEs: B2 - AC 0 Difference Quotients Formulas for 1st, 2d Derivatives Heart of Sequential C Program u[j+1][i] = 2.0*(1.0-L)*u[j][i] + L*(u[j][i+1] + u[j][i-1]) - u[j-1][i]; Parallel Program Design Associate primitive task with each element of matrix Examine communication pattern Agglomerate tasks in same column Static number of identical tasks Regular communication pattern Strategy: agglomerate columns, assign one block of columns to each task Result of Agglomerat