#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 /*<<< orphan*/ assert(int);
int find_xb(double, double const *, int, int);
int index_2D_to_1D(int, int, int, int);
double lagrange_interp_2D(double, double, double *, int, double *, int,
                          double *);
double lagrange_interp_2Dfo(double x, double y, int order_x, int order_y,
                            const double x_values[], int num_x,
                            const double y_values[], int num_y,
                            const double z_values[]) {
  assert(order_x % 2 == 0);
  assert(order_y % 2 == 0);
  assert(num_x > order_x);
  assert(num_y > order_y);
  double mx[order_x + 1];
  double my[order_y + 1];
  double mz[(order_x + 1) * (order_y + 1)];
  int ix = find_xb(x, x_values, num_x, order_x / 2);
  int iy = find_xb(y, y_values, num_y, order_y / 2);
  int iz_small, iz_big;
  for (int i = 0; i < order_x + 1; i++) {
    mx[i] = x_values[ix + i - order_x / 2];
    for (int j = 0; j < order_y + 1; j++) {
      my[j] = y_values[iy + j - order_y / 2];
      iz_small = index_2D_to_1D(i, order_x + 1, j, order_y + 1);
      iz_big = index_2D_to_1D(ix + i - order_x / 2, num_x, iy + j - order_y / 2,
                              num_y);
      mz[iz_small] = z_values[iz_big];
    }
  }
  return lagrange_interp_2D(x, y, mx, order_x + 1, my, order_y + 1, mz);
}


