#pragma omp parallel for
  for (i=0; i<N; ++i) {
    D[i] = sqrt(A[i]) ;
  }


  #pragma omp parallel for private(dx, dy)
  for (i=0; i<N; ++i) {
    dx = AX[i] + BX[i] ;
    dy = AY[i] + BY[i] ;
    D[i] = sqrt(dx*dx + dy*dy) ;
  }

  dp = 0.0 ;
  #pragma omp parallel for reduction(+:dp)
  for (i=0; i<N; ++i) {
    dp = dp + A[i] * B[i] ;
  }