void polint(xa, ya, n, x, y, dy) float xa[], ya[], x, *y, *dy;
int n;
{
  int i, m, ns = 1;
  float den, dif, dift, ho, hp, w;
  float *c, *d, *vector();
  void nrerror(), free_vector();
  dif = fabs(x - xa[1]);
  c = vector(1, n);
  d = vector(1, n);
  for (i = 1; 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 = 1; m < n; m++) {
    for (i = 1; i <= n - m; i++) {
      ho = xa[i] - x;
      hp = xa[i + m] - x;
      w = c[i + 1] - d[i];
      if ((den = ho - hp) == 0.0)
        nrerror("Error in routine POLINT");
      den = w / den;
      d[i] = hp * den;
      c[i] = ho * den;
    }
    *y += (*dy = (2 * ns < (n - m) ? c[ns + 1] : d[ns--]));
  }
  free_vector(d, 1, n);
  free_vector(c, 1, n);
}


x86 assembly (clang, -O0):
<polint>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
cvtsd2ss %xmm0,%xmm0
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movss  %xmm0,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movl   $0x1,-0x34(%rbp)
movss  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
subss  0x4(%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movaps 0xe8c(%rip),%xmm1
pand   %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x14(%rbp),%esi
mov    $0x1,%edi
mov    $0x0,%al
call   1050 <vector@plt>
mov    %rax,-0x58(%rbp)
mov    -0x14(%rbp),%esi
mov    $0x1,%edi
mov    $0x0,%al
call   1050 <vector@plt>
mov    %rax,-0x60(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1244 <polint+0x114>
movss  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movaps 0xe29(%rip),%xmm1
pand   %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm1
movss  %xmm1,-0x40(%rbp)
movss  -0x3c(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    1202 <polint+0xd2>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
movss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x58(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x60(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11ae <polint+0x7e>
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    %ecx,%edx
add    $0xffffffff,%edx
mov    %edx,-0x34(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    13d7 <polint+0x2a7>
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x30(%rbp),%ecx
cmp    %ecx,%eax
jg     1364 <polint+0x234>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x58(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x60(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x44(%rbp),%xmm0
subss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1319 <polint+0x1e9>
jp     1319 <polint+0x1e9>
lea    0xcfe(%rip),%rdi
mov    $0x0,%al
call   1030 <nrerror@plt>
movss  -0x4c(%rbp),%xmm0
divss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x48(%rbp),%xmm0
mulss  -0x38(%rbp),%xmm0
mov    -0x60(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x44(%rbp),%xmm0
mulss  -0x38(%rbp),%xmm0
mov    -0x58(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    127d <polint+0x14d>
mov    -0x34(%rbp),%eax
shl    $0x1,%eax
mov    -0x14(%rbp),%ecx
sub    -0x30(%rbp),%ecx
cmp    %ecx,%eax
jge    1394 <polint+0x264>
mov    -0x58(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x64(%rbp)
jmp    13b0 <polint+0x280>
mov    -0x60(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    %ecx,%edx
add    $0xffffffff,%edx
mov    %edx,-0x34(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x64(%rbp)
movss  -0x64(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
addss  (%rax),%xmm0
movss  %xmm0,(%rax)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    126a <polint+0x13a>
mov    -0x60(%rbp),%rdi
mov    -0x14(%rbp),%edx
mov    $0x1,%esi
mov    $0x0,%al
call   1040 <free_vector@plt>
mov    -0x58(%rbp),%rdi
mov    -0x14(%rbp),%edx
mov    $0x1,%esi
mov    $0x0,%al
call   1040 <free_vector@plt>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<polint>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %r8,0x48(%rsp)
mov    %rcx,%r12
movsd  %xmm0,0x10(%rsp)
mov    %edx,%r13d
mov    %rsi,%r15
mov    %rdi,0x28(%rsp)
movss  0x4(%rdi),%xmm0
movss  %xmm0,0x8(%rsp)
mov    $0x1,%r14d
mov    $0x1,%edi
mov    %edx,%esi
xor    %eax,%eax
call   1050 <vector@plt>
mov    %rax,%rbx
mov    $0x1,%edi
mov    %r13d,%esi
xor    %eax,%eax
call   1050 <vector@plt>
mov    %rax,%rbp
test   %r13d,%r13d
jle    138a <polint+0x25a>
movsd  0x10(%rsp),%xmm0
cvtsd2ss %xmm0,%xmm3
movaps %xmm3,%xmm0
subss  0x8(%rsp),%xmm0
andps  0xe52(%rip),%xmm0
lea    0x1(%r13),%eax
mov    $0x1,%ecx
movaps 0xe42(%rip),%xmm1
mov    0x28(%rsp),%rsi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movaps %xmm3,%xmm2
subss  (%rsi,%rcx,4),%xmm2
andps  %xmm1,%xmm2
ucomiss %xmm2,%xmm0
cmova  %ecx,%r14d
minss  %xmm0,%xmm2
movss  (%r15,%rcx,4),%xmm0
movss  %xmm0,(%rbx,%rcx,4)
movss  (%r15,%rcx,4),%xmm0
movss  %xmm0,0x0(%rbp,%rcx,4)
add    $0x1,%rcx
movaps %xmm2,%xmm0
cmp    %rcx,%rax
jne    11d0 <polint+0xa0>
movslq %r14d,%rax
movss  (%r15,%rax,4),%xmm0
mov    %r12,0x40(%rsp)
movss  %xmm0,(%r12)
cmp    $0x2,%r13d
jl     1396 <polint+0x266>
add    $0xffffffff,%r14d
mov    %r13d,%eax
mov    %rax,0x30(%rsp)
lea    0x8(%rsi),%r12
mov    $0x1,%eax
mov    %r13d,%ecx
mov    %r13d,0xc(%rsp)
mov    %r13,0x38(%rsp)
movss  %xmm3,0x1c(%rsp)
jmp    12bd <polint+0x18d>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    (%r14,%r14,1),%eax
xor    %ecx,%ecx
cmp    0x18(%rsp),%eax
setge  %cl
movslq %r14d,%r14
lea    0x4(%rbx,%r14,4),%rax
lea    0x0(%rbp,%r14,4),%rdx
cmovl  %rax,%rdx
sub    %ecx,%r14d
movss  (%rdx),%xmm0
mov    0x48(%rsp),%rax
movss  %xmm0,(%rax)
mov    0x40(%rsp),%rax
addss  (%rax),%xmm0
movss  %xmm0,(%rax)
mov    0x50(%rsp),%rax
add    $0x1,%rax
addl   $0xffffffff,0xc(%rsp)
add    $0x4,%r12
cmp    0x30(%rsp),%rax
mov    0x38(%rsp),%r13
je     1396 <polint+0x266>
mov    %r13d,%ecx
mov    %rax,0x50(%rsp)
sub    %eax,%ecx
mov    %ecx,0x18(%rsp)
jle    1260 <polint+0x130>
mov    0xc(%rsp),%r15d
add    $0xffffffffffffffff,%r15
xor    %r13d,%r13d
jmp    1307 <polint+0x1d7>
nopl   0x0(%rax,%rax,1)
divss  %xmm4,%xmm2
mulss  %xmm2,%xmm1
movss  %xmm1,0x4(%rbp,%r13,4)
mulss  %xmm2,%xmm0
movss  %xmm0,0x4(%rbx,%r13,4)
add    $0x1,%r13
cmp    %r13,%r15
je     1260 <polint+0x130>
movss  0x4(%rsi,%r13,4),%xmm0
subss  %xmm3,%xmm0
movss  (%r12,%r13,4),%xmm1
movss  0x8(%rbx,%r13,4),%xmm2
subss  0x4(%rbp,%r13,4),%xmm2
subss  %xmm3,%xmm1
movaps %xmm0,%xmm4
subss  %xmm1,%xmm4
ucomiss 0xcd8(%rip),%xmm4
jne    12e0 <polint+0x1b0>
jp     12e0 <polint+0x1b0>
lea    0xcd1(%rip),%rdi
xor    %eax,%eax
movss  %xmm0,0x10(%rsp)
movss  %xmm1,0x8(%rsp)
movss  %xmm2,0x24(%rsp)
movss  %xmm4,0x20(%rsp)
call   1030 <nrerror@plt>
movss  0x20(%rsp),%xmm4
movss  0x24(%rsp),%xmm2
movss  0x8(%rsp),%xmm1
movss  0x10(%rsp),%xmm0
movss  0x1c(%rsp),%xmm3
mov    0x28(%rsp),%rsi
jmp    12e0 <polint+0x1b0>
movss  0x4(%r15),%xmm0
movss  %xmm0,(%r12)
mov    %rbp,%rdi
mov    $0x1,%esi
mov    %r13d,%edx
xor    %eax,%eax
call   1040 <free_vector@plt>
mov    %rbx,%rdi
mov    $0x1,%esi
mov    %r13d,%edx
xor    %eax,%eax
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <free_vector@plt>

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

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