#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 sample(int /*<<< orphan*/, double *, int);
int samplebinsearch(int /*<<< orphan*/, double *, int);
void csamplestatepathH(int T, int *obs, double **fwds, int R, int N,
                       double ****mats, double *****rmats, int *states,
                       int *storedpos, int *storedfpos, double ***sneak) {
  int j, s, e, p, md;
  int pos, cs, js, je, fpos, si;
  double *distribution;
  md = T % R;
  if (md == 0) {
    md = R;
    p = T / R;
  } else {
    p = T / R + 1;
  }
  e = T;
  s = T - md;
  je = T - 2;
  double tmp[N];
  tmp[0] = fwds[p][0];
  int i;
  for (i = 1; i < N; i++)
    tmp[i] = tmp[i - 1] + fwds[p][i];
  states[e - 1] = sample(0, tmp, N);
  while (s >= 0) {
    p--;
    cs = states[e - 1];
    if (s > 0) {
      js = s - 1;
      pos = storedpos[s];
      fpos = storedfpos[s];
      states[js] = samplebinsearch(0, sneak[p + 1][cs], N);
    } else {
      js = 0;
      pos = storedpos[1];
      fpos = storedfpos[1];
      states[js] = samplebinsearch(0, sneak[p + 1][cs], N);
    }
    for (; js < je; js++) {
      distribution = rmats[pos][fpos][states[js]][cs];
      states[js + 1] = samplebinsearch(0, distribution, N);
      pos = storedpos[js + 2];
      fpos = storedfpos[js + 2];
    }
    s -= R;
    e = s + R;
    je = R + s - 2;
  }
}


