CowFreedom
Joined 5 years ago
Karma
86
Submissions
8
Member since
Nov 17
Recent submissions
Comment · 1 year ago
I wonder how the shared memory model is mapped to this approach. Looking at it now
Comment · 1 year ago
The gist of the post is that optimizations and interpretations thereof must always be made with respect to the underlying hardware.
Comment · 1 year ago
The BLAS GEMM routines I have seen use normal blocked algorithms.
Comment · 1 year ago
Even the integrated Intel HD Graphics would be an interesting comparison.
Comment · 5 years ago
Variations of the Navier-Stokes equations are ubiquitious in many domains in scientific computing.
[Chapter 38. Fast Fluid Dynamics Simulation on the GPU](https://developer.download.nvidia.com/books/HTML/gpugems/gpu...) in GPU Germs
is concerned with solving these equations efficiently on the GPU.
There, the numerical version of Poisson's equation involves solving a massive but sparse linear system,
<img src="https://latex.codecogs.com/svg.image?Ax=b" title="Ax=b"/>
In GPU Germs this system is solved with Jacobi iterations. These work well for many applications but might converge slowly, especially for large matrices. In GPU Germs it is mentioned in passing, that multigrid schemes are a possible way to solve such matrices more efficiently.
It took me some time to understand the concepts but I now have a working GPU accelerated multigrid formulation for the Poisson problem. In some situations, this has accelerated my simulations at times more than twofold. In order to help you avoid making the same mistakes that I did, I published my notes and source code on
[keinefirma.xyz/documents/fluid_tutorial](keinefirma.xyz/documents/fluid_tutorial)
For now, I continuously update these notes and rectify errors and more.
I hope to help somebody with this!