620 likes | 635 Views
Performance analysis tools applied to a finite adaptive mesh free boundary seepage parallel algorithm. S. Boeriu 1 and J.C. Bruch, Jr. 2 1 Center for Computational Science and Engineering 2 Department of Mechanical and Environmental Engineering and Department of Mathematics
E N D
Performance analysis tools applied to a finite adaptive mesh free boundary seepage parallel algorithm S. Boeriu1 and J.C. Bruch, Jr.2 1Center for Computational Science and Engineering 2Department of Mechanical and Environmental Engineering and Department of Mathematics University of California, Santa Barbara http://www.engineering.ucsb.edu/~hpscicom
Acknowledgements Thismaterial is based upon work supported by the National Science Foundation under Grant #0086262. This research was supported in part by NSF cooperative agreement ACI-9619020 through computing resources provided by the National Partnership for Advanced Computational Infrastructure at the San Diego Supercomputer Center. http://www.npaci.edu/Horizon/guide_linked/bh_tools_txt.html
Outline of Presentation • Introduction (Physical problem) • Problem formulation • Fixed domain formulation • Numerical algorithm • Test case • Performance tools and considerations a. VAMPIR b. PARAVER • Diagnostic example • Conclusions
Physical problem Figure 1. Seepage through a rectangular dam.
Simplifying assumptions • The soil in the flowfield is homogeneous and isotropic • Capillary and evaporation effects are neglected • The flow obeys Darcy’s Law • Two-dimensional • Steady state
Mathematical formulation • Darcy’s Law: • Potential Function: Velocity Components: Continuity Equation: Irrotationality Condition: Cauchy-Riemann Equations: Laplace’s Equations:
Problem formulation Figure 2. Mathematical formulation of physical problem.
Extension of solution domain The solution domain is extended to the known region Then extend continuously to be defined on by setting
This yields in the sense of distributions where
Fixed domain formulation Figure 3. Fixed domain mathematical formulation.
Numerical Algorithm A minimization problem can be formulated in terms of the functional where ais a bilinear form, continuous, symmetric, positive definite on R and i.e.,
The functional J has one and only one minimum on a closed convex set. The minimum is found using the following algorithm:
Finite Element Error Analysis Adaptive Mesh Finite Element Analysis (FEA) General Equation for FEA:
Error Analysis Error Definition: where is the approximation of the exact solution ; is the calculated of an element (constant); is the shape function; and
Averaging Technique: Error Estimate in an Element:
Error Norm of the Whole Computation Domain: Percentage Error:
Local Mesh Refinement Desired Criteria: Desired Local Error Criteria: Error Ratio: New Element Size:
Performance tools and considerations The parallel program is monitored while it is executed. Monitoring produces performance data that is interpreted in order to reveal areas of poor performance. The program is then altered and the process is repeated until an acceptable level of performance is reached.
VAMPIR(Visualization and Analysis of MPI Resources – 2.0) • VAMPIR 2.0 is a post-mortem trace visualization tool from Pallas GmbH http://www.pallas.com It uses the profile extensions to MPI and permits analysis of the message events where data is transmitted between processors during execution of a parallel program. It has a convenient user-interface and an excellent zooming and filtering. Global displays show all selected processes.
Global Timeline: detailed application execution over time axis • Activity Chart: presents per-process profiling information • Summaric Chart: aggregated profiling information • Communication Statistics: message statistics for each process pair • Global Communication Statistics: collective operations statistics • I/O Statistics: MPI I/O operation statistics • Calling Tree: global dynamic calling tree
PARAVER(Parallel Program Visualization and Analysis Tool) • PARAVER is a flexible parallel program visualization and analysis tool based on an easy-to-use Motif GUI (graphical user interface) PARAVER was developed to respond to the basic need to have a qualitative perception of the application behavior by visual inspection and then to be able to focus on the detailed quantitative analysis of the problems.
Paraver (Parallel Program Visualization and Analysis Tool) • Powerful flexible parallel program visualization tool based on an easy-to-use Motif GUI (graphical user interface) • Developed by : European Center for Parallelism of Barcelona (CEPBA) Universitat Politecnica de Catalunya http://www.cepba.upc.es/
Paraver is designed to visualize and analyze - Communication and load balance - Combining OpenMP and MPI - Hardware performance and counters • Usage - Compile programs with special libraries - Run programs to produce trace files - View and analyze traces - Designed to help in program understanding and optimization
Inefficient programming example • Load imbalance (inefficient memory use) • TLB (translation lookaside buffer) misses