Integration of Trilinos Into The Cactus Code Framework

Size: px
Start display at page:

Download "Integration of Trilinos Into The Cactus Code Framework"

Transcription

1 Integration of Trilinos Into The Cactus Code Framework Josh Abadie Research programmer Center for Computation & Technology Louisiana State University

2 Summary Motivation Objectives The Cactus Code Trilinos Process of Integration Further Work Acknowledgments

3 Motivation If in other science we should arrive at certainty without doubt and truth without error it behooves us to place the foundations of knowledge in mathematics Roger Bacon Mathematics can be used to describe accurately amazingly complex and otherwise incomprehensible things Complex math problems can be solved with computers Scientific Computing Programs which solve these problems can be very complex Software reuse is therefore important in Scientific Computing because of this complexity Frameworks, like Cactus and libraries, like Trilinos facilitate reuse of code, save time, leave less room for error, and encourage cooperation in parallel code development

4 Motivation: Parallel Code Development Parallel computing is an attempt to maximize the infinite but seemingly scarce commodity called time. Traditionally, software has been written for serial computation In the simplest sense, parallel computing is the simultaneous use of multiple compute resources to solve a computational problem Source:

5 Motivation: Parallel Programming Complexity Parallel applications usually much more complex than corresponding serial applications, perhaps an order of magnitude Not only do you have multiple instruction streams executing at the same time, but you also have data flowing between them So complex is parallel programming that frameworks like Cactus exist to make the job easier Source:

6 Objectives Cactus is not the answer to all our prayers: other frameworks exist that contain specialized code (solvers, preconditioners ) that we need to use in order to solve complex problems Since Cactus is extensible, we want to leverage the resource that is Trilinos for use in Cactus The integration of the Trilinos library in Cactus should be seamless, the code should have a clear interface that is easy to use

7 Cactus Code Framework Open source application framework designed primarily for scientists and engineers Simple and easy to extend in both C/C++ and Fortran The Flesh Glue for the framework The Thorns - actual computation Modules of Cactus

8 Structure of Cactus Plug-In Thorns (modules) driver input/output extensible APIs ANSI C CCL Interpreter scheduling Core Flesh remote steering Fortran/C/C++ equations of state interpolation SOR solver error handling parameters grid variables make system black holes boundary conditions wave evolvers multigrid coordinates Source: Tom Goodale

9 Standard Thorns HDF5 Toolkit PETSc Toolkit Web Browser Toolkit CFD Toolkit I/O using HDF5 data format Library of Elliptic equation solvers Web steering of computations Solving problems in CFD

10 Trilinos Set of Libraries solvers, preconditioners vector, matrix and sparse matrix classes with parallel capability Developed at Sandia National labs Easy to use, well documented common set of Data structures - Epetra

11 Trilinos Libraries LAPACK Algebraic preconditioner IFPACK Algebraic preconditioner ML Multilevel preconditioner Aztec Linear Solver NOX Collection of nonlinear solvers Anasazi Collection on eigen solvers Komplex Complex Solver Belos Krylov linear solver

12 Process of Integration Create Trilinos Libraries Documentation concise and very helpful Easy to make multiple builds with different configurations Create a Cactus thorn Create a thorn which tells Cactus about Trilinos and where the libraries are located To use Trilinos, just require the thorn, set some environment variables, write your thorn

13 Using Trilinos In Cactus Plug-In Thorns (module) extensible APIs ANSI C parameters Trilinos-Test scheduling Core Flesh Dependency Checking Trilinos error handling make system grid variables

14 Future Work Utility library A set of functions to simplify the process of using Trilinos in Cactus MPI transparency Data translation Epetra Driver Cactus Thorn which stores all data does all memory management with Trilinos Epetra classes Ultimate integration of Trilinos and Cactus Benchmarking

15 Acknowledgements Yaakoub El Khamra supervisor and mentor, Frameworks at CCT Gabrielle Allen Assistant Director of the Center for Computation and Technology Tom Goodale Chief Architect of the Cactus Code Framework Kathryn Traxler Education and Outreach at CCT The Trilinos Project Team Development of such a great library The Computer Science Department and Louisiana State University Assistance with funding

Fourteen years of Cactus Community

Fourteen years of Cactus Community Fourteen years of Cactus Community Frank Löffler Center for Computation and Technology Louisiana State University, Baton Rouge, LA September 6th 2012 Outline Motivation scenario from Astrophysics Cactus

More information

Cactus Tutorial. Introduction to Cactus. Yaakoub El Khamra. Cactus Developer, Frameworks Group CCT 27 March, 2007

Cactus Tutorial. Introduction to Cactus. Yaakoub El Khamra. Cactus Developer, Frameworks Group CCT 27 March, 2007 Cactus Tutorial Introduction to Cactus Yaakoub El Khamra Cactus Developer, Frameworks Group CCT 27 March, 2007 Agenda Introduction to Cactus What is Cactus Flesh and thorns Cactus Computational Toolkit

More information

What is Cactus? Cactus is a framework for developing portable, modular applications

