#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 scalar_t__ uint;

/* Variables and functions */
int /*<<< orphan*/ DC(int, double *, scalar_t__ *, scalar_t__, double, double,
                      double);
int /*<<< orphan*/ drand49();
double fabs(double);
double log(int /*<<< orphan*/);
double pow(double, double);
scalar_t__ round(double);
void POISSON_SHOT1(double amplitude, double frequency, double width,
                   double *output, uint *index, uint Ni, double expon,
                   double srate, double dt) {
  double tmp;
  uint i, o, m, j, k;
  if (frequency == 0.) {
    DC(0., output, index, Ni, expon, srate, dt);
    return;
  }
  if ((expon == -1) || (expon == 0))
    amplitude = fabs(amplitude);
  else if (expon != 1)
    amplitude = pow(amplitude, expon);
  m = (uint)(width * srate / 1000.);
  tmp = srate / frequency;
  if (frequency > 0.) {
    j = 0;
    while (j < Ni) {
      k = (uint)(-log(drand49()) * tmp);
      o = j + k;
      if ((o + m) < Ni) {
        for (i = 0; i < k; i++)
          output[(*index)++] = 0.;
        for (i = 0; i < m; i++)
          output[(*index)++] = amplitude;
        j = o + m;
      } else {
        (*index) = (*index) + (Ni - j);
        j = Ni;
      }
    }
  } else {
    k = (uint)round(-srate / frequency);
    i = -1;
    for (j = 0; j < Ni; j++, (*index)++) {
      if (j % k == 0)
        i++;
      if ((j - i * k) % m == (j - i * k))
        output[*index] = amplitude;
    }
  }
  return;
}


