#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 */
double cost(int *, int);
double exp(double);
int integerrng(int, int);
int /*<<< orphan*/ self_adj_index_search(int *, double, double *, int);
double uniformrng();
void Mutation_int_Gibbs(int *z, double *fz, int N_dimension, double *theta,
                        double *grid_points, int grid_size, double temp,
                        int Krep, int *z_new) {
  int k_old;
  int k_new;
  int i;
  double Pr;
  int lab;
  int gg;
  double fz_new;
  double rat;
  double un;
  for (i = 1; i <= N_dimension; i++)
    z_new[i] = z[i];
  for (gg = 1; gg <= Krep; gg++) {
    lab = integerrng(1, N_dimension);
    self_adj_index_search(&k_old, *fz, grid_points, grid_size);
    z_new[lab] = 1 - z_new[lab];
    fz_new = cost(z_new, N_dimension);
    self_adj_index_search(&k_new, fz_new, grid_points, grid_size);
    rat = theta[k_new] + fz_new / temp - theta[k_old] - *fz / temp;
    Pr = 1.0 / (1.0 + exp(rat));
    un = uniformrng();
    if (Pr >= un) {
      *fz = fz_new;
      z[lab] = z_new[lab];
    } else {
      z_new[lab] = 1 - z_new[lab];
    }
  }
}