What is Cactus? Cactus is a framework for developing portable, modular applications What is Cactus? Cactus is a framework for developing portable, modular applications What is Cactus? Cactus is a framework for developing portable, modular applications focusing, although not exclusively,

More information

Parallel PDE Solvers in Python

Parallel PDE Solvers in Python Parallel PDE Solvers in Python Bill Spotz Sandia National Laboratories Scientific Python 2006 August 18, 2006 Computational Sciences at Sandia Chemically reacting flows Climate modeling Combustion Compressible

More information

Teko: A Package for Multiphysics Preconditioners

Teko: A Package for Multiphysics Preconditioners SAND 2009-7242P Teko: A Package for Multiphysics Preconditioners Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energyʼs

More information

High Performance and Grid Computing Applications with the Cactus Framework. HPCC Program Grand Challenges (1995)

High Performance and Grid Computing Applications with the Cactus Framework. HPCC Program Grand Challenges (1995) High Performance and Grid Computing Applications with the Cactus Framework Gabrielle Allen Department of Computer Science Center for Computation & Technology Louisiana State University HPCC Program Grand

More information

The Cactus Framework: Design, Applications and Future Directions. Cactus Code

The Cactus Framework: Design, Applications and Future Directions. Cactus Code The Cactus Framework: Design, Applications and Future Directions Gabrielle Allen gallen@cct.lsu.edu Center for Computation & Technology Departments of Computer Science & Physics Louisiana State University

More information

The Cactus Framework. Erik Schnetter September 2006

The Cactus Framework. Erik Schnetter September 2006 The Cactus Framework Erik Schnetter September 2006 Outline History The Cactus User Community Cactus Usage Patterns Bird s eye view Cactus is a freely available, portable, and manageable environment for

More information

Cactus Framework: Scaling and Lessons Learnt

Cactus Framework: Scaling and Lessons Learnt Cactus Framework: Scaling and Lessons Learnt Gabrielle Allen, Erik Schnetter, Jian Tao Center for Computation & Technology Departments of Computer Science & Physics Louisiana State University Also: Christian

More information

MPI Related Software

MPI Related Software 1 MPI Related Software Profiling Libraries and Tools Visualizing Program Behavior Timing Performance Measurement and Tuning High Level Libraries Profiling Libraries MPI provides mechanism to intercept

More information

MPI Related Software. Profiling Libraries. Performance Visualization with Jumpshot

MPI Related Software. Profiling Libraries. Performance Visualization with Jumpshot 1 MPI Related Software Profiling Libraries and Tools Visualizing Program Behavior Timing Performance Measurement and Tuning High Level Libraries Performance Visualization with Jumpshot For detailed analysis

More information

Overview of Trilinos and PT-Scotch

Overview of Trilinos and PT-Scotch 29.03.2012 Outline PT-Scotch 1 PT-Scotch The Dual Recursive Bipartitioning Algorithm Parallel Graph Bipartitioning Methods 2 Overview of the Trilinos Packages Examples on using Trilinos PT-Scotch The Scotch

More information

PyTrilinos: A Python Interface to Trilinos, a Set of Object-Oriented Solver Packages

PyTrilinos: A Python Interface to Trilinos, a Set of Object-Oriented Solver Packages PyTrilinos: A Python Interface to Trilinos, a Set of Object-Oriented Solver Packages Bill Spotz Sandia National Laboratories SciPy 2005 Pasadena, CA 22 Sep 2005 With special thanks to Marzio Sala, Eric

More information

ACCELERATING CFD AND RESERVOIR SIMULATIONS WITH ALGEBRAIC MULTI GRID Chris Gottbrath, Nov 2016

ACCELERATING CFD AND RESERVOIR SIMULATIONS WITH ALGEBRAIC MULTI GRID Chris Gottbrath, Nov 2016 ACCELERATING CFD AND RESERVOIR SIMULATIONS WITH ALGEBRAIC MULTI GRID Chris Gottbrath, Nov 2016 Challenges What is Algebraic Multi-Grid (AMG)? AGENDA Why use AMG? When to use AMG? NVIDIA AmgX Results 2

More information

Linear Algebra libraries in Debian. DebConf 10 New York 05/08/2010 Sylvestre

Linear Algebra libraries in Debian. DebConf 10 New York 05/08/2010 Sylvestre Linear Algebra libraries in Debian Who I am? Core developer of Scilab (daily job) Debian Developer Involved in Debian mainly in Science and Java aspects sylvestre.ledru@scilab.org / sylvestre@debian.org

More information

An Overview of Trilinos

An Overview of Trilinos SAND REPORT SAND2003-2927 Unlimited Release Printed August 2003 An Overview of Trilinos Michael Heroux, Roscoe Bartlett, Vicki Howle Robert Hoekstra, Jonathan Hu, Tamara Kolda, Richard Lehoucq, Kevin Long,

More information

Algorithms, System and Data Centre Optimisation for Energy Efficient HPC

Algorithms, System and Data Centre Optimisation for Energy Efficient HPC 2015-09-14 Algorithms, System and Data Centre Optimisation for Energy Efficient HPC Vincent Heuveline URZ Computing Centre of Heidelberg University EMCL Engineering Mathematics and Computing Lab 1 Energy

