As part of a HECToR distributed CSE support project, a general-purpose 2D decomposition (also known as 'pencil' or 'drawer' decomposition) communication library has been developed. This Fortran library provides a powerful and flexible framework to build applications based on 3D Cartesian data structures and spatially implicit numerical schemes (such as the compact finite difference method or spectral method). The library also supports shared-memory architecture which becomes increasingly popular. A user-friendly FFT interface has been built on top of the communication library to perform distributed three-dimensional FFTs. Both the decomposition library and the FFT interface scale well to tens of thousands of cores on Cray XT systems. The library has been applied to Incompact3D, a CFD application performing large-scale Direct Numerical Simulations of turbulence, enabling exciting scientific studies to be conducted.
[1]
Warren E. Ferguson.
A simple derivation of Glassman's general N fast fourier transform☆
,
1982
.
[2]
P. Swarztrauber.
THE METHODS OF CYCLIC REDUCTION, FOURIER ANALYSIS AND THE FACR ALGORITHM FOR THE DISCRETE SOLUTION OF POISSON'S EQUATION ON A RECTANGLE*
,
1977
.
[3]
Sylvain Laizet,et al.
High-order compact schemes for incompressible flows: A simple and efficient method with quasi-spectral accuracy
,
2009,
J. Comput. Phys..
[4]
John Christos Vassilicos,et al.
Dissipation and decay of fractal-generated turbulence
,
2007
.
[5]
John Christos Vassilicos,et al.
Scalings and decay of fractal-generated turbulence
,
2007
.
[6]
Robert B. Wilhelmson,et al.
Direct solutions for Poisson's equation in three dimensions
,
1977
.