#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 */
typedef int int64_t;

/* Variables and functions */
int omp_get_max_threads();
int omp_get_thread_num();
void voodoo(double *a, double *b, double *c, int nelem, int device,
            int offload) {
  int mthreads;
#pragma offload target(mic : device) if (offload)
  { mthreads = omp_get_max_threads(); }
  const int64_t nworkers = nelem > mthreads ? mthreads : 1;
  const int64_t work_split = nelem / nworkers;
  const int64_t work_spill = nelem % nworkers;
#pragma offload target(mic                                                     \
                       : device) in(a                                          \
                                    : length(0) alloc_if(0) free_if(0)),       \
    in(b                                                                       \
       : length(0) alloc_if(0) free_if(0)),                                    \
    in(c                                                                       \
       : length(0) alloc_if(0)                                                 \
           free_if(0))in(work_split)in(work_spill) if (offload)
  {
#pragma omp parallel num_threads(nworkers)
    {
      for (int i = 0; i < 100; ++i) {
        const int tid = omp_get_thread_num();
        int64_t work = 0, work_offset = 0, work_end = 0;
        if (tid < work_spill) {
          work = work_split + 1;
          work_offset = tid * work;
        } else {
          work = work_split;
          work_offset = tid * work + work_spill;
        }
        work_end = work_offset + work;
        if (work) {
          for (int eidx = work_offset; eidx < work_end; ++eidx) {
            c[eidx] = ((a[eidx] + b[eidx]) / 100.0);
          }
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<voodoo>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
mov    $0x0,%al
call   1040 <omp_get_max_threads@plt>
mov    %eax,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x28(%rbp),%eax
jle    1160 <voodoo+0x40>
mov    -0x28(%rbp),%eax
mov    %eax,-0x50(%rbp)
jmp    116d <voodoo+0x4d>
mov    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    116d <voodoo+0x4d>
mov    -0x50(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x1c(%rbp),%eax
cltd
idivl  -0x2c(%rbp)
mov    %eax,-0x30(%rbp)
mov    -0x1c(%rbp),%eax
cltd
idivl  -0x2c(%rbp)
mov    %edx,-0x34(%rbp)
movl   $0x0,-0x38(%rbp)
cmpl   $0x64,-0x38(%rbp)
jge    126c <voodoo+0x14c>
mov    $0x0,%al
call   1030 <omp_get_thread_num@plt>
mov    %eax,-0x3c(%rbp)
movl   $0x0,-0x40(%rbp)
movl   $0x0,-0x44(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    11db <voodoo+0xbb>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x3c(%rbp),%eax
imul   -0x40(%rbp),%eax
mov    %eax,-0x44(%rbp)
jmp    11ee <voodoo+0xce>
mov    -0x30(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x3c(%rbp),%eax
imul   -0x40(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x48(%rbp)
cmpl   $0x0,-0x40(%rbp)
je     1259 <voodoo+0x139>
mov    -0x44(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x48(%rbp),%eax
jge    1254 <voodoo+0x134>
mov    -0x8(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  0xdcb(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1207 <voodoo+0xe7>
jmp    1259 <voodoo+0x139>
jmp    125e <voodoo+0x13e>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    118e <voodoo+0x6e>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<voodoo>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %ecx,%ebp
mov    %rdx,%r14
mov    %rsi,%r15
mov    %rdi,%r13
xor    %eax,%eax
call   1040 <omp_get_max_threads@plt>
cmp    %ebp,%eax
mov    $0x1,%ecx
cmovl  %eax,%ecx
mov    %ebp,%eax
cltd
idiv   %ecx
mov    %edx,%ebx
mov    %eax,%ebp
lea    0x1(%rbp),%eax
mov    %eax,0xc(%rsp)
lea    0x10(%r14),%rax
mov    %rax,0x20(%rsp)
lea    0x10(%r15),%rax
mov    %rax,0x18(%rsp)
lea    0x10(%r13),%rax
mov    %rax,0x10(%rsp)
xor    %r12d,%r12d
mov    %rbp,(%rsp)
jmp    118e <voodoo+0x6e>
xchg   %ax,%ax
add    $0x1,%r12d
cmp    $0x64,%r12d
je     133f <voodoo+0x21f>
xor    %eax,%eax
call   1030 <omp_get_thread_num@plt>
movsd  0xe73(%rip),%xmm2
mov    %eax,%esi
mov    0xc(%rsp),%edi
imul   %edi,%esi
mov    %eax,%ecx
imul   %ebp,%ecx
add    %ebx,%ecx
cmp    %ebx,%eax
mov    %ebp,%edx
cmovl  %edi,%edx
cmovl  %esi,%ecx
test   %edx,%edx
jle    1180 <voodoo+0x60>
add    %ecx,%edx
movslq %ecx,%rax
movslq %edx,%rcx
lea    0x1(%rax),%rdx
cmp    %rcx,%rdx
mov    %rcx,%r10
cmovg  %rdx,%r10
sub    %rax,%r10
cmp    $0x2,%r10
jb     1310 <voodoo+0x1f0>
mov    %ebx,%r11d
lea    (%r14,%rax,8),%rsi
cmp    %rcx,%rdx
cmovle %rcx,%rdx
lea    (%r14,%rdx,8),%rdi
lea    0x0(,%rax,8),%rbp
add    %r13,%rbp
lea    0x0(,%rdx,8),%rbx
add    %r13,%rbx
lea    (%r15,%rax,8),%r8
lea    (%r15,%rdx,8),%rdx
cmp    %rbx,%rsi
setb   %r9b
cmp    %rdi,%rbp
setb   %bl
cmp    %rdx,%rsi
setb   %dl
cmp    %rdi,%r8
setb   %sil
test   %bl,%r9b
movapd 0xdcd(%rip),%xmm3
jne    1309 <voodoo+0x1e9>
and    %sil,%dl
jne    1309 <voodoo+0x1e9>
mov    %r10,%r8
and    $0xfffffffffffffffe,%r8
lea    -0x2(%r8),%rdx
mov    %rdx,%r9
shr    %r9
add    $0x1,%r9
test   %rdx,%rdx
je     1335 <voodoo+0x215>
mov    %r9,%rdi
and    $0xfffffffffffffffe,%rdi
mov    0x20(%rsp),%rdx
lea    (%rdx,%rax,8),%rbp
mov    0x18(%rsp),%rdx
lea    (%rdx,%rax,8),%rdx
mov    0x10(%rsp),%rsi
lea    (%rsi,%rax,8),%rbx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd -0x10(%rbx,%rsi,8),%xmm0
movupd -0x10(%rdx,%rsi,8),%xmm1
addpd  %xmm0,%xmm1
divpd  %xmm3,%xmm1
movupd %xmm1,-0x10(%rbp,%rsi,8)
movupd (%rbx,%rsi,8),%xmm0
movupd (%rdx,%rsi,8),%xmm1
addpd  %xmm0,%xmm1
divpd  %xmm3,%xmm1
movupd %xmm1,0x0(%rbp,%rsi,8)
add    $0x4,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1290 <voodoo+0x170>
test   $0x1,%r9b
je     12f4 <voodoo+0x1d4>
add    %rax,%rsi
movupd 0x0(%r13,%rsi,8),%xmm0
movupd (%r15,%rsi,8),%xmm1
addpd  %xmm0,%xmm1
divpd  0xd12(%rip),%xmm1
movupd %xmm1,(%r14,%rsi,8)
cmp    %r8,%r10
mov    %r11d,%ebx
mov    (%rsp),%rbp
je     1180 <voodoo+0x60>
add    %r8,%rax
jmp    1310 <voodoo+0x1f0>
mov    %r11d,%ebx
mov    (%rsp),%rbp
movsd  0x0(%r13,%rax,8),%xmm0
addsd  (%r15,%rax,8),%xmm0
divsd  %xmm2,%xmm0
movsd  %xmm0,(%r14,%rax,8)
add    $0x1,%rax
cmp    %rcx,%rax
jl     1310 <voodoo+0x1f0>
jmp    1180 <voodoo+0x60>
xor    %esi,%esi
test   $0x1,%r9b
jne    12d2 <voodoo+0x1b2>
jmp    12f4 <voodoo+0x1d4>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @voodoo(double* %0, double* %1, double* %2, i32 %3, i32 %4, i32 %5) {
%7 = alloca double*
%8 = alloca double*
%9 = alloca double*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca i32
store double* %0, double** %7
store double* %1, double** %8
store double* %2, double** %9
store i32 %3, i32* %10
store i32 %4, i32* %11
store i32 %5, i32* %12
%23 = call i32 (...) @omp_get_max_threads()
store i32 %23, i32* %13
%24 = load i32, i32* %10
%25 = load i32, i32* %13
%26 = icmp sgt i32 %24, %25
br i1 %26, label %27, label %29
27: ; preds = %6
%28 = load i32, i32* %13
br label %30
29: ; preds = %6
br label %30
30: ; preds = %29, %27
%31 = phi i32 [ %28, %27 ], [ 1, %29 ]
store i32 %31, i32* %14
%32 = load i32, i32* %10
%33 = load i32, i32* %14
%34 = sdiv i32 %32, %33
store i32 %34, i32* %15
%35 = load i32, i32* %10
%36 = load i32, i32* %14
%37 = srem i32 %35, %36
store i32 %37, i32* %16
store i32 0, i32* %17
br label %38
38: ; preds = %93, %30
%39 = load i32, i32* %17
%40 = icmp slt i32 %39, 100
br i1 %40, label %41, label %96
41: ; preds = %38
%42 = call i32 (...) @omp_get_thread_num()
store i32 %42, i32* %18
store i32 0, i32* %19
store i32 0, i32* %20
store i32 0, i32* %21
%43 = load i32, i32* %18
%44 = load i32, i32* %16
%45 = icmp slt i32 %43, %44
br i1 %45, label %46, label %52
46: ; preds = %41
%47 = load i32, i32* %15
%48 = add nsw i32 %47, 1
store i32 %48, i32* %19
%49 = load i32, i32* %18
%50 = load i32, i32* %19
%51 = mul nsw i32 %49, %50
store i32 %51, i32* %20
br label %59
52: ; preds = %41
%53 = load i32, i32* %15
store i32 %53, i32* %19
%54 = load i32, i32* %18
%55 = load i32, i32* %19
%56 = mul nsw i32 %54, %55
%57 = load i32, i32* %16
%58 = add nsw i32 %56, %57
store i32 %58, i32* %20
br label %59
59: ; preds = %52, %46
%60 = load i32, i32* %20
%61 = load i32, i32* %19
%62 = add nsw i32 %60, %61
store i32 %62, i32* %21
%63 = load i32, i32* %19
%64 = icmp ne i32 %63, 0
br i1 %64, label %65, label %92
65: ; preds = %59
%66 = load i32, i32* %20
store i32 %66, i32* %22
br label %67
67: ; preds = %88, %65
%68 = load i32, i32* %22
%69 = load i32, i32* %21
%70 = icmp slt i32 %68, %69
br i1 %70, label %71, label %91
71: ; preds = %67
%72 = load double*, double** %7
%73 = load i32, i32* %22
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds double, double* %72, i64 %74
%76 = load double, double* %75
%77 = load double*, double** %8
%78 = load i32, i32* %22
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds double, double* %77, i64 %79
%81 = load double, double* %80
%82 = fadd double %76, %81
%83 = fdiv double %82, 1.000000e+02
%84 = load double*, double** %9
%85 = load i32, i32* %22
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %84, i64 %86
store double %83, double* %87
br label %88
88: ; preds = %71
%89 = load i32, i32* %22
%90 = add nsw i32 %89, 1
store i32 %90, i32* %22
br label %67
91: ; preds = %67
br label %92
92: ; preds = %91, %59
br label %93
93: ; preds = %92
%94 = load i32, i32* %17
%95 = add nsw i32 %94, 1
store i32 %95, i32* %17
br label %38
96: ; preds = %38
ret void
}

LLVM IR (-O2):
define void @voodoo(double* nocapture %0, double* nocapture %1, double* nocapture writeonly %2, i32 %3, i32 %4, i32 %5) {
%7 = tail call i32 (...) @omp_get_max_threads()
%8 = icmp slt i32 %7, %3
%9 = select i1 %8, i32 %7, i32 1
%10 = sdiv i32 %3, %9
%11 = srem i32 %3, %9
%12 = add nsw i32 %10, 1
br label %14
13: ; preds = %117
ret void
14: ; preds = %6, %117
%15 = phi i32 [ 0, %6 ], [ %118, %117 ]
%16 = tail call i32 (...) @omp_get_thread_num()
%17 = icmp slt i32 %16, %11
%18 = mul nsw i32 %16, %12
%19 = mul nsw i32 %16, %10
%20 = add nsw i32 %19, %11
%21 = select i1 %17, i32 %12, i32 %10
%22 = select i1 %17, i32 %18, i32 %20
%23 = icmp sgt i32 %21, 0
br i1 %23, label %24, label %117
24: ; preds = %14
%25 = add nsw i32 %22, %21
%26 = sext i32 %22 to i64
%27 = sext i32 %25 to i64
%28 = add nsw i64 %26, 1
%29 = call i64 @llvm.smax.i64(i64 %28, i64 %27)
%30 = sub i64 %29, %26
%31 = icmp ult i64 %30, 2
br i1 %31, label %104, label %32
32: ; preds = %24
%33 = getelementptr double, double* %2, i64 %26
%34 = add nsw i64 %26, 1
%35 = call i64 @llvm.smax.i64(i64 %34, i64 %27)
%36 = getelementptr double, double* %2, i64 %35
%37 = getelementptr double, double* %0, i64 %26
%38 = getelementptr double, double* %0, i64 %35
%39 = getelementptr double, double* %1, i64 %26
%40 = getelementptr double, double* %1, i64 %35
%41 = icmp ult double* %33, %38
%42 = icmp ult double* %37, %36
%43 = and i1 %41, %42
%44 = icmp ult double* %33, %40
%45 = icmp ult double* %39, %36
%46 = and i1 %44, %45
%47 = or i1 %43, %46
br i1 %47, label %104, label %48
48: ; preds = %32
%49 = and i64 %30, -2
%50 = add i64 %49, %26
%51 = add i64 %49, -2
%52 = lshr exact i64 %51, 1
%53 = add nuw i64 %52, 1
%54 = and i64 %53, 1
%55 = icmp eq i64 %51, 0
br i1 %55, label %87, label %56
56: ; preds = %48
%57 = and i64 %53, -2
br label %58
58: ; preds = %58, %56
%59 = phi i64 [ 0, %56 ], [ %84, %58 ]
%60 = phi i64 [ 0, %56 ], [ %85, %58 ]
%61 = add i64 %59, %26
%62 = getelementptr inbounds double, double* %0, i64 %61
%63 = bitcast double* %62 to <2 x double>*
%64 = load <2 x double>, <2 x double>* %63
%65 = getelementptr inbounds double, double* %1, i64 %61
%66 = bitcast double* %65 to <2 x double>*
%67 = load <2 x double>, <2 x double>* %66
%68 = fadd <2 x double> %64, %67
%69 = fdiv <2 x double> %68, <double 1.000000e+02, double 1.000000e+02>
%70 = getelementptr inbounds double, double* %2, i64 %61
%71 = bitcast double* %70 to <2 x double>*
store <2 x double> %69, <2 x double>* %71
%72 = or i64 %59, 2
%73 = add i64 %72, %26
%74 = getelementptr inbounds double, double* %0, i64 %73
%75 = bitcast double* %74 to <2 x double>*
%76 = load <2 x double>, <2 x double>* %75
%77 = getelementptr inbounds double, double* %1, i64 %73
%78 = bitcast double* %77 to <2 x double>*
%79 = load <2 x double>, <2 x double>* %78
%80 = fadd <2 x double> %76, %79
%81 = fdiv <2 x double> %80, <double 1.000000e+02, double 1.000000e+02>
%82 = getelementptr inbounds double, double* %2, i64 %73
%83 = bitcast double* %82 to <2 x double>*
store <2 x double> %81, <2 x double>* %83
%84 = add nuw i64 %59, 4
%85 = add nuw i64 %60, 2
%86 = icmp eq i64 %85, %57
br i1 %86, label %87, label %58
87: ; preds = %58, %48
%88 = phi i64 [ 0, %48 ], [ %84, %58 ]
%89 = icmp eq i64 %54, 0
br i1 %89, label %102, label %90
90: ; preds = %87
%91 = add i64 %88, %26
%92 = getelementptr inbounds double, double* %0, i64 %91
%93 = bitcast double* %92 to <2 x double>*
%94 = load <2 x double>, <2 x double>* %93
%95 = getelementptr inbounds double, double* %1, i64 %91
%96 = bitcast double* %95 to <2 x double>*
%97 = load <2 x double>, <2 x double>* %96
%98 = fadd <2 x double> %94, %97
%99 = fdiv <2 x double> %98, <double 1.000000e+02, double 1.000000e+02>
%100 = getelementptr inbounds double, double* %2, i64 %91
%101 = bitcast double* %100 to <2 x double>*
store <2 x double> %99, <2 x double>* %101
br label %102
102: ; preds = %87, %90
%103 = icmp eq i64 %30, %49
br i1 %103, label %117, label %104
104: ; preds = %32, %24, %102
%105 = phi i64 [ %26, %32 ], [ %26, %24 ], [ %50, %102 ]
br label %106
106: ; preds = %104, %106
%107 = phi i64 [ %115, %106 ], [ %105, %104 ]
%108 = getelementptr inbounds double, double* %0, i64 %107
%109 = load double, double* %108
%110 = getelementptr inbounds double, double* %1, i64 %107
%111 = load double, double* %110
%112 = fadd double %109, %111
%113 = fdiv double %112, 1.000000e+02
%114 = getelementptr inbounds double, double* %2, i64 %107
store double %113, double* %114
%115 = add nsw i64 %107, 1
%116 = icmp slt i64 %115, %27
br i1 %116, label %106, label %117
117: ; preds = %106, %102, %14
%118 = add nuw nsw i32 %15, 1
%119 = icmp eq i32 %118, 100
br i1 %119, label %13, label %14
}