More information

Ocean Modeling. Infrastructure (COMI) at LSU

Ocean Modeling. Infrastructure (COMI) at LSU Development of Coastal Ocean Modeling Infrastructure (COMI) at LSU Q. Jim Chen Department of Civil and Environmental Engineering & Center for Computation and Technology Louisiana State University Acknowledgements

More information

ML 3.1 Smoothed Aggregation User s Guide

ML 3.1 Smoothed Aggregation User s Guide SAND2004 4819 Unlimited Release Printed September 2004 ML 3.1 Smoothed Aggregation User s Guide Marzio Sala Computational Math & Algorithms Sandia National Laboratories P.O. Box 5800 Albuquerque, NM 87185-1110

More information

Preparation of Codes for Trinity

Preparation of Codes for Trinity Preparation of Codes for Trinity Courtenay T. Vaughan, Mahesh Rajan, Dennis C. Dinge, Clark R. Dohrmann, Micheal W. Glass, Kenneth J. Franko, Kendall H. Pierson, and Michael R. Tupek Sandia National Laboratories

More information

Web Interface to Materials Simulations

Web Interface to Materials Simulations Web Interface to Materials Simulations Web Interface Generator and Legacy Application Façade Portal Development Team Funding Akos J. Czikmantory (JPL - Wiglaf) DARPA-PROM Hook Hua (JPL - Wiglaf ) JPL SRRF

More information

Grid Computing in Numerical Relativity and Astrophysics

Grid Computing in Numerical Relativity and Astrophysics Grid Computing in Numerical Relativity and Astrophysics Gabrielle Allen: gallen@cct.lsu.edu Depts Computer Science & Physics Center for Computation & Technology (CCT) Louisiana State University Challenge

More information

Epetra Performance Optimization Guide

Epetra Performance Optimization Guide SAND2005-1668 Unlimited elease Printed March 2005 Updated for Trilinos 9.0 February 2009 Epetra Performance Optimization Guide Michael A. Heroux Scalable Algorithms Department Sandia National Laboratories

More information

Using Java for Scientific Computing. Mark Bul EPCC, University of Edinburgh

Using Java for Scientific Computing. Mark Bul EPCC, University of Edinburgh Using Java for Scientific Computing Mark Bul EPCC, University of Edinburgh markb@epcc.ed.ac.uk Java and Scientific Computing? Benefits of Java for Scientific Computing Portability Network centricity Software

More information

computational Fluid Dynamics - Prof. V. Esfahanian

computational Fluid Dynamics - Prof. V. Esfahanian Three boards categories: Experimental Theoretical Computational Crucial to know all three: Each has their advantages and disadvantages. Require validation and verification. School of Mechanical Engineering

More information

Research Collection. WebParFE A web interface for the high performance parallel finite element solver ParFE. Report. ETH Library

Research Collection. WebParFE A web interface for the high performance parallel finite element solver ParFE. Report. ETH Library Research Collection Report WebParFE A web interface for the high performance parallel finite element solver ParFE Author(s): Paranjape, Sumit; Kaufmann, Martin; Arbenz, Peter Publication Date: 2009 Permanent

More information

Intel Performance Libraries

Intel Performance Libraries Intel Performance Libraries Powerful Mathematical Library Intel Math Kernel Library (Intel MKL) Energy Science & Research Engineering Design Financial Analytics Signal Processing Digital Content Creation

More information

Dynamic Selection of Auto-tuned Kernels to the Numerical Libraries in the DOE ACTS Collection

Dynamic Selection of Auto-tuned Kernels to the Numerical Libraries in the DOE ACTS Collection Numerical Libraries in the DOE ACTS Collection The DOE ACTS Collection SIAM Parallel Processing for Scientific Computing, Savannah, Georgia Feb 15, 2012 Tony Drummond Computational Research Division Lawrence

More information

Ilya Lashuk, Merico Argentati, Evgenii Ovtchinnikov, Andrew Knyazev (speaker)

Ilya Lashuk, Merico Argentati, Evgenii Ovtchinnikov, Andrew Knyazev (speaker) Ilya Lashuk, Merico Argentati, Evgenii Ovtchinnikov, Andrew Knyazev (speaker) Department of Mathematics and Center for Computational Mathematics University of Colorado at Denver SIAM Conference on Parallel

More information

Anna Morajko.

Anna Morajko. Performance analysis and tuning of parallel/distributed applications Anna Morajko Anna.Morajko@uab.es 26 05 2008 Introduction Main research projects Develop techniques and tools for application performance

More information

Numerical Methods in Scientific Computation

Numerical Methods in Scientific Computation Numerical Methods in Scientific Computation Programming and Software Introduction to error analysis 1 Packages vs. Programming Packages MATLAB Excel Mathematica Maple Packages do the work for you Most

More information

Introduction to Parallel Computing

