A collection of software packages for the solution and analysis of large-scale, parallel simulations have been developed at Sandia National Labs. With these tools, it is possible for the sophisticated user to rapidly build a world-class modeling capability, whether it be for partial differential equations or molecular theories. The open-source Trilinos framework includes tools for creating the application, including an Automatic Differentiation capability (no more deriving and debugging Jacobians!) and distributed-memory parallel data structures. Even desktop computers now have 4-8 cores that can be efficiently utilized with only modest effort. The solution capabilities in Trilinos include iterative and direct linear solvers, algebraic and multi-level preconditioners, and a wide range of nonlinear solver algorithms. Since the purpose of computational work is not to run a single simulation, but to perform analysis and design, a considerable effort has been put into the high-level analysis capabilities, including eigensolvers for stability analysis, continuation and bifurcation algorithms, transient and periodic solution tracking, and efficient optimization algorithms.

In this talk, I will show examples of these numerous packages working together to perform analysis of some complex PDE systems, including a multi-phase reacting flow application and 4D stability problems.