Последовательный и параллельный алгоритм вычисления произведения квадратной матрицы на вектор

Умножение матрицы на вектор
// (поcледовательный алгоритм)
void main(int argc, char* argv[]) {
 double* pMatrix; // Исходная матрица
 double* pVector; // Исходный вектор
 double* pResult; // Вектор-результат
 int Size; // Размер матрицы и векторов
 
 // Инициализация данных
 ProcessInitialization(pMatrix,pVector,pResult,Size);
 // Матрично-векторное умножение
SerialResultCalculation(pMatrix,pVector,pResult,Size);
 // Завершение вычислений
 ProcessTermination(pMatrix,pVector,pResult);
}
// Функция выделения памяти и инициализации данных
void ProcessInitialization (double* &pMatrix,
 double* &pVector, double* &pResult, int &Size) {
 int i;
 do {
 printf("\nВведите размер матрицы и векторов: ");
 scanf("%d", &Size);
 if (Size <= 0) {
 printf("Размер должен быть больше 0! \n ");
}
 } while (Size <= 0);
 pMatrix = new double [Size*Size];
 pVector = new double [Size];
pResult = new double [Size];
for (i=0; i<Size; i++)
 pResult[i] = 0;
 RandomDataInitialization(pMatrix, pVector, Size);
}
// Функция матрично-векторного умножения
void SerialResultCalculation(double* pMatrix,
 double* pVector, double* pResult, int Size) {
 int i, j;
 for (i=0; i<Size; i++) {
 for (j=0; j<Size; j++)
 pResult[i] += pMatrix[i*Size+j]*pVector[j];
 }
}