Introduction to Parallel Computing Introduction to Parallel Computing W. P. Petersen Seminar for Applied Mathematics Department of Mathematics, ETHZ, Zurich wpp@math. ethz.ch P. Arbenz Institute for Scientific Computing Department Informatik,

More information

A Python HPC framework: PyTrilinos, ODIN, and Seamless

A Python HPC framework: PyTrilinos, ODIN, and Seamless A Python HPC framework: PyTrilinos, ODIN, and Seamless K.W. Smith Enthought, Inc. 515 Congress Ave. Austin, TX 78701 ksmith@enthought.com W.F. Spotz Sandia National Laboratories P.O. Box 5800 Albuquerque,

More information

Intel Visual Fortran Compiler Professional Edition 11.0 for Windows* In-Depth

Intel Visual Fortran Compiler Professional Edition 11.0 for Windows* In-Depth Intel Visual Fortran Compiler Professional Edition 11.0 for Windows* In-Depth Contents Intel Visual Fortran Compiler Professional Edition for Windows*........................ 3 Features...3 New in This

More information

An algebraic multi-grid implementation in FEniCS for solving 3d fracture problems using XFEM

An algebraic multi-grid implementation in FEniCS for solving 3d fracture problems using XFEM An algebraic multi-grid implementation in FEniCS for solving 3d fracture problems using XFEM Axel Gerstenberger, Garth N. Wells, Chris Richardson, David Bernstein, and Ray Tuminaro FEniCS 14, Paris, France

More information

Mathematical Libraries and Application Software on JUQUEEN and JURECA

Mathematical Libraries and Application Software on JUQUEEN and JURECA Mitglied der Helmholtz-Gemeinschaft Mathematical Libraries and Application Software on JUQUEEN and JURECA JSC Training Course May 2017 I.Gutheil Outline General Informations Sequential Libraries Parallel

More information

Mathematical Libraries and Application Software on JUQUEEN and JURECA

Mathematical Libraries and Application Software on JUQUEEN and JURECA Mitglied der Helmholtz-Gemeinschaft Mathematical Libraries and Application Software on JUQUEEN and JURECA JSC Training Course November 2015 I.Gutheil Outline General Informations Sequential Libraries Parallel

More information

Recent Updates to the CFD General Notation System (CGNS)

Recent Updates to the CFD General Notation System (CGNS) Recent Updates to the CFD General Notation System (CGNS) C. L. Rumsey NASA Langley Research Center B. Wedan Computational Engineering Solutions T. Hauser University of Colorado M. Poinot ONERA AIAA-2012-1264,

More information

AmgX 2.0: Scaling toward CORAL Joe Eaton, November 19, 2015

AmgX 2.0: Scaling toward CORAL Joe Eaton, November 19, 2015 AmgX 2.0: Scaling toward CORAL Joe Eaton, November 19, 2015 Agenda Introduction to AmgX Current Capabilities Scaling V2.0 Roadmap for the future 2 AmgX Fast, scalable linear solvers, emphasis on iterative

More information

AllScale Pilots Applications AmDaDos Adaptive Meshing and Data Assimilation for the Deepwater Horizon Oil Spill

AllScale Pilots Applications AmDaDos Adaptive Meshing and Data Assimilation for the Deepwater Horizon Oil Spill This project has received funding from the European Union s Horizon 2020 research and innovation programme under grant agreement No. 671603 An Exascale Programming, Multi-objective Optimisation and Resilience

More information

Performance of deal.ii on a node

Performance of deal.ii on a node Performance of deal.ii on a node Bruno Turcksin Texas A&M University, Dept. of Mathematics Bruno Turcksin Deal.II on a node 1/37 Outline 1 Introduction 2 Architecture 3 Paralution 4 Other Libraries 5 Conclusions

More information

Presentation Outline. Some RSI Customers

Presentation Outline. Some RSI Customers Research Systems, Inc. Software Vision Presentation Outline Introducing Research Systems IDL overview ENVI - remote sensing application Visible Human - anatomical CD reference VIP RiverTools NeoSys The

More information

Lubuntu Linux Virtual Machine

Lubuntu Linux Virtual Machine Lubuntu Linux 18.04 Virtual Machine About Us Slide / 01 Founded in 2015, Sourcery Institute is a California nonprofit public-benefit corporation engaged in research, education, and advisory services in

More information

Teaching numerical methods : a first experience. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai.

Teaching numerical methods : a first experience. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai. Teaching numerical methods : a first experience Ronojoy Adhikari The Institute of Mathematical Sciences Chennai. Numerical Methods : the first class 14 students, spread across second and third years of

More information

Day 2 August 06, 2004 (Friday)

Day 2 August 06, 2004 (Friday) An Overview of Grid Computing Day 2 August 06, 2004 (Friday) By CDAC Experts Contact :vcvrao@cdacindia.com; betatest@cdacindia.com URL : http://www.cs.umn.edu/~vcvrao 1 Betatesting Group,NPSF, C-DAC,Pune

More information

An Object Oriented Finite Element Library

