void polint(double xa[], double ya[], int n, double x, double *y, double *dy) {
  char *func = {"polint(): "};
  int i, m, ns = 1;
  double den, dif, dift, ho, hp, w;
  double *c, *d;
  dif = fabs(x - xa[0]);
  c = (double *)mem_alloc((n) * sizeof(double), func, "c[]");
  d = (double *)mem_alloc((n) * sizeof(double), func, "d[]");
  for (i = 0; i < n; i++) {
    if ((dift = fabs(x - xa[i])) < dif) {
      ns = i;
      dif = dift;
    }
    c[i] = ya[i];
    d[i] = ya[i];
  }
  *y = ya[ns--];
  for (m = 0; m < (n - 1); m++) {
    for (i = 0; i < n - m; i++) {
      ho = xa[i] - x;
      hp = xa[i + m] - x;
      w = c[i + 1] - d[i];
      if ((den = ho - hp) < 1.0E-10) {
        printf("\n\n Error in function polint(): ");
        printf("\nden = ho - hp = %4.1E -- too small\n", den);
        exit(1);
      }
      den = w / den;
      d[i] = hp * den;
      c[i] = ho * den;
    }
    *y += (*dy = (2 * ns < (n - m) ? c[ns + 1] : d[ns--]));
  }
  free(d);
  free(c);
}


