#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int omp_get_num_threads();
int omp_get_thread_num();
double omp_get_wtime();
int /*<<< orphan*/ printf(char *, double);
int main() {
  double startTime = omp_get_wtime();
  long num_steps = 100000;
  double step;
  double pi = 0.0;
  int totalNthreads = omp_get_num_threads();
  step = 1.0 / (double)num_steps;
#pragma omp parallel
  {
    int id = omp_get_thread_num();
    double x = 0.0;
    double sum = 0.0;
    int istart = id * num_steps / totalNthreads;
    int iend = (id + 1) * num_steps / totalNthreads;
    if (id - 1 == totalNthreads)
      iend = num_steps;
    for (int i = istart; i < iend; i++) {
      x = (i + 0.5) * step;
      sum += 4.0 / (1.0 + x * x);
    }
    sum = sum * step;
#pragma omp critical
    pi += sum;
  }
  double endTime = omp_get_wtime();
  printf("PIGRECO => %.16g\n", pi);
  printf("TIME :%.16g\n ", endTime - startTime);
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
movl   $0x0,-0x4(%rbp)
mov    $0x0,%al
call   1040 <omp_get_wtime@plt>
movsd  %xmm0,-0x10(%rbp)
movq   $0x186a0,-0x18(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    $0x0,%al
call   1060 <omp_get_num_threads@plt>
mov    %eax,-0x2c(%rbp)
cvtsi2sdq -0x18(%rbp),%xmm1
movsd  0xe7d(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    $0x0,%al
call   1030 <omp_get_thread_num@plt>
mov    %eax,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movslq -0x30(%rbp),%rax
imul   -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cqto
idiv   %rcx
mov    %eax,-0x44(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
cltq
imul   -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cqto
idiv   %rcx
mov    %eax,-0x48(%rbp)
mov    -0x30(%rbp),%eax
sub    $0x1,%eax
cmp    -0x2c(%rbp),%eax
jne    11ea <main+0xaa>
mov    -0x18(%rbp),%rax
mov    %eax,-0x48(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x48(%rbp),%eax
jge    1255 <main+0x115>
cvtsi2sdl -0x4c(%rbp),%xmm0
movsd  0xe07(%rip),%xmm1
addsd  %xmm1,%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm0
mulsd  %xmm0,%xmm1
movsd  0xdd3(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  0xdcf(%rip),%xmm0
divsd  %xmm1,%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11f0 <main+0xb0>
movsd  -0x40(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
addsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    $0x0,%al
call   1040 <omp_get_wtime@plt>
movsd  %xmm0,-0x58(%rbp)
movsd  -0x28(%rbp),%xmm0
lea    0xd8d(%rip),%rdi
call   1050 <printf@plt>
movsd  -0x58(%rbp),%xmm0
subsd  -0x10(%rbp),%xmm0
lea    0xd89(%rip),%rdi
call   1050 <printf@plt>
xor    %eax,%eax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbx
sub    $0x20,%rsp
xor    %eax,%eax
call   1040 <omp_get_wtime@plt>
movsd  %xmm0,0x18(%rsp)
xor    %eax,%eax
call   1060 <omp_get_num_threads@plt>
mov    %eax,%ebx
xor    %eax,%eax
call   1030 <omp_get_thread_num@plt>
mov    %eax,%esi
cltq
imul   $0x186a0,%rax,%rax
movslq %ebx,%rdi
mov    %rax,%rcx
or     %rdi,%rcx
shr    $0x20,%rcx
je     1186 <main+0x46>
cqto
idiv   %rdi
mov    %rax,%rcx
jmp    118c <main+0x4c>
xor    %edx,%edx
div    %ebx
mov    %eax,%ecx
lea    0x1(%rsi),%eax
cltq
imul   $0x186a0,%rax,%rax
mov    %rax,%rdx
or     %rdi,%rdx
shr    $0x20,%rdx
je     11ab <main+0x6b>
cqto
idiv   %rdi
jmp    11af <main+0x6f>
xor    %edx,%edx
div    %ebx
add    $0xffffffff,%esi
cmp    %ebx,%esi
mov    $0x186a0,%edx
cmovne %eax,%edx
xorpd  %xmm1,%xmm1
mov    %edx,%esi
sub    %ecx,%esi
jle    12a0 <main+0x160>
lea    0x1(%rcx),%eax
test   $0x1,%sil
jne    11e0 <main+0xa0>
xorpd  %xmm0,%xmm0
cmp    %eax,%edx
jne    121d <main+0xdd>
jmp    1290 <main+0x150>
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
addsd  0xe11(%rip),%xmm0
mulsd  0xe11(%rip),%xmm0
mulsd  %xmm0,%xmm0
addsd  0xe0d(%rip),%xmm0
movsd  0xe0d(%rip),%xmm1
divsd  %xmm0,%xmm1
xorpd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
mov    %eax,%ecx
cmp    %eax,%edx
je     1290 <main+0x150>
movsd  0xddb(%rip),%xmm1
movsd  0xddb(%rip),%xmm2
movsd  0xddb(%rip),%xmm3
movsd  0xddb(%rip),%xmm4
nopl   (%rax)
xorps  %xmm5,%xmm5
cvtsi2sd %ecx,%xmm5
addsd  %xmm1,%xmm5
mulsd  %xmm2,%xmm5
mulsd  %xmm5,%xmm5
addsd  %xmm3,%xmm5
movapd %xmm4,%xmm6
divsd  %xmm5,%xmm6
addsd  %xmm0,%xmm6
lea    0x1(%rcx),%eax
xorps  %xmm5,%xmm5
cvtsi2sd %eax,%xmm5
addsd  %xmm1,%xmm5
mulsd  %xmm2,%xmm5
mulsd  %xmm5,%xmm5
addsd  %xmm3,%xmm5
movapd %xmm4,%xmm0
divsd  %xmm5,%xmm0
addsd  %xmm6,%xmm0
add    $0x2,%ecx
cmp    %edx,%ecx
jne    1240 <main+0x100>
mulsd  0xd70(%rip),%xmm0
xorpd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
xor    %eax,%eax
call   1040 <omp_get_wtime@plt>
movsd  %xmm0,0x10(%rsp)
lea    0xd66(%rip),%rdi
movsd  0x8(%rsp),%xmm0
call   1050 <printf@plt>
movsd  0x10(%rsp),%xmm0
subsd  0x18(%rsp),%xmm0
lea    0xd5a(%rip),%rdi
call   1050 <printf@plt>
xor    %eax,%eax
add    $0x20,%rsp
pop    %rbx
ret

LLVM IR (-O0):
@.str = [18 x i8] c"PIGRECO => %.16g\0A\00"
@.str.1 = [14 x i8] c"TIME :%.16g\0A \00"
define i32 @main() {
%1 = alloca i32
%2 = alloca double
%3 = alloca i64
%4 = alloca double
%5 = alloca double
%6 = alloca i32
%7 = alloca i32
%8 = alloca double
%9 = alloca double
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
store i32 0, i32* %1
%14 = call double (...) @omp_get_wtime()
store double %14, double* %2
store i64 100000, i64* %3
store double 0.000000e+00, double* %5
%15 = call i32 (...) @omp_get_num_threads()
store i32 %15, i32* %6
%16 = load i64, i64* %3
%17 = sitofp i64 %16 to double
%18 = fdiv double 1.000000e+00, %17
store double %18, double* %4
%19 = call i32 (...) @omp_get_thread_num()
store i32 %19, i32* %7
store double 0.000000e+00, double* %8
store double 0.000000e+00, double* %9
%20 = load i32, i32* %7
%21 = sext i32 %20 to i64
%22 = load i64, i64* %3
%23 = mul nsw i64 %21, %22
%24 = load i32, i32* %6
%25 = sext i32 %24 to i64
%26 = sdiv i64 %23, %25
%27 = trunc i64 %26 to i32
store i32 %27, i32* %10
%28 = load i32, i32* %7
%29 = add nsw i32 %28, 1
%30 = sext i32 %29 to i64
%31 = load i64, i64* %3
%32 = mul nsw i64 %30, %31
%33 = load i32, i32* %6
%34 = sext i32 %33 to i64
%35 = sdiv i64 %32, %34
%36 = trunc i64 %35 to i32
store i32 %36, i32* %11
%37 = load i32, i32* %7
%38 = sub nsw i32 %37, 1
%39 = load i32, i32* %6
%40 = icmp eq i32 %38, %39
br i1 %40, label %41, label %44
41: ; preds = %0
%42 = load i64, i64* %3
%43 = trunc i64 %42 to i32
store i32 %43, i32* %11
br label %44
44: ; preds = %41, %0
%45 = load i32, i32* %10
store i32 %45, i32* %12
br label %46
46: ; preds = %62, %44
%47 = load i32, i32* %12
%48 = load i32, i32* %11
%49 = icmp slt i32 %47, %48
br i1 %49, label %50, label %65
50: ; preds = %46
%51 = load i32, i32* %12
%52 = sitofp i32 %51 to double
%53 = fadd double %52, 5.000000e-01
%54 = load double, double* %4
%55 = fmul double %53, %54
store double %55, double* %8
%56 = load double, double* %8
%57 = load double, double* %8
%58 = call double @llvm.fmuladd.f64(double %56, double %57, double 1.000000e+00)
%59 = fdiv double 4.000000e+00, %58
%60 = load double, double* %9
%61 = fadd double %60, %59
store double %61, double* %9
br label %62
62: ; preds = %50
%63 = load i32, i32* %12
%64 = add nsw i32 %63, 1
store i32 %64, i32* %12
br label %46
65: ; preds = %46
%66 = load double, double* %9
%67 = load double, double* %4
%68 = fmul double %66, %67
store double %68, double* %9
%69 = load double, double* %9
%70 = load double, double* %5
%71 = fadd double %70, %69
store double %71, double* %5
%72 = call double (...) @omp_get_wtime()
store double %72, double* %13
%73 = load double, double* %5
%74 = call i32 @printf(i8* getelementptr inbounds ([18 x i8], [18 x i8]* @.str, i64 0, i64 0), double %73)
%75 = load double, double* %13
%76 = load double, double* %2
%77 = fsub double %75, %76
%78 = call i32 @printf(i8* getelementptr inbounds ([14 x i8], [14 x i8]* @.str.1, i64 0, i64 0), double %77)
ret i32 0
}

LLVM IR (-O2):
@.str = [18 x i8] c"PIGRECO => %.16g\0A\00"
@.str.1 = [14 x i8] c"TIME :%.16g\0A \00"
define i32 @main() {
%1 = tail call double (...) @omp_get_wtime()
%2 = tail call i32 (...) @omp_get_num_threads()
%3 = tail call i32 (...) @omp_get_thread_num()
%4 = sext i32 %3 to i64
%5 = mul nsw i64 %4, 100000
%6 = sext i32 %2 to i64
%7 = sdiv i64 %5, %6
%8 = trunc i64 %7 to i32
%9 = add nsw i32 %3, 1
%10 = sext i32 %9 to i64
%11 = mul nsw i64 %10, 100000
%12 = sdiv i64 %11, %6
%13 = trunc i64 %12 to i32
%14 = add nsw i32 %3, -1
%15 = icmp eq i32 %14, %2
%16 = select i1 %15, i32 100000, i32 %13
%17 = icmp sgt i32 %16, %8
br i1 %17, label %18, label %40
18: ; preds = %0
%19 = sub i32 %16, %8
%20 = add i32 %8, 1
%21 = and i32 %19, 1
%22 = icmp eq i32 %21, 0
br i1 %22, label %31, label %23
23: ; preds = %18
%24 = sitofp i32 %8 to double
%25 = fadd double %24, 5.000000e-01
%26 = fmul double %25, 1.000000e-05
%27 = tail call double @llvm.fmuladd.f64(double %26, double %26, double 1.000000e+00)
%28 = fdiv double 4.000000e+00, %27
%29 = fadd double %28, 0.000000e+00
%30 = add nsw i32 %8, 1
br label %31
31: ; preds = %23, %18
%32 = phi double [ undef, %18 ], [ %29, %23 ]
%33 = phi i32 [ %8, %18 ], [ %30, %23 ]
%34 = phi double [ 0.000000e+00, %18 ], [ %29, %23 ]
%35 = icmp eq i32 %16, %20
br i1 %35, label %36, label %46
36: ; preds = %46, %31
%37 = phi double [ %32, %31 ], [ %61, %46 ]
%38 = fmul double %37, 1.000000e-05
%39 = fadd double %38, 0.000000e+00
br label %40
40: ; preds = %36, %0
%41 = phi double [ 0.000000e+00, %0 ], [ %39, %36 ]
%42 = tail call double (...) @omp_get_wtime()
%43 = tail call i32 @printf(i8* getelementptr inbounds ([18 x i8], [18 x i8]* @.str, i64 0, i64 0), double %41)
%44 = fsub double %42, %1
%45 = tail call i32 @printf(i8* getelementptr inbounds ([14 x i8], [14 x i8]* @.str.1, i64 0, i64 0), double %44)
ret i32 0
46: ; preds = %31, %46
%47 = phi i32 [ %62, %46 ], [ %33, %31 ]
%48 = phi double [ %61, %46 ], [ %34, %31 ]
%49 = sitofp i32 %47 to double
%50 = fadd double %49, 5.000000e-01
%51 = fmul double %50, 1.000000e-05
%52 = tail call double @llvm.fmuladd.f64(double %51, double %51, double 1.000000e+00)
%53 = fdiv double 4.000000e+00, %52
%54 = fadd double %48, %53
%55 = add nsw i32 %47, 1
%56 = sitofp i32 %55 to double
%57 = fadd double %56, 5.000000e-01
%58 = fmul double %57, 1.000000e-05
%59 = tail call double @llvm.fmuladd.f64(double %58, double %58, double 1.000000e+00)
%60 = fdiv double 4.000000e+00, %59
%61 = fadd double %54, %60
%62 = add nsw i32 %47, 2
%63 = icmp eq i32 %62, %16
br i1 %63, label %36, label %46
}