An Object Oriented Finite Element Library An Object Oriented Finite Element Library Release 3.1.0 Rachid Touzani Laboratoire de Mathématiques Blaise Pascal Université Clermont Auvergne 63177 Aubière, France e-mail: Rachid.Touzani@univ-bpclermont.fr

More information

Supercomputing and Science An Introduction to High Performance Computing

Supercomputing and Science An Introduction to High Performance Computing Supercomputing and Science An Introduction to High Performance Computing Part VII: Scientific Computing Henry Neeman, Director OU Supercomputing Center for Education & Research Outline Scientific Computing

More information

Parallelism. Wolfgang Kastaun. May 9, 2008

Parallelism. Wolfgang Kastaun. May 9, 2008 Parallelism Wolfgang Kastaun May 9, 2008 Outline Parallel computing Frameworks MPI and the batch system Running MPI code at TAT The CACTUS framework Overview Mesh refinement Writing Cactus modules Links

More information

Mathematical Libraries and Application Software on JUROPA, JUGENE, and JUQUEEN. JSC Training Course

Mathematical Libraries and Application Software on JUROPA, JUGENE, and JUQUEEN. JSC Training Course Mitglied der Helmholtz-Gemeinschaft Mathematical Libraries and Application Software on JUROPA, JUGENE, and JUQUEEN JSC Training Course May 22, 2012 Outline General Informations Sequential Libraries Parallel

More information

Module 4. Computer-Aided Design (CAD) systems

Module 4. Computer-Aided Design (CAD) systems Module 4. Computer-Aided Design (CAD) systems Nowadays the design of complex systems is unconceivable without computers. The fast computers, the sophisticated developing environments and the well elaborated

More information

On Level Scheduling for Incomplete LU Factorization Preconditioners on Accelerators

On Level Scheduling for Incomplete LU Factorization Preconditioners on Accelerators On Level Scheduling for Incomplete LU Factorization Preconditioners on Accelerators Karl Rupp, Barry Smith rupp@mcs.anl.gov Mathematics and Computer Science Division Argonne National Laboratory FEMTEC

More information

A Brief History of Numerical Libraries. Sven Hammarling NAG Ltd, Oxford & University of Manchester

A Brief History of Numerical Libraries. Sven Hammarling NAG Ltd, Oxford & University of Manchester A Brief History of Numerical Libraries Sven Hammarling NAG Ltd, Oxford & University of Manchester John Reid at Gatlinburg VII (1977) Solution of large finite element systems of linear equations out of

More information

Performances and Tuning for Designing a Fast Parallel Hemodynamic Simulator. Bilel Hadri

Performances and Tuning for Designing a Fast Parallel Hemodynamic Simulator. Bilel Hadri Performances and Tuning for Designing a Fast Parallel Hemodynamic Simulator Bilel Hadri University of Tennessee Innovative Computing Laboratory Collaboration: Dr Marc Garbey, University of Houston, Department

More information

An HPC Implementation of the Finite Element Method

An HPC Implementation of the Finite Element Method An HPC Implementation of the Finite Element Method John Rugis Interdisciplinary research group: David Yule Physiology James Sneyd, Shawn Means, Di Zhu Mathematics John Rugis Computer Science Project funding:

More information

On a Future Software Platform for Demanding Multi-Scale and Multi-Physics Problems

On a Future Software Platform for Demanding Multi-Scale and Multi-Physics Problems On a Future Software Platform for Demanding Multi-Scale and Multi-Physics Problems H. P. Langtangen X. Cai Simula Research Laboratory, Oslo Department of Informatics, Univ. of Oslo SIAM-CSE07, February

More information

Moderator: Edward Seidel, Director, Center for Computation & Technology, Louisiana State University

Moderator: Edward Seidel, Director, Center for Computation & Technology, Louisiana State University Ask A Grid Expert Panel & Interaction Thursday January 6, 2005, 4:00-5:30PM Moderator: Edward Seidel, Director, Center for Computation & Technology, Louisiana State University Panelists: Gabrielle Allen,

More information

Introduction to Multigrid and its Parallelization

Introduction to Multigrid and its Parallelization Introduction to Multigrid and its Parallelization! Thomas D. Economon Lecture 14a May 28, 2014 Announcements 2 HW 1 & 2 have been returned. Any questions? Final projects are due June 11, 5 pm. If you are

More information

Speedup Altair RADIOSS Solvers Using NVIDIA GPU

Speedup Altair RADIOSS Solvers Using NVIDIA GPU Innovation Intelligence Speedup Altair RADIOSS Solvers Using NVIDIA GPU Eric LEQUINIOU, HPC Director Hongwei Zhou, Senior Software Developer May 16, 2012 Innovation Intelligence ALTAIR OVERVIEW Altair

More information

An Overview of the Trilinos Project

An Overview of the Trilinos Project SAND REPORT SAND2003-xxxx Unlimited Release February 2003 An Overview of the Trilinos Project Michael A. Heroux, Sandia National Laboratories Prepared by Sandia National Laboratories Albuquerque, New Mexico

More information

Large-scale Gas Turbine Simulations on GPU clusters

