%----------------------------------------------------------------------------------------------
% 			UE - PP - Matrix Parallel Programming
%----------------------------------------------------------------------------------------------
\documentclass[a4paper,11pt]{article}

%\usepackage[portuguese]{babel}
\usepackage[utf8]{inputenc}
\usepackage[T1]{fontenc}
\usepackage{graphicx}
\usepackage{float}
\usepackage{amsmath}
\usepackage{geometry}
\usepackage{listings}
\usepackage{xcolor}

\definecolor{codegreen}{rgb}{0,0.6,0}
\definecolor{codegray}{rgb}{0.5,0.5,0.5}
\definecolor{codepurple}{rgb}{0.58,0,0.82}
\definecolor{backcolour}{rgb}{0.95,0.95,0.92}

\lstdefinestyle{mystyle}{
    backgroundcolor=\color{backcolour},   
    commentstyle=\color{codegreen},
    keywordstyle=\color{magenta},
    numberstyle=\tiny\color{codegray},
    stringstyle=\color{codepurple},
    basicstyle=\ttfamily\footnotesize,
    breakatwhitespace=false,         
    breaklines=true,                 
    captionpos=b,                    
    keepspaces=true,                 
    numbers=left,                    
    numbersep=2pt,                  
    showspaces=false,                
    showstringspaces=false,
    showtabs=false,                  
    tabsize=1
}

\lstset{style=mystyle}
\geometry{top=2.5cm, bottom=2.5cm, left=2.5cm, right=2.5cm}
\usepackage{fancyhdr}
\pagestyle{fancy}
\fancyhf{} 

\lhead{Programação Paralela}
\rhead{Matrix Search Parallelization}
\fancypagestyle{noheadrule}{
  \fancyhf{}
  \renewcommand{\headrulewidth}{0pt}
}
\lfoot{Universidade de Évora}
\rfoot{\thepage}


\title{Matrix Search Parallelization}
\author{Vitor Gonçalo Costa, m70323}
\date{Universidade de Évora\\ \today}
\begin{document}
\maketitle

\abstract{This study aims to investigate the performance gain of a program when applying parallel programming techniques to divide operations between several threads using POSIX framework. }

\section*{Presenting the Problem}
We want to write a program that finds the biggest element in a \textbf{square matrix} with thousands of rows and columns, built during run time. We want to evaluate the metrics \underline{speedup}, \underline{efficiency}, and \underline{scalability} when converting the program from sequencial to parallel computation and we should analyse how does the amount of processors and threads influences the program execution time.\\\\The Matrix is built with the logic: $matrix[r][c] = r*m\_cols+c$, inside a nested for loop for rows and columns, where r and c represent the iteration variable for row and column respectively. This way, the output of $matrix[nrows][ncolumns]$ would always be the highest element in the matrix.

\subsection*{Resources Available}
The machine where we will be compiling and executing the program is a MacBook Pro with an i7 dual-core processor of 2,5 GHz - This cpu has available 2 physical cores with 4 virtual cores each, which means we might be able to invoke up to 8 threads of computation.

\subsection*{The Sequencial Program}
\begin{figure}[h]
\begin{minipage}[c]{0.45\textwidth}
	\centering
	\includegraphics[width=0.925\textwidth]{imagens/matrixsearchsequencial.png}  
	%\lstinputlisting[language=C]{SequencialMatrixSearch2.c}
\end{minipage}
\begin{minipage}[c]{0.5\textwidth}
\small

The execution of the program on the right relies on a lot of memory available in order to compute such big matrix. If we pick the value of 40000 for the amount of rows and columns, and knowing that an integer requires 4 bytes of memory, we need around 6,4 GB to compute such matrix. This machine has 16 GB of LPDDR3 installed.\\
 
The current implementation takes around 5,45 seconds to be executed in this machine, and throws a bus error if m\_rows is updated to 40000, which is not expected at first glance, since the system has 16GB of memory available. \\
	
This might happen due to OS imposed limits, so we allocate resources with POSIX library. 
\end{minipage}
\end{figure}
	
We update the above program to make use of posix\_memalign() function in order to align the matrix memory with the cpu's 64 bytes cache lines. This technique is called \textit{dynamic memory assignment}, is usually applied to very large matrices and a requirement for \textit{High-Performance Computation} programs.

The program was executed three times using the time command, and the average execution time was computed:

