Legacy codes are in ubiquitous use in scientific simulations; they are well-tested and there is significant time investment in their use. However, one challenge is the adoption of new, sometimes incompatible computing paradigms, such as GPU hardware. In this paper, we explore using automated code translation to enable execution of legacy multigrid solver code on GPUs without significant time investment and while avoiding intrusive changes to the codebase. We developed a thin, reusable translation layer that parses Fortran 2003 at compile time, interfacing with the existing library Loopy (Klöckner, 2014) to transpile to C++/GPU code, which is then managed by a custom MPI runtime system that we created. With this low-effort approach, we are able to achieve a payoff of an approximately 2–3× speedup over a full CPU socket, and 6× in multi-node settings.
ChristenMSchenkOBurkhartH (2011) Patus: a code generation and autotuning framework for parallel iterative stencil computations on modern microarchitectures. 2011 IEEE International Parallel & Distributed Processing Symposium: 676–687. https://ieeexplore.ieee.org/document/6012879/
HagedornBStoltzfusLSteuwerM, et al. (2018) High performance stencil code generation with lift. New York: ACM, 100–112.
7.
HerdmanJAGaudinWPPerksO, et al. (2014) Achieving portability and performance through OpenACC. In: 2014 First Workshop on Accelerator Programming Using Directives, 19–26. https://ieeexplore.ieee.org/document/7081674
8.
HolewinskiJPouchetLNSadayappanP (2012) High-performance code generation for stencil computations on gpu architectures. New York: ACM, 311–320.
9.
KlöcknerA (2014) Loo.py: transformation-based code generation for GPUs and CPUs. In: Proceedings of ACM SIGPLAN International Workshop on Libraries, Languages, and Compilers for Array Programming, ARRAY’14. New York, NY, USA: Association for Computing Machinery, 82–87.
10.
KlöcknerA (2015) Loo.py: from Fortran to performance via transformation and substitution rules. In: Proceedings of the 2nd ACM SIGPLAN International Workshop on Libraries, Languages, and Compilers for Array Programming, ARRAY 2015. New York, NY, USA: Association for Computing Machinery, 1–6.
KulkarniK (2023) Domain-Specific Code Transformations for Computational Science Based on the Polyhedral Model. PhD thesis. Urbana-Champaign, Urbana, IL: University of Illinois. https://www.ideals.illinois.edu/items/129117
SteuwerMFenschCLindleyS, et al. (2015) Generating performance portable code using rewrite rules: from high-level functional expressions to high-performance opencl code. ACM SIGPLAN Notices50: 205–217.
20.
VerdoolaegeS (2010) isl: an integer set library for the polyhedral model. In: FukudaKJvdHJoswigM, et al. (eds) Mathematical Software – ICMS 2010. Berlin, Heidelberg: Springer Berlin Heidelberg, 299–302.