x86 assembly (clang, -O0):
<POISSON_SHOT1>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
mov    %rdi,-0x20(%rbp)
mov    %rsi,-0x28(%rbp)
mov    %rdx,-0x30(%rbp)
movsd  %xmm3,-0x38(%rbp)
movsd  %xmm4,-0x40(%rbp)
movsd  %xmm5,-0x48(%rbp)
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11c4 <POISSON_SHOT1+0x74>
jp     11c4 <POISSON_SHOT1+0x74>
mov    -0x20(%rbp),%rsi
mov    -0x28(%rbp),%rdx
mov    -0x30(%rbp),%rcx
movsd  -0x38(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
movsd  -0x48(%rbp),%xmm2
xor    %edi,%edi
call   1030 <DC@plt>
jmp    149a <POISSON_SHOT1+0x34a>
movsd  -0x38(%rbp),%xmm0
movsd  0xe2f(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    11e6 <POISSON_SHOT1+0x96>
jp     11e6 <POISSON_SHOT1+0x96>
jmp    11fe <POISSON_SHOT1+0xae>
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1218 <POISSON_SHOT1+0xc8>
jp     1218 <POISSON_SHOT1+0xc8>
movsd  -0x8(%rbp),%xmm0
movaps 0xe16(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1253 <POISSON_SHOT1+0x103>
movsd  -0x38(%rbp),%xmm0
movsd  0xde3(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    123a <POISSON_SHOT1+0xea>
jp     123a <POISSON_SHOT1+0xea>
jmp    124e <POISSON_SHOT1+0xfe>
movsd  -0x8(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm1
call   1050 <pow@plt>
movsd  %xmm0,-0x8(%rbp)
jmp    1253 <POISSON_SHOT1+0x103>
movsd  -0x18(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  0xdab(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rax
mov    %rax,-0x68(%rbp)
movsd  -0x40(%rbp),%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    13c0 <POISSON_SHOT1+0x270>
movq   $0x0,-0x70(%rbp)
mov    -0x70(%rbp),%rax
cmp    -0x30(%rbp),%rax
jge    13bb <POISSON_SHOT1+0x26b>
mov    $0x0,%al
call   1070 <drand49@plt>
mov    %eax,%edi
call   1040 <log@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x50(%rbp),%xmm0
cvttsd2si %xmm0,%rax
mov    %rax,-0x78(%rbp)
mov    -0x70(%rbp),%rax
add    -0x78(%rbp),%rax
mov    %rax,-0x60(%rbp)
mov    -0x60(%rbp),%rax
add    -0x68(%rbp),%rax
cmp    -0x30(%rbp),%rax
jge    1395 <POISSON_SHOT1+0x245>
movq   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%rax
cmp    -0x78(%rbp),%rax
jge    133e <POISSON_SHOT1+0x1ee>
mov    -0x20(%rbp),%rax
mov    -0x28(%rbp),%rdx
mov    (%rdx),%rcx
mov    %rcx,%rsi
add    $0x1,%rsi
mov    %rsi,(%rdx)
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x58(%rbp)
jmp    1302 <POISSON_SHOT1+0x1b2>
movq   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%rax
cmp    -0x68(%rbp),%rax
jge    1384 <POISSON_SHOT1+0x234>
movsd  -0x8(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x28(%rbp),%rdx
mov    (%rdx),%rcx
mov    %rcx,%rsi
add    $0x1,%rsi
mov    %rsi,(%rdx)
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x58(%rbp)
jmp    1346 <POISSON_SHOT1+0x1f6>
mov    -0x60(%rbp),%rax
add    -0x68(%rbp),%rax
mov    %rax,-0x70(%rbp)
jmp    13b6 <POISSON_SHOT1+0x266>
mov    -0x28(%rbp),%rax
mov    (%rax),%rcx
mov    -0x30(%rbp),%rax
sub    -0x70(%rbp),%rax
add    %rax,%rcx
mov    -0x28(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x30(%rbp),%rax
mov    %rax,-0x70(%rbp)
jmp    129b <POISSON_SHOT1+0x14b>
jmp    1495 <POISSON_SHOT1+0x345>
movsd  -0x40(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
divsd  -0x10(%rbp),%xmm0
call   1060 <round@plt>
mov    %rax,-0x78(%rbp)
movq   $0xffffffffffffffff,-0x58(%rbp)
ff
movq   $0x0,-0x70(%rbp)
mov    -0x70(%rbp),%rax
cmp    -0x30(%rbp),%rax
jge    1490 <POISSON_SHOT1+0x340>
mov    -0x70(%rbp),%rax
cqto
idivq  -0x78(%rbp)
cmp    $0x0,%rdx
jne    1428 <POISSON_SHOT1+0x2d8>
mov    -0x58(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x58(%rbp)
mov    -0x70(%rbp),%rax
mov    -0x58(%rbp),%rcx
imul   -0x78(%rbp),%rcx
sub    %rcx,%rax
cqto
idivq  -0x68(%rbp)
mov    -0x70(%rbp),%rax
mov    -0x58(%rbp),%rcx
imul   -0x78(%rbp),%rcx
sub    %rcx,%rax
cmp    %rax,%rdx
jne    146c <POISSON_SHOT1+0x31c>
movsd  -0x8(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x28(%rbp),%rcx
mov    (%rcx),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1471 <POISSON_SHOT1+0x321>
mov    -0x70(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x70(%rbp)
mov    -0x28(%rbp),%rax
mov    (%rax),%rcx
add    $0x1,%rcx
mov    %rcx,(%rax)
jmp    13fa <POISSON_SHOT1+0x2aa>
jmp    1495 <POISSON_SHOT1+0x345>
jmp    149a <POISSON_SHOT1+0x34a>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<POISSON_SHOT1>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x98,%rsp
mov    %rdx,%rcx
mov    %rsi,%rdx
mov    %rdi,%rbp
movapd %xmm1,%xmm6
movapd %xmm0,%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm6
jne    11b2 <POISSON_SHOT1+0x52>
jp     11b2 <POISSON_SHOT1+0x52>
xor    %edi,%edi
mov    %rbp,%rsi
movaps %xmm3,%xmm0
movaps %xmm4,%xmm1
movaps %xmm5,%xmm2
add    $0x98,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <DC@plt>
ucomisd 0xe46(%rip),%xmm3
mov    %rcx,0x10(%rsp)
mov    %rdx,0x18(%rsp)
jne    11c8 <POISSON_SHOT1+0x68>
jnp    121b <POISSON_SHOT1+0xbb>
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm3
jne    11d4 <POISSON_SHOT1+0x74>
jnp    121b <POISSON_SHOT1+0xbb>
ucomisd 0xe2c(%rip),%xmm3
jne    11e0 <POISSON_SHOT1+0x80>
jnp    1223 <POISSON_SHOT1+0xc3>
movapd %xmm1,%xmm0
movapd %xmm3,%xmm1
movaps %xmm4,0x30(%rsp)
movsd  %xmm6,0x8(%rsp)
movsd  %xmm2,0x20(%rsp)
call   1060 <pow@plt>
movsd  0x20(%rsp),%xmm2
movsd  0x8(%rsp),%xmm6
movapd 0x30(%rsp),%xmm4
mov    0x10(%rsp),%rcx
movapd %xmm0,%xmm1
jmp    1223 <POISSON_SHOT1+0xc3>
andpd  0xdfd(%rip),%xmm1
mov    %rbp,0x8(%rsp)
mulsd  %xmm4,%xmm2
divsd  0xddc(%rip),%xmm2
cvttsd2si %xmm2,%rbp
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm6
jbe    1445 <POISSON_SHOT1+0x2e5>
test   %rcx,%rcx
jle    1526 <POISSON_SHOT1+0x3c6>
divsd  %xmm6,%xmm4
lea    -0x4(%rbp),%rax
mov    %rax,0x70(%rsp)
shr    $0x2,%rax
add    $0x1,%rax
mov    %rbp,%rcx
and    $0xfffffffffffffffc,%rcx
mov    %rcx,0x68(%rsp)
movapd %xmm1,%xmm0
unpcklpd %xmm1,%xmm0
mov    %eax,%ecx
and    $0x3,%ecx
and    $0xfffffffffffffffc,%rax
mov    %rax,0x50(%rsp)
mov    0x8(%rsp),%rax
lea    0x70(%rax),%rdx
mov    %rdx,0x48(%rsp)
add    $0x10,%rax
mov    %rax,0x58(%rsp)
mov    %rcx,0x60(%rsp)
mov    %rcx,%rbx
shl    $0x5,%rbx
xor    %r14d,%r14d
movapd %xmm1,0x20(%rsp)
movapd %xmm4,0x30(%rsp)
movapd %xmm0,0x80(%rsp)
mov    %rbx,0x78(%rsp)
jmp    12df <POISSON_SHOT1+0x17f>
nopl   0x0(%rax)
mov    %rax,(%r8)
mov    %r13,%r14
cmp    %rdi,%r13
jge    1526 <POISSON_SHOT1+0x3c6>
xor    %eax,%eax
call   1080 <drand49@plt>
mov    %eax,%edi
call   1050 <log@plt>
movapd 0x30(%rsp),%xmm4
mov    0x10(%rsp),%rdi
xorpd  0xd30(%rip),%xmm0
mulsd  %xmm4,%xmm0
cvttsd2si %xmm0,%r15
lea    (%r14,%rbp,1),%r13
add    %r15,%r13
cmp    %rdi,%r13
jge    151b <POISSON_SHOT1+0x3bb>
mov    %rbp,%r12
test   %r15,%r15
mov    0x8(%rsp),%r14
mov    0x18(%rsp),%r8
jle    1352 <POISSON_SHOT1+0x1f2>
mov    %r8,%rbx
mov    (%r8),%rbp
lea    (%r14,%rbp,8),%rdi
lea    0x0(,%r15,8),%rdx
xor    %esi,%esi
call   1030 <memset@plt>
mov    0x10(%rsp),%rdi
add    %r15,%rbp
mov    %rbp,(%rbx)
mov    %rbx,%r8
mov    %r12,%rbp
test   %r12,%r12
movapd 0x20(%rsp),%xmm1
movapd 0x80(%rsp),%xmm0
mov    0x78(%rsp),%rbx
jle    12d3 <POISSON_SHOT1+0x173>
mov    (%r8),%rax
cmp    $0x4,%rbp
jae    1390 <POISSON_SHOT1+0x230>
xor    %edx,%edx
jmp    1428 <POISSON_SHOT1+0x2c8>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
cmpq   $0xc,0x70(%rsp)
jae    139c <POISSON_SHOT1+0x23c>
xor    %ecx,%ecx
jmp    13e9 <POISSON_SHOT1+0x289>
mov    0x48(%rsp),%rcx
lea    (%rcx,%rax,8),%rdx
mov    0x50(%rsp),%rsi
xor    %ecx,%ecx
nopl   0x0(%rax)
movupd %xmm0,-0x70(%rdx,%rcx,8)
movupd %xmm0,-0x60(%rdx,%rcx,8)
movupd %xmm0,-0x50(%rdx,%rcx,8)
movupd %xmm0,-0x40(%rdx,%rcx,8)
movupd %xmm0,-0x30(%rdx,%rcx,8)
movupd %xmm0,-0x20(%rdx,%rcx,8)
movupd %xmm0,-0x10(%rdx,%rcx,8)
movupd %xmm0,(%rdx,%rcx,8)
add    $0x10,%rcx
add    $0xfffffffffffffffc,%rsi
jne    13b0 <POISSON_SHOT1+0x250>
cmpq   $0x0,0x60(%rsp)
je     1414 <POISSON_SHOT1+0x2b4>
add    %rax,%rcx
mov    0x58(%rsp),%rdx
lea    (%rdx,%rcx,8),%rcx
xor    %edx,%edx
nop
movupd %xmm0,-0x10(%rcx,%rdx,1)
movupd %xmm0,(%rcx,%rdx,1)
add    $0x20,%rdx
cmp    %rdx,%rbx
jne    1400 <POISSON_SHOT1+0x2a0>
mov    0x68(%rsp),%rcx
add    %rcx,%rax
mov    %rcx,%rdx
cmp    %rbp,%rcx
je     12d0 <POISSON_SHOT1+0x170>
mov    %rbp,%rcx
sub    %rdx,%rcx
xchg   %ax,%ax
movsd  %xmm1,(%r14,%rax,8)
add    $0x1,%rax
add    $0xffffffffffffffff,%rcx
jne    1430 <POISSON_SHOT1+0x2d0>
jmp    12d0 <POISSON_SHOT1+0x170>
xorpd  0xbe3(%rip),%xmm4
divsd  %xmm6,%xmm4
movapd %xmm4,%xmm0
movapd %xmm1,0x20(%rsp)
call   1070 <round@plt>
movapd 0x20(%rsp),%xmm0
mov    0x10(%rsp),%r11
test   %r11,%r11
jle    1526 <POISSON_SHOT1+0x3c6>
mov    %rax,%rcx
mov    0x18(%rsp),%rax
mov    %rax,%r10
mov    (%rax),%r8
mov    0x8(%rsp),%rax
lea    (%rax,%r8,8),%r9
mov    %rcx,%r14
neg    %r14
mov    $0xffffffffffffffff,%rbx
xor    %esi,%esi
jmp    14a9 <POISSON_SHOT1+0x349>
nopl   0x0(%rax)
add    $0x1,%rsi
cmp    %rsi,%r11
je     1513 <POISSON_SHOT1+0x3b3>
mov    %rsi,%rax
or     %rcx,%rax
shr    $0x20,%rax
je     14c0 <POISSON_SHOT1+0x360>
mov    %rsi,%rax
cqto
idiv   %rcx
jmp    14c6 <POISSON_SHOT1+0x366>
nop
mov    %esi,%eax
xor    %edx,%edx
div    %ecx
cmp    $0x1,%rdx
adc    $0x0,%rbx
mov    %r14,%rdi
imul   %rbx,%rdi
add    %rsi,%rdi
mov    %rdi,%rax
or     %rbp,%rax
shr    $0x20,%rax
je     1500 <POISSON_SHOT1+0x3a0>
mov    %rdi,%rax
cqto
idiv   %rbp
cmp    %rdx,%rdi
jne    14a0 <POISSON_SHOT1+0x340>
jmp    150b <POISSON_SHOT1+0x3ab>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %edi,%eax
xor    %edx,%edx
div    %ebp
cmp    %rdx,%rdi
jne    14a0 <POISSON_SHOT1+0x340>
movsd  %xmm0,(%r9,%rsi,8)
jmp    14a0 <POISSON_SHOT1+0x340>
add    %rsi,%r8
mov    %r8,(%r10)
jmp    1526 <POISSON_SHOT1+0x3c6>
sub    %r14,%rdi
mov    0x18(%rsp),%rax
add    %rdi,(%rax)
add    $0x98,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @POISSON_SHOT1(double %0, double %1, double %2, double* %3, i64* %4, i64 %5, double %6, double %7, double %8) {
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double*
%14 = alloca i64*
%15 = alloca i64
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca i64
%21 = alloca i64
%22 = alloca i64
%23 = alloca i64
%24 = alloca i64
store double %0, double* %10
store double %1, double* %11
store double %2, double* %12
store double* %3, double** %13
store i64* %4, i64** %14
store i64 %5, i64* %15
store double %6, double* %16
store double %7, double* %17
store double %8, double* %18
%25 = load double, double* %11
%26 = fcmp oeq double %25, 0.000000e+00
br i1 %26, label %27, label %35
27: ; preds = %9
%28 = load double*, double** %13
%29 = load i64*, i64** %14
%30 = load i64, i64* %15
%31 = load double, double* %16
%32 = load double, double* %17
%33 = load double, double* %18
%34 = call i32 @DC(i32 0, double* %28, i64* %29, i64 %30, double %31, double %32, double %33)
br label %174
35: ; preds = %9
%36 = load double, double* %16
%37 = fcmp oeq double %36, -1.000000e+00
br i1 %37, label %41, label %38
38: ; preds = %35
%39 = load double, double* %16
%40 = fcmp oeq double %39, 0.000000e+00
br i1 %40, label %41, label %44
41: ; preds = %38, %35
%42 = load double, double* %10
%43 = call double @llvm.fabs.f64(double %42)
store double %43, double* %10
br label %52
44: ; preds = %38
%45 = load double, double* %16
%46 = fcmp une double %45, 1.000000e+00
br i1 %46, label %47, label %51
47: ; preds = %44
%48 = load double, double* %10
%49 = load double, double* %16
%50 = call double @pow(double %48, double %49)
store double %50, double* %10
br label %51
51: ; preds = %47, %44
br label %52
52: ; preds = %51, %41
%53 = load double, double* %12
%54 = load double, double* %17
%55 = fmul double %53, %54
%56 = fdiv double %55, 1.000000e+03
%57 = fptosi double %56 to i64
store i64 %57, i64* %22
%58 = load double, double* %17
%59 = load double, double* %11
%60 = fdiv double %58, %59
store double %60, double* %19
%61 = load double, double* %11
%62 = fcmp ogt double %61, 0.000000e+00
br i1 %62, label %63, label %127
63: ; preds = %52
store i64 0, i64* %23
br label %64
64: ; preds = %125, %63
%65 = load i64, i64* %23
%66 = load i64, i64* %15
%67 = icmp slt i64 %65, %66
br i1 %67, label %68, label %126
68: ; preds = %64
%69 = call i32 (...) @drand49()
%70 = call double @log(i32 %69)
%71 = fneg double %70
%72 = load double, double* %19
%73 = fmul double %71, %72
%74 = fptosi double %73 to i64
store i64 %74, i64* %24
%75 = load i64, i64* %23
%76 = load i64, i64* %24
%77 = add nsw i64 %75, %76
store i64 %77, i64* %21
%78 = load i64, i64* %21
%79 = load i64, i64* %22
%80 = add nsw i64 %78, %79
%81 = load i64, i64* %15
%82 = icmp slt i64 %80, %81
br i1 %82, label %83, label %116
83: ; preds = %68
store i64 0, i64* %20
br label %84
84: ; preds = %94, %83
%85 = load i64, i64* %20
%86 = load i64, i64* %24
%87 = icmp slt i64 %85, %86
br i1 %87, label %88, label %97
88: ; preds = %84
%89 = load double*, double** %13
%90 = load i64*, i64** %14
%91 = load i64, i64* %90
%92 = add nsw i64 %91, 1
store i64 %92, i64* %90
%93 = getelementptr inbounds double, double* %89, i64 %91
store double 0.000000e+00, double* %93
br label %94
94: ; preds = %88
%95 = load i64, i64* %20
%96 = add nsw i64 %95, 1
store i64 %96, i64* %20
br label %84
97: ; preds = %84
store i64 0, i64* %20
br label %98
98: ; preds = %109, %97
%99 = load i64, i64* %20
%100 = load i64, i64* %22
%101 = icmp slt i64 %99, %100
br i1 %101, label %102, label %112
102: ; preds = %98
%103 = load double, double* %10
%104 = load double*, double** %13
%105 = load i64*, i64** %14
%106 = load i64, i64* %105
%107 = add nsw i64 %106, 1
store i64 %107, i64* %105
%108 = getelementptr inbounds double, double* %104, i64 %106
store double %103, double* %108
br label %109
109: ; preds = %102
%110 = load i64, i64* %20
%111 = add nsw i64 %110, 1
store i64 %111, i64* %20
br label %98
112: ; preds = %98
%113 = load i64, i64* %21
%114 = load i64, i64* %22
%115 = add nsw i64 %113, %114
store i64 %115, i64* %23
br label %125
116: ; preds = %68
%117 = load i64*, i64** %14
%118 = load i64, i64* %117
%119 = load i64, i64* %15
%120 = load i64, i64* %23
%121 = sub nsw i64 %119, %120
%122 = add nsw i64 %118, %121
%123 = load i64*, i64** %14
store i64 %122, i64* %123
%124 = load i64, i64* %15
store i64 %124, i64* %23
br label %125
125: ; preds = %116, %112
br label %64
126: ; preds = %64
br label %173
127: ; preds = %52
%128 = load double, double* %17
%129 = fneg double %128
%130 = load double, double* %11
%131 = fdiv double %129, %130
%132 = call i64 @round(double %131)
store i64 %132, i64* %24
store i64 -1, i64* %20
store i64 0, i64* %23
br label %133
133: ; preds = %166, %127
%134 = load i64, i64* %23
%135 = load i64, i64* %15
%136 = icmp slt i64 %134, %135
br i1 %136, label %137, label %172
137: ; preds = %133
%138 = load i64, i64* %23
%139 = load i64, i64* %24
%140 = srem i64 %138, %139
%141 = icmp eq i64 %140, 0
br i1 %141, label %142, label %145
142: ; preds = %137
%143 = load i64, i64* %20
%144 = add nsw i64 %143, 1
store i64 %144, i64* %20
br label %145
145: ; preds = %142, %137
%146 = load i64, i64* %23
%147 = load i64, i64* %20
%148 = load i64, i64* %24
%149 = mul nsw i64 %147, %148
%150 = sub nsw i64 %146, %149
%151 = load i64, i64* %22
%152 = srem i64 %150, %151
%153 = load i64, i64* %23
%154 = load i64, i64* %20
%155 = load i64, i64* %24
%156 = mul nsw i64 %154, %155
%157 = sub nsw i64 %153, %156
%158 = icmp eq i64 %152, %157
br i1 %158, label %159, label %165
159: ; preds = %145
%160 = load double, double* %10
%161 = load double*, double** %13
%162 = load i64*, i64** %14
%163 = load i64, i64* %162
%164 = getelementptr inbounds double, double* %161, i64 %163
store double %160, double* %164
br label %165
165: ; preds = %159, %145
br label %166
166: ; preds = %165
%167 = load i64, i64* %23
%168 = add nsw i64 %167, 1
store i64 %168, i64* %23
%169 = load i64*, i64** %14
%170 = load i64, i64* %169
%171 = add nsw i64 %170, 1
store i64 %171, i64* %169
br label %133
172: ; preds = %133
br label %173
173: ; preds = %172, %126
br label %174
174: ; preds = %173, %27
ret void
}

LLVM IR (-O2):
define void @POISSON_SHOT1(double %0, double %1, double %2, double* %3, i64* %4, i64 %5, double %6, double %7, double %8) {
%10 = fcmp oeq double %1, 0.000000e+00
br i1 %10, label %11, label %13
11: ; preds = %9
%12 = tail call i32 @DC(i32 0, double* %3, i64* %4, i64 %5, double %6, double %7, double %8)
br label %158
13: ; preds = %9
%14 = fcmp oeq double %6, -1.000000e+00
%15 = fcmp oeq double %6, 0.000000e+00
%16 = or i1 %14, %15
br i1 %16, label %17, label %19
17: ; preds = %13
%18 = tail call double @llvm.fabs.f64(double %0)
br label %23
19: ; preds = %13
%20 = fcmp une double %6, 1.000000e+00
br i1 %20, label %21, label %23
21: ; preds = %19
%22 = tail call double @pow(double %0, double %6)
br label %23
23: ; preds = %19, %21, %17
%24 = phi double [ %18, %17 ], [ %22, %21 ], [ %0, %19 ]
%25 = fmul double %2, %7
%26 = fdiv double %25, 1.000000e+03
%27 = fptosi double %26 to i64
%28 = fdiv double %7, %1
%29 = fcmp ogt double %1, 0.000000e+00
br i1 %29, label %30, label %132
30: ; preds = %23
%31 = icmp sgt i64 %5, 0
br i1 %31, label %32, label %158
32: ; preds = %30
%33 = icmp sgt i64 %27, 0
%34 = add i64 %27, -4
%35 = lshr i64 %34, 2
%36 = add nuw nsw i64 %35, 1
%37 = icmp ult i64 %27, 4
%38 = and i64 %27, -4
%39 = insertelement <2 x double> poison, double %24, i64 0
%40 = shufflevector <2 x double> %39, <2 x double> poison, <2 x i32> zeroinitializer
%41 = insertelement <2 x double> poison, double %24, i64 0
%42 = shufflevector <2 x double> %41, <2 x double> poison, <2 x i32> zeroinitializer
%43 = and i64 %36, 3
%44 = icmp ult i64 %34, 12
%45 = and i64 %36, 9223372036854775804
%46 = icmp eq i64 %43, 0
%47 = icmp eq i64 %38, %27
br label %48
48: ; preds = %32, %130
%49 = phi i64 [ 0, %32 ], [ %56, %130 ]
%50 = tail call i32 (...) @drand49()
%51 = tail call double @log(i32 %50)
%52 = fneg double %51
%53 = fmul double %28, %52
%54 = fptosi double %53 to i64
%55 = add i64 %49, %27
%56 = add i64 %55, %54
%57 = icmp slt i64 %56, %5
br i1 %57, label %58, label %124
58: ; preds = %48
%59 = icmp sgt i64 %54, 0
br i1 %59, label %60, label %66
60: ; preds = %58
%61 = load i64, i64* %4
%62 = getelementptr double, double* %3, i64 %61
%63 = bitcast double* %62 to i8*
%64 = shl nuw i64 %54, 3
%65 = add i64 %61, %54
store i64 %65, i64* %4
br label %66
66: ; preds = %60, %58
br i1 %33, label %67, label %130
67: ; preds = %66
%68 = load i64, i64* %4
br i1 %37, label %114, label %69
69: ; preds = %67
%70 = add i64 %68, %38
br i1 %44, label %100, label %71
71: ; preds = %69, %71
%72 = phi i64 [ %97, %71 ], [ 0, %69 ]
%73 = phi i64 [ %98, %71 ], [ 0, %69 ]
%74 = add i64 %68, %72
%75 = getelementptr inbounds double, double* %3, i64 %74
%76 = bitcast double* %75 to <2 x double>*
store <2 x double> %40, <2 x double>* %76
%77 = getelementptr inbounds double, double* %75, i64 2
%78 = bitcast double* %77 to <2 x double>*
store <2 x double> %42, <2 x double>* %78
%79 = or i64 %72, 4
%80 = add i64 %68, %79
%81 = getelementptr inbounds double, double* %3, i64 %80
%82 = bitcast double* %81 to <2 x double>*
store <2 x double> %40, <2 x double>* %82
%83 = getelementptr inbounds double, double* %81, i64 2
%84 = bitcast double* %83 to <2 x double>*
store <2 x double> %42, <2 x double>* %84
%85 = or i64 %72, 8
%86 = add i64 %68, %85
%87 = getelementptr inbounds double, double* %3, i64 %86
%88 = bitcast double* %87 to <2 x double>*
store <2 x double> %40, <2 x double>* %88
%89 = getelementptr inbounds double, double* %87, i64 2
%90 = bitcast double* %89 to <2 x double>*
store <2 x double> %42, <2 x double>* %90
%91 = or i64 %72, 12
%92 = add i64 %68, %91
%93 = getelementptr inbounds double, double* %3, i64 %92
%94 = bitcast double* %93 to <2 x double>*
store <2 x double> %40, <2 x double>* %94
%95 = getelementptr inbounds double, double* %93, i64 2
%96 = bitcast double* %95 to <2 x double>*
store <2 x double> %42, <2 x double>* %96
%97 = add nuw i64 %72, 16
%98 = add nuw i64 %73, 4
%99 = icmp eq i64 %98, %45
br i1 %99, label %100, label %71
100: ; preds = %71, %69
%101 = phi i64 [ 0, %69 ], [ %97, %71 ]
br i1 %46, label %113, label %102
102: ; preds = %100, %102
%103 = phi i64 [ %110, %102 ], [ %101, %100 ]
%104 = phi i64 [ %111, %102 ], [ 0, %100 ]
%105 = add i64 %68, %103
%106 = getelementptr inbounds double, double* %3, i64 %105
%107 = bitcast double* %106 to <2 x double>*
store <2 x double> %40, <2 x double>* %107
%108 = getelementptr inbounds double, double* %106, i64 2
%109 = bitcast double* %108 to <2 x double>*
store <2 x double> %42, <2 x double>* %109
%110 = add nuw i64 %103, 4
%111 = add i64 %104, 1
%112 = icmp eq i64 %111, %43
br i1 %112, label %113, label %102
113: ; preds = %102, %100
br i1 %47, label %128, label %114
114: ; preds = %67, %113
%115 = phi i64 [ %68, %67 ], [ %70, %113 ]
%116 = phi i64 [ 0, %67 ], [ %38, %113 ]
br label %117
117: ; preds = %114, %117
%118 = phi i64 [ %120, %117 ], [ %115, %114 ]
%119 = phi i64 [ %122, %117 ], [ %116, %114 ]
%120 = add nsw i64 %118, 1
%121 = getelementptr inbounds double, double* %3, i64 %118
store double %24, double* %121
%122 = add nuw nsw i64 %119, 1
%123 = icmp eq i64 %122, %27
br i1 %123, label %128, label %117
124: ; preds = %48
%125 = load i64, i64* %4
%126 = sub i64 %5, %49
%127 = add nsw i64 %126, %125
store i64 %127, i64* %4
br label %158
128: ; preds = %117, %113
%129 = phi i64 [ %70, %113 ], [ %120, %117 ]
store i64 %129, i64* %4
br label %130
130: ; preds = %66, %128
%131 = icmp slt i64 %56, %5
br i1 %131, label %48, label %158
132: ; preds = %23
%133 = fneg double %7
%134 = fdiv double %133, %1
%135 = tail call i64 @round(double %134)
%136 = icmp sgt i64 %5, 0
br i1 %136, label %137, label %158
137: ; preds = %132
%138 = load i64, i64* %4
br label %139
139: ; preds = %137, %153
%140 = phi i64 [ %138, %137 ], [ %155, %153 ]
%141 = phi i64 [ 0, %137 ], [ %154, %153 ]
%142 = phi i64 [ -1, %137 ], [ %146, %153 ]
%143 = srem i64 %141, %135
%144 = icmp eq i64 %143, 0
%145 = zext i1 %144 to i64
%146 = add nsw i64 %142, %145
%147 = mul nsw i64 %146, %135
%148 = sub nsw i64 %141, %147
%149 = srem i64 %148, %27
%150 = icmp eq i64 %149, %148
br i1 %150, label %151, label %153
151: ; preds = %139
%152 = getelementptr inbounds double, double* %3, i64 %140
store double %24, double* %152
br label %153
153: ; preds = %139, %151
%154 = add nuw nsw i64 %141, 1
%155 = add nsw i64 %140, 1
%156 = icmp eq i64 %154, %5
br i1 %156, label %157, label %139
157: ; preds = %153
store i64 %155, i64* %4
br label %158
158: ; preds = %130, %124, %132, %157, %30, %11
ret void
}