параллельный алгоритм умножения матриц С++

[Kurvochka]

Добрый день всем! Надеюсь что не откажете в помощи) Нужно написать программу параллельного умножения матриц. Ниже приведен примерный код выполнения только он не компилируется. Может кто-то глянет к этому коду или буду рада другому рабочему варианту параллельной программы.
Заранее благодарна!
Листинг программы с++

void main(int argc, char* argv[])
{
double* pAMatrix; // Матрица А
double* pBMatrix; // Матрица B
double* pCMatrix; // Матрица С
int Size; // Размер матриц
// Инициализация данных
ProcessInitialization(pAMatrix, pBMatrix, pCMatrix, Size);
void ProcessInitialization (double* &pAMatrix, double* &pBMatrix, double* &pCMatrix, int &Size)
{
int i, j;
do
{
printf("\nВведите размер матриц: ");
scanf("%d", &Size);
if (Size <= 0)
{
printf("Размер должен быть больше 0! \n ");
} } while (Size <= 0);
pAMatrix = new double [Size*Size];
pBMatrix = new double [Size*Size];
pCMatrix = new double [Size*Size];
for (i=0; i<Size; i++) for (j=0; j<Size; j++) pCMatrix[i*Size+j] = 0;
RandomDataInitialization(pAMatrix, pBMatrix, Size);
}
// Умножение матриц
ParallelResultCalculation(pAMatrix, pBMatrix, pCMatrix, Size);
void ParallelResultCalculation (double* pAMatrix, double* pBMatrix, double* pCMatrix, int Size)
{
int i, j, k;
#pragma omp paralell for private (j, k) for (i=0; i<Size; i++) for (j=0;
j<Size; j++) for (k=0; k<Size; k++) pCMatrix[i*Size+j] += pAMatrix[i*Size+k]* pBMatrix[k*Size+j];
}
// Завершение вычислений
ProcessTermination(pAMatrix, pBMatrix, pCMatrix, Size);
}