#include <stdio.h>
#include <stdlib.h>

#include "omp.h"

#define NELEMENTS 1024*16

int global_NumThreads;

float vecA[NELEMENTS];
float vecB[NELEMENTS];
float vecC[NELEMENTS];

float dp = 0.0;

int main ()
{

  #pragma omp parallel for schedule(static)
  for(int i=0; i < NELEMENTS; i ++) {
    vecA[i] = 1.0 ;
    vecB[i] = 2.0 ;
  }

  #pragma omp parallel for schedule(static)
  for(int i=0; i < NELEMENTS; i ++) {
    vecC[i] = vecA[i] * vecB[i];
  }

  // main nato sesteje vse delne produkte:

  for(int i=0; i < NELEMENTS; i ++) {
    dp += vecC[i];
  }

  printf("Dot product = %f \n", dp);
  
  return 0;
}