x86 assembly (clang, -O0):
<Mutation_int_Gibbs>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
movsd  %xmm0,-0x38(%rbp)
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     11b2 <Mutation_int_Gibbs+0x62>
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1182 <Mutation_int_Gibbs+0x32>
movl   $0x1,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    0x10(%rbp),%eax
jg     1305 <Mutation_int_Gibbs+0x1b5>
mov    -0x14(%rbp),%esi
mov    $0x1,%edi
call   1060 <integerrng@plt>
mov    %eax,-0x54(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x28(%rbp),%rsi
mov    -0x2c(%rbp),%edx
lea    -0x3c(%rbp),%rdi
call   1050 <self_adj_index_search@plt>
mov    0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    $0x1,%edx
sub    (%rax,%rcx,4),%edx
mov    0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    0x18(%rbp),%rdi
mov    -0x14(%rbp),%esi
call   1040 <cost@plt>
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
mov    -0x28(%rbp),%rsi
mov    -0x2c(%rbp),%edx
lea    -0x40(%rbp),%rdi
call   1050 <self_adj_index_search@plt>
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x60(%rbp),%xmm1
divsd  -0x38(%rbp),%xmm1
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm1
divsd  -0x38(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
call   1070 <exp@plt>
movsd  0xd82(%rip),%xmm1
addsd  %xmm0,%xmm1
movsd  0xd76(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    $0x0,%al
call   1030 <uniformrng@plt>
movsd  %xmm0,-0x70(%rbp)
movsd  -0x50(%rbp),%xmm0
ucomisd -0x70(%rbp),%xmm0
jb     12d7 <Mutation_int_Gibbs+0x187>
movsd  -0x60(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
mov    0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x8(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
jmp    12f2 <Mutation_int_Gibbs+0x1a2>
mov    0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    $0x1,%edx
sub    (%rax,%rcx,4),%edx
mov    0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
jmp    12f7 <Mutation_int_Gibbs+0x1a7>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    11b9 <Mutation_int_Gibbs+0x69>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Mutation_int_Gibbs>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
movsd  %xmm0,0x38(%rsp)
mov    %r9d,0xc(%rsp)
mov    %r8,0x30(%rsp)
mov    %rcx,0x28(%rsp)
mov    %rsi,%r12
mov    0x88(%rsp),%r9
mov    0x80(%rsp),%eax
mov    %edx,0x4(%rsp)
test   %edx,%edx
jle    12fb <Mutation_int_Gibbs+0x1ab>
mov    0x4(%rsp),%esi
mov    %esi,%edx
lea    0x1(%rdx),%r8
mov    $0x1,%ecx
cmp    $0x8,%esi
jb     1297 <Mutation_int_Gibbs+0x147>
lea    0x4(%r9),%rsi
lea    (%rdi,%r8,4),%rbp
cmp    %rbp,%rsi
jae    11c7 <Mutation_int_Gibbs+0x77>
lea    (%r9,%r8,4),%rsi
lea    0x4(%rdi),%rbp
cmp    %rsi,%rbp
jb     1297 <Mutation_int_Gibbs+0x147>
mov    %edx,%ecx
and    $0xfffffff8,%ecx
lea    -0x8(%rcx),%rbx
mov    %rbx,%rsi
shr    $0x3,%rsi
add    $0x1,%rsi
mov    %esi,%ebp
and    $0x3,%ebp
cmp    $0x18,%rbx
jae    11ea <Mutation_int_Gibbs+0x9a>
xor    %ebx,%ebx
jmp    1256 <Mutation_int_Gibbs+0x106>
and    $0xfffffffffffffffc,%rsi
xor    %ebx,%ebx
movups 0x4(%rdi,%rbx,4),%xmm0
movups 0x14(%rdi,%rbx,4),%xmm1
movups %xmm0,0x4(%r9,%rbx,4)
movups %xmm1,0x14(%r9,%rbx,4)
movups 0x24(%rdi,%rbx,4),%xmm0
movups 0x34(%rdi,%rbx,4),%xmm1
movups %xmm0,0x24(%r9,%rbx,4)
movups %xmm1,0x34(%r9,%rbx,4)
movups 0x44(%rdi,%rbx,4),%xmm0
movups 0x54(%rdi,%rbx,4),%xmm1
movups %xmm0,0x44(%r9,%rbx,4)
movups %xmm1,0x54(%r9,%rbx,4)
movupd 0x64(%rdi,%rbx,4),%xmm0
movupd 0x74(%rdi,%rbx,4),%xmm1
movupd %xmm0,0x64(%r9,%rbx,4)
movupd %xmm1,0x74(%r9,%rbx,4)
add    $0x20,%rbx
add    $0xfffffffffffffffc,%rsi
jne    11f0 <Mutation_int_Gibbs+0xa0>
test   %rbp,%rbp
je     128e <Mutation_int_Gibbs+0x13e>
nopl   0x0(%rax,%rax,1)
lea    0x0(,%rbx,4),%rsi
or     $0x4,%rsi
movupd (%rdi,%rsi,1),%xmm0
movupd 0x10(%rdi,%rsi,1),%xmm1
movupd %xmm0,(%r9,%rsi,1)
movupd %xmm1,0x10(%r9,%rsi,1)
add    $0x8,%rbx
add    $0xffffffffffffffff,%rbp
jne    1260 <Mutation_int_Gibbs+0x110>
cmp    %rcx,%rdx
je     12fb <Mutation_int_Gibbs+0x1ab>
or     $0x1,%rcx
mov    %r8d,%esi
sub    %ecx,%esi
mov    %rcx,%rdx
not    %rdx
add    %r8,%rdx
and    $0x3,%rsi
je     12c1 <Mutation_int_Gibbs+0x171>
nopl   0x0(%rax,%rax,1)
mov    (%rdi,%rcx,4),%ebp
mov    %ebp,(%r9,%rcx,4)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    12b0 <Mutation_int_Gibbs+0x160>
cmp    $0x3,%rdx
jb     12fb <Mutation_int_Gibbs+0x1ab>
nopw   0x0(%rax,%rax,1)
mov    (%rdi,%rcx,4),%edx
mov    %edx,(%r9,%rcx,4)
mov    0x4(%rdi,%rcx,4),%edx
mov    %edx,0x4(%r9,%rcx,4)
mov    0x8(%rdi,%rcx,4),%edx
mov    %edx,0x8(%r9,%rcx,4)
mov    0xc(%rdi,%rcx,4),%edx
mov    %edx,0xc(%r9,%rcx,4)
add    $0x4,%rcx
cmp    %rcx,%r8
jne    12d0 <Mutation_int_Gibbs+0x180>
test   %eax,%eax
jle    1421 <Mutation_int_Gibbs+0x2d1>
mov    %rdi,0x20(%rsp)
mov    0x4(%rsp),%ebx
mov    0x30(%rsp),%r14
mov    0xc(%rsp),%r13d
jmp    133a <Mutation_int_Gibbs+0x1ea>
nopl   0x0(%rax,%rax,1)
mov    $0x1,%eax
sub    (%r15,%rbp,4),%eax
mov    %eax,(%r15,%rbp,4)
mov    0x8(%rsp),%eax
add    $0xffffffff,%eax
je     1421 <Mutation_int_Gibbs+0x2d1>
mov    %eax,0x8(%rsp)
mov    $0x1,%edi
mov    %ebx,%esi
call   1060 <integerrng@plt>
mov    %eax,%r15d
movsd  (%r12),%xmm0
lea    0x14(%rsp),%rdi
mov    %r14,%rsi
mov    %r13d,%edx
call   1050 <self_adj_index_search@plt>
movslq %r15d,%rbp
mov    $0x1,%eax
mov    0x88(%rsp),%r15
sub    (%r15,%rbp,4),%eax
mov    %eax,(%r15,%rbp,4)
mov    %r15,%rdi
mov    %ebx,%esi
call   1040 <cost@plt>
movsd  %xmm0,0x18(%rsp)
lea    0x10(%rsp),%rdi
mov    %r14,%rsi
mov    %r13d,%edx
call   1050 <self_adj_index_search@plt>
movslq 0x10(%rsp),%rax
movsd  0x18(%rsp),%xmm0
movsd  0x38(%rsp),%xmm2
divsd  %xmm2,%xmm0
mov    0x28(%rsp),%rcx
addsd  (%rcx,%rax,8),%xmm0
movslq 0x14(%rsp),%rax
subsd  (%rcx,%rax,8),%xmm0
movsd  (%r12),%xmm1
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
call   1070 <exp@plt>
movsd  0xc21(%rip),%xmm1
addsd  %xmm1,%xmm0
divsd  %xmm0,%xmm1
movsd  %xmm1,0x40(%rsp)
xor    %eax,%eax
call   1030 <uniformrng@plt>
movsd  0x40(%rsp),%xmm1
ucomisd %xmm0,%xmm1
jb     1320 <Mutation_int_Gibbs+0x1d0>
movsd  0x18(%rsp),%xmm0
movsd  %xmm0,(%r12)
mov    (%r15,%rbp,4),%eax
mov    0x20(%rsp),%rdi
mov    %eax,(%rdi,%rbp,4)
jmp    132d <Mutation_int_Gibbs+0x1dd>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @Mutation_int_Gibbs(i32* %0, double* %1, i32 %2, double* %3, double* %4, i32 %5, double %6, i32 %7, i32* %8) {
%10 = alloca i32*
%11 = alloca double*
%12 = alloca i32
%13 = alloca double*
%14 = alloca double*
%15 = alloca i32
%16 = alloca double
%17 = alloca i32
%18 = alloca i32*
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca double
%23 = alloca i32
%24 = alloca i32
%25 = alloca double
%26 = alloca double
%27 = alloca double
store i32* %0, i32** %10
store double* %1, double** %11
store i32 %2, i32* %12
store double* %3, double** %13
store double* %4, double** %14
store i32 %5, i32* %15
store double %6, double* %16
store i32 %7, i32* %17
store i32* %8, i32** %18
store i32 1, i32* %21
br label %28
28: ; preds = %42, %9
%29 = load i32, i32* %21
%30 = load i32, i32* %12
%31 = icmp sle i32 %29, %30
br i1 %31, label %32, label %45
32: ; preds = %28
%33 = load i32*, i32** %10
%34 = load i32, i32* %21
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds i32, i32* %33, i64 %35
%37 = load i32, i32* %36
%38 = load i32*, i32** %18
%39 = load i32, i32* %21
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds i32, i32* %38, i64 %40
store i32 %37, i32* %41
br label %42
42: ; preds = %32
%43 = load i32, i32* %21
%44 = add nsw i32 %43, 1
store i32 %44, i32* %21
br label %28
45: ; preds = %28
store i32 1, i32* %24
br label %46
46: ; preds = %127, %45
%47 = load i32, i32* %24
%48 = load i32, i32* %17
%49 = icmp sle i32 %47, %48
br i1 %49, label %50, label %130
50: ; preds = %46
%51 = load i32, i32* %12
%52 = call i32 @integerrng(i32 1, i32 %51)
store i32 %52, i32* %23
%53 = load double*, double** %11
%54 = load double, double* %53
%55 = load double*, double** %14
%56 = load i32, i32* %15
%57 = call i32 @self_adj_index_search(i32* %19, double %54, double* %55, i32 %56)
%58 = load i32*, i32** %18
%59 = load i32, i32* %23
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds i32, i32* %58, i64 %60
%62 = load i32, i32* %61
%63 = sub nsw i32 1, %62
%64 = load i32*, i32** %18
%65 = load i32, i32* %23
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds i32, i32* %64, i64 %66
store i32 %63, i32* %67
%68 = load i32*, i32** %18
%69 = load i32, i32* %12
%70 = call double @cost(i32* %68, i32 %69)
store double %70, double* %25
%71 = load double, double* %25
%72 = load double*, double** %14
%73 = load i32, i32* %15
%74 = call i32 @self_adj_index_search(i32* %20, double %71, double* %72, i32 %73)
%75 = load double*, double** %13
%76 = load i32, i32* %20
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %75, i64 %77
%79 = load double, double* %78
%80 = load double, double* %25
%81 = load double, double* %16
%82 = fdiv double %80, %81
%83 = fadd double %79, %82
%84 = load double*, double** %13
%85 = load i32, i32* %19
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %84, i64 %86
%88 = load double, double* %87
%89 = fsub double %83, %88
%90 = load double*, double** %11
%91 = load double, double* %90
%92 = load double, double* %16
%93 = fdiv double %91, %92
%94 = fsub double %89, %93
store double %94, double* %26
%95 = load double, double* %26
%96 = call double @exp(double %95)
%97 = fadd double 1.000000e+00, %96
%98 = fdiv double 1.000000e+00, %97
store double %98, double* %22
%99 = call double (...) @uniformrng()
store double %99, double* %27
%100 = load double, double* %22
%101 = load double, double* %27
%102 = fcmp oge double %100, %101
br i1 %102, label %103, label %115
103: ; preds = %50
%104 = load double, double* %25
%105 = load double*, double** %11
store double %104, double* %105
%106 = load i32*, i32** %18
%107 = load i32, i32* %23
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds i32, i32* %106, i64 %108
%110 = load i32, i32* %109
%111 = load i32*, i32** %10
%112 = load i32, i32* %23
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds i32, i32* %111, i64 %113
store i32 %110, i32* %114
br label %126
115: ; preds = %50
%116 = load i32*, i32** %18
%117 = load i32, i32* %23
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds i32, i32* %116, i64 %118
%120 = load i32, i32* %119
%121 = sub nsw i32 1, %120
%122 = load i32*, i32** %18
%123 = load i32, i32* %23
%124 = sext i32 %123 to i64
%125 = getelementptr inbounds i32, i32* %122, i64 %124
store i32 %121, i32* %125
br label %126
126: ; preds = %115, %103
br label %127
127: ; preds = %126
%128 = load i32, i32* %24
%129 = add nsw i32 %128, 1
store i32 %129, i32* %24
br label %46
130: ; preds = %46
ret void
}

LLVM IR (-O2):
define void @Mutation_int_Gibbs(i32* nocapture %0, double* nocapture %1, i32 %2, double* nocapture %3, double* %4, i32 %5, double %6, i32 %7, i32* %8) {
%10 = alloca i32
%11 = alloca i32
%12 = bitcast i32* %10 to i8*
%13 = bitcast i32* %11 to i8*
%14 = icmp slt i32 %2, 1
br i1 %14, label %129, label %15
15: ; preds = %9
%16 = add nuw i32 %2, 1
%17 = zext i32 %16 to i64
%18 = add nsw i64 %17, -1
%19 = icmp ult i64 %18, 8
br i1 %19, label %110, label %20
20: ; preds = %15
%21 = getelementptr i32, i32* %8, i64 1
%22 = getelementptr i32, i32* %8, i64 %17
%23 = getelementptr i32, i32* %0, i64 1
%24 = getelementptr i32, i32* %0, i64 %17
%25 = icmp ult i32* %21, %24
%26 = icmp ult i32* %23, %22
%27 = and i1 %25, %26
br i1 %27, label %110, label %28
28: ; preds = %20
%29 = and i64 %18, -8
%30 = or i64 %29, 1
%31 = add nsw i64 %29, -8
%32 = lshr exact i64 %31, 3
%33 = add nuw nsw i64 %32, 1
%34 = and i64 %33, 3
%35 = icmp ult i64 %31, 24
br i1 %35, label %88, label %36
36: ; preds = %28
%37 = and i64 %33, 4611686018427387900
br label %38
38: ; preds = %38, %36
%39 = phi i64 [ 0, %36 ], [ %85, %38 ]
%40 = phi i64 [ 0, %36 ], [ %86, %38 ]
%41 = or i64 %39, 1
%42 = getelementptr inbounds i32, i32* %0, i64 %41
%43 = bitcast i32* %42 to <4 x i32>*
%44 = load <4 x i32>, <4 x i32>* %43
%45 = getelementptr inbounds i32, i32* %42, i64 4
%46 = bitcast i32* %45 to <4 x i32>*
%47 = load <4 x i32>, <4 x i32>* %46
%48 = getelementptr inbounds i32, i32* %8, i64 %41
%49 = bitcast i32* %48 to <4 x i32>*
store <4 x i32> %44, <4 x i32>* %49
%50 = getelementptr inbounds i32, i32* %48, i64 4
%51 = bitcast i32* %50 to <4 x i32>*
store <4 x i32> %47, <4 x i32>* %51
%52 = or i64 %39, 9
%53 = getelementptr inbounds i32, i32* %0, i64 %52
%54 = bitcast i32* %53 to <4 x i32>*
%55 = load <4 x i32>, <4 x i32>* %54
%56 = getelementptr inbounds i32, i32* %53, i64 4
%57 = bitcast i32* %56 to <4 x i32>*
%58 = load <4 x i32>, <4 x i32>* %57
%59 = getelementptr inbounds i32, i32* %8, i64 %52
%60 = bitcast i32* %59 to <4 x i32>*
store <4 x i32> %55, <4 x i32>* %60
%61 = getelementptr inbounds i32, i32* %59, i64 4
%62 = bitcast i32* %61 to <4 x i32>*
store <4 x i32> %58, <4 x i32>* %62
%63 = or i64 %39, 17
%64 = getelementptr inbounds i32, i32* %0, i64 %63
%65 = bitcast i32* %64 to <4 x i32>*
%66 = load <4 x i32>, <4 x i32>* %65
%67 = getelementptr inbounds i32, i32* %64, i64 4
%68 = bitcast i32* %67 to <4 x i32>*
%69 = load <4 x i32>, <4 x i32>* %68
%70 = getelementptr inbounds i32, i32* %8, i64 %63
%71 = bitcast i32* %70 to <4 x i32>*
store <4 x i32> %66, <4 x i32>* %71
%72 = getelementptr inbounds i32, i32* %70, i64 4
%73 = bitcast i32* %72 to <4 x i32>*
store <4 x i32> %69, <4 x i32>* %73
%74 = or i64 %39, 25
%75 = getelementptr inbounds i32, i32* %0, i64 %74
%76 = bitcast i32* %75 to <4 x i32>*
%77 = load <4 x i32>, <4 x i32>* %76
%78 = getelementptr inbounds i32, i32* %75, i64 4
%79 = bitcast i32* %78 to <4 x i32>*
%80 = load <4 x i32>, <4 x i32>* %79
%81 = getelementptr inbounds i32, i32* %8, i64 %74
%82 = bitcast i32* %81 to <4 x i32>*
store <4 x i32> %77, <4 x i32>* %82
%83 = getelementptr inbounds i32, i32* %81, i64 4
%84 = bitcast i32* %83 to <4 x i32>*
store <4 x i32> %80, <4 x i32>* %84
%85 = add nuw i64 %39, 32
%86 = add i64 %40, 4
%87 = icmp eq i64 %86, %37
br i1 %87, label %88, label %38
88: ; preds = %38, %28
%89 = phi i64 [ 0, %28 ], [ %85, %38 ]
%90 = icmp eq i64 %34, 0
br i1 %90, label %108, label %91
91: ; preds = %88, %91
%92 = phi i64 [ %105, %91 ], [ %89, %88 ]
%93 = phi i64 [ %106, %91 ], [ 0, %88 ]
%94 = or i64 %92, 1
%95 = getelementptr inbounds i32, i32* %0, i64 %94
%96 = bitcast i32* %95 to <4 x i32>*
%97 = load <4 x i32>, <4 x i32>* %96
%98 = getelementptr inbounds i32, i32* %95, i64 4
%99 = bitcast i32* %98 to <4 x i32>*
%100 = load <4 x i32>, <4 x i32>* %99
%101 = getelementptr inbounds i32, i32* %8, i64 %94
%102 = bitcast i32* %101 to <4 x i32>*
store <4 x i32> %97, <4 x i32>* %102
%103 = getelementptr inbounds i32, i32* %101, i64 4
%104 = bitcast i32* %103 to <4 x i32>*
store <4 x i32> %100, <4 x i32>* %104
%105 = add nuw i64 %92, 8
%106 = add i64 %93, 1
%107 = icmp eq i64 %106, %34
br i1 %107, label %108, label %91
108: ; preds = %91, %88
%109 = icmp eq i64 %18, %29
br i1 %109, label %129, label %110
110: ; preds = %20, %15, %108
%111 = phi i64 [ 1, %20 ], [ 1, %15 ], [ %30, %108 ]
%112 = sub nsw i64 %17, %111
%113 = xor i64 %111, -1
%114 = add nsw i64 %113, %17
%115 = and i64 %112, 3
%116 = icmp eq i64 %115, 0
br i1 %116, label %126, label %117
117: ; preds = %110, %117
%118 = phi i64 [ %123, %117 ], [ %111, %110 ]
%119 = phi i64 [ %124, %117 ], [ 0, %110 ]
%120 = getelementptr inbounds i32, i32* %0, i64 %118
%121 = load i32, i32* %120
%122 = getelementptr inbounds i32, i32* %8, i64 %118
store i32 %121, i32* %122
%123 = add nuw nsw i64 %118, 1
%124 = add i64 %119, 1
%125 = icmp eq i64 %124, %115
br i1 %125, label %126, label %117
126: ; preds = %117, %110
%127 = phi i64 [ %111, %110 ], [ %123, %117 ]
%128 = icmp ult i64 %114, 3
br i1 %128, label %129, label %131
129: ; preds = %126, %131, %108, %9
%130 = icmp slt i32 %7, 1
br i1 %130, label %189, label %150
131: ; preds = %126, %131
%132 = phi i64 [ %148, %131 ], [ %127, %126 ]
%133 = getelementptr inbounds i32, i32* %0, i64 %132
%134 = load i32, i32* %133
%135 = getelementptr inbounds i32, i32* %8, i64 %132
store i32 %134, i32* %135
%136 = add nuw nsw i64 %132, 1
%137 = getelementptr inbounds i32, i32* %0, i64 %136
%138 = load i32, i32* %137
%139 = getelementptr inbounds i32, i32* %8, i64 %136
store i32 %138, i32* %139
%140 = add nuw nsw i64 %132, 2
%141 = getelementptr inbounds i32, i32* %0, i64 %140
%142 = load i32, i32* %141
%143 = getelementptr inbounds i32, i32* %8, i64 %140
store i32 %142, i32* %143
%144 = add nuw nsw i64 %132, 3
%145 = getelementptr inbounds i32, i32* %0, i64 %144
%146 = load i32, i32* %145
%147 = getelementptr inbounds i32, i32* %8, i64 %144
store i32 %146, i32* %147
%148 = add nuw nsw i64 %132, 4
%149 = icmp eq i64 %148, %17
br i1 %149, label %129, label %131
150: ; preds = %129, %186
%151 = phi i32 [ %187, %186 ], [ 1, %129 ]
%152 = call i32 @integerrng(i32 1, i32 %2)
%153 = load double, double* %1
%154 = call i32 @self_adj_index_search(i32* %10, double %153, double* %4, i32 %5)
%155 = sext i32 %152 to i64
%156 = getelementptr inbounds i32, i32* %8, i64 %155
%157 = load i32, i32* %156
%158 = sub nsw i32 1, %157
store i32 %158, i32* %156
%159 = call double @cost(i32* %8, i32 %2)
%160 = call i32 @self_adj_index_search(i32* %11, double %159, double* %4, i32 %5)
%161 = load i32, i32* %11
%162 = sext i32 %161 to i64
%163 = getelementptr inbounds double, double* %3, i64 %162
%164 = load double, double* %163
%165 = fdiv double %159, %6
%166 = fadd double %165, %164
%167 = load i32, i32* %10
%168 = sext i32 %167 to i64
%169 = getelementptr inbounds double, double* %3, i64 %168
%170 = load double, double* %169
%171 = fsub double %166, %170
%172 = load double, double* %1
%173 = fdiv double %172, %6
%174 = fsub double %171, %173
%175 = call double @exp(double %174)
%176 = fadd double %175, 1.000000e+00
%177 = fdiv double 1.000000e+00, %176
%178 = call double (...) @uniformrng()
%179 = fcmp ult double %177, %178
br i1 %179, label %183, label %180
180: ; preds = %150
store double %159, double* %1
%181 = load i32, i32* %156
%182 = getelementptr inbounds i32, i32* %0, i64 %155
store i32 %181, i32* %182
br label %186
183: ; preds = %150
%184 = load i32, i32* %156
%185 = sub nsw i32 1, %184
store i32 %185, i32* %156
br label %186
186: ; preds = %180, %183
%187 = add nuw i32 %151, 1
%188 = icmp eq i32 %151, %7
br i1 %188, label %189, label %150
189: ; preds = %186, %129
ret void
}