Large-scale Gas Turbine Simulations on GPU clusters Large-scale Gas Turbine Simulations on GPU clusters Tobias Brandvik and Graham Pullan Whittle Laboratory University of Cambridge A large-scale simulation Overview PART I: Turbomachinery PART II: Stencil-based

More information

Using Linear Programming for Management Decisions

Using Linear Programming for Management Decisions Using Linear Programming for Management Decisions By Tim Wright Linear programming creates mathematical models from real-world business problems to maximize profits, reduce costs and allocate resources.

More information

An Example of Porting PETSc Applications to Heterogeneous Platforms with OpenACC

An Example of Porting PETSc Applications to Heterogeneous Platforms with OpenACC An Example of Porting PETSc Applications to Heterogeneous Platforms with OpenACC Pi-Yueh Chuang The George Washington University Fernanda S. Foertter Oak Ridge National Laboratory Goal Develop an OpenACC

More information

Numerical Libraries for Petascale Computing. Performance Estimate. Why Use Libraries? Faster (Better Code)

Numerical Libraries for Petascale Computing. Performance Estimate. Why Use Libraries? Faster (Better Code) Why Use Libraries? Numerical Libraries for Petascale Computing William Gropp www.cs.uiuc.edu/homes/wgropp There are many reasons to use libraries: Faster - Code tricks Faster Better Algorithms Correct

More information

Preconditioning Linear Systems Arising from Graph Laplacians of Complex Networks

Preconditioning Linear Systems Arising from Graph Laplacians of Complex Networks Preconditioning Linear Systems Arising from Graph Laplacians of Complex Networks Kevin Deweese 1 Erik Boman 2 1 Department of Computer Science University of California, Santa Barbara 2 Scalable Algorithms

More information

Numlua: a numerical package for Lua

Numlua: a numerical package for Lua Numlua: a numerical package for Lua Luis Carvalho Brown University carvalho@dam.brown.edu Lua Workshop July, 2008 Introduction A lot of problems in scientific computing need high performance numerical

More information

Computational Steering. Nate Woody Drew Dolgert

Computational Steering. Nate Woody Drew Dolgert Computational Steering Nate Woody Drew Dolgert Lab Materials In with the other labs. compsteer/simple compsteer/gauss_steer 12/9/2010 www.cac.cornell.edu 2 What is Computational Steering? Interactivity

More information

A Parallel Implementation of the BDDC Method for Linear Elasticity

A Parallel Implementation of the BDDC Method for Linear Elasticity A Parallel Implementation of the BDDC Method for Linear Elasticity Jakub Šístek joint work with P. Burda, M. Čertíková, J. Mandel, J. Novotný, B. Sousedík Institute of Mathematics of the AS CR, Prague

More information

Development of an Integrated Modeling Framework for Simulations of Coastal Processes in Deltaic Environments Using High-Performance Computing

Development of an Integrated Modeling Framework for Simulations of Coastal Processes in Deltaic Environments Using High-Performance Computing Development of an Integrated Modeling Framework for Simulations of Coastal Processes in Deltaic Environments Using High-Performance Computing Q. Jim Chen Department of Civil and Environmental Engineering

More information

Review of previous examinations TMA4280 Introduction to Supercomputing

Review of previous examinations TMA4280 Introduction to Supercomputing Review of previous examinations TMA4280 Introduction to Supercomputing NTNU, IMF April 24. 2017 1 Examination The examination is usually comprised of: one problem related to linear algebra operations with

More information

MS6021 Scientific Computing. MatLab and Python for Mathematical Modelling. Aimed at the absolute beginner.

MS6021 Scientific Computing. MatLab and Python for Mathematical Modelling. Aimed at the absolute beginner. MS6021 Scientific Computing MatLab and Python for Mathematical Modelling. Aimed at the absolute beginner. Natalia Kopteva Email: natalia.kopteva@ul.ie Web: http://www.staff.ul.ie/natalia/ Room: B2037 Office

More information

Computational Steering

Computational Steering Computational Steering Nate Woody 10/13/2009 www.cac.cornell.edu 1 Lab Materials I ve placed some sample code in ~train100 that performs the operations that I ll demonstrate during this talk. We ll walk

More information

Scalable Algorithms in Optimization: Computational Experiments

Scalable Algorithms in Optimization: Computational Experiments Scalable Algorithms in Optimization: Computational Experiments Steven J. Benson, Lois McInnes, Jorge J. Moré, and Jason Sarich Mathematics and Computer Science Division, Argonne National Laboratory, Argonne,

More information

Scientific Components and Frameworks

Scientific Components and Frameworks High Performance Computing: Concepts, Methods & Means Scientific Components and Frameworks Prof. Daniel S. Katz Department of Electrical and Computer Engineering Louisiana State University April 24 th,

More information

The GridPACK Toolkit for Developing Power Grid Simulations on High Performance Computing Platforms

The GridPACK Toolkit for Developing Power Grid Simulations on High Performance Computing Platforms The GridPACK Toolkit for Developing Power Grid Simulations on High Performance Computing Platforms Bruce Palmer, Bill Perkins, Kevin Glass, Yousu Chen, Shuangshuang Jin, Ruisheng Diao, Mark Rice, David