\begin{table}[htp]
\begin{center}
\begin{tabular}{ccc}
m\_rows & no HPC (sec) & with HPC (sec)\\
\hline
30k & 5,482 & 5,395 \\
40k & error & 9.286 \\
\end{tabular}
\end{center}
\caption{program's total time execution per amount of matrix rows.}
\label{sequencial program time results}
\end{table}%

We could now use the \textbf{time.h} library to evaluate the time spent on each operation, 1st allocate memory, 2nd building the matrix and 3rd search for the biggest value, we perform this exercise in next section. 

\section*{Parallel Strategies}
From the three tasks presented previously, the search for the biggest value is the one that we will parallelise first, therefore, a certain number of threads are spawned to split the task among them with the usage of pthread\_create() and pthread\_join() functions from \textbf{pthread.h} library.

With the current implementation of our program, the global variable hv, used to store the highest value found during the search operation, will be shared between threads by the time we execute a parallelised version of it. We need to take into account this \textit{race condition} and decide on a strategy to overcome it. The usage of a pre variable lock() and after unlock() \textbf{muttex} pattern will stop the execution of threads waiting for the access to that variable, due to the nature of our problem, we update this variable to be \textbf{atomic int} instead, but this change alone actually increases the computation time up to 18,35 seconds with $m\_rows = 40k$, which again, suggests the use of \textbf{time.h} library to better understand the contribution of each parallel strategy presented in the following sections of this document. The next picture displays the time our program toke on each task after applying the implementations mentioned before:

 \begin{figure}[h]
	\centering
	\includegraphics[width=0.6\textwidth]{imagens/matrixsearchsequencial-02.png}  
\end{figure}

From the previous image we can be sure the 3rd task is the one to parallelise. We should also notice that the changes performed to the program actually increase it's time of execution, the view operation performed on this type of variable is more resource consuming than when performed on a variable of type integer. This can be proved by analysing the time for each task on the first version of our program, populating the matrix toke around 4.5 seconds while the search matrix task toke less than 4 seconds for 40k rows. The wall-clock time of the above figure is 18,931 seconds, this representes the real time a users waits for a program to finish,

\subsection*{Columns Split}
Because the matrix in the spotlight is a square one, at first approach, one might think it doesn't really matter if a split is performed by columns or rows, but it does. Since C approaches 2D arrays per row, by the time the cpu access the 32 or 64 bytes of the matrix memory stored in ram to make it available closer, at LP2 or LP3, it brings the closest row values also, making it faster to loop though rows instead of columns. 

Again, one might thinks that what matters to the performance of such program, is the speed one of the about to be spawned threads takes to find the highest value present in the matrix. But the amount of memory access during a program execution will definitely affect it's performance. For such big matrix, with this many elements, for sure the way we loop through it will affect performance. %This way, threads will only perform read on hv global var and not write, which is a more time consuming computer operation.

A new function, parallel\_search\_cols(void *args), was defined in order to share the work between threads, it expects the start and end column indices as arguments in order to split the matrix in as many slots as needed. A second function, run\_parallel\_search(int num\_threads), was also defined to make the program more useful by allowing to define the amount of threads to be created. 

The best execution times obtained on this scenario were with the creation of 600 to 800 threads, presenting a total cpu computing time of around 8,5 seconds on the 3rd task. The next figure computes the \underline{sum up time taken by all threads} per its cardinality for both strategies presented in this document. 

\begin{figure}[h]
	\centering
	\includegraphics[width=0.6\textwidth]{imagens/matrixsearchsequencial-03.png}
	\caption{Time taken by all threads on search matrix task.}
\end{figure}

A performance analysis is computed later this document, but we shall advance that the wall-clock time for 800 threads is around 8,21 seconds (around 1 millisecond per thread), which is better than the sequencial program in table 1. The next figure plots the amount of user time taken by the columns parallel solution implemented, it's more comum to observe then the total cpu time.

\begin{figure}[h]
	\centering
	\includegraphics[width=0.6\textwidth]{imagens/matrixsearchsequencial-04.png}
	\caption{Total time taken by program with columns parallelisation.}
\end{figure}


\subsection*{Blocks Split}
We updated the function run\_parallel\_search (int num\_threads), in order to divide the matrix in blocks of the same size. First it splits the matrix in two and calculate the amount of blocks for each side by the square root of the amount of blocks to be created. If we pick 400 threads for example, means we have 40 blocks to be worked out during the 3rd task. Each of this blocks has 2000 rows and columns, $\sqrt{400} = 20$,which means 20 blocks for each side of the first split, and then we can validate the amount of items in the blocks with: $20*(2000*2000)*2=160k$. 

We also defined a new struct object to carry the start/end row and column for each block assigned to a created thread, this helps to define the block's rows and columns boundaries. The function parallel\_search\_block() receives this list of arguments by the time is invoked by pthread\_create() and performs the comparison of matrix index with the atomic int hv variable.

From figure 1 we can advance that, in general, both strategies have similar performances, but the second solution appears to scale better than the first one. Some of the spikes present in the graph for this solution could be explained by the size of sub-matrixes to assign on threads. The amount of threads used in those cases do not divide the matrix in perfectly square sub-matrixes, all threads access the atomic hv var on each iteration, therefor some weird behaviours could happen when trying to access it's address.

Next we plot the total time taken by the program with this parallel solution.

\begin{figure}[h]
	\centering
	\includegraphics[width=0.6\textwidth]{imagens/matrixsearchsequencial-05.png}
	\caption{Total time taken by program with blocks parallelisation.}
\end{figure}

\subsubsection*{Performance Analysis}
In order to evaluate the performance of our program when parallelised with both of previous strategies we compute the wall-clock cpu time obtained per threads, available in table 2. 

The \textbf{speedup} metric is determined by the quotient between the best time obtained by a sequential program and a parallelised version of it: $Speedup = \frac{t^*}{t_p} $. 

The \textbf{efficiency} of a parallel version of a program can then be calculated by the fraction $efficiency = \frac{speedup}{p}$.

Finally, the \textbf{scalability} of a program measures how well a parallel program handles increasing levels of parallelism. %It is calculated by the formula: $E = \frac{S_p}{p}$.

\begin{table}[htp]
\begin{center}
\begin{tabular}{ccc}
\# Threads & Columns & Blocks \\
\hline
4 & 16.627 & 16.626\\
8 & 12.812 & 17.125\\
16 & 10.572 & 10.372\\
32 & 9.265 & 9.466 \\
49 & 8.664 & 8.422 \\
64 & 8.28 &  8.520 \\
400 & 7.749 & 10.884 \\
800 & 7.916 & 10.613 \\
1000 & 8.126 & 10.560 
\end{tabular}
\end{center}
\caption{Program user time for each strategy per threads, for matrix with 40k rows and columns.}
\label{wall-cpu time of parallel strategies}
\end{table}%

Combining values from tables 1 and 2, and using the formulas presented before we can compute the performance metrics for both strategies:

\begin{table}[htp]
\begin{center}
\begin{tabular}{cccc}
Strategy & \#Threads & Speedup & Efficiency (\%)\\
\hline
Columns & 400 & 1.1467 & 0.007\\
Blocks & 400 & 0.8164 & 0.2 \\
Columns & 49 & 1.026 & 2.1 \\
Blocks & 49 & 1.0526 & 2.1
\end{tabular}
\end{center}
\caption{Performance of parallel strategies.}
\label{performance of parallel strategies}
\end{table}%

We can conclude that the parallelisation of our initial program didn't got a significant performance gain. Probably due to the nature of our problem, since the biggest value is always in the last entry of the matrix, threads will have to loop through the entire batch of items to compare values. But also because of the implementation applied, as mentioned before, the function run\_parallel\_search (int num\_threads) could be improved in order to access the atomic variable only once per thread. 

The metric scalability can be divided in two: \underline{strong scaling}, when the objective is to evaluate how faster a program gets as resources (or threads) are added to a fixed problem size - we already computed this analysis in table 2. There, we can see that the strategy to parallelise the search matrix task by columns has a better performance with the increase of resources. The blocks strategy algorithm loses performance sooner, and if we look at figure 1, we can notice that it needs much more resources to achieve similar total cpu time as the columns split. The \underline{weak scalability} is measured to evaluate if a program can handle larger problems in the same amount of time by adding more resources, ideally, the execution time should remain constant as both work and threads are scaled. From the next plot we can conclude that the columns split is a better strategy for the peculiar situation in study.

\begin{figure}[h]
	\centering
	\includegraphics[width=0.6\textwidth]{imagens/matrixsearchsequencial-06.png}
	\caption{Weak scalability evaluation from 1 to 8 threads.}
\end{figure}

Since there wasn't a big gain when parallelising the search task we could either try to develop a new parallel strategy for the search task or parallelise the matrix initialisation instead. 

\section*{Improving performance}
tbd

\section*{Conclusions}
Modern CPUs can spawn thousands of virtual threads to perform calculations.\\

The action of populating the matrix could also update the value of global variable hv instead of the 3rd task. Search the matrix could have a local variable to store the highest value found and compare with the global var hv at the end.\\

The clock() function from time.h library is useless to compute wall-clock time and therefore should be ignored on metrics evaluations.


TODO:
- strategy to split and continue threads;
- evaluate strong/weak scalability with plot per number of threads 


\end{document}