GPGPU
From CFD-Wiki
This article needs good pictures and diagrams.
Contents |
Introduction
GPGPU is an acronym for General Purpose Graphic Processor Unit. A GPGPU is any graphics processor being used for general computing beyond graphics. GPUsare widely available, and often targeted at the computer gaming industry. Graphics workloads are very parallel, and so GPUs developed as large-scale parallel computation machines. Originally GPGPU processing was done by tricking the GPU by disguising computation loads as graphic loads. In recent years, GPU manufacturers have been actively encouraging GPGPU computing with the release of specialized languages which support GPGPU commands. GPUs incorporate many more computational cores than their equivalent CPUs, and so the performance of parallel operations can be greatly enhanced. Programming in parallel on a GPU has the same justification given for parallel computing in general.
Application to CFD
GPGPU computing offers large amounts of compute power, which can be tapped for parallel components of CFD algorithms, while the CPU performs the serial portions of the algorithm. GPGPU languages also support data-parallel computation, similar to vector processors. In short, modern GPUs provide raw computational power orders of magnitude larger than a CPU and can fit inside a single computer case.
Graphics Architecture
A GPU has a main memory, many stages and parallel processors. Traditional GPUs have a linear pipeline with several distinct stages: application, command, geometry, rasterization, fragment, and display. Intel's project Larrabee promises a reconfigurable graphics pipeline, with many of the traditional steps being handled in software. Such a development would expose even more of the GPUs compute power to parallel programmers.
Traditional Pipeline
A traditional pipeline will have three main computation stages: geometry, rasterization, and fragment. Graphics is traditionally done with triangles, and a GPU will operate on a batch of triangle verticies to first create fragments, which will help create the pixels that end up on the monitor.
Geometry
Vertex processing is handled in the geometry step. Geometry from the CPU is transformed based on the vertex shaders (programs) written to the GPU. These processors specialized in matrix transformations. Common operations include projecting 3D coordinates onto 2D screen coordinates. The closest analogue would be a vector or quaternion processor since each vector operation takes a series of components which represent a triangle vertex. Lagrangian frame computations might be well suited to vertex shaders.
Rasterization
Rasterization takes the transformed vectors from the geometry step and creates fragments from the geometry. The easiest way to think of rasterization, is "chunking" a large triangle into many fragments. This stage is typically done with fixed function specialized hardware.
Fragment
Fragment processing requires floating point math, as the fragments are colored and filtered to become pixels. This stage is where muchof the interesting compute for CFD can happen, as the parallel floating point processors can be repurposed with either fragment shaders or special purpose languages to do non-graphics floating point math.
Performance
NVidia has been actively promoting GPGPU computing in recent years. They introduced CUDA in 2007, giving programmer direct access to the GPU. NVidia maintains a CUDA community showcase on their CUDA website showing the performance boost to a variety of applications when making use of the GPU.
When programming GPGPU systems, it is important to remember that there is significant overhead involved in transferring data between the CPU and the GPU. Serial programs with memory requirements beyond those available on the GPU will not see the same dramatic performance enhancements as programs with low memory transfer overhead.
Languages
There exist several languages which support direct control of the GPU. OpenCL, Microsoft's DirectCompute, and Nvidia's CUDA are good examples of these. For more information about how to program on the GPU, see the Wikipedia site on GPGPU computing. Coding syntax is similar to C/C++ programming syntax.
Additionally, if a programmer wants to use the graphics card to compute, and then display the results of their CFD equations. OpenGL hooks into C++ to enable that, and Microsoft's DirectX performs a similar role.
All of these languages have active communities participating in their further development. There are many code samples and tools on the internet to help a programmer get started with GPGPU computing. Go for it!
References
Coding Resources
Academic Resources
Resources with a * require academic subscriptions for papers, or book purchases, you can often use your local university library.
Brandvik, Tobias and Graham Pullan (2008), "Acceleration of a 3D Euler Solver using Commodity Graphics Hardware", Proceedings of the 46th AIAA Aerospace Sciences Meeting (view).
Cohen, J. and Molemaker, M.J. (1997), "A fast double precision CFD code using CUDA", Journal of Physics Soc. Japan vol. 66 pp. 2237 - 2341, 1997(view).
Christen, M., O. Schenk, E. Neufeld, P. Messmer, and H. Burkhart (2009), "Parallel Data-Locality Aware Stencil Computations on Modern Micro-architectures", Proceedings of the 2009 IEEE International Symposium on Parallel & Distributed Processing-Volume(view).
Elsen, E. and Legresley, P. and Darve, E. (2008), "Large calculation of the flow over a hypersonic vehicle using a GPU", Journal of Computational Physics vol 227 num. 24 pp. 10148 - 10161(view).
Hagen, Trond Runar, Knut-Andreas Lie and Jostein R. Natvig (2006), Solving the Euler Equations on Graphics Processing Units, in collection Proceedings of the 6th International Conference on Computational Science, pp. 220 - 227 (view*).
Harris, Mark (2004), Fast Fluid Dynamics Simulation on the GPU, Ed. Fernando Randima, in collection GPU Gems, Ch. 38 (view*).
Kolb, Andreas and Nicolas Cuntz (2005), "Dynamic Particle Coupling for GPU-based Fluid Simulation", Proceedings of the 18th Symposium on Simulation Technique, pp. 722 - 727 (view).
Phillips, Everett H., Yao Zhang, Roger L. Davis, and John D. Owens (2009), "Rapid Aerodynamic Performance Prediction on a Cluster of Graphics Processing Units", Proceedings of the 47th AIAA Aerospace Sciences Meeting (view).
Liu, Youguan, Xuehui Liu and Enhua Wu (2004), "Real-Time 3D Fluid Simulation on the GPU with complex obstacles", Proceedings of Pacific Graphics 2004, pp. 247 - 256 (view*).