More information

WORHP Lab The Graphical User Interface for Optimisation and Optimal Control

WORHP Lab The Graphical User Interface for Optimisation and Optimal Control WORHP Lab The Graphical User Interface for Optimisation and Optimal Control M. Knauer, C. Büskens Zentrum für Universität Bremen 3rd European Optimisation in Space Engineering 17 th - 18 th September 2015

More information

The Cactus Framework and Toolkit: Design and Applications

The Cactus Framework and Toolkit: Design and Applications The Cactus Framework and Toolkit: Design and Applications Tom Goodale 1, Gabrielle Allen 1, Gerd Lanfermann 1, Joan Massó 2, Thomas Radke 1, Edward Seidel 1, and John Shalf 3 1 Max-Planck-Institut für

More information

Introduction to Parallel Programming & Cluster Computing

Introduction to Parallel Programming & Cluster Computing Introduction to Parallel Programming & Cluster Computing Scientific Libraries & I/O Libraries Joshua Alexander, U Oklahoma Ivan Babic, Earlham College Michial Green, Contra Costa College Mobeen Ludin,

More information

AMS526: Numerical Analysis I (Numerical Linear Algebra)

AMS526: Numerical Analysis I (Numerical Linear Algebra) AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 20: Sparse Linear Systems; Direct Methods vs. Iterative Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 26

More information

Java Performance Analysis for Scientific Computing

Java Performance Analysis for Scientific Computing Java Performance Analysis for Scientific Computing Roldan Pozo Leader, Mathematical Software Group National Institute of Standards and Technology USA UKHEC: Java for High End Computing Nov. 20th, 2000

More information

Open Compute Stack (OpenCS) Overview. D.D. Nikolić Updated: 20 August 2018 DAE Tools Project,

Open Compute Stack (OpenCS) Overview. D.D. Nikolić Updated: 20 August 2018 DAE Tools Project, Open Compute Stack (OpenCS) Overview D.D. Nikolić Updated: 20 August 2018 DAE Tools Project, http://www.daetools.com/opencs What is OpenCS? A framework for: Platform-independent model specification 1.

More information

PETSc Satish Balay, Kris Buschelman, Bill Gropp, Dinesh Kaushik, Lois McInnes, Barry Smith

PETSc   Satish Balay, Kris Buschelman, Bill Gropp, Dinesh Kaushik, Lois McInnes, Barry Smith PETSc http://www.mcs.anl.gov/petsc Satish Balay, Kris Buschelman, Bill Gropp, Dinesh Kaushik, Lois McInnes, Barry Smith PDE Application Codes PETSc PDE Application Codes! ODE Integrators! Nonlinear Solvers,!

More information

Parallel Libraries And ToolBoxes for PDEs Luca Heltai

Parallel Libraries And ToolBoxes for PDEs Luca Heltai The 2nd Workshop on High Performance Computing Parallel Libraries And ToolBoxes for PDEs Luca Heltai SISSA/eLAB - Trieste Shahid Beheshti University, Institute for Studies in Theoretical Physics and Mathematics

More information

CMSC 714 Lecture 6 MPI vs. OpenMP and OpenACC. Guest Lecturer: Sukhyun Song (original slides by Alan Sussman)

CMSC 714 Lecture 6 MPI vs. OpenMP and OpenACC. Guest Lecturer: Sukhyun Song (original slides by Alan Sussman) CMSC 714 Lecture 6 MPI vs. OpenMP and OpenACC Guest Lecturer: Sukhyun Song (original slides by Alan Sussman) Parallel Programming with Message Passing and Directives 2 MPI + OpenMP Some applications can

More information

PySparse and PyFemax: A Python framework for large scale sparse linear algebra

PySparse and PyFemax: A Python framework for large scale sparse linear algebra PySparse and PyFemax: A Python framework for large scale sparse linear algebra Roman Geus and Peter Arbenz Swiss Federal Institute of Technology Zurich Institute of Scientific Computing mailto:geus@inf.ethz.ch

More information

Computational Steering

Computational Steering Computational Steering Nate Woody 10/23/2008 www.cac.cornell.edu 1 What is computational steering? Generally, computational steering can be thought of as a method (or set of methods) for providing interactivity

More information

What is Multigrid? They have been extended to solve a wide variety of other problems, linear and nonlinear.

What is Multigrid? They have been extended to solve a wide variety of other problems, linear and nonlinear. AMSC 600/CMSC 760 Fall 2007 Solution of Sparse Linear Systems Multigrid, Part 1 Dianne P. O Leary c 2006, 2007 What is Multigrid? Originally, multigrid algorithms were proposed as an iterative method to

More information

Performance of Implicit Solver Strategies on GPUs

Performance of Implicit Solver Strategies on GPUs 9. LS-DYNA Forum, Bamberg 2010 IT / Performance Performance of Implicit Solver Strategies on GPUs Prof. Dr. Uli Göhner DYNAmore GmbH Stuttgart, Germany Abstract: The increasing power of GPUs can be used