x86 assembly (clang, -O0):
<polint>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
lea    0xecb(%rip),%rax
mov    %rax,-0x38(%rbp)
movl   $0x1,-0x44(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
subsd  (%rax),%xmm0
movaps 0xe77(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movslq -0x14(%rbp),%rdi
shl    $0x3,%rdi
mov    -0x38(%rbp),%rsi
lea    0xe73(%rip),%rdx
mov    $0x0,%al
call   1050 <mem_alloc@plt>
cltq
mov    %rax,-0x80(%rbp)
movslq -0x14(%rbp),%rdi
shl    $0x3,%rdi
mov    -0x38(%rbp),%rsi
lea    0xe57(%rip),%rdx
mov    $0x0,%al
call   1050 <mem_alloc@plt>
cltq
mov    %rax,-0x88(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    126e <polint+0x12e>
movsd  -0x20(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
movaps 0xdff(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  %xmm1,-0x60(%rbp)
movsd  -0x58(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1229 <polint+0xe9>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x44(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x80(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x88(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11dc <polint+0x9c>
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    %ecx,%edx
add    $0xffffffff,%edx
mov    %edx,-0x44(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1430 <polint+0x2f0>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x40(%rbp),%ecx
cmp    %ecx,%eax
jge    13b1 <polint+0x271>
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x3c(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x80(%rbp),%rax
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x88(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x68(%rbp),%xmm1
subsd  -0x70(%rbp),%xmm1
movsd  %xmm1,-0x50(%rbp)
movsd  0xce2(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1363 <polint+0x223>
lea    0xce1(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
movsd  -0x50(%rbp),%xmm0
lea    0xcee(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
mov    $0x1,%edi
call   1060 <exit@plt>
movsd  -0x78(%rbp),%xmm0
divsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x70(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
mov    -0x88(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x68(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
mov    -0x80(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    12ac <polint+0x16c>
mov    -0x44(%rbp),%eax
shl    $0x1,%eax
mov    -0x14(%rbp),%ecx
sub    -0x40(%rbp),%ecx
cmp    %ecx,%eax
jge    13e4 <polint+0x2a4>
mov    -0x80(%rbp),%rax
mov    -0x44(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
jmp    1406 <polint+0x2c6>
mov    -0x88(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    %ecx,%edx
add    $0xffffffff,%edx
mov    %edx,-0x44(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1294 <polint+0x154>
mov    -0x88(%rbp),%rdi
call   1030 <free@plt>
mov    -0x80(%rbp),%rdi
call   1030 <free@plt>
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<polint>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r8,0x20(%rsp)
mov    %rcx,%r12
movsd  %xmm0,0x18(%rsp)
mov    %edx,%r15d
mov    %rsi,%rbp
mov    %rdi,%r13
movsd  (%rdi),%xmm0
movsd  %xmm0,0x10(%rsp)
movslq %edx,%rax
mov    %rax,0x8(%rsp)
lea    0x0(,%rax,8),%r14
lea    0xeaf(%rip),%rsi
lea    0xe8b(%rip),%rdx
mov    %r14,%rdi
xor    %eax,%eax
call   1050 <mem_alloc@plt>
movslq %eax,%rbx
lea    0xe7b(%rip),%rdx
mov    %r14,%rdi
lea    0xe8a(%rip),%rsi
xor    %eax,%eax
call   1050 <mem_alloc@plt>
movslq %eax,%rdi
cmpl   $0x0,0x8(%rsp)
jle    1314 <polint+0x1d4>
movsd  0x18(%rsp),%xmm4
movapd %xmm4,%xmm0
subsd  0x10(%rsp),%xmm0
andpd  0xe28(%rip),%xmm0
mov    %r15d,%ecx
mov    $0x1,%eax
xor    %edx,%edx
movapd 0xe16(%rip),%xmm1
nopw   0x0(%rax,%rax,1)
movapd %xmm4,%xmm2
subsd  0x0(%r13,%rdx,8),%xmm2
andpd  %xmm1,%xmm2
ucomisd %xmm2,%xmm0
cmova  %edx,%eax
minsd  %xmm0,%xmm2
movsd  0x0(%rbp,%rdx,8),%xmm0
movsd  %xmm0,(%rbx,%rdx,8)
movsd  0x0(%rbp,%rdx,8),%xmm0
movsd  %xmm0,(%rdi,%rdx,8)
add    $0x1,%rdx
movapd %xmm2,%xmm0
cmp    %rdx,%rcx
jne    11f0 <polint+0xb0>
movslq %eax,%rdx
movsd  0x0(%rbp,%rdx,8),%xmm0
movsd  %xmm0,(%r12)
cmp    $0x2,%ecx
mov    0x20(%rsp),%r11
jl     131f <polint+0x1df>
lea    -0x1(%r15),%r8d
add    $0xffffffff,%eax
xor    %r10d,%r10d
movsd  0xdb4(%rip),%xmm0
mov    %r13,%rsi
mov    %r15d,%ebp
jmp    12b1 <polint+0x171>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
lea    (%rax,%rax,1),%ecx
xor    %edx,%edx
cmp    %r9d,%ecx
setge  %dl
cltq
lea    0x8(%rbx,%rax,8),%r9
lea    (%rdi,%rax,8),%rcx
cmovl  %r9,%rcx
sub    %edx,%eax
movsd  (%rcx),%xmm1
movsd  %xmm1,(%r11)
addsd  (%r12),%xmm1
movsd  %xmm1,(%r12)
add    $0x1,%r10
add    $0xffffffff,%ebp
add    $0x8,%rsi
cmp    %r8,%r10
je     131f <polint+0x1df>
mov    %ebp,%ebp
mov    %r15d,%r9d
sub    %r10d,%r9d
jle    1270 <polint+0x130>
xor    %ecx,%ecx
nopl   (%rax)
movsd  0x0(%r13,%rcx,8),%xmm1
subsd  %xmm4,%xmm1
movsd  (%rsi,%rcx,8),%xmm2
subsd  %xmm4,%xmm2
movapd %xmm1,%xmm5
subsd  %xmm2,%xmm5
ucomisd %xmm5,%xmm0
ja     133a <polint+0x1fa>
movsd  0x8(%rbx,%rcx,8),%xmm3
subsd  (%rdi,%rcx,8),%xmm3
divsd  %xmm5,%xmm3
mulsd  %xmm3,%xmm2
movsd  %xmm2,(%rdi,%rcx,8)
mulsd  %xmm3,%xmm1
movsd  %xmm1,(%rbx,%rcx,8)
lea    0x1(%rcx),%rdx
mov    %rdx,%rcx
cmp    %rdx,%rbp
jne    12c0 <polint+0x180>
jmp    1270 <polint+0x130>
movsd  0x8(%rbp),%xmm0
movsd  %xmm0,(%r12)
call   1030 <free@plt>
mov    %rbx,%rdi
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>
lea    0xcdf(%rip),%rdi
xor    %eax,%eax
movsd  %xmm5,0x8(%rsp)
call   1040 <printf@plt>
lea    0xceb(%rip),%rdi
movsd  0x8(%rsp),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
mov    $0x1,%edi
call   1060 <exit@plt>

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

LLVM IR (-O2):
@.str = [11 x i8] c"polint(): \00"
@.str.1 = [4 x i8] c"c[]\00"
@.str.2 = [4 x i8] c"d[]\00"
@.str.3 = [32 x i8] c"\0A\0A Error in function polint(): \00"
@.str.4 = [37 x i8] c"\0Aden = ho - hp = %4.1E -- too small\0A\00"
define void @polint(double* nocapture %0, double* nocapture %1, i32 %2, double %3, double* nocapture %4, double* nocapture writeonly %5) {
%7 = load double, double* %0
%8 = sext i32 %2 to i64
%9 = shl nsw i64 %8, 3
%10 = tail call i32 (i64, i8*, i8*, ...) bitcast (i32 (...)* @mem_alloc to i32 (i64, i8*, i8*, ...)*)(i64 %9, i8* getelementptr inbounds ([11 x i8], [11 x i8]* @.str, i64 0, i64 0), i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.1, i64 0, i64 0))
%11 = sext i32 %10 to i64
%12 = inttoptr i64 %11 to double*
%13 = tail call i32 (i64, i8*, i8*, ...) bitcast (i32 (...)* @mem_alloc to i32 (i64, i8*, i8*, ...)*)(i64 %9, i8* getelementptr inbounds ([11 x i8], [11 x i8]* @.str, i64 0, i64 0), i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0))
%14 = sext i32 %13 to i64
%15 = inttoptr i64 %14 to double*
%16 = icmp sgt i32 %2, 0
br i1 %16, label %20, label %17
17: ; preds = %6
%18 = getelementptr inbounds double, double* %1, i64 1
%19 = load double, double* %18
store double %19, double* %4
br label %105
20: ; preds = %6
%21 = fsub double %3, %7
%22 = tail call double @llvm.fabs.f64(double %21)
%23 = zext i32 %2 to i64
br label %24
24: ; preds = %20, %24
%25 = phi i64 [ 0, %20 ], [ %41, %24 ]
%26 = phi i32 [ 1, %20 ], [ %35, %24 ]
%27 = phi double [ %22, %20 ], [ %33, %24 ]
%28 = getelementptr inbounds double, double* %0, i64 %25
%29 = load double, double* %28
%30 = fsub double %3, %29
%31 = tail call double @llvm.fabs.f64(double %30)
%32 = fcmp olt double %31, %27
%33 = select i1 %32, double %31, double %27
%34 = trunc i64 %25 to i32
%35 = select i1 %32, i32 %34, i32 %26
%36 = getelementptr inbounds double, double* %1, i64 %25
%37 = load double, double* %36
%38 = getelementptr inbounds double, double* %12, i64 %25
store double %37, double* %38
%39 = load double, double* %36
%40 = getelementptr inbounds double, double* %15, i64 %25
store double %39, double* %40
%41 = add nuw nsw i64 %25, 1
%42 = icmp eq i64 %41, %23
br i1 %42, label %43, label %24
43: ; preds = %24
%44 = sext i32 %35 to i64
%45 = getelementptr inbounds double, double* %1, i64 %44
%46 = load double, double* %45
store double %46, double* %4
%47 = icmp sgt i32 %2, 1
br i1 %47, label %48, label %105
48: ; preds = %43
%49 = add nsw i32 %2, -1
%50 = add nsw i32 %35, -1
%51 = zext i32 %49 to i64
br label %52
52: ; preds = %48, %87
%53 = phi i64 [ 0, %48 ], [ %102, %87 ]
%54 = phi i32 [ %2, %48 ], [ %103, %87 ]
%55 = phi i32 [ %50, %48 ], [ %97, %87 ]
%56 = trunc i64 %53 to i32
%57 = sub nsw i32 %2, %56
%58 = icmp sgt i32 %57, 0
br i1 %58, label %59, label %87
59: ; preds = %52
%60 = zext i32 %54 to i64
br label %61
61: ; preds = %59, %77
%62 = phi i64 [ 0, %59 ], [ %70, %77 ]
%63 = getelementptr inbounds double, double* %0, i64 %62
%64 = load double, double* %63
%65 = fsub double %64, %3
%66 = add nuw nsw i64 %62, %53
%67 = getelementptr inbounds double, double* %0, i64 %66
%68 = load double, double* %67
%69 = fsub double %68, %3
%70 = add nuw nsw i64 %62, 1
%71 = getelementptr inbounds double, double* %15, i64 %62
%72 = fsub double %65, %69
%73 = fcmp olt double %72, 1.000000e-10
br i1 %73, label %74, label %77
74: ; preds = %61
%75 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([32 x i8], [32 x i8]* @.str.3, i64 0, i64 0))
%76 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([37 x i8], [37 x i8]* @.str.4, i64 0, i64 0), double %72)
tail call void @exit(i32 1)
unreachable
77: ; preds = %61
%78 = getelementptr inbounds double, double* %12, i64 %70
%79 = load double, double* %78
%80 = load double, double* %71
%81 = fsub double %79, %80
%82 = fdiv double %81, %72
%83 = fmul double %69, %82
store double %83, double* %71
%84 = fmul double %65, %82
%85 = getelementptr inbounds double, double* %12, i64 %62
store double %84, double* %85
%86 = icmp eq i64 %70, %60
br i1 %86, label %87, label %61
87: ; preds = %77, %52
%88 = shl nsw i32 %55, 1
%89 = icmp slt i32 %88, %57
%90 = add nsw i32 %55, 1
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double, double* %12, i64 %91
%93 = sext i32 %55 to i64
%94 = getelementptr inbounds double, double* %15, i64 %93
%95 = xor i1 %89, true
%96 = sext i1 %95 to i32
%97 = add nsw i32 %55, %96
%98 = select i1 %89, double* %92, double* %94
%99 = load double, double* %98
store double %99, double* %5
%100 = load double, double* %4
%101 = fadd double %99, %100
store double %101, double* %4
%102 = add nuw nsw i64 %53, 1
%103 = add i32 %54, -1
%104 = icmp eq i64 %102, %51
br i1 %104, label %105, label %52
105: ; preds = %87, %17, %43
%106 = inttoptr i64 %14 to i8*
tail call void @free(i8* %106)
%107 = inttoptr i64 %11 to i8*
tail call void @free(i8* %107)
ret void
}