x86 assembly (clang, -O0):
<csamplestatepathH>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    -0x4(%rbp),%eax
cltd
idivl  -0x1c(%rbp)
mov    %edx,-0x3c(%rbp)
cmpl   $0x0,-0x3c(%rbp)
jne    117e <csamplestatepathH+0x5e>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x4(%rbp),%eax
cltd
idivl  -0x1c(%rbp)
mov    %eax,-0x38(%rbp)
jmp    118b <csamplestatepathH+0x6b>
mov    -0x4(%rbp),%eax
cltd
idivl  -0x1c(%rbp)
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x4(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x4(%rbp),%eax
mov    -0x3c(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,-0x30(%rbp)
mov    -0x4(%rbp),%eax
add    $0xfffffffe,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,%ecx
mov    %rsp,%rax
mov    %rax,-0x68(%rbp)
lea    0xf(,%rcx,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rax
sub    %rdx,%rax
mov    %rax,-0x80(%rbp)
mov    %rax,%rsp
mov    %rcx,-0x70(%rbp)
mov    -0x18(%rbp),%rcx
movslq -0x38(%rbp),%rdx
mov    (%rcx,%rdx,8),%rcx
movsd  (%rcx),%xmm0
movsd  %xmm0,(%rax)
movl   $0x1,-0x74(%rbp)
mov    -0x74(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1233 <csamplestatepathH+0x113>
mov    -0x80(%rbp),%rax
mov    -0x74(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rcx
movslq -0x38(%rbp),%rdx
mov    (%rcx,%rdx,8),%rcx
movslq -0x74(%rbp),%rdx
addsd  (%rcx,%rdx,8),%xmm0
movslq -0x74(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    11e9 <csamplestatepathH+0xc9>
mov    -0x80(%rbp),%rsi
mov    -0x20(%rbp),%edx
xor    %edi,%edi
call   1040 <sample@plt>
mov    %eax,%edx
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
cmpl   $0x0,-0x30(%rbp)
jl     13e2 <csamplestatepathH+0x2c2>
mov    -0x38(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x38(%rbp)
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x44(%rbp)
cmpl   $0x0,-0x30(%rbp)
jle    12dd <csamplestatepathH+0x1bd>
mov    -0x30(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x48(%rbp)
mov    0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
mov    0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x50(%rbp)
mov    0x30(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rsi
mov    -0x20(%rbp),%edx
xor    %edi,%edi
call   1030 <samplebinsearch@plt>
mov    %eax,%edx
mov    0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
jmp    1328 <csamplestatepathH+0x208>
movl   $0x0,-0x48(%rbp)
mov    0x20(%rbp),%rax
mov    0x4(%rax),%eax
mov    %eax,-0x40(%rbp)
mov    0x28(%rbp),%rax
mov    0x4(%rax),%eax
mov    %eax,-0x50(%rbp)
mov    0x30(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rsi
mov    -0x20(%rbp),%edx
xor    %edi,%edi
call   1030 <samplebinsearch@plt>
mov    %eax,%edx
mov    0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
jmp    132d <csamplestatepathH+0x20d>
mov    -0x48(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jge    13bd <csamplestatepathH+0x29d>
mov    0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    0x18(%rbp),%rcx
movslq -0x48(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x60(%rbp)
mov    -0x60(%rbp),%rsi
mov    -0x20(%rbp),%edx
xor    %edi,%edi
call   1030 <samplebinsearch@plt>
mov    %eax,%edx
mov    0x18(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    0x20(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
mov    0x28(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x50(%rbp)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    132d <csamplestatepathH+0x20d>
mov    -0x1c(%rbp),%ecx
mov    -0x30(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
add    -0x1c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x1c(%rbp),%eax
add    -0x30(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x4c(%rbp)
jmp    1253 <csamplestatepathH+0x133>
mov    -0x68(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<csamplestatepathH>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rdx,%r9
mov    %edi,%r12d
mov    %edi,%eax
cltd
idiv   %ecx
xor    %edi,%edi
test   %edx,%edx
mov    %ecx,-0x3c(%rbp)
cmove  %ecx,%edx
setne  %dil
add    %eax,%edi
mov    %r12d,%r14d
sub    %edx,%r14d
mov    %r8d,%r10d
mov    %rsp,%rax
lea    0xf(,%r10,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rax,%rsi
sub    %rdx,%rsi
mov    %rsi,%rsp
neg    %rdx
movslq %edi,%rcx
mov    %rcx,%r15
mov    (%r9,%rcx,8),%rcx
movsd  (%rcx),%xmm0
movsd  %xmm0,(%rax,%rdx,1)
cmp    $0x2,%r8d
jl     1213 <csamplestatepathH+0xf3>
lea    -0x1(%r10),%rdi
add    $0xfffffffffffffffe,%r10
mov    %edi,%eax
and    $0x3,%eax
mov    $0x1,%edx
cmp    $0x3,%r10
jb     11ed <csamplestatepathH+0xcd>
and    $0xfffffffffffffffc,%rdi
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
addsd  0x8(%rcx,%rdx,8),%xmm0
movsd  %xmm0,0x8(%rsi,%rdx,8)
addsd  0x10(%rcx,%rdx,8),%xmm0
movsd  %xmm0,0x10(%rsi,%rdx,8)
addsd  0x18(%rcx,%rdx,8),%xmm0
movsd  %xmm0,0x18(%rsi,%rdx,8)
addsd  0x20(%rcx,%rdx,8),%xmm0
movsd  %xmm0,0x20(%rsi,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%rdi
jne    11b0 <csamplestatepathH+0x90>
add    $0x1,%rdx
test   %rax,%rax
je     1213 <csamplestatepathH+0xf3>
lea    (%rsi,%rdx,8),%rdi
lea    (%rcx,%rdx,8),%rcx
xor    %edx,%edx
nopl   0x0(%rax)
addsd  (%rcx,%rdx,8),%xmm0
movsd  %xmm0,(%rdi,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1200 <csamplestatepathH+0xe0>
xor    %edi,%edi
mov    %r8d,-0x2c(%rbp)
mov    %r8d,%edx
call   1040 <sample@plt>
movslq %r12d,%rcx
mov    0x18(%rbp),%rdx
mov    %eax,-0x4(%rdx,%rcx,4)
test   %r14d,%r14d
js     138e <csamplestatepathH+0x26e>
mov    %r15,%rsi
mov    0x28(%rbp),%r8
mov    0x20(%rbp),%rbx
jmp    126c <csamplestatepathH+0x14c>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    -0x48(%rbp),%rsi
add    $0xffffffffffffffff,%rsi
mov    -0x50(%rbp),%r12
mov    %r12d,%r14d
sub    -0x3c(%rbp),%r14d
mov    %rdi,%rbx
js     138e <csamplestatepathH+0x26e>
mov    %r12d,%r15d
mov    %r14d,%ecx
cmp    $0x2,%r14d
mov    $0x1,%edi
cmovae %r14d,%edi
movslq %r12d,%rax
add    $0xfffffffe,%r15d
movslq -0x4(%rdx,%rax,4),%rdx
test   %r14d,%r14d
mov    %rcx,-0x50(%rbp)
mov    %rsi,-0x48(%rbp)
mov    %edi,-0x38(%rbp)
je     12d0 <csamplestatepathH+0x1b0>
lea    -0x1(%rcx),%r12d
mov    %ecx,%eax
mov    (%rbx,%rax,4),%r14d
mov    (%r8,%rax,4),%r13d
mov    0x30(%rbp),%rax
mov    (%rax,%rsi,8),%rax
mov    %rdx,%rbx
mov    (%rax,%rdx,8),%rsi
xor    %edi,%edi
mov    -0x2c(%rbp),%edx
call   1030 <samplebinsearch@plt>
mov    0x10(%rbp),%rsi
mov    0x18(%rbp),%rdx
mov    %eax,(%rdx,%r12,4)
jmp    12fe <csamplestatepathH+0x1de>
mov    0x4(%rbx),%r14d
mov    0x4(%r8),%r13d
mov    0x30(%rbp),%rax
mov    (%rax,%rsi,8),%rax
mov    %rdx,%rbx
mov    (%rax,%rdx,8),%rsi
xor    %r12d,%r12d
xor    %edi,%edi
mov    -0x2c(%rbp),%edx
call   1030 <samplebinsearch@plt>
mov    0x10(%rbp),%rsi
mov    0x18(%rbp),%rdx
mov    %eax,(%rdx)
cmp    %r15d,%r12d
mov    -0x2c(%rbp),%r9d
mov    0x28(%rbp),%r8
mov    0x20(%rbp),%rdi
jge    1250 <csamplestatepathH+0x130>
mov    -0x38(%rbp),%eax
add    $0xffffffff,%eax
movslq %eax,%r12
movslq %r15d,%rax
mov    %rax,-0x38(%rbp)
mov    (%rdx,%r12,4),%eax
nopw   0x0(%rax,%rax,1)
movslq %r14d,%rcx
mov    (%rsi,%rcx,8),%rcx
movslq %r13d,%rdx
mov    (%rcx,%rdx,8),%rcx
cltq
mov    (%rcx,%rax,8),%rax
mov    (%rax,%rbx,8),%rsi
mov    %rdi,%r13
xor    %edi,%edi
mov    %r9d,%edx
mov    %r8,%r14
mov    %r9d,%r15d
call   1030 <samplebinsearch@plt>
mov    %r15d,%r9d
mov    %r13,%rdi
mov    %r14,%r8
mov    0x10(%rbp),%rsi
mov    0x18(%rbp),%rdx
mov    %eax,0x4(%rdx,%r12,4)
mov    0x8(%r13,%r12,4),%r14d
mov    0x8(%r8,%r12,4),%r13d
lea    0x1(%r12),%rcx
mov    %rcx,%r12
cmp    -0x38(%rbp),%rcx
jl     1330 <csamplestatepathH+0x210>
jmp    1250 <csamplestatepathH+0x130>
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @csamplestatepathH(i32 %0, i32* %1, double** %2, i32 %3, i32 %4, double**** %5, double***** %6, i32* %7, i32* %8, i32* %9, double*** %10) {
%12 = alloca i32
%13 = alloca i32*
%14 = alloca double**
%15 = alloca i32
%16 = alloca i32
%17 = alloca double****
%18 = alloca double*****
%19 = alloca i32*
%20 = alloca i32*
%21 = alloca i32*
%22 = alloca double***
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca i32
%27 = alloca i32
%28 = alloca i32
%29 = alloca i32
%30 = alloca i32
%31 = alloca i32
%32 = alloca i32
%33 = alloca i32
%34 = alloca double*
%35 = alloca i8*
%36 = alloca i64
%37 = alloca i32
store i32 %0, i32* %12
store i32* %1, i32** %13
store double** %2, double*** %14
store i32 %3, i32* %15
store i32 %4, i32* %16
store double**** %5, double***** %17
store double***** %6, double****** %18
store i32* %7, i32** %19
store i32* %8, i32** %20
store i32* %9, i32** %21
store double*** %10, double**** %22
%38 = load i32, i32* %12
%39 = load i32, i32* %15
%40 = srem i32 %38, %39
store i32 %40, i32* %27
%41 = load i32, i32* %27
%42 = icmp eq i32 %41, 0
br i1 %42, label %43, label %48
43: ; preds = %11
%44 = load i32, i32* %15
store i32 %44, i32* %27
%45 = load i32, i32* %12
%46 = load i32, i32* %15
%47 = sdiv i32 %45, %46
store i32 %47, i32* %26
br label %53
48: ; preds = %11
%49 = load i32, i32* %12
%50 = load i32, i32* %15
%51 = sdiv i32 %49, %50
%52 = add nsw i32 %51, 1
store i32 %52, i32* %26
br label %53
53: ; preds = %48, %43
%54 = load i32, i32* %12
store i32 %54, i32* %25
%55 = load i32, i32* %12
%56 = load i32, i32* %27
%57 = sub nsw i32 %55, %56
store i32 %57, i32* %24
%58 = load i32, i32* %12
%59 = sub nsw i32 %58, 2
store i32 %59, i32* %31
%60 = load i32, i32* %16
%61 = zext i32 %60 to i64
%62 = call i8* @llvm.stacksave()
store i8* %62, i8** %35
%63 = alloca double, i64 %61
store i64 %61, i64* %36
%64 = load double**, double*** %14
%65 = load i32, i32* %26
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double*, double** %64, i64 %66
%68 = load double*, double** %67
%69 = getelementptr inbounds double, double* %68, i64 0
%70 = load double, double* %69
%71 = getelementptr inbounds double, double* %63, i64 0
store double %70, double* %71
store i32 1, i32* %37
br label %72
72: ; preds = %95, %53
%73 = load i32, i32* %37
%74 = load i32, i32* %16
%75 = icmp slt i32 %73, %74
br i1 %75, label %76, label %98
76: ; preds = %72
%77 = load i32, i32* %37
%78 = sub nsw i32 %77, 1
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds double, double* %63, i64 %79
%81 = load double, double* %80
%82 = load double**, double*** %14
%83 = load i32, i32* %26
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double*, double** %82, i64 %84
%86 = load double*, double** %85
%87 = load i32, i32* %37
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds double, double* %86, i64 %88
%90 = load double, double* %89
%91 = fadd double %81, %90
%92 = load i32, i32* %37
%93 = sext i32 %92 to i64
%94 = getelementptr inbounds double, double* %63, i64 %93
store double %91, double* %94
br label %95
95: ; preds = %76
%96 = load i32, i32* %37
%97 = add nsw i32 %96, 1
store i32 %97, i32* %37
br label %72
98: ; preds = %72
%99 = load i32, i32* %16
%100 = call i32 @sample(i32 0, double* %63, i32 %99)
%101 = load i32*, i32** %19
%102 = load i32, i32* %25
%103 = sub nsw i32 %102, 1
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds i32, i32* %101, i64 %104
store i32 %100, i32* %105
br label %106
106: ; preds = %222, %98
%107 = load i32, i32* %24
%108 = icmp sge i32 %107, 0
br i1 %108, label %109, label %233
109: ; preds = %106
%110 = load i32, i32* %26
%111 = add nsw i32 %110, -1
store i32 %111, i32* %26
%112 = load i32*, i32** %19
%113 = load i32, i32* %25
%114 = sub nsw i32 %113, 1
%115 = sext i32 %114 to i64
%116 = getelementptr inbounds i32, i32* %112, i64 %115
%117 = load i32, i32* %116
store i32 %117, i32* %29
%118 = load i32, i32* %24
%119 = icmp sgt i32 %118, 0
br i1 %119, label %120, label %149
120: ; preds = %109
%121 = load i32, i32* %24
%122 = sub nsw i32 %121, 1
store i32 %122, i32* %30
%123 = load i32*, i32** %20
%124 = load i32, i32* %24
%125 = sext i32 %124 to i64
%126 = getelementptr inbounds i32, i32* %123, i64 %125
%127 = load i32, i32* %126
store i32 %127, i32* %28
%128 = load i32*, i32** %21
%129 = load i32, i32* %24
%130 = sext i32 %129 to i64
%131 = getelementptr inbounds i32, i32* %128, i64 %130
%132 = load i32, i32* %131
store i32 %132, i32* %32
%133 = load double***, double**** %22
%134 = load i32, i32* %26
%135 = add nsw i32 %134, 1
%136 = sext i32 %135 to i64
%137 = getelementptr inbounds double**, double*** %133, i64 %136
%138 = load double**, double*** %137
%139 = load i32, i32* %29
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds double*, double** %138, i64 %140
%142 = load double*, double** %141
%143 = load i32, i32* %16
%144 = call i32 @samplebinsearch(i32 0, double* %142, i32 %143)
%145 = load i32*, i32** %19
%146 = load i32, i32* %30
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds i32, i32* %145, i64 %147
store i32 %144, i32* %148
br label %172
149: ; preds = %109
store i32 0, i32* %30
%150 = load i32*, i32** %20
%151 = getelementptr inbounds i32, i32* %150, i64 1
%152 = load i32, i32* %151
store i32 %152, i32* %28
%153 = load i32*, i32** %21
%154 = getelementptr inbounds i32, i32* %153, i64 1
%155 = load i32, i32* %154
store i32 %155, i32* %32
%156 = load double***, double**** %22
%157 = load i32, i32* %26
%158 = add nsw i32 %157, 1
%159 = sext i32 %158 to i64
%160 = getelementptr inbounds double**, double*** %156, i64 %159
%161 = load double**, double*** %160
%162 = load i32, i32* %29
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds double*, double** %161, i64 %163
%165 = load double*, double** %164
%166 = load i32, i32* %16
%167 = call i32 @samplebinsearch(i32 0, double* %165, i32 %166)
%168 = load i32*, i32** %19
%169 = load i32, i32* %30
%170 = sext i32 %169 to i64
%171 = getelementptr inbounds i32, i32* %168, i64 %170
store i32 %167, i32* %171
br label %172
172: ; preds = %149, %120
br label %173
173: ; preds = %219, %172
%174 = load i32, i32* %30
%175 = load i32, i32* %31
%176 = icmp slt i32 %174, %175
br i1 %176, label %177, label %222
177: ; preds = %173
%178 = load double*****, double****** %18
%179 = load i32, i32* %28
%180 = sext i32 %179 to i64
%181 = getelementptr inbounds double****, double***** %178, i64 %180
%182 = load double****, double***** %181
%183 = load i32, i32* %32
%184 = sext i32 %183 to i64
%185 = getelementptr inbounds double***, double**** %182, i64 %184
%186 = load double***, double**** %185
%187 = load i32*, i32** %19
%188 = load i32, i32* %30
%189 = sext i32 %188 to i64
%190 = getelementptr inbounds i32, i32* %187, i64 %189
%191 = load i32, i32* %190
%192 = sext i32 %191 to i64
%193 = getelementptr inbounds double**, double*** %186, i64 %192
%194 = load double**, double*** %193
%195 = load i32, i32* %29
%196 = sext i32 %195 to i64
%197 = getelementptr inbounds double*, double** %194, i64 %196
%198 = load double*, double** %197
store double* %198, double** %34
%199 = load double*, double** %34
%200 = load i32, i32* %16
%201 = call i32 @samplebinsearch(i32 0, double* %199, i32 %200)
%202 = load i32*, i32** %19
%203 = load i32, i32* %30
%204 = add nsw i32 %203, 1
%205 = sext i32 %204 to i64
%206 = getelementptr inbounds i32, i32* %202, i64 %205
store i32 %201, i32* %206
%207 = load i32*, i32** %20
%208 = load i32, i32* %30
%209 = add nsw i32 %208, 2
%210 = sext i32 %209 to i64
%211 = getelementptr inbounds i32, i32* %207, i64 %210
%212 = load i32, i32* %211
store i32 %212, i32* %28
%213 = load i32*, i32** %21
%214 = load i32, i32* %30
%215 = add nsw i32 %214, 2
%216 = sext i32 %215 to i64
%217 = getelementptr inbounds i32, i32* %213, i64 %216
%218 = load i32, i32* %217
store i32 %218, i32* %32
br label %219
219: ; preds = %177
%220 = load i32, i32* %30
%221 = add nsw i32 %220, 1
store i32 %221, i32* %30
br label %173
222: ; preds = %173
%223 = load i32, i32* %15
%224 = load i32, i32* %24
%225 = sub nsw i32 %224, %223
store i32 %225, i32* %24
%226 = load i32, i32* %24
%227 = load i32, i32* %15
%228 = add nsw i32 %226, %227
store i32 %228, i32* %25
%229 = load i32, i32* %15
%230 = load i32, i32* %24
%231 = add nsw i32 %229, %230
%232 = sub nsw i32 %231, 2
store i32 %232, i32* %31
br label %106
233: ; preds = %106
%234 = load i8*, i8** %35
ret void
}

LLVM IR (-O2):
define void @csamplestatepathH(i32 %0, i32* nocapture readnone %1, double** nocapture %2, i32 %3, i32 %4, double**** nocapture readnone %5, double***** nocapture %6, i32* nocapture %7, i32* nocapture %8, i32* nocapture %9, double*** nocapture %10) {
%12 = srem i32 %0, %3
%13 = icmp eq i32 %12, 0
%14 = sdiv i32 %0, %3
%15 = select i1 %13, i32 %3, i32 %12
%16 = xor i1 %13, true
%17 = zext i1 %16 to i32
%18 = add nsw i32 %14, %17
%19 = sub nsw i32 %0, %15
%20 = zext i32 %4 to i64
%21 = alloca double, i64 %20
%22 = sext i32 %18 to i64
%23 = getelementptr inbounds double*, double** %2, i64 %22
%24 = load double*, double** %23
%25 = load double, double* %24
store double %25, double* %21
%26 = icmp sgt i32 %4, 1
br i1 %26, label %27, label %75
27: ; preds = %11
%28 = add nsw i64 %20, -1
%29 = add nsw i64 %20, -2
%30 = and i64 %28, 3
%31 = icmp ult i64 %29, 3
br i1 %31, label %60, label %32
32: ; preds = %27
%33 = and i64 %28, -4
br label %34
34: ; preds = %34, %32
%35 = phi double [ %25, %32 ], [ %55, %34 ]
%36 = phi i64 [ 1, %32 ], [ %57, %34 ]
%37 = phi i64 [ 0, %32 ], [ %58, %34 ]
%38 = getelementptr inbounds double, double* %24, i64 %36
%39 = load double, double* %38
%40 = fadd double %35, %39
%41 = getelementptr inbounds double, double* %21, i64 %36
store double %40, double* %41
%42 = add nuw nsw i64 %36, 1
%43 = getelementptr inbounds double, double* %24, i64 %42
%44 = load double, double* %43
%45 = fadd double %40, %44
%46 = getelementptr inbounds double, double* %21, i64 %42
store double %45, double* %46
%47 = add nuw nsw i64 %36, 2
%48 = getelementptr inbounds double, double* %24, i64 %47
%49 = load double, double* %48
%50 = fadd double %45, %49
%51 = getelementptr inbounds double, double* %21, i64 %47
store double %50, double* %51
%52 = add nuw nsw i64 %36, 3
%53 = getelementptr inbounds double, double* %24, i64 %52
%54 = load double, double* %53
%55 = fadd double %50, %54
%56 = getelementptr inbounds double, double* %21, i64 %52
store double %55, double* %56
%57 = add nuw nsw i64 %36, 4
%58 = add i64 %37, 4
%59 = icmp eq i64 %58, %33
br i1 %59, label %60, label %34
60: ; preds = %34, %27
%61 = phi double [ %25, %27 ], [ %55, %34 ]
%62 = phi i64 [ 1, %27 ], [ %57, %34 ]
%63 = icmp eq i64 %30, 0
br i1 %63, label %75, label %64
64: ; preds = %60, %64
%65 = phi double [ %70, %64 ], [ %61, %60 ]
%66 = phi i64 [ %72, %64 ], [ %62, %60 ]
%67 = phi i64 [ %73, %64 ], [ 0, %60 ]
%68 = getelementptr inbounds double, double* %24, i64 %66
%69 = load double, double* %68
%70 = fadd double %65, %69
%71 = getelementptr inbounds double, double* %21, i64 %66
store double %70, double* %71
%72 = add nuw nsw i64 %66, 1
%73 = add i64 %67, 1
%74 = icmp eq i64 %73, %30
br i1 %74, label %75, label %64
75: ; preds = %60, %64, %11
%76 = call i32 @sample(i32 0, double* %21, i32 %4)
%77 = add nsw i32 %0, -1
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds i32, i32* %7, i64 %78
store i32 %76, i32* %79
%80 = icmp sgt i32 %19, -1
br i1 %80, label %81, label %160
81: ; preds = %75
%82 = getelementptr inbounds i32, i32* %8, i64 1
%83 = getelementptr inbounds i32, i32* %9, i64 1
br label %84
84: ; preds = %157, %81
%85 = phi i64 [ %22, %81 ], [ %92, %157 ]
%86 = phi i32 [ %19, %81 ], [ %158, %157 ]
%87 = phi i32 [ %0, %81 ], [ %86, %157 ]
%88 = call i32 @llvm.umax.i32(i32 %86, i32 1)
%89 = add i32 %88, -1
%90 = sext i32 %89 to i64
%91 = add nsw i32 %87, -2
%92 = add i64 %85, -1
%93 = add nsw i32 %87, -1
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds i32, i32* %7, i64 %94
%96 = load i32, i32* %95
%97 = icmp eq i32 %86, 0
br i1 %97, label %113, label %98
98: ; preds = %84
%99 = add nsw i32 %86, -1
%100 = zext i32 %86 to i64
%101 = getelementptr inbounds i32, i32* %8, i64 %100
%102 = load i32, i32* %101
%103 = getelementptr inbounds i32, i32* %9, i64 %100
%104 = load i32, i32* %103
%105 = getelementptr inbounds double**, double*** %10, i64 %85
%106 = load double**, double*** %105
%107 = sext i32 %96 to i64
%108 = getelementptr inbounds double*, double** %106, i64 %107
%109 = load double*, double** %108
%110 = call i32 @samplebinsearch(i32 0, double* %109, i32 %4)
%111 = zext i32 %99 to i64
%112 = getelementptr inbounds i32, i32* %7, i64 %111
store i32 %110, i32* %112
br label %122
113: ; preds = %84
%114 = load i32, i32* %82
%115 = load i32, i32* %83
%116 = getelementptr inbounds double**, double*** %10, i64 %85
%117 = load double**, double*** %116
%118 = sext i32 %96 to i64
%119 = getelementptr inbounds double*, double** %117, i64 %118
%120 = load double*, double** %119
%121 = call i32 @samplebinsearch(i32 0, double* %120, i32 %4)
store i32 %121, i32* %7
br label %122
122: ; preds = %113, %98
%123 = phi i32 [ %102, %98 ], [ %114, %113 ]
%124 = phi i32 [ %99, %98 ], [ 0, %113 ]
%125 = phi i32 [ %104, %98 ], [ %115, %113 ]
%126 = icmp slt i32 %124, %91
br i1 %126, label %127, label %157
127: ; preds = %122
%128 = sext i32 %96 to i64
%129 = sext i32 %91 to i64
%130 = getelementptr inbounds i32, i32* %7, i64 %90
%131 = load i32, i32* %130
br label %132
132: ; preds = %127, %132
%133 = phi i32 [ %131, %127 ], [ %148, %132 ]
%134 = phi i64 [ %90, %127 ], [ %149, %132 ]
%135 = phi i32 [ %125, %127 ], [ %155, %132 ]
%136 = phi i32 [ %123, %127 ], [ %153, %132 ]
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds double****, double***** %6, i64 %137
%139 = load double****, double***** %138
%140 = sext i32 %135 to i64
%141 = getelementptr inbounds double***, double**** %139, i64 %140
%142 = load double***, double**** %141
%143 = sext i32 %133 to i64
%144 = getelementptr inbounds double**, double*** %142, i64 %143
%145 = load double**, double*** %144
%146 = getelementptr inbounds double*, double** %145, i64 %128
%147 = load double*, double** %146
%148 = call i32 @samplebinsearch(i32 0, double* %147, i32 %4)
%149 = add nsw i64 %134, 1
%150 = getelementptr inbounds i32, i32* %7, i64 %149
store i32 %148, i32* %150
%151 = add nsw i64 %134, 2
%152 = getelementptr inbounds i32, i32* %8, i64 %151
%153 = load i32, i32* %152
%154 = getelementptr inbounds i32, i32* %9, i64 %151
%155 = load i32, i32* %154
%156 = icmp slt i64 %149, %129
br i1 %156, label %132, label %157
157: ; preds = %132, %122
%158 = sub nsw i32 %86, %3
%159 = icmp sgt i32 %158, -1
br i1 %159, label %84, label %160
160: ; preds = %157, %75
ret void
}