More information

Dynamic Deployment of a Component Framework with the Ubiqis System

Dynamic Deployment of a Component Framework with the Ubiqis System Dynamic Deployment of a Component Framework with the Ubiqis System Steven Brandt (1) Gabrielle Allen (1,2) Matthew Eastman (1) Matthew Kemp (1) Erik Schnetter (1,3) (1) Center for Computation & Technology,

More information

Massimiliano Lupo Pasini

Massimiliano Lupo Pasini Massimiliano Lupo Pasini 400 Dowman Drive, W401 - Atlanta, GA, 30322, USA +1 (678) 536 5533 massimiliano.lupo.pasini@emory.edu http://www.mathcs.emory.edu/ mlupopa/ Curriculum Vitae Citizenship USA VISA

More information

Advanced High Performance Computing CSCI 580

Advanced High Performance Computing CSCI 580 Advanced High Performance Computing CSCI 580 2:00 pm - 3:15 pm Tue & Thu Marquez Hall 322 Timothy H. Kaiser, Ph.D. tkaiser@mines.edu CTLM 241A http://inside.mines.edu/~tkaiser/csci580fall13/ 1 Two Similar

More information

Multithreaded Programming in Cilk. Matteo Frigo

Multithreaded Programming in Cilk. Matteo Frigo Multithreaded Programming in Cilk Matteo Frigo Multicore challanges Development time: Will you get your product out in time? Where will you find enough parallel-programming talent? Will you be forced to

More information

Speeding up MATLAB Applications Sean de Wolski Application Engineer

Speeding up MATLAB Applications Sean de Wolski Application Engineer Speeding up MATLAB Applications Sean de Wolski Application Engineer 2014 The MathWorks, Inc. 1 Non-rigid Displacement Vector Fields 2 Agenda Leveraging the power of vector and matrix operations Addressing

More information

Introduction to MATLAB

Introduction to MATLAB Introduction to MATLAB Contents 1.1 Objectives... 1 1.2 Lab Requirement... 1 1.3 Background of MATLAB... 1 1.4 The MATLAB System... 1 1.5 Start of MATLAB... 3 1.6 Working Modes of MATLAB... 4 1.7 Basic

More information

Cactus Tutorial. Thorn Writing 101. Yaakoub Y El Khamra. Frameworks Group, CCT 17 Feb, 2006

Cactus Tutorial. Thorn Writing 101. Yaakoub Y El Khamra. Frameworks Group, CCT 17 Feb, 2006 Cactus Tutorial Thorn Writing 101 Yaakoub Y El Khamra Frameworks Group, CCT 17 Feb, 2006 Agenda Downloading Cactus The GetCactus Script CVS Thorn writing Making a new thorn Editing the.ccl files Adding

More information

High-level Abstraction for Block Structured Applications: A lattice Boltzmann Exploration

High-level Abstraction for Block Structured Applications: A lattice Boltzmann Exploration High-level Abstraction for Block Structured Applications: A lattice Boltzmann Exploration Jianping Meng, Xiao-Jun Gu, David R. Emerson, Gihan Mudalige, István Reguly and Mike B Giles Scientific Computing

More information

XSEDE Visualization Use Cases

XSEDE Visualization Use Cases XSEDE Visualization Use Cases July 24, 2014 Version 1.4 XSEDE Visualization Use Cases Page i Table of Contents A. Document History iii B. Document Scope iv XSEDE Visualization Use Cases Page ii A. Document

More information

Introduction to Linear Programming. Algorithmic and Geometric Foundations of Optimization

Introduction to Linear Programming. Algorithmic and Geometric Foundations of Optimization Introduction to Linear Programming Algorithmic and Geometric Foundations of Optimization Optimization and Linear Programming Mathematical programming is a class of methods for solving problems which ask

More information

Parallel High-Order Geometric Multigrid Methods on Adaptive Meshes for Highly Heterogeneous Nonlinear Stokes Flow Simulations of Earth s Mantle

Parallel High-Order Geometric Multigrid Methods on Adaptive Meshes for Highly Heterogeneous Nonlinear Stokes Flow Simulations of Earth s Mantle ICES Student Forum The University of Texas at Austin, USA November 4, 204 Parallel High-Order Geometric Multigrid Methods on Adaptive Meshes for Highly Heterogeneous Nonlinear Stokes Flow Simulations of

More information

Optimization and Scalability

Optimization and Scalability Optimization and Scalability Drew Dolgert CAC 29 May 2009 Intro to Parallel Computing 5/29/2009 www.cac.cornell.edu 1 Great Little Program What happens when I run it on the cluster? How can I make it faster?

More information

CCA-LISI: On Designing A CCA Parallel Sparse Linear Solver Interface

CCA-LISI: On Designing A CCA Parallel Sparse Linear Solver Interface CCA-LISI: On Designing A CCA Parallel Sparse Linear Solver Interface Fang (Cherry) Liu and Randall Bramley Indiana University Bloomington, Indiana {fangliu,bramley}@indiana.edu Abstract Sparse linear solvers

More information