Abstract:Developers of Molecular Dynamics (MD) codes face significant challenges when adapting existing simulation packages to new hardware. In a continuously diversifying hardware landscape it becomes increasingly difficult for scientists to be experts both in their own domain (physics/chemistry/biology) and specialists in the low level parallelisation and optimisation of their codes. To address this challenge, we describe a "Separation of Concerns" approach for the development of parallel and optimised MD codes: the science specialist writes code at a high abstraction level in a domain specific language (DSL), which is then translated into efficient computer code by a scientific programmer. In a related context, an abstraction for the solution of partial differential equations with grid based methods has recently been implemented in the (Py)OP2 library. Inspired by this approach, we develop a Python code generation system for molecular dynamics simulations on different parallel architectures, including massively parallel distributed memory systems and GPUs. We demonstrate the efficiency of the auto-generated code by studying its performance and scalability on different hardware and compare it to other state-of-the-art simulation packages. With growing data volumes the extraction of physically meaningful information from the simulation becomes increasingly challenging and requires equally efficient implementations. A particular advantage of our approach is the easy expression of such analysis algorithms. We consider two popular methods for deducing the crystalline structure of a material from the local environment of each atom, show how they can be expressed in our abstraction and implement them in the code generation framework.

Devito: Towards a generic Finite Difference DSL using Symbolic Python

Automated MPI-X code generation for scalable finite-difference solvers

A shared compilation stack for distributed-memory parallelism in stencil DSLs

Scaling through abstractions -- high-performance vectorial wave simulations for seismic inversion with Devito

Yet Another Tensor Toolbox for discontinuous Galerkin methods and other applications

One DSL to Rule Them All: IDE-Assisted Code Generation for Agile Data Analysis

Domain-specific languages: from design to implementation application to video device drivers generation

An exploration to visualize finite element data with a DSL

Design and implementation of DeepDSL: A DSL for deep learning

Mat2Stencil: A Modular Matrix-Based DSL for Explicit and Implicit Matrix-Free PDE Solvers on Structured Grid

DeepDSL: A Compilation-based Domain-Specific Language for Deep Learning

A Domain Specific Language for Performance Portable Molecular Dynamics Algorithms

GPU Support for Automatic Generation of Finite-Differences Stencil Kernels

Developing Web-based Geographic Information Systems with a DSL: Proposal and Case Study

Domain-Specific Language Techniques for Visual Computing: A Comprehensive Study

Generating Finite Element Codes combining Adaptive Octrees with Complex Geometries

Applying the Polyhedral Model to Tile Time Loops in Devito

A Domain-Specific Language and Editor for Parallel Particle Methods

Building-Blocks for Performance Oriented DSLs

Automatic Code Generation and Optimization of Large-scale Stencil Computation on Many-core Processors.

Python-based DSL for generating Verilog model of Synchronous Digital Circuits