…
|
||
---|---|---|
.. | ||
experiments | ||
src | ||
README.md |
README for the matrix-vector multiplication demo code
Synopsis
This program implements the multiplication of a matrix and a vector. It is written in C and has been parallelized using the Pthreads parallel programming model. Each thread gets assigned a contiguous set of rows of the matrix to work on and the results are stored in the output vector.
The code initializes the data, executes the matrix-vector multiplication, and checks the correctness of the results. In case of an error, a message to this extent is printed and the program aborts. Otherwise it prints a one line message on the screen.
About this code
This is a standalone code, not a library. It is meant as a simple example to experiment with gprofng.
Directory structure
There are four directories:
-
bindir
- after the build, it contains the executable. -
experiments
- after the installation, it contains the executable and also has an example profiling script calledprofile.sh
. -
objects
- after the build, it contains the object files. -
src
- contains the source code and the make file to build, install, and check correct functioning of the executable.
Code internals
This is the main execution flow:
- Parse the user options.
- Compute the internal settings for the algorithm.
- Initialize the data and compute the reference results needed for the correctness check.
- Create and execute the threads. Each thread performs the matrix-vector multiplication on a pre-determined set of rows.
- Verify the results are correct.
- Print statistics and release the allocated memory.
Installation
The Makefile in the src
subdirectory can be used to build, install and check the
code.
Use make
at the command line to (re)build the executable called mxv-pthreads
. It will be
stored in the directory bindir
:
$ make
gcc -o ../objects/main.o -c -g -O -Wall -Werror=undef -Wstrict-prototypes main.c
gcc -o ../objects/manage_data.o -c -g -O -Wall -Werror=undef -Wstrict-prototypes manage_data.c
gcc -o ../objects/workload.o -c -g -O -Wall -Werror=undef -Wstrict-prototypes workload.c
gcc -o ../objects/mxv.o -c -g -O -Wall -Werror=undef -Wstrict-prototypes mxv.c
gcc -o ../bindir/mxv-pthreads ../objects/main.o ../objects/manage_data.o ../objects/workload.o ../objects/mxv.o -lm -lpthread
ldd ../bindir/mxv-pthreads
linux-vdso.so.1 (0x0000ffff9ea8b000)
libm.so.6 => /lib64/libm.so.6 (0x0000ffff9e9ad000)
libc.so.6 => /lib64/libc.so.6 (0x0000ffff9e7ff000)
/lib/ld-linux-aarch64.so.1 (0x0000ffff9ea4e000)
$
The make install
command installs the executable in directory experiments
.
$ make install
Installed mxv-pthreads in ../experiments
$
The make check
command may be used to verify the program works as expected:
$ make check
Running mxv-pthreads in ../experiments
mxv: error check passed - rows = 1000 columns = 1500 threads = 2
$
The make clean
comand removes the object files from the objects
directory
and the executable from the bindir
directory.
The make veryclean
command implies make clean
, but also removes the
executable from directory experiments
.
Usage
The code takes several options, but all have a default value. If the code is
executed without any options, these defaults will be used. To get an overview of
all the options supported, and the defaults, use the -h
option:
$ ./mxv-pthreads -h
Usage: ./mxv-pthreads [-m <number of rows>] [-n <number of columns] [-r <repeat count>] [-t <number of threads] [-v] [-h]
-m - number of rows, default = 2000
-n - number of columns, default = 3000
-r - the number of times the algorithm is repeatedly executed, default = 200
-t - the number of threads used, default = 1
-v - enable verbose mode, off by default
-h - print this usage overview and exit
$
For more extensive run time diagnostic messages use the -v
option.
As an example, these are the options to compute the product of a 2000x1000 matrix with a vector of length 1000 and use 4 threads. Verbose mode has been enabled:
$ ./mxv-pthreads -m 2000 -n 1000 -t 4 -v
Verbose mode enabled
Allocated data structures
Initialized matrix and vectors
Defined workload distribution
Assigned work to threads
Thread 0 has been created
Thread 1 has been created
Thread 2 has been created
Thread 3 has been created
Matrix vector multiplication has completed
Verify correctness of result
Error check passed
mxv: error check passed - rows = 2000 columns = 1000 threads = 4
$
Executing the examples
Directory experiments
contains the profile.sh
script. This script
checks if gprofng can be found and for the executable to be installed.
The script will then run a data collection experiment, followed by a series
of invocations of gprofng display text
to show various views. The results
are printed on stdout.
To include the commands executed in the output of the script, and store the
results in a file called LOG
, execute the script as follows:
$ bash -x ./profile.sh >& LOG
Additional comments
-
The reason that compiler based inlining is disabled is to make the call tree look more interesting. For the same reason, the core multiplication function
mxv_core
has inlining disabled through thevoid __attribute__ ((noinline))
attribute. Of course you're free to change this. It certainly does not affect the workings of the code. -
This distribution includes a script called
profile.sh
. It is in theexperiments
directory and meant as an example for (new) users of gprofng. It can be used to produce profiles at the command line. It is also suitable as a starting point to develop your own profiling script(s).