#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 fabs(double);
double stub1(double);
double stub2(double);
double stub3(double);
double romberg(double (*g)(double), double a, double b, double eps) {
  double t[12 + 1][12 + 1], h, fourj, sum, integ;
  int m, k, i, j;
  h = b - a;
  m = 1;
  t[1][1] = h * ((*g)(a) + (*g)(b)) / 2.;
  for (k = 2; k <= 12; k++) {
    h /= 2.;
    m *= 2;
    for (i = 1, sum = 0.; i <= m; i += 2)
      sum += (*g)(a + i * h);
    t[k][1] = t[k - 1][1] * .5 + sum * h;
    for (j = 2, fourj = 1.; j <= k; j++) {
      fourj *= 4.;
      t[k][j] = t[k][j - 1] + (t[k][j - 1] - t[k - 1][j - 1]) / (fourj - 1.);
    }
    if (fabs((t[k][k] - t[k - 1][k - 1]) / t[k][k]) <= eps) {
      integ = t[k][k];
      return (integ);
    }
  }
  integ = t[12][12];
  return (integ);
}


x86 assembly (clang, -O0):
<romberg>:
push   %rbp
mov    %rsp,%rbp
sub    $0x5b0,%rsp
mov    %rdi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
movsd  %xmm2,-0x28(%rbp)
movsd  -0x20(%rbp),%xmm0
subsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x578(%rbp)
ff
movl   $0x1,-0x594(%rbp)
movsd  -0x578(%rbp),%xmm0
ff
movsd  %xmm0,-0x5a8(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  -0x18(%rbp),%xmm0
call   *%rax
movsd  %xmm0,-0x5b0(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  -0x20(%rbp),%xmm0
call   *%rax
movsd  -0x5b0(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x5a8(%rbp),%xmm0
ff
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xe75(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x500(%rbp)
ff
movl   $0x2,-0x598(%rbp)
cmpl   $0xc,-0x598(%rbp)
jg     1487 <romberg+0x387>
movsd  -0x578(%rbp),%xmm0
ff
movsd  0xe42(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x578(%rbp)
ff
mov    -0x594(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x594(%rbp)
movl   $0x1,-0x59c(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x588(%rbp)
ff
mov    -0x59c(%rbp),%eax
cmp    -0x594(%rbp),%eax
jg     1247 <romberg+0x147>
mov    -0x10(%rbp),%rax
movsd  -0x18(%rbp),%xmm1
cvtsi2sdl -0x59c(%rbp),%xmm0
ff
movsd  -0x578(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   *%rax
addsd  -0x588(%rbp),%xmm0
ff
movsd  %xmm0,-0x588(%rbp)
ff
mov    -0x59c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x59c(%rbp)
jmp    11ee <romberg+0xee>
mov    -0x598(%rbp),%eax
sub    $0x1,%eax
movslq %eax,%rcx
lea    -0x570(%rbp),%rax
imul   $0x68,%rcx,%rcx
add    %rcx,%rax
movsd  0x8(%rax),%xmm0
movsd  -0x588(%rbp),%xmm1
ff
mulsd  -0x578(%rbp),%xmm1
ff
movsd  0xd92(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movslq -0x598(%rbp),%rcx
lea    -0x570(%rbp),%rax
imul   $0x68,%rcx,%rcx
add    %rcx,%rax
movsd  %xmm0,0x8(%rax)
movl   $0x2,-0x5a0(%rbp)
movsd  0xd56(%rip),%xmm0
movsd  %xmm0,-0x580(%rbp)
ff
mov    -0x5a0(%rbp),%eax
cmp    -0x598(%rbp),%eax
jg     13ac <romberg+0x2ac>
movsd  0xd44(%rip),%xmm0
mulsd  -0x580(%rbp),%xmm0
ff
movsd  %xmm0,-0x580(%rbp)
ff
movslq -0x598(%rbp),%rcx
lea    -0x570(%rbp),%rax
imul   $0x68,%rcx,%rcx
add    %rcx,%rax
mov    -0x5a0(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movslq -0x598(%rbp),%rcx
lea    -0x570(%rbp),%rax
imul   $0x68,%rcx,%rcx
add    %rcx,%rax
mov    -0x5a0(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x598(%rbp),%eax
sub    $0x1,%eax
movslq %eax,%rcx
lea    -0x570(%rbp),%rax
imul   $0x68,%rcx,%rcx
add    %rcx,%rax
mov    -0x5a0(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm1
movsd  -0x580(%rbp),%xmm2
ff
movsd  0xc9d(%rip),%xmm3
subsd  %xmm3,%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movslq -0x598(%rbp),%rcx
lea    -0x570(%rbp),%rax
imul   $0x68,%rcx,%rcx
add    %rcx,%rax
movslq -0x5a0(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x5a0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5a0(%rbp)
jmp    12ba <romberg+0x1ba>
movslq -0x598(%rbp),%rcx
lea    -0x570(%rbp),%rax
imul   $0x68,%rcx,%rcx
add    %rcx,%rax
movslq -0x598(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x598(%rbp),%eax
sub    $0x1,%eax
movslq %eax,%rcx
lea    -0x570(%rbp),%rax
imul   $0x68,%rcx,%rcx
add    %rcx,%rax
mov    -0x598(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm1
movslq -0x598(%rbp),%rcx
lea    -0x570(%rbp),%rax
imul   $0x68,%rcx,%rcx
add    %rcx,%rax
movslq -0x598(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm1
movaps 0xc00(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  -0x28(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     146e <romberg+0x36e>
movslq -0x598(%rbp),%rcx
lea    -0x570(%rbp),%rax
imul   $0x68,%rcx,%rcx
add    %rcx,%rax
movslq -0x598(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x590(%rbp)
ff
movsd  -0x590(%rbp),%xmm0
ff
movsd  %xmm0,-0x8(%rbp)
jmp    14a1 <romberg+0x3a1>
jmp    1473 <romberg+0x373>
mov    -0x598(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x598(%rbp)
jmp    11a1 <romberg+0xa1>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x590(%rbp)
ff
movsd  -0x590(%rbp),%xmm0
ff
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x5b0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<romberg>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x578,%rsp
movsd  %xmm2,0x18(%rsp)
movsd  %xmm1,(%rsp)
mov    %rdi,%r14
movapd %xmm1,%xmm2
subsd  %xmm0,%xmm2
movsd  %xmm2,0x8(%rsp)
movsd  %xmm0,0x28(%rsp)
call   *%rdi
movsd  %xmm0,0x10(%rsp)
movsd  (%rsp),%xmm0
call   *%r14
movsd  0x8(%rsp),%xmm4
addsd  0x10(%rsp),%xmm0
mulsd  %xmm4,%xmm0
mulsd  0xea5(%rip),%xmm0
movsd  %xmm0,0xa0(%rsp)
lea    0x118(%rsp),%rbx
mov    $0x1,%ebp
mov    $0x3,%esi
mov    $0x2,%r12d
xor    %r13d,%r13d
movsd  0xe79(%rip),%xmm5
nopw   0x0(%rax,%rax,1)
mov    %rsi,0x20(%rsp)
mulsd  %xmm5,%xmm4
add    %ebp,%ebp
xorps  %xmm3,%xmm3
movsd  0xe62(%rip),%xmm0
mov    $0x1,%r15d
movsd  %xmm4,0x8(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  %xmm0,(%rsp)
movsd  %xmm3,0x10(%rsp)
movsd  (%rsp),%xmm0
mulsd  %xmm4,%xmm0
addsd  0x28(%rsp),%xmm0
call   *%r14
movsd  0x10(%rsp),%xmm3
movsd  0xe25(%rip),%xmm1
movsd  0x8(%rsp),%xmm4
addsd  %xmm0,%xmm3
movsd  (%rsp),%xmm0
addsd  %xmm1,%xmm0
add    $0x2,%r15d
cmp    %ebp,%r15d
jbe    11c0 <romberg+0xc0>
lea    -0x1(%r12),%rax
imul   $0x68,%r12,%rcx
movsd  -0x30(%rsp,%rcx,1),%xmm1
mulsd  %xmm4,%xmm3
movapd %xmm1,%xmm2
movsd  0xdda(%rip),%xmm5
mulsd  %xmm5,%xmm2
addsd  %xmm3,%xmm2
movsd  %xmm2,0x38(%rsp,%rcx,1)
movapd %xmm2,%xmm0
subsd  %xmm1,%xmm0
divsd  0xdd4(%rip),%xmm0
addsd  %xmm2,%xmm0
movsd  %xmm0,0x40(%rsp,%rcx,1)
mov    0x20(%rsp),%rsi
cmp    $0x3,%rsi
movsd  0xdc1(%rip),%xmm6
movsd  0xdc1(%rip),%xmm7
je     129c <romberg+0x19c>
xor    %edx,%edx
movapd %xmm6,%xmm1
nop
mulsd  %xmm6,%xmm1
movapd %xmm0,%xmm2
subsd  -0x70(%rbx,%rdx,1),%xmm2
movapd %xmm1,%xmm3
addsd  %xmm7,%xmm3
divsd  %xmm3,%xmm2
addsd  %xmm2,%xmm0
movsd  %xmm0,(%rbx,%rdx,1)
add    $0x8,%rdx
cmp    %rdx,%r13
jne    1270 <romberg+0x170>
add    %rsp,%rcx
add    $0x30,%rcx
movsd  (%rcx,%r12,8),%xmm0
imul   $0x68,%rax,%rcx
add    %rsp,%rcx
add    $0x30,%rcx
movapd %xmm0,%xmm1
subsd  (%rcx,%rax,8),%xmm1
divsd  %xmm0,%xmm1
andpd  0xd67(%rip),%xmm1
movsd  0x18(%rsp),%xmm2
ucomisd %xmm1,%xmm2
jae    12f8 <romberg+0x1f8>
add    $0x1,%r12
add    $0x1,%rsi
add    $0x68,%rbx
add    $0x8,%r13
cmp    $0xd,%r12
jne    1190 <romberg+0x90>
movsd  0x570(%rsp),%xmm0
add    $0x578,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define double @romberg(double (double)* nocapture %0, double %1, double %2, double %3) {
%5 = alloca [13 x [13 x double]]
%6 = bitcast [13 x [13 x double]]* %5 to i8*
%7 = fsub double %2, %1
%8 = tail call double %0(double %1)
%9 = tail call double %0(double %2)
%10 = fadd double %8, %9
%11 = fmul double %7, %10
%12 = fmul double %11, 5.000000e-01
%13 = getelementptr inbounds [13 x [13 x double]], [13 x [13 x double]]* %5, i64 0, i64 1, i64 1
store double %12, double* %13
br label %18
14: ; preds = %61
%15 = add nuw nsw i64 %19, 1
%16 = add nuw nsw i64 %20, 1
%17 = icmp eq i64 %15, 13
br i1 %17, label %70, label %18
18: ; preds = %14, %4
%19 = phi i64 [ 2, %4 ], [ %15, %14 ]
%20 = phi i64 [ 3, %4 ], [ %16, %14 ]
%21 = phi i32 [ 1, %4 ], [ %24, %14 ]
%22 = phi double [ %7, %4 ], [ %23, %14 ]
%23 = fmul double %22, 5.000000e-01
%24 = shl nsw i32 %21, 1
br label %25
25: ; preds = %18, %25
%26 = phi i32 [ %32, %25 ], [ 1, %18 ]
%27 = phi double [ %31, %25 ], [ 0.000000e+00, %18 ]
%28 = sitofp i32 %26 to double
%29 = tail call double @llvm.fmuladd.f64(double %28, double %23, double %1)
%30 = tail call double %0(double %29)
%31 = fadd double %27, %30
%32 = add nuw nsw i32 %26, 2
%33 = icmp ugt i32 %32, %24
br i1 %33, label %34, label %25
34: ; preds = %25
%35 = add nsw i64 %19, -1
%36 = getelementptr inbounds [13 x [13 x double]], [13 x [13 x double]]* %5, i64 0, i64 %35, i64 1
%37 = load double, double* %36
%38 = fmul double %23, %31
%39 = tail call double @llvm.fmuladd.f64(double %37, double 5.000000e-01, double %38)
%40 = getelementptr inbounds [13 x [13 x double]], [13 x [13 x double]]* %5, i64 0, i64 %19, i64 1
store double %39, double* %40
%41 = fsub double %39, %37
%42 = fdiv double %41, 3.000000e+00
%43 = fadd double %39, %42
%44 = getelementptr inbounds [13 x [13 x double]], [13 x [13 x double]]* %5, i64 0, i64 %19, i64 2
store double %43, double* %44
%45 = icmp eq i64 %20, 3
br i1 %45, label %61, label %46
46: ; preds = %34, %46
%47 = phi i64 [ %59, %46 ], [ 3, %34 ]
%48 = phi double [ %57, %46 ], [ %43, %34 ]
%49 = phi double [ %53, %46 ], [ 4.000000e+00, %34 ]
%50 = phi i64 [ %47, %46 ], [ 2, %34 ]
%51 = getelementptr inbounds [13 x [13 x double]], [13 x [13 x double]]* %5, i64 0, i64 %35, i64 %50
%52 = load double, double* %51
%53 = fmul double %49, 4.000000e+00
%54 = fsub double %48, %52
%55 = fadd double %53, -1.000000e+00
%56 = fdiv double %54, %55
%57 = fadd double %48, %56
%58 = getelementptr inbounds [13 x [13 x double]], [13 x [13 x double]]* %5, i64 0, i64 %19, i64 %47
store double %57, double* %58
%59 = add nuw nsw i64 %47, 1
%60 = icmp eq i64 %59, %20
br i1 %60, label %61, label %46
61: ; preds = %46, %34
%62 = getelementptr inbounds [13 x [13 x double]], [13 x [13 x double]]* %5, i64 0, i64 %19, i64 %19
%63 = load double, double* %62
%64 = getelementptr inbounds [13 x [13 x double]], [13 x [13 x double]]* %5, i64 0, i64 %35, i64 %35
%65 = load double, double* %64
%66 = fsub double %63, %65
%67 = fdiv double %66, %63
%68 = tail call double @llvm.fabs.f64(double %67)
%69 = fcmp ugt double %68, %3
br i1 %69, label %14, label %73
70: ; preds = %14
%71 = getelementptr inbounds [13 x [13 x double]], [13 x [13 x double]]* %5, i64 0, i64 12, i64 12
%72 = load double, double* %71
br label %73
73: ; preds = %61, %70
%74 = phi double [ %72, %70 ], [ %63, %61 ]
ret double %74
}