Permalink
Find file
Fetching contributors…
Cannot retrieve contributors at this time
26 lines (20 sloc) 1.24 KB
Magnetic Pendulum, OpenCL version
=================================
This is an implementation of a magnetic pendulum in OpenCL. A spherical
pendulum is influenced by magnets. For every point of origin the program
finds the magnet the pendulum's body will be nearest when it comes to rest
(friction is taken into account). It will create an image, where every pixel
(points in the $\phi$-$\theta$-plane) is mapped to a color representing
each magnet.
The problem is ideally suited for parallel computing as every initial
condition can be solved independently from the others. I did a Python version
some time ago that utilized Scipy. Although it was quite optimized with C code,
it performed very bad. To create a map of 1600x1600 pixels, it took about 2700
minutes! Leveraging the power hidden in GPUs, I wanted to know what speed-up
would be reached.
The code works quite good so far. I began coding a CPU version in C that is a
lot faster than the original Python version, but even my decent GPU out-
performs that by a factor of ten. The map from above can now be done in 284
seconds. That is seconds, not minutes! :-)
However, I think there is still plenty of room for optimization and a
real-time visualization of the mapping process would be really cool.