The purpose of this article is to show how to speed up some matrix operations. This improvement is obtained by using assembly language with SSE (Streaming SIMD Extension) technology.
In this article, we are going to see how to improve the speed of the mathematical operation: Matrix Vector multiplication. I have chosen this operation because it is often used in programming (image processing, 2D, 3D etc.).
Since the Pentium III Processor, Intel encloses the SSE technology. This technology was included in order to improve image, video, audio processing. The SSE technology uses SIMD (Single Instruction Multiple Data) instructions. For example, this allows to perform four multiplications (on single precision float) with only one instruction. By using these characteristics, I tried to speed up the matrix/vector multiplication.
In order to use the SSE instruction set, you need a computer with a Pentium III processor or newer. Next, Visual C++ can't compile the SSE instructions without the Processor Pack. (Instructions to install it in the Zip file.)
Using the code
The program provided by the link on the top performs a matrix/vector multiplication. It displays the time spent in the C++ function and the time spent in the assembly function. It also displays the matrix and the two vectors (multiplication and result).
You can find two ways to proceed this operation (one in C++ and another in assembler). The code shows you the time (in CPU clock cycles) spent in each function.
First in the C++ function, the code is optimized when you compile in Release mode. So the matching must be done in the Release mode.
Here is the
main of the project:
int main(int argc, char* argv)
printf("Enter the size of the matrix:\n");
float* matrix=(float*) malloc (size*size*sizeof(float));
float* vector=(float*) malloc (size*sizeof(float));
float* result=(float*) malloc (size*sizeof(float));
float* matrix1=(float*) malloc (size*size*sizeof(float));
for(i=0; i<size*size; i++)
MatrixVectorWritting(matrix, vector, size);
for(i=0; i<100; i++)
printf("Time spend en C++ fonction: %d clock cycles.\n",
for(i=0; i<100; i++)
result=MatrixVector_SSE(matrix1, vector, size);
printf("Time spend en Asm SSE fonction:"
" %d clock cycles.\n",time_SSE);
MatrixVectorDisplay(matrix, vector, result, size);
We can see in the
main function that we make a "
for" loop (100 times) on each function we compare. At the beginning of each loop we read the CPU clock. And at the end of each loop we read again the CPU timer clock. Then we take the difference between the two times read. And we obtain the number of clock cycles spent in each function (independent of the CPU used for the test).
Here is the code used to read the CPU timer called RDTSC:
<PRE lang=mc++>/// Read the computer's timer RDTSC
rdtsc // Resad the RDTSC Timer
mov dword ptr[clock], eax // Store the value in EAX and EDX Registers
mov dword ptr[clock+4], edx
Points of Interest
As you can see with the program, there are some matrix dimensions where the assembly function is not very efficient. This happens when the size of the matrix is lower than 3x3. For very low sizes of the matrix the code is not very optimized. But it's possible to make a specific function for the matrix/vector multiplication of 3x3 or 2x2 matrixes. You can find a white paper on the Intel web site.
Some sections of the code can be improved. So if I find anything else I'll update the code. If you have some ideas, you can propose.