The DiamondCandy Algorithm for Maximum Performance Vectorized Cross-Stencil Computation

Abstract

Алгоритм DiamondCandyV для повышения производительности конечноразностных вычислений c применением векторизации На основе поиска оптимальной пространственно-временной декомпозиции 3D1T пространства операций построен новый алгоритм DiamondCandyV для эффективной реализации конечно-разностных вычислений с использованием аппаратной векторизации. Базовый элемент разбиения пространства операций DiamondCandy получен пересечением областей зависимостей и влияния (конусоидов) для схем с шаблоном типа крест. Благодаря этому алгоритм характеризуется высокой вычислительной интенсивностью и локализует обрабатываемые данные на верхних уровнях иерархии памяти современных компьютеров. Ключевой особенностью алгоритма является легко реализуемая в программном коде поддержка двух основных средств повышения производительности современных процессоров, а именно аппаратной векторизации (с использованием SIMD расширения AVX) и вычислительных потоков с общей памятью (many-core CPU). Обсуждаются детали программной реализации поддержки параллельности различных уровней на примере численного решения волнового уравнения. Результаты тестирования реализации алгоритма показывают повышение производительности на порядок по сравнению с традиционными алгоритмами с пошаговой синхронизацией. Также, в отличие от традиционного подхода, с увеличением размера обрабатываемых данных производительность не деградирует.

Numerical simulation of increasing initial perturbations of a bubble in the bubble–shock interaction problem

“…A special algorithm of implementation called 'DiamondTorre' (Perepelkina et al 2014, Korneev and) is used for this purpose. It uses the concepts of spatio-temporal tiling (Wonnacott 2002, Bandishti et al 2012 and locally recursive non-local asynchronicity (Levchenko 2005, Levchenko et al 2011, Perepelkina and Levchenko 2015. The algorithm is based on the topological properties of the stencil.…”