x86 assembly (clang, -O0):
<lagrange_interp_2Dfo>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    0x10(%rbp),%rax
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
mov    %edi,-0x14(%rbp)
mov    %esi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %ecx,-0x24(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9d,-0x34(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,%edx
shr    $0x1f,%edx
mov    %eax,%ecx
add    %edx,%ecx
and    $0xfffffffe,%ecx
sub    %ecx,%eax
sete   %al
movzbl %al,%edi
call   1050 <assert@plt>
mov    -0x18(%rbp),%eax
mov    %eax,%edx
shr    $0x1f,%edx
mov    %eax,%ecx
add    %edx,%ecx
and    $0xfffffffe,%ecx
sub    %ecx,%eax
sete   %al
movzbl %al,%edi
call   1050 <assert@plt>
mov    -0x24(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    %ecx,%eax
setg   %al
movzbl %al,%edi
call   1050 <assert@plt>
mov    -0x34(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    %ecx,%eax
setg   %al
movzbl %al,%edi
call   1050 <assert@plt>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %rsp,%rcx
mov    %rcx,-0x40(%rbp)
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x88(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x48(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x80(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x50(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    -0x18(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x78(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x58(%rbp)
movsd  -0x8(%rbp),%xmm0
mov    -0x20(%rbp),%rdi
mov    -0x24(%rbp),%esi
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%edx
call   1040 <find_xb@plt>
mov    %eax,-0x5c(%rbp)
movsd  -0x10(%rbp),%xmm0
mov    -0x30(%rbp),%rdi
mov    -0x34(%rbp),%esi
mov    -0x18(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%edx
call   1040 <find_xb@plt>
mov    %eax,-0x60(%rbp)
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    13c8 <lagrange_interp_2Dfo+0x288>
mov    -0x20(%rbp),%rcx
mov    -0x5c(%rbp),%eax
add    -0x6c(%rbp),%eax
mov    %eax,-0x8c(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x2,%esi
cltd
idiv   %esi
mov    -0x8c(%rbp),%edx
mov    %eax,%esi
mov    -0x88(%rbp),%rax
sub    %esi,%edx
movslq %edx,%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x6c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x70(%rbp)
mov    -0x70(%rbp),%eax
mov    -0x18(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    13b5 <lagrange_interp_2Dfo+0x275>
mov    -0x30(%rbp),%rcx
mov    -0x60(%rbp),%eax
add    -0x70(%rbp),%eax
mov    %eax,-0x94(%rbp)
mov    -0x18(%rbp),%eax
mov    $0x2,%esi
cltd
idiv   %esi
mov    -0x94(%rbp),%edx
mov    %eax,%esi
mov    -0x80(%rbp),%rax
sub    %esi,%edx
movslq %edx,%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x70(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x6c(%rbp),%edi
mov    -0x14(%rbp),%esi
add    $0x1,%esi
mov    -0x70(%rbp),%edx
mov    -0x18(%rbp),%ecx
add    $0x1,%ecx
call   1060 <index_2D_to_1D@plt>
mov    %eax,-0x64(%rbp)
mov    -0x5c(%rbp),%edi
add    -0x6c(%rbp),%edi
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
sub    %eax,%edi
mov    -0x24(%rbp),%esi
mov    -0x60(%rbp),%eax
add    -0x70(%rbp),%eax
mov    %eax,-0x90(%rbp)
mov    -0x18(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    -0x90(%rbp),%edx
sub    %eax,%edx
mov    -0x34(%rbp),%ecx
call   1060 <index_2D_to_1D@plt>
mov    %eax,%ecx
mov    -0x78(%rbp),%rax
mov    %ecx,-0x68(%rbp)
mov    0x10(%rbp),%rcx
movslq -0x68(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x64(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x70(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x70(%rbp)
jmp    12e6 <lagrange_interp_2Dfo+0x1a6>
jmp    13ba <lagrange_interp_2Dfo+0x27a>
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    1291 <lagrange_interp_2Dfo+0x151>
mov    -0x78(%rbp),%r8
mov    -0x80(%rbp),%rdx
mov    -0x88(%rbp),%rdi
movsd  -0x8(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
mov    -0x14(%rbp),%esi
add    $0x1,%esi
mov    -0x18(%rbp),%ecx
add    $0x1,%ecx
call   1030 <lagrange_interp_2D@plt>
mov    -0x40(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<lagrange_interp_2Dfo>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %r9d,-0x2c(%rbp)
mov    %r8,%r13
mov    %ecx,%r14d
mov    %ecx,-0x30(%rbp)
mov    %rdx,-0x70(%rbp)
mov    %esi,%r12d
mov    %edi,%r15d
movsd  %xmm1,-0x60(%rbp)
movsd  %xmm0,-0x58(%rbp)
mov    %esi,%ebx
not    %ebx
not    %edi
and    $0x1,%edi
call   1050 <assert@plt>
and    $0x1,%ebx
mov    %ebx,%edi
mov    -0x2c(%rbp),%ebx
call   1050 <assert@plt>
xor    %edi,%edi
cmp    %r15d,%r14d
setg   %dil
call   1050 <assert@plt>
xor    %edi,%edi
cmp    %r12d,%ebx
setg   %dil
call   1050 <assert@plt>
lea    0x1(%r15),%edx
mov    %rsp,%rcx
lea    0xf(,%rdx,8),%rax
mov    %rdx,-0x80(%rbp)
and    $0xfffffffffffffff0,%rax
sub    %rax,%rcx
mov    %rcx,-0x68(%rbp)
mov    %rcx,%rsp
lea    0x1(%r12),%esi
mov    %rsi,-0x88(%rbp)
mov    %rsp,%rcx
lea    0xf(,%rsi,8),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%rcx
mov    %rcx,-0x40(%rbp)
mov    %rcx,%rsp
mov    %esi,%eax
imul   %edx,%eax
lea    0xf(,%rax,8),%rax
and    $0xfffffffffffffff0,%rax
mov    %rsp,%rcx
sub    %rax,%rcx
mov    %rcx,-0x78(%rbp)
mov    %rcx,%rsp
mov    %r15,-0x48(%rbp)
mov    %r15d,%r14d
shr    $0x1f,%r14d
add    %r15d,%r14d
sar    %r14d
movsd  -0x58(%rbp),%xmm0
mov    -0x70(%rbp),%rdi
mov    -0x30(%rbp),%esi
mov    %r14d,%edx
call   1040 <find_xb@plt>
mov    %rax,-0x38(%rbp)
mov    %r12d,%r15d
shr    $0x1f,%r15d
mov    %r12,-0x90(%rbp)
add    %r12d,%r15d
sar    %r15d
movsd  -0x60(%rbp),%xmm0
mov    %r13,%rdi
mov    %ebx,%esi
mov    %r15d,%edx
call   1040 <find_xb@plt>
mov    -0x88(%rbp),%r12
mov    -0x80(%rbp),%rsi
cmpl   $0x0,-0x48(%rbp)
js     1356 <lagrange_interp_2Dfo+0x216>
mov    %eax,%edx
mov    -0x38(%rbp),%rax
sub    %r14d,%eax
mov    %rax,-0x38(%rbp)
sub    %r15d,%edx
xor    %edi,%edi
mov    %r13,-0xa0(%rbp)
mov    %rdx,-0x98(%rbp)
jmp    12ad <lagrange_interp_2Dfo+0x16d>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rdi
cmp    %rsi,%rdi
je     1356 <lagrange_interp_2Dfo+0x216>
mov    -0x38(%rbp),%rax
add    %edi,%eax
mov    %eax,-0x4c(%rbp)
cltq
mov    -0x70(%rbp),%rcx
movsd  (%rcx,%rax,8),%xmm0
mov    -0x68(%rbp),%rax
movsd  %xmm0,(%rax,%rdi,8)
cmpl   $0x0,-0x90(%rbp)
mov    -0x40(%rbp),%rcx
js     12a0 <lagrange_interp_2Dfo+0x160>
xor    %r15d,%r15d
nopw   0x0(%rax,%rax,1)
lea    (%rdx,%r15,1),%rax
movslq %eax,%rbx
movsd  0x0(%r13,%rbx,8),%xmm0
movsd  %xmm0,(%rcx,%r15,8)
mov    %rdi,%r13
mov    %rsi,%r14
mov    %r15d,%edx
mov    %r12d,%ecx
call   1060 <index_2D_to_1D@plt>
mov    %eax,-0x48(%rbp)
mov    -0x4c(%rbp),%edi
mov    -0x30(%rbp),%esi
mov    %ebx,%edx
mov    -0x2c(%rbp),%ecx
call   1060 <index_2D_to_1D@plt>
mov    %r13,%rdi
mov    0x10(%rbp),%rcx
mov    -0x98(%rbp),%rdx
mov    %r14,%rsi
cltq
movsd  (%rcx,%rax,8),%xmm0
mov    -0x40(%rbp),%rcx
movslq -0x48(%rbp),%rax
mov    -0xa0(%rbp),%r13
mov    -0x78(%rbp),%rbx
movsd  %xmm0,(%rbx,%rax,8)
add    $0x1,%r15
cmp    %r15,%r12
jne    12e0 <lagrange_interp_2Dfo+0x1a0>
jmp    12a0 <lagrange_interp_2Dfo+0x160>
movsd  -0x58(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm1
mov    -0x68(%rbp),%rdi
mov    -0x40(%rbp),%rdx
mov    %r12d,%ecx
mov    -0x78(%rbp),%r8
call   1030 <lagrange_interp_2D@plt>
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define double @lagrange_interp_2Dfo(double %0, double %1, i32 %2, i32 %3, double* %4, i32 %5, double* %6, i32 %7, double* nocapture %8) {
%10 = and i32 %2, 1
%11 = xor i32 %10, 1
%12 = tail call i32 @assert(i32 %11)
%13 = and i32 %3, 1
%14 = xor i32 %13, 1
%15 = tail call i32 @assert(i32 %14)
%16 = icmp sgt i32 %5, %2
%17 = zext i1 %16 to i32
%18 = tail call i32 @assert(i32 %17)
%19 = icmp sgt i32 %7, %3
%20 = zext i1 %19 to i32
%21 = tail call i32 @assert(i32 %20)
%22 = add i32 %2, 1
%23 = zext i32 %22 to i64
%24 = alloca double, i64 %23
%25 = add i32 %3, 1
%26 = zext i32 %25 to i64
%27 = alloca double, i64 %26
%28 = mul nsw i32 %25, %22
%29 = zext i32 %28 to i64
%30 = alloca double, i64 %29
%31 = sdiv i32 %2, 2
%32 = tail call i32 @find_xb(double %0, double* %4, i32 %5, i32 %31)
%33 = sdiv i32 %3, 2
%34 = tail call i32 @find_xb(double %1, double* %6, i32 %7, i32 %33)
%35 = icmp slt i32 %2, 0
br i1 %35, label %40, label %36
36: ; preds = %9
%37 = sub i32 %32, %31
%38 = icmp slt i32 %3, 0
%39 = sub i32 %34, %33
br label %42
40: ; preds = %50, %9
%41 = call double @lagrange_interp_2D(double %0, double %1, double* %24, i32 %22, double* %27, i32 %25, double* %30)
ret double %41
42: ; preds = %36, %50
%43 = phi i64 [ 0, %36 ], [ %51, %50 ]
%44 = trunc i64 %43 to i32
%45 = add i32 %37, %44
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds double, double* %4, i64 %46
%48 = load double, double* %47
%49 = getelementptr inbounds double, double* %24, i64 %43
store double %48, double* %49
br i1 %38, label %50, label %53
50: ; preds = %53, %42
%51 = add nuw nsw i64 %43, 1
%52 = icmp eq i64 %51, %23
br i1 %52, label %40, label %42
53: ; preds = %42, %53
%54 = phi i64 [ %68, %53 ], [ 0, %42 ]
%55 = trunc i64 %54 to i32
%56 = add i32 %39, %55
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds double, double* %6, i64 %57
%59 = load double, double* %58
%60 = getelementptr inbounds double, double* %27, i64 %54
store double %59, double* %60
%61 = tail call i32 @index_2D_to_1D(i32 %44, i32 %22, i32 %55, i32 %25)
%62 = tail call i32 @index_2D_to_1D(i32 %45, i32 %5, i32 %56, i32 %7)
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds double, double* %8, i64 %63
%65 = load double, double* %64
%66 = sext i32 %61 to i64
%67 = getelementptr inbounds double, double* %30, i64 %66
store double %65, double* %67
%68 = add nuw nsw i64 %54, 1
%69 = icmp eq i64 %68, %26
br i1 %69, label %50, label %53
}