void multiply_matrices_IKJ()
{
	for (int i = 0; i < ROWS; i++)
		for (int k = 0; k < COLUMNS; k++)
			for (int j = 0; j < COLUMNS; j++)
				matrix_r[i][j] += matrix_a[i][k] * matrix_b[k][j];

}

void multiply_matrices_IKJ_ATOMIC()
{
#pragma omp parallel for
	for (int i = 0; i < ROWS; i++)
		for (int k = 0; k < COLUMNS; k++)
			for (int j = 0; j < COLUMNS; j++)
			{
#pragma omp atomic
				matrix_r[i][j] += matrix_a[i][k] * matrix_b[k][j];
			}
}