/* 

sum2.c 

compute sum = atan(1) + atan(2) + ... + atan(n)

 */

#include <stdio.h> 
#include <assert.h> 
#include <math.h> 

#include <omp.h> 


int main(int argc, char** argv) {

  long n  = 100000000; 
  double  sum; 
  long  i; 
  double t1, t2; 
  
  t1 = omp_get_wtime(); 

  //initialize sum 
  sum = 0; 

#pragma omp parallel 
  { 
    int thread_id = omp_get_thread_num(); 
    printf("Start(%d)\n", thread_id); 
    
    for (i=0; i<n; i++) {
      sum = sum + atan(i); 
    }
    printf("Done(%d)\n", thread_id); 
  } //end omp 

  t2 = omp_get_wtime();   
  
  printf("sum=%.1f ", sum);
  printf("Done. time = %.1f\n",  t2-t1);
  
  return 0; 
}