Thus I am glad to hear that a new conference is going to take place this year for the first time in Utah, June 6-9, 2010:  The First International Meeting on Image Formationin X-Ray Computed Tomography.

OpenCL is being created by the Khronos Group with the participation of many industry-leading companies and institutions including AMD, Apple, Broadcom, IBM, Intel,  NVIDIA and many more. OpenCL aims to be the first royalty-free standard for general-purpose parallel programming of heterogeneous systems. It provides a uniform programming environment for software developers to write efficient, portable code for high-performance compute servers, desktop computer systems and handheld devices using a diverse mix of multi-core CPUs, GPUs, Cell-type architectures and other parallel processors such as DSPs.

Developers can apply to become a GPU Computing Registered Developer at:  www.nvidia.com/opencl. I registered for it. So look forward to some interesting OpenCL code in a future post.

Intel announced already some time ago a new GPU code-named Larrabee. It is likely that Intel targets also the scientific programming community, as the platform has some nice benefits:

1. It is x86 compatible and your code will run with a simple recompilation.
2. Full support for OpenMP or Intel TBB.
3. Enhanced 512-bit vector processing units.

However, yet no hardware prototypes are available but writing code is already possible as Intel  provided some C++ implementation of the new Larrabee instruction set extension called LRBni. First prototypes are expected for the end of 2009.

First Medical Experiences at the Fully3D?

At this years Fully3D conference a high-performance workshop will take place. According to the topic list, Larrabee seems to be of some importance. Maybe we will see first prototype applications in Beijing 2009.

For those who don’t want to wait until then, check out these links for more information:

• First Look at the Larrabee New Instructions (LRBni)
• Intel fact sheet about Larrabee
• Intel’s SIGGRAPH 2008 paper on Larrabee

So use the additional time to fine tune your papers …

• Given: Set of discrete delements E(i). Likelihood P(i) for the occurence of an element E(i).
• Goal: A function which returns a random element from E based on the distribution described by P.

What is this useful for? Well, I wanted to implement some kind of stochastic gradient descent algorithm for image registration without gradients but given a certain likelihood for a point to increase the cost function value. However, I think this kind of function comes in handy for many situations.

Matlab Code & Example

I wrote a Matlab function randelement() which solves exactly the problem. It is based on the theory of inverse transform sampling. The link to the download is given below. It is well-documented, that’s why I will only provide a little usage example:

The histogram will then look somehow similar to this:

So … this function seems to do what we wanted - I like it!
The latest version of randelement.m can be downloaded here.

However, NVIDIA seems to realize that this situation is not amusing to any Laptop-Programmer (like me). So yesterday I found on the download site the following item:

Beta Notebook Driver for Developers (181.22)

So downloaded, installed, and same error as with the other drivers. Well, at least it is a step into the right direction.  Maybe it works for others.

Standard alternative

However, there is still the option to install a modified driver for your notebook. This worked for me but is not the clean solution. Those drivers can be found at http://www.laptopvideo2go.com/

So the quest of today was: Find a fast way for cubic B-spline evaluation using CUDA!

Of course Google helped me out and pointed meto the following papers:

Daniel Ruijters, Bart M. ter Haar Romeny, and Paul Suetens, “ Accuracy of GPU-based B-Spline Evaluation,” In Proc. Tenth IASTED International Conference on Computer Graphics and Imaging (CGIM), Innsbruck, Austria, pp. 117-122, February 13-15, 2008.

Christian Sigg and Markus Hadwiger, “ Fast Third-Order Texture Filtering,” In GPU Gems 2: Programming Techniques for High-Performance Graphics and General-Purpose Computation, Matt Pharr (ed.), Addison-Wesley; chapter 20, pp. 313-329, 2005.

They showed that a fast and accurate evaluation of B-splines is possible by replacing the nearest neighbor lookups during the B-spline evaulation by linear interpolation which is hard-wired on the GPU.

Finally the best of all: Daniel Ruijters provides some clean CUDA code for download on a website. It helped me out very well. In overall I gained a speed-up factor of 7 in comparison to my first naive implementation!