#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(double *, int);
double exp(double);
double normalrng();
double pow(int, double);
int /*<<< orphan*/ self_adj_index_search(int *, double, double *, int);
double uniformrng();
void Mutation_Metropolis(double *z, double *fz, int N_dimension, double *theta,
                         double *grid_points, int grid_size, double temp,
                         double scl, double *accpr_pop, double *z_new) {
  int k_old;
  int k_new;
  int i;
  double fz_new;
  double rat;
  double un;
  self_adj_index_search(&k_old, *fz, grid_points, grid_size);
  for (i = 1; i <= N_dimension; i++)
    z_new[i] = z[i] + normalrng() * scl * pow(k_old / (grid_size + 1.0), 0.5);
  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;
  *accpr_pop = ((rat > 0.0) ? 1.0 : exp(rat));
  un = uniformrng();
  if (*accpr_pop > un) {
    *fz = fz_new;
    for (i = 1; i <= N_dimension; i++)
      z[i] = z_new[i];
  }
}


x86 assembly (clang, -O0):
<Mutation_Metropolis>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
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)
movsd  %xmm1,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x28(%rbp),%rsi
mov    -0x2c(%rbp),%edx
lea    -0x44(%rbp),%rdi
call   1050 <self_adj_index_search@plt>
movl   $0x1,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     123d <Mutation_Metropolis+0xdd>
mov    -0x8(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    $0x0,%al
call   1070 <normalrng@plt>
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
cvtsi2sdl -0x44(%rbp),%xmm0
cvtsi2sdl -0x2c(%rbp),%xmm1
movsd  0xe0c(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
movsd  0xe00(%rip),%xmm0
call   1080 <pow@plt>
movsd  -0x78(%rbp),%xmm1
movaps %xmm0,%xmm2
movsd  -0x70(%rbp),%xmm0
mulsd  %xmm2,%xmm0
addsd  %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    11b3 <Mutation_Metropolis+0x53>
mov    0x18(%rbp),%rdi
mov    -0x14(%rbp),%esi
call   1040 <cost@plt>
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    -0x28(%rbp),%rsi
mov    -0x2c(%rbp),%edx
lea    -0x48(%rbp),%rdi
call   1050 <self_adj_index_search@plt>
mov    -0x20(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x58(%rbp),%xmm1
divsd  -0x38(%rbp),%xmm1
subsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm1
divsd  -0x38(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    12dc <Mutation_Metropolis+0x17c>
movsd  0xd2e(%rip),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    12eb <Mutation_Metropolis+0x18b>
movsd  -0x60(%rbp),%xmm0
call   1060 <exp@plt>
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
mov    0x10(%rbp),%rax
movsd  %xmm0,(%rax)
mov    $0x0,%al
call   1030 <uniformrng@plt>
movsd  %xmm0,-0x68(%rbp)
mov    0x10(%rbp),%rax
movsd  (%rax),%xmm0
ucomisd -0x68(%rbp),%xmm0
jbe    1364 <Mutation_Metropolis+0x204>
movsd  -0x58(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     135f <Mutation_Metropolis+0x1ff>
mov    0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    132b <Mutation_Metropolis+0x1cb>
jmp    1364 <Mutation_Metropolis+0x204>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Mutation_Metropolis>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
movsd  %xmm1,0x48(%rsp)
movsd  %xmm0,0x30(%rsp)
mov    %r9d,%r12d
mov    %rcx,0x28(%rsp)
mov    %edx,%ebp
mov    %rdi,%rbx
mov    0x98(%rsp),%r14
mov    %rsi,0x38(%rsp)
movsd  (%rsi),%xmm0
lea    0x4(%rsp),%rdi
mov    %r8,0x20(%rsp)
mov    %r8,%rsi
mov    %r9d,%edx
call   1050 <self_adj_index_search@plt>
mov    %rbp,0x18(%rsp)
test   %ebp,%ebp
jle    123a <Mutation_Metropolis+0xda>
xorps  %xmm0,%xmm0
cvtsi2sd %r12d,%xmm0
addsd  0xe36(%rip),%xmm0
movsd  %xmm0,0x40(%rsp)
mov    0x18(%rsp),%rax
lea    0x1(%rax),%r13d
mov    $0x1,%r15d
nop
movsd  (%rbx,%r15,8),%xmm0
movsd  %xmm0,0x8(%rsp)
xor    %eax,%eax
call   1070 <normalrng@plt>
xorps  %xmm1,%xmm1
cvtsi2sdl 0x4(%rsp),%xmm1
mulsd  0x48(%rsp),%xmm0
movsd  %xmm0,0x50(%rsp)
divsd  0x40(%rsp),%xmm1
cvttsd2si %xmm1,%edi
movsd  0xdee(%rip),%xmm0
call   1080 <pow@plt>
mulsd  0x50(%rsp),%xmm0
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,(%r14,%r15,8)
add    $0x1,%r15
cmp    %r15,%r13
jne    11e0 <Mutation_Metropolis+0x80>
mov    0x90(%rsp),%r13
mov    %r14,%rdi
mov    0x18(%rsp),%r15
mov    %r15d,%esi
call   1040 <cost@plt>
movsd  %xmm0,0x8(%rsp)
lea    0x14(%rsp),%rdi
mov    0x20(%rsp),%rsi
mov    %r12d,%edx
call   1050 <self_adj_index_search@plt>
movslq 0x14(%rsp),%rax
mov    0x28(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm1
xorpd  0xd8f(%rip),%xmm1
movsd  0x8(%rsp),%xmm0
movsd  0x30(%rsp),%xmm2
divsd  %xmm2,%xmm0
subsd  %xmm0,%xmm1
movslq 0x4(%rsp),%rax
addsd  (%rcx,%rax,8),%xmm1
mov    0x38(%rsp),%rbp
movsd  0x0(%rbp),%xmm0
divsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    12c5 <Mutation_Metropolis+0x165>
movsd  0xd3d(%rip),%xmm0
jmp    12ca <Mutation_Metropolis+0x16a>
call   1060 <exp@plt>
movsd  %xmm0,0x0(%r13)
xor    %eax,%eax
call   1030 <uniformrng@plt>
movsd  0x0(%r13),%xmm1
ucomisd %xmm0,%xmm1
jbe    147b <Mutation_Metropolis+0x31b>
movsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x0(%rbp)
test   %r15d,%r15d
jle    147b <Mutation_Metropolis+0x31b>
mov    %r15d,%edx
lea    0x1(%rdx),%rax
mov    $0x1,%ecx
cmp    $0x4,%r15d
jb     1407 <Mutation_Metropolis+0x2a7>
lea    0x8(%rbx),%rsi
lea    (%r14,%rax,8),%rdi
cmp    %rdi,%rsi
jae    132f <Mutation_Metropolis+0x1cf>
lea    (%rbx,%rax,8),%rsi
lea    0x8(%r14),%rdi
cmp    %rsi,%rdi
jb     1407 <Mutation_Metropolis+0x2a7>
mov    %edx,%ecx
and    $0xfffffffc,%ecx
lea    -0x4(%rcx),%rdi
mov    %rdi,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %esi,%ebp
and    $0x3,%ebp
cmp    $0xc,%rdi
jae    1352 <Mutation_Metropolis+0x1f2>
xor    %edi,%edi
jmp    13c6 <Mutation_Metropolis+0x266>
and    $0xfffffffffffffffc,%rsi
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
movups 0x8(%r14,%rdi,8),%xmm0
movups 0x18(%r14,%rdi,8),%xmm1
movups %xmm0,0x8(%rbx,%rdi,8)
movups %xmm1,0x18(%rbx,%rdi,8)
movups 0x28(%r14,%rdi,8),%xmm0
movups 0x38(%r14,%rdi,8),%xmm1
movups %xmm0,0x28(%rbx,%rdi,8)
movups %xmm1,0x38(%rbx,%rdi,8)
movups 0x48(%r14,%rdi,8),%xmm0
movups 0x58(%r14,%rdi,8),%xmm1
movups %xmm0,0x48(%rbx,%rdi,8)
movups %xmm1,0x58(%rbx,%rdi,8)
movupd 0x68(%r14,%rdi,8),%xmm0
movupd 0x78(%r14,%rdi,8),%xmm1
movupd %xmm0,0x68(%rbx,%rdi,8)
movupd %xmm1,0x78(%rbx,%rdi,8)
add    $0x10,%rdi
add    $0xfffffffffffffffc,%rsi
jne    1360 <Mutation_Metropolis+0x200>
test   %rbp,%rbp
je     13fe <Mutation_Metropolis+0x29e>
nopl   0x0(%rax,%rax,1)
lea    0x0(,%rdi,8),%rsi
or     $0x8,%rsi
movupd (%r14,%rsi,1),%xmm0
movupd 0x10(%r14,%rsi,1),%xmm1
movupd %xmm0,(%rbx,%rsi,1)
movupd %xmm1,0x10(%rbx,%rsi,1)
add    $0x4,%rdi
add    $0xffffffffffffffff,%rbp
jne    13d0 <Mutation_Metropolis+0x270>
cmp    %rcx,%rdx
je     147b <Mutation_Metropolis+0x31b>
or     $0x1,%rcx
mov    %eax,%esi
sub    %ecx,%esi
mov    %rcx,%rdx
not    %rdx
add    %rax,%rdx
and    $0x3,%rsi
je     1435 <Mutation_Metropolis+0x2d5>
nopw   0x0(%rax,%rax,1)
movsd  (%r14,%rcx,8),%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    1420 <Mutation_Metropolis+0x2c0>
cmp    $0x3,%rdx
jb     147b <Mutation_Metropolis+0x31b>
nopl   0x0(%rax,%rax,1)
movsd  (%r14,%rcx,8),%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
movsd  0x8(%r14,%rcx,8),%xmm0
movsd  %xmm0,0x8(%rbx,%rcx,8)
movsd  0x10(%r14,%rcx,8),%xmm0
movsd  %xmm0,0x10(%rbx,%rcx,8)
movsd  0x18(%r14,%rcx,8),%xmm0
movsd  %xmm0,0x18(%rbx,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    1440 <Mutation_Metropolis+0x2e0>
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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