= Matrix Multiply on GPU =
We have implemented single/double precision matrix multiply program for RV770/Cypress. In our implementation, we use two input streams for computing C=AB. One is transposed input matrix A (i.e. column major) and other is input matrix B in normal format (i.e. row major). Output matrix C is also row major. We adopted 8x8 block for single precision and 4x4 for double precision. Here is benchmark result for each case. Note only kernel execution time is measured.

Add double-double precision performance. We used 2x2 block in this case. On Cypress architecture GPU, we take advantage of FMA_64 instruction.

== Peformance Summary ==
|| board  || Pmax || Nmax || prec  || reg. usage || MAD peak ||
|| HD4850 || 736    || 3328 || SP || 25 || 1040 ||
|| HD5870 || 2140   || 7424 || SP || 25 || 2720 ||
|| HD4850 || 177    || 1408 || DP || 19 || 208 ||
|| HD5870 || 475    || 2048 || DP || 19 || 544 ||
|| HD4850 || 7.5    || 768  || DDP || 21 ||  ||
|| HD5870 || 20     || 1024 || DDP || 21 ||  ||
|| HD5870 FMA || 31     || 1024 || DDP || 18 ||   ||

Pmax & MAD in GFLOPS

== Source code ==
Will be posted later.

== Single precision ==
[[Image(SMM.png)]]

== Double precision ==
[[Image(DMM.png)]]

== Double-Double precision ==
[[Image(DDMM.png)]]


= Useful forum discussions =
== Discussion on a highly optimized MM kernel ==
http://forum.beyond3d.com/showthread.php?t=54842

== Discussion on MM kernels in OpenCL ==
http://forums.amd.com/devforum/messageview.cfm?catid=390&threadid=127963

=== IL code generator in C++ ===
CAL++ http://sourceforge.net/projects/calpp/

Meta-programing works in reality. Impressive work!


[[Image(MM1.png)]]