#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 */
float sqrt(float);
void nrfit(float *x, float *y, int ndata, float *sig, int mwt, float *a,
           float *b, float *siga, float *sigb, float *chi2, float *q) {
  int i;
  float wt, t, sxoss, sx = 0.0, sy = 0.0, st2 = 0.0, ss, sigdat;
  *b = 0.0;
  if (mwt) {
    ss = 0.0;
    for (i = 0; i < ndata; i++) {
      wt = 1.0 / (sig[i] * sig[i]);
      ss += wt;
      sx += x[i] * wt;
      sy += y[i] * wt;
    }
  } else {
    for (i = 0; i < ndata; i++) {
      sx += x[i];
      sy += y[i];
    }
    ss = ndata;
  }
  sxoss = sx / ss;
  if (mwt) {
    for (i = 0; i < ndata; i++) {
      t = (x[i] - sxoss) / sig[i];
      st2 += t * t;
      *b += t * y[i] / sig[i];
    }
  } else {
    for (i = 0; i < ndata; i++) {
      t = x[i] - sxoss;
      st2 += t * t;
      *b += t * y[i];
    }
  }
  *b /= st2;
  *a = (sy - sx * (*b)) / ss;
  *siga = sqrt((1.0 + sx * sx / (ss * st2)) / ss);
  *sigb = sqrt(1.0 / st2);
  *chi2 = 0.0;
  if (mwt == 0) {
    for (i = 0; i < ndata; i++)
      *chi2 += (y[i] - (*a) - (*b) * x[i]) * (y[i] - (*a) - (*b) * x[i]);
    *q = 1.0;
    sigdat = sqrt((*chi2) / (ndata - 2));
    *siga *= sigdat;
    *sigb *= sigdat;
  } else {
    for (i = 0; i < ndata; i++)
      *chi2 += ((y[i] - (*a) - (*b) * x[i]) / sig[i]) *
               ((y[i] - (*a) - (*b) * x[i]) / sig[i]);
    *q = 1.0;
  }
}


x86 assembly (clang, -O0):
<nrfit>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
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    %r8d,-0x24(%rbp)
mov    %r9,-0x30(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x48(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4c(%rbp)
mov    0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
cmpl   $0x0,-0x24(%rbp)
je     1228 <nrfit+0x118>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1223 <nrfit+0x113>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe4f(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x38(%rbp),%xmm0
addss  -0x50(%rbp),%xmm0
movss  %xmm0,-0x50(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x38(%rbp),%xmm2
movss  -0x44(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x38(%rbp),%xmm2
movss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    117f <nrfit+0x6f>
jmp    1281 <nrfit+0x171>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1277 <nrfit+0x167>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    122f <nrfit+0x11f>
cvtsi2ssl -0x14(%rbp),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x44(%rbp),%xmm0
divss  -0x50(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
cmpl   $0x0,-0x24(%rbp)
je     132b <nrfit+0x21b>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1326 <nrfit+0x216>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x40(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x3c(%rbp),%xmm0
movss  -0x3c(%rbp),%xmm2
movss  -0x4c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x3c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm0
mov    0x10(%rbp),%rax
addss  (%rax),%xmm0
movss  %xmm0,(%rax)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    12a1 <nrfit+0x191>
jmp    13aa <nrfit+0x29a>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    13a5 <nrfit+0x295>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x3c(%rbp),%xmm0
movss  -0x3c(%rbp),%xmm2
movss  -0x4c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x3c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    0x10(%rbp),%rax
movss  (%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1332 <nrfit+0x222>
jmp    13aa <nrfit+0x29a>
movss  -0x4c(%rbp),%xmm1
mov    0x10(%rbp),%rax
movss  (%rax),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax)
movss  -0x48(%rbp),%xmm1
movss  -0x44(%rbp),%xmm0
mov    0x10(%rbp),%rax
movss  (%rax),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
divss  -0x50(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x44(%rbp),%xmm0
mulss  -0x44(%rbp),%xmm0
movss  -0x50(%rbp),%xmm1
mulss  -0x4c(%rbp),%xmm1
divss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xbe9(%rip),%xmm0
addsd  %xmm1,%xmm0
movss  -0x50(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1030 <sqrt@plt>
mov    0x18(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x4c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xbb6(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1030 <sqrt@plt>
mov    0x20(%rbp),%rax
movss  %xmm0,(%rax)
mov    0x28(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
cmpl   $0x0,-0x24(%rbp)
jne    157e <nrfit+0x46e>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1527 <nrfit+0x417>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x30(%rbp),%rax
subss  (%rax),%xmm1
mov    0x10(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x30(%rbp),%rax
subss  (%rax),%xmm1
mov    0x10(%rbp),%rax
movss  (%rax),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm3
movd   %xmm2,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm2
mulss  %xmm3,%xmm2
addss  %xmm1,%xmm2
mov    0x28(%rbp),%rax
movss  (%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    147b <nrfit+0x36b>
mov    0x30(%rbp),%rax
movss  0xad5(%rip),%xmm0
movss  %xmm0,(%rax)
mov    0x28(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x14(%rbp),%eax
sub    $0x2,%eax
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
call   1030 <sqrt@plt>
movss  %xmm0,-0x54(%rbp)
movss  -0x54(%rbp),%xmm0
mov    0x18(%rbp),%rax
mulss  (%rax),%xmm0
movss  %xmm0,(%rax)
movss  -0x54(%rbp),%xmm0
mov    0x20(%rbp),%rax
mulss  (%rax),%xmm0
movss  %xmm0,(%rax)
jmp    165b <nrfit+0x54b>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    164b <nrfit+0x53b>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x30(%rbp),%rax
subss  (%rax),%xmm1
mov    0x10(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x30(%rbp),%rax
subss  (%rax),%xmm1
mov    0x10(%rbp),%rax
movss  (%rax),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm3
movd   %xmm2,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm2
mulss  %xmm3,%xmm2
addss  %xmm1,%xmm2
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm2
mov    0x28(%rbp),%rax
movss  (%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1585 <nrfit+0x475>
mov    0x30(%rbp),%rax
movss  0x9b1(%rip),%xmm0
movss  %xmm0,(%rax)
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<nrfit>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r9,%r14
mov    %rcx,%rbp
mov    %rsi,%rbx
mov    %rdi,%r15
mov    0x60(%rsp),%r9
movl   $0x0,(%r9)
test   %r8d,%r8d
je     11a4 <nrfit+0x94>
test   %edx,%edx
jle    11c0 <nrfit+0xb0>
mov    %edx,%eax
xorps  %xmm4,%xmm4
xorps  %xmm0,%xmm0
xor    %ecx,%ecx
movss  0xeab(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
nop
movaps %xmm4,%xmm3
movss  0x0(%rbp,%rcx,4),%xmm4
mulss  %xmm4,%xmm4
movaps %xmm2,%xmm1
divss  %xmm4,%xmm1
addss  %xmm1,%xmm0
movss  (%r15,%rcx,4),%xmm4
movss  (%rbx,%rcx,4),%xmm5
unpcklps %xmm4,%xmm5
shufps $0x0,%xmm1,%xmm1
mulps  %xmm5,%xmm1
movaps %xmm1,%xmm4
addps  %xmm3,%xmm4
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1160 <nrfit+0x50>
addps  %xmm3,%xmm1
jmp    1270 <nrfit+0x160>
test   %edx,%edx
jle    11cb <nrfit+0xbb>
mov    %edx,%esi
lea    -0x1(%rsi),%rcx
mov    %esi,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    11d3 <nrfit+0xc3>
xorps  %xmm1,%xmm1
xor    %ecx,%ecx
jmp    1233 <nrfit+0x123>
xorps  %xmm1,%xmm1
xorps  %xmm0,%xmm0
jmp    12e9 <nrfit+0x1d9>
xorps  %xmm1,%xmm1
jmp    1269 <nrfit+0x159>
and    $0xfffffffc,%esi
xorps  %xmm1,%xmm1
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movss  (%r15,%rcx,4),%xmm0
movss  0x4(%r15,%rcx,4),%xmm2
movss  0x8(%r15,%rcx,4),%xmm3
movss  0xc(%r15,%rcx,4),%xmm4
movss  (%rbx,%rcx,4),%xmm5
unpcklps %xmm0,%xmm5
addps  %xmm1,%xmm5
movss  0x4(%rbx,%rcx,4),%xmm0
unpcklps %xmm2,%xmm0
addps  %xmm5,%xmm0
movss  0x8(%rbx,%rcx,4),%xmm2
unpcklps %xmm3,%xmm2
addps  %xmm0,%xmm2
movss  0xc(%rbx,%rcx,4),%xmm1
unpcklps %xmm4,%xmm1
addps  %xmm2,%xmm1
add    $0x4,%rcx
cmp    %rcx,%rsi
jne    11e0 <nrfit+0xd0>
test   %rax,%rax
je     1269 <nrfit+0x159>
lea    (%rbx,%rcx,4),%rdi
lea    (%r15,%rcx,4),%rcx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%rcx,%rsi,4),%xmm0
movss  (%rdi,%rsi,4),%xmm2
unpcklps %xmm0,%xmm2
addps  %xmm2,%xmm1
add    $0x1,%rsi
cmp    %rsi,%rax
jne    1250 <nrfit+0x140>
xorps  %xmm0,%xmm0
cvtsi2ss %edx,%xmm0
movaps %xmm1,%xmm3
shufps $0x55,%xmm1,%xmm3
divss  %xmm0,%xmm3
test   %r8d,%r8d
je     12cc <nrfit+0x1bc>
test   %edx,%edx
jle    12e9 <nrfit+0x1d9>
mov    %edx,%eax
xorps  %xmm2,%xmm2
xor    %ecx,%ecx
xorps  %xmm7,%xmm7
xchg   %ax,%ax
movss  (%r15,%rcx,4),%xmm4
subss  %xmm3,%xmm4
movss  0x0(%rbp,%rcx,4),%xmm5
divss  %xmm5,%xmm4
movaps %xmm4,%xmm6
mulss  %xmm4,%xmm6
addss  %xmm6,%xmm7
mulss  (%rbx,%rcx,4),%xmm4
divss  %xmm5,%xmm4
addss  %xmm4,%xmm2
movss  %xmm2,(%r9)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1290 <nrfit+0x180>
jmp    12ef <nrfit+0x1df>
test   %edx,%edx
jle    12e9 <nrfit+0x1d9>
mov    %edx,%eax
cmp    $0x1,%edx
jne    1522 <nrfit+0x412>
xorps  %xmm2,%xmm2
xor    %ecx,%ecx
test   $0x1,%al
je     12ef <nrfit+0x1df>
jmp    158c <nrfit+0x47c>
xorps  %xmm2,%xmm2
xorps  %xmm7,%xmm7
divss  %xmm7,%xmm2
movss  %xmm2,(%r9)
movaps %xmm1,%xmm3
shufps $0x55,%xmm1,%xmm3
mulss  %xmm3,%xmm2
subss  %xmm2,%xmm1
divss  %xmm0,%xmm1
mulss  %xmm3,%xmm3
xorps  %xmm2,%xmm2
cvtss2sd %xmm0,%xmm2
mulss  %xmm7,%xmm0
divss  %xmm0,%xmm3
xorps  %xmm0,%xmm0
cvtss2sd %xmm3,%xmm0
addsd  0xcdb(%rip),%xmm0
mov    0x68(%rsp),%rax
divsd  %xmm2,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm1,(%r14)
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
mov    %r8d,0xc(%rsp)
mov    %edx,%r13d
jb     1355 <nrfit+0x245>
sqrtss %xmm0,%xmm0
jmp    136d <nrfit+0x25d>
mov    %rax,%r12
movaps %xmm7,0x10(%rsp)
call   1030 <sqrt@plt>
xorps  %xmm1,%xmm1
movaps 0x10(%rsp),%xmm7
mov    %r12,%rax
mov    0x78(%rsp),%r12
mov    0x70(%rsp),%r9
movss  %xmm0,(%rax)
movss  0xc7d(%rip),%xmm0
divss  %xmm7,%xmm0
ucomiss %xmm1,%xmm0
jb     1392 <nrfit+0x282>
sqrtss %xmm0,%xmm0
jmp    139c <nrfit+0x28c>
call   1030 <sqrt@plt>
mov    0x70(%rsp),%r9
mov    0x60(%rsp),%rax
mov    0xc(%rsp),%ecx
mov    0x80(%rsp),%r8
movss  %xmm0,(%r9)
movl   $0x0,(%r12)
test   %ecx,%ecx
je     1411 <nrfit+0x301>
test   %r13d,%r13d
jle    1405 <nrfit+0x2f5>
mov    %r13d,%ecx
xorps  %xmm0,%xmm0
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movss  (%rbx,%rdx,4),%xmm1
subss  (%r14),%xmm1
movss  (%rax),%xmm2
mulss  (%r15,%rdx,4),%xmm2
subss  %xmm2,%xmm1
divss  0x0(%rbp,%rdx,4),%xmm1
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm0
movss  %xmm0,(%r12)
add    $0x1,%rdx
cmp    %rdx,%rcx
jne    13d0 <nrfit+0x2c0>
movl   $0x3f800000,(%r8)
jmp    1513 <nrfit+0x403>
test   %r13d,%r13d
jle    14c2 <nrfit+0x3b2>
mov    %r13d,%ecx
cmp    $0x1,%r13d
jne    142a <nrfit+0x31a>
xorps  %xmm0,%xmm0
xor    %edx,%edx
jmp    1497 <nrfit+0x387>
mov    %ecx,%esi
and    $0xfffffffe,%esi
xorps  %xmm0,%xmm0
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rbx,%rdx,4),%xmm1
subss  (%r14),%xmm1
movss  (%rax),%xmm2
mulss  (%r15,%rdx,4),%xmm2
subss  %xmm2,%xmm1
mulss  %xmm1,%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%r12)
movss  0x4(%rbx,%rdx,4),%xmm0
subss  (%r14),%xmm0
movss  (%rax),%xmm2
mulss  0x4(%r15,%rdx,4),%xmm2
subss  %xmm2,%xmm0
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%r12)
add    $0x2,%rdx
cmp    %rdx,%rsi
jne    1440 <nrfit+0x330>
test   $0x1,%cl
je     14c2 <nrfit+0x3b2>
movss  (%rbx,%rdx,4),%xmm1
subss  (%r14),%xmm1
movss  (%rax),%xmm2
mulss  (%r15,%rdx,4),%xmm2
subss  %xmm2,%xmm1
mulss  %xmm1,%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%r12)
movl   $0x3f800000,(%r8)
movss  (%r12),%xmm0
add    $0xfffffffe,%r13d
xorps  %xmm1,%xmm1
cvtsi2ss %r13d,%xmm1
divss  %xmm1,%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jb     14ed <nrfit+0x3dd>
sqrtss %xmm0,%xmm0
jmp    14f8 <nrfit+0x3e8>
mov    %r9,%rbx
call   1030 <sqrt@plt>
mov    %rbx,%r9
mov    0x68(%rsp),%rax
movss  (%rax),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%rax)
mulss  (%r9),%xmm0
movss  %xmm0,(%r9)
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %eax,%esi
and    $0xfffffffe,%esi
xorps  %xmm2,%xmm2
xor    %ecx,%ecx
nopl   0x0(%rax)
movss  (%r15,%rcx,4),%xmm4
subss  %xmm3,%xmm4
movss  (%rbx,%rcx,4),%xmm5
unpcklps %xmm4,%xmm5
shufps $0x0,%xmm4,%xmm4
mulps  %xmm5,%xmm4
addps  %xmm2,%xmm4
movss  %xmm4,(%r9)
movss  0x4(%r15,%rcx,4),%xmm2
subss  %xmm3,%xmm2
movss  0x4(%rbx,%rcx,4),%xmm5
unpcklps %xmm2,%xmm5
shufps $0x0,%xmm2,%xmm2
mulps  %xmm5,%xmm2
addps  %xmm4,%xmm2
movss  %xmm2,(%r9)
add    $0x2,%rcx
cmp    %rcx,%rsi
jne    1530 <nrfit+0x420>
movaps %xmm2,%xmm7
shufps $0x55,%xmm2,%xmm7
test   $0x1,%al
je     12ef <nrfit+0x1df>
movss  (%r15,%rcx,4),%xmm4
subss  %xmm3,%xmm4
movss  (%rbx,%rcx,4),%xmm3
unpcklps %xmm4,%xmm3
shufps $0x0,%xmm4,%xmm4
mulps  %xmm3,%xmm4
addps  %xmm4,%xmm2
movss  %xmm2,(%r9)
movaps %xmm2,%xmm7
shufps $0x55,%xmm2,%xmm7
jmp    12ef <nrfit+0x1df>

LLVM IR (-O0):
define void @nrfit(float* %0, float* %1, i32 %2, float* %3, i32 %4, float* %5, float* %6, float* %7, float* %8, float* %9, float* %10) {
%12 = alloca float*
%13 = alloca float*
%14 = alloca i32
%15 = alloca float*
%16 = alloca i32
%17 = alloca float*
%18 = alloca float*
%19 = alloca float*
%20 = alloca float*
%21 = alloca float*
%22 = alloca float*
%23 = alloca i32
%24 = alloca float
%25 = alloca float
%26 = alloca float
%27 = alloca float
%28 = alloca float
%29 = alloca float
%30 = alloca float
%31 = alloca float
store float* %0, float** %12
store float* %1, float** %13
store i32 %2, i32* %14
store float* %3, float** %15
store i32 %4, i32* %16
store float* %5, float** %17
store float* %6, float** %18
store float* %7, float** %19
store float* %8, float** %20
store float* %9, float** %21
store float* %10, float** %22
store float 0.000000e+00, float* %27
store float 0.000000e+00, float* %28
store float 0.000000e+00, float* %29
%32 = load float*, float** %18
store float 0.000000e+00, float* %32
%33 = load i32, i32* %16
%34 = icmp ne i32 %33, 0
br i1 %34, label %35, label %78
35: ; preds = %11
store float 0.000000e+00, float* %30
store i32 0, i32* %23
br label %36
36: ; preds = %74, %35
%37 = load i32, i32* %23
%38 = load i32, i32* %14
%39 = icmp slt i32 %37, %38
br i1 %39, label %40, label %77
40: ; preds = %36
%41 = load float*, float** %15
%42 = load i32, i32* %23
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds float, float* %41, i64 %43
%45 = load float, float* %44
%46 = load float*, float** %15
%47 = load i32, i32* %23
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds float, float* %46, i64 %48
%50 = load float, float* %49
%51 = fmul float %45, %50
%52 = fpext float %51 to double
%53 = fdiv double 1.000000e+00, %52
%54 = fptrunc double %53 to float
store float %54, float* %24
%55 = load float, float* %24
%56 = load float, float* %30
%57 = fadd float %56, %55
store float %57, float* %30
%58 = load float*, float** %12
%59 = load i32, i32* %23
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds float, float* %58, i64 %60
%62 = load float, float* %61
%63 = load float, float* %24
%64 = load float, float* %27
%65 = call float @llvm.fmuladd.f32(float %62, float %63, float %64)
store float %65, float* %27
%66 = load float*, float** %13
%67 = load i32, i32* %23
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds float, float* %66, i64 %68
%70 = load float, float* %69
%71 = load float, float* %24
%72 = load float, float* %28
%73 = call float @llvm.fmuladd.f32(float %70, float %71, float %72)
store float %73, float* %28
br label %74
74: ; preds = %40
%75 = load i32, i32* %23
%76 = add nsw i32 %75, 1
store i32 %76, i32* %23
br label %36
77: ; preds = %36
br label %104
78: ; preds = %11
store i32 0, i32* %23
br label %79
79: ; preds = %98, %78
%80 = load i32, i32* %23
%81 = load i32, i32* %14
%82 = icmp slt i32 %80, %81
br i1 %82, label %83, label %101
83: ; preds = %79
%84 = load float*, float** %12
%85 = load i32, i32* %23
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds float, float* %84, i64 %86
%88 = load float, float* %87
%89 = load float, float* %27
%90 = fadd float %89, %88
store float %90, float* %27
%91 = load float*, float** %13
%92 = load i32, i32* %23
%93 = sext i32 %92 to i64
%94 = getelementptr inbounds float, float* %91, i64 %93
%95 = load float, float* %94
%96 = load float, float* %28
%97 = fadd float %96, %95
store float %97, float* %28
br label %98
98: ; preds = %83
%99 = load i32, i32* %23
%100 = add nsw i32 %99, 1
store i32 %100, i32* %23
br label %79
101: ; preds = %79
%102 = load i32, i32* %14
%103 = sitofp i32 %102 to float
store float %103, float* %30
br label %104
104: ; preds = %101, %77
%105 = load float, float* %27
%106 = load float, float* %30
%107 = fdiv float %105, %106
store float %107, float* %26
%108 = load i32, i32* %16
%109 = icmp ne i32 %108, 0
br i1 %109, label %110, label %153
110: ; preds = %104
store i32 0, i32* %23
br label %111
111: ; preds = %149, %110
%112 = load i32, i32* %23
%113 = load i32, i32* %14
%114 = icmp slt i32 %112, %113
br i1 %114, label %115, label %152
115: ; preds = %111
%116 = load float*, float** %12
%117 = load i32, i32* %23
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds float, float* %116, i64 %118
%120 = load float, float* %119
%121 = load float, float* %26
%122 = fsub float %120, %121
%123 = load float*, float** %15
%124 = load i32, i32* %23
%125 = sext i32 %124 to i64
%126 = getelementptr inbounds float, float* %123, i64 %125
%127 = load float, float* %126
%128 = fdiv float %122, %127
store float %128, float* %25
%129 = load float, float* %25
%130 = load float, float* %25
%131 = load float, float* %29
%132 = call float @llvm.fmuladd.f32(float %129, float %130, float %131)
store float %132, float* %29
%133 = load float, float* %25
%134 = load float*, float** %13
%135 = load i32, i32* %23
%136 = sext i32 %135 to i64
%137 = getelementptr inbounds float, float* %134, i64 %136
%138 = load float, float* %137
%139 = fmul float %133, %138
%140 = load float*, float** %15
%141 = load i32, i32* %23
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds float, float* %140, i64 %142
%144 = load float, float* %143
%145 = fdiv float %139, %144
%146 = load float*, float** %18
%147 = load float, float* %146
%148 = fadd float %147, %145
store float %148, float* %146
br label %149
149: ; preds = %115
%150 = load i32, i32* %23
%151 = add nsw i32 %150, 1
store i32 %151, i32* %23
br label %111
152: ; preds = %111
br label %183
153: ; preds = %104
store i32 0, i32* %23
br label %154
154: ; preds = %179, %153
%155 = load i32, i32* %23
%156 = load i32, i32* %14
%157 = icmp slt i32 %155, %156
br i1 %157, label %158, label %182
158: ; preds = %154
%159 = load float*, float** %12
%160 = load i32, i32* %23
%161 = sext i32 %160 to i64
%162 = getelementptr inbounds float, float* %159, i64 %161
%163 = load float, float* %162
%164 = load float, float* %26
%165 = fsub float %163, %164
store float %165, float* %25
%166 = load float, float* %25
%167 = load float, float* %25
%168 = load float, float* %29
%169 = call float @llvm.fmuladd.f32(float %166, float %167, float %168)
store float %169, float* %29
%170 = load float, float* %25
%171 = load float*, float** %13
%172 = load i32, i32* %23
%173 = sext i32 %172 to i64
%174 = getelementptr inbounds float, float* %171, i64 %173
%175 = load float, float* %174
%176 = load float*, float** %18
%177 = load float, float* %176
%178 = call float @llvm.fmuladd.f32(float %170, float %175, float %177)
store float %178, float* %176
br label %179
179: ; preds = %158
%180 = load i32, i32* %23
%181 = add nsw i32 %180, 1
store i32 %181, i32* %23
br label %154
182: ; preds = %154
br label %183
183: ; preds = %182, %152
%184 = load float, float* %29
%185 = load float*, float** %18
%186 = load float, float* %185
%187 = fdiv float %186, %184
store float %187, float* %185
%188 = load float, float* %28
%189 = load float, float* %27
%190 = load float*, float** %18
%191 = load float, float* %190
%192 = fneg float %189
%193 = call float @llvm.fmuladd.f32(float %192, float %191, float %188)
%194 = load float, float* %30
%195 = fdiv float %193, %194
%196 = load float*, float** %17
store float %195, float* %196
%197 = load float, float* %27
%198 = load float, float* %27
%199 = fmul float %197, %198
%200 = load float, float* %30
%201 = load float, float* %29
%202 = fmul float %200, %201
%203 = fdiv float %199, %202
%204 = fpext float %203 to double
%205 = fadd double 1.000000e+00, %204
%206 = load float, float* %30
%207 = fpext float %206 to double
%208 = fdiv double %205, %207
%209 = fptrunc double %208 to float
%210 = call float @sqrt(float %209)
%211 = load float*, float** %19
store float %210, float* %211
%212 = load float, float* %29
%213 = fpext float %212 to double
%214 = fdiv double 1.000000e+00, %213
%215 = fptrunc double %214 to float
%216 = call float @sqrt(float %215)
%217 = load float*, float** %20
store float %216, float* %217
%218 = load float*, float** %21
store float 0.000000e+00, float* %218
%219 = load i32, i32* %16
%220 = icmp eq i32 %219, 0
br i1 %220, label %221, label %284
221: ; preds = %183
store i32 0, i32* %23
br label %222
222: ; preds = %264, %221
%223 = load i32, i32* %23
%224 = load i32, i32* %14
%225 = icmp slt i32 %223, %224
br i1 %225, label %226, label %267
226: ; preds = %222
%227 = load float*, float** %13
%228 = load i32, i32* %23
%229 = sext i32 %228 to i64
%230 = getelementptr inbounds float, float* %227, i64 %229
%231 = load float, float* %230
%232 = load float*, float** %17
%233 = load float, float* %232
%234 = fsub float %231, %233
%235 = load float*, float** %18
%236 = load float, float* %235
%237 = load float*, float** %12
%238 = load i32, i32* %23
%239 = sext i32 %238 to i64
%240 = getelementptr inbounds float, float* %237, i64 %239
%241 = load float, float* %240
%242 = fneg float %236
%243 = call float @llvm.fmuladd.f32(float %242, float %241, float %234)
%244 = load float*, float** %13
%245 = load i32, i32* %23
%246 = sext i32 %245 to i64
%247 = getelementptr inbounds float, float* %244, i64 %246
%248 = load float, float* %247
%249 = load float*, float** %17
%250 = load float, float* %249
%251 = fsub float %248, %250
%252 = load float*, float** %18
%253 = load float, float* %252
%254 = load float*, float** %12
%255 = load i32, i32* %23
%256 = sext i32 %255 to i64
%257 = getelementptr inbounds float, float* %254, i64 %256
%258 = load float, float* %257
%259 = fneg float %253
%260 = call float @llvm.fmuladd.f32(float %259, float %258, float %251)
%261 = load float*, float** %21
%262 = load float, float* %261
%263 = call float @llvm.fmuladd.f32(float %243, float %260, float %262)
store float %263, float* %261
br label %264
264: ; preds = %226
%265 = load i32, i32* %23
%266 = add nsw i32 %265, 1
store i32 %266, i32* %23
br label %222
267: ; preds = %222
%268 = load float*, float** %22
store float 1.000000e+00, float* %268
%269 = load float*, float** %21
%270 = load float, float* %269
%271 = load i32, i32* %14
%272 = sub nsw i32 %271, 2
%273 = sitofp i32 %272 to float
%274 = fdiv float %270, %273
%275 = call float @sqrt(float %274)
store float %275, float* %31
%276 = load float, float* %31
%277 = load float*, float** %19
%278 = load float, float* %277
%279 = fmul float %278, %276
store float %279, float* %277
%280 = load float, float* %31
%281 = load float*, float** %20
%282 = load float, float* %281
%283 = fmul float %282, %280
store float %283, float* %281
br label %344
284: ; preds = %183
store i32 0, i32* %23
br label %285
285: ; preds = %339, %284
%286 = load i32, i32* %23
%287 = load i32, i32* %14
%288 = icmp slt i32 %286, %287
br i1 %288, label %289, label %342
289: ; preds = %285
%290 = load float*, float** %13
%291 = load i32, i32* %23
%292 = sext i32 %291 to i64
%293 = getelementptr inbounds float, float* %290, i64 %292
%294 = load float, float* %293
%295 = load float*, float** %17
%296 = load float, float* %295
%297 = fsub float %294, %296
%298 = load float*, float** %18
%299 = load float, float* %298
%300 = load float*, float** %12
%301 = load i32, i32* %23
%302 = sext i32 %301 to i64
%303 = getelementptr inbounds float, float* %300, i64 %302
%304 = load float, float* %303
%305 = fneg float %299
%306 = call float @llvm.fmuladd.f32(float %305, float %304, float %297)
%307 = load float*, float** %15
%308 = load i32, i32* %23
%309 = sext i32 %308 to i64
%310 = getelementptr inbounds float, float* %307, i64 %309
%311 = load float, float* %310
%312 = fdiv float %306, %311
%313 = load float*, float** %13
%314 = load i32, i32* %23
%315 = sext i32 %314 to i64
%316 = getelementptr inbounds float, float* %313, i64 %315
%317 = load float, float* %316
%318 = load float*, float** %17
%319 = load float, float* %318
%320 = fsub float %317, %319
%321 = load float*, float** %18
%322 = load float, float* %321
%323 = load float*, float** %12
%324 = load i32, i32* %23
%325 = sext i32 %324 to i64
%326 = getelementptr inbounds float, float* %323, i64 %325
%327 = load float, float* %326
%328 = fneg float %322
%329 = call float @llvm.fmuladd.f32(float %328, float %327, float %320)
%330 = load float*, float** %15
%331 = load i32, i32* %23
%332 = sext i32 %331 to i64
%333 = getelementptr inbounds float, float* %330, i64 %332
%334 = load float, float* %333
%335 = fdiv float %329, %334
%336 = load float*, float** %21
%337 = load float, float* %336
%338 = call float @llvm.fmuladd.f32(float %312, float %335, float %337)
store float %338, float* %336
br label %339
339: ; preds = %289
%340 = load i32, i32* %23
%341 = add nsw i32 %340, 1
store i32 %341, i32* %23
br label %285
342: ; preds = %285
%343 = load float*, float** %22
store float 1.000000e+00, float* %343
br label %344
344: ; preds = %342, %267
ret void
}

LLVM IR (-O2):
define void @nrfit(float* nocapture %0, float* nocapture %1, i32 %2, float* nocapture %3, i32 %4, float* nocapture %5, float* nocapture %6, float* nocapture %7, float* nocapture %8, float* nocapture %9, float* nocapture writeonly %10) {
store float 0.000000e+00, float* %6
%12 = icmp eq i32 %4, 0
%13 = icmp sgt i32 %2, 0
br i1 %12, label %17, label %14
14: ; preds = %11
br i1 %13, label %15, label %193
15: ; preds = %14
%16 = zext i32 %2 to i64
br label %25
17: ; preds = %11
br i1 %13, label %18, label %102
18: ; preds = %17
%19 = zext i32 %2 to i64
%20 = add nsw i64 %19, -1
%21 = and i64 %19, 3
%22 = icmp ult i64 %20, 3
br i1 %22, label %83, label %23
23: ; preds = %18
%24 = and i64 %19, 4294967292
br label %45
25: ; preds = %15, %25
%26 = phi i64 [ 0, %15 ], [ %43, %25 ]
%27 = phi float [ 0.000000e+00, %15 ], [ %33, %25 ]
%28 = phi <2 x float> [ zeroinitializer, %15 ], [ %42, %25 ]
%29 = getelementptr inbounds float, float* %3, i64 %26
%30 = load float, float* %29
%31 = fmul float %30, %30
%32 = fdiv float 1.000000e+00, %31
%33 = fadd float %27, %32
%34 = getelementptr inbounds float, float* %0, i64 %26
%35 = load float, float* %34
%36 = getelementptr inbounds float, float* %1, i64 %26
%37 = load float, float* %36
%38 = insertelement <2 x float> poison, float %37, i64 0
%39 = insertelement <2 x float> %38, float %35, i64 1
%40 = insertelement <2 x float> poison, float %32, i64 0
%41 = shufflevector <2 x float> %40, <2 x float> poison, <2 x i32> zeroinitializer
%42 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %39, <2 x float> %41, <2 x float> %28)
%43 = add nuw nsw i64 %26, 1
%44 = icmp eq i64 %43, %16
br i1 %44, label %105, label %25
45: ; preds = %45, %23
%46 = phi i64 [ 0, %23 ], [ %80, %45 ]
%47 = phi <2 x float> [ zeroinitializer, %23 ], [ %79, %45 ]
%48 = phi i64 [ 0, %23 ], [ %81, %45 ]
%49 = getelementptr inbounds float, float* %0, i64 %46
%50 = load float, float* %49
%51 = getelementptr inbounds float, float* %1, i64 %46
%52 = load float, float* %51
%53 = insertelement <2 x float> poison, float %52, i64 0
%54 = insertelement <2 x float> %53, float %50, i64 1
%55 = fadd <2 x float> %47, %54
%56 = or i64 %46, 1
%57 = getelementptr inbounds float, float* %0, i64 %56
%58 = load float, float* %57
%59 = getelementptr inbounds float, float* %1, i64 %56
%60 = load float, float* %59
%61 = insertelement <2 x float> poison, float %60, i64 0
%62 = insertelement <2 x float> %61, float %58, i64 1
%63 = fadd <2 x float> %55, %62
%64 = or i64 %46, 2
%65 = getelementptr inbounds float, float* %0, i64 %64
%66 = load float, float* %65
%67 = getelementptr inbounds float, float* %1, i64 %64
%68 = load float, float* %67
%69 = insertelement <2 x float> poison, float %68, i64 0
%70 = insertelement <2 x float> %69, float %66, i64 1
%71 = fadd <2 x float> %63, %70
%72 = or i64 %46, 3
%73 = getelementptr inbounds float, float* %0, i64 %72
%74 = load float, float* %73
%75 = getelementptr inbounds float, float* %1, i64 %72
%76 = load float, float* %75
%77 = insertelement <2 x float> poison, float %76, i64 0
%78 = insertelement <2 x float> %77, float %74, i64 1
%79 = fadd <2 x float> %71, %78
%80 = add nuw nsw i64 %46, 4
%81 = add i64 %48, 4
%82 = icmp eq i64 %81, %24
br i1 %82, label %83, label %45
83: ; preds = %45, %18
%84 = phi <2 x float> [ undef, %18 ], [ %79, %45 ]
%85 = phi i64 [ 0, %18 ], [ %80, %45 ]
%86 = phi <2 x float> [ zeroinitializer, %18 ], [ %79, %45 ]
%87 = icmp eq i64 %21, 0
br i1 %87, label %102, label %88
88: ; preds = %83, %88
%89 = phi i64 [ %99, %88 ], [ %85, %83 ]
%90 = phi <2 x float> [ %98, %88 ], [ %86, %83 ]
%91 = phi i64 [ %100, %88 ], [ 0, %83 ]
%92 = getelementptr inbounds float, float* %0, i64 %89
%93 = load float, float* %92
%94 = getelementptr inbounds float, float* %1, i64 %89
%95 = load float, float* %94
%96 = insertelement <2 x float> poison, float %95, i64 0
%97 = insertelement <2 x float> %96, float %93, i64 1
%98 = fadd <2 x float> %90, %97
%99 = add nuw nsw i64 %89, 1
%100 = add i64 %91, 1
%101 = icmp eq i64 %100, %21
br i1 %101, label %102, label %88
102: ; preds = %83, %88, %17
%103 = phi <2 x float> [ zeroinitializer, %17 ], [ %84, %83 ], [ %98, %88 ]
%104 = sitofp i32 %2 to float
br label %107
105: ; preds = %25
%106 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %39, <2 x float> %41, <2 x float> %28)
br label %107
107: ; preds = %105, %102
%108 = phi float [ %104, %102 ], [ %33, %105 ]
%109 = phi <2 x float> [ %103, %102 ], [ %106, %105 ]
%110 = extractelement <2 x float> %109, i64 1
%111 = fdiv float %110, %108
%112 = icmp sgt i32 %2, 0
br i1 %12, label %116, label %113
113: ; preds = %107
br i1 %112, label %114, label %193
114: ; preds = %113
%115 = zext i32 %2 to i64
br label %123
116: ; preds = %107
br i1 %112, label %117, label %193
117: ; preds = %116
%118 = zext i32 %2 to i64
%119 = and i64 %118, 1
%120 = icmp eq i32 %2, 1
br i1 %120, label %174, label %121
121: ; preds = %117
%122 = and i64 %118, 4294967294
br label %141
123: ; preds = %114, %123
%124 = phi float [ 0.000000e+00, %114 ], [ %138, %123 ]
%125 = phi i64 [ 0, %114 ], [ %139, %123 ]
%126 = phi float [ 0.000000e+00, %114 ], [ %133, %123 ]
%127 = getelementptr inbounds float, float* %0, i64 %125
%128 = load float, float* %127
%129 = fsub float %128, %111
%130 = getelementptr inbounds float, float* %3, i64 %125
%131 = load float, float* %130
%132 = fdiv float %129, %131
%133 = tail call float @llvm.fmuladd.f32(float %132, float %132, float %126)
%134 = getelementptr inbounds float, float* %1, i64 %125
%135 = load float, float* %134
%136 = fmul float %132, %135
%137 = fdiv float %136, %131
%138 = fadd float %124, %137
store float %138, float* %6
%139 = add nuw nsw i64 %125, 1
%140 = icmp eq i64 %139, %115
br i1 %140, label %193, label %123
141: ; preds = %141, %121
%142 = phi i64 [ 0, %121 ], [ %168, %141 ]
%143 = phi <2 x float> [ zeroinitializer, %121 ], [ %166, %141 ]
%144 = phi i64 [ 0, %121 ], [ %169, %141 ]
%145 = getelementptr inbounds float, float* %0, i64 %142
%146 = load float, float* %145
%147 = fsub float %146, %111
%148 = getelementptr inbounds float, float* %1, i64 %142
%149 = load float, float* %148
%150 = insertelement <2 x float> poison, float %147, i64 0
%151 = shufflevector <2 x float> %150, <2 x float> poison, <2 x i32> zeroinitializer
%152 = insertelement <2 x float> poison, float %149, i64 0
%153 = insertelement <2 x float> %152, float %147, i64 1
%154 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %151, <2 x float> %153, <2 x float> %143)
%155 = extractelement <2 x float> %154, i64 0
store float %155, float* %6
%156 = or i64 %142, 1
%157 = getelementptr inbounds float, float* %0, i64 %156
%158 = load float, float* %157
%159 = fsub float %158, %111
%160 = getelementptr inbounds float, float* %1, i64 %156
%161 = load float, float* %160
%162 = insertelement <2 x float> poison, float %159, i64 0
%163 = shufflevector <2 x float> %162, <2 x float> poison, <2 x i32> zeroinitializer
%164 = insertelement <2 x float> poison, float %161, i64 0
%165 = insertelement <2 x float> %164, float %159, i64 1
%166 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %163, <2 x float> %165, <2 x float> %154)
%167 = extractelement <2 x float> %166, i64 0
store float %167, float* %6
%168 = add nuw nsw i64 %142, 2
%169 = add i64 %144, 2
%170 = icmp eq i64 %169, %122
br i1 %170, label %171, label %141
171: ; preds = %141
%172 = extractelement <2 x float> %166, i64 0
%173 = extractelement <2 x float> %166, i64 1
br label %174
174: ; preds = %171, %117
%175 = phi float [ undef, %117 ], [ %172, %171 ]
%176 = phi float [ undef, %117 ], [ %173, %171 ]
%177 = phi i64 [ 0, %117 ], [ %168, %171 ]
%178 = phi <2 x float> [ zeroinitializer, %117 ], [ %166, %171 ]
%179 = icmp eq i64 %119, 0
br i1 %179, label %193, label %180
180: ; preds = %174
%181 = getelementptr inbounds float, float* %0, i64 %177
%182 = load float, float* %181
%183 = fsub float %182, %111
%184 = getelementptr inbounds float, float* %1, i64 %177
%185 = load float, float* %184
%186 = insertelement <2 x float> poison, float %183, i64 0
%187 = shufflevector <2 x float> %186, <2 x float> poison, <2 x i32> zeroinitializer
%188 = insertelement <2 x float> poison, float %185, i64 0
%189 = insertelement <2 x float> %188, float %183, i64 1
%190 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %187, <2 x float> %189, <2 x float> %178)
%191 = extractelement <2 x float> %190, i64 0
store float %191, float* %6
%192 = extractelement <2 x float> %190, i64 1
br label %193
193: ; preds = %123, %180, %174, %14, %113, %116
%194 = phi float [ %108, %116 ], [ %108, %113 ], [ 0.000000e+00, %14 ], [ %108, %174 ], [ %108, %180 ], [ %108, %123 ]
%195 = phi float [ 0.000000e+00, %116 ], [ 0.000000e+00, %113 ], [ 0.000000e+00, %14 ], [ %175, %174 ], [ %191, %180 ], [ %138, %123 ]
%196 = phi float [ 0.000000e+00, %116 ], [ 0.000000e+00, %113 ], [ 0.000000e+00, %14 ], [ %176, %174 ], [ %192, %180 ], [ %133, %123 ]
%197 = phi <2 x float> [ %109, %116 ], [ %109, %113 ], [ zeroinitializer, %14 ], [ %109, %174 ], [ %109, %180 ], [ %109, %123 ]
%198 = fdiv float %195, %196
store float %198, float* %6
%199 = extractelement <2 x float> %197, i64 1
%200 = fneg float %199
%201 = extractelement <2 x float> %197, i64 0
%202 = tail call float @llvm.fmuladd.f32(float %200, float %198, float %201)
%203 = fdiv float %202, %194
store float %203, float* %5
%204 = fmul float %199, %199
%205 = fmul float %194, %196
%206 = fdiv float %204, %205
%207 = fpext float %206 to double
%208 = fadd double %207, 1.000000e+00
%209 = fpext float %194 to double
%210 = fdiv double %208, %209
%211 = fptrunc double %210 to float
%212 = tail call float @sqrt(float %211)
store float %212, float* %7
%213 = fdiv float 1.000000e+00, %196
%214 = tail call float @sqrt(float %213)
store float %214, float* %8
store float 0.000000e+00, float* %9
%215 = icmp sgt i32 %2, 0
br i1 %12, label %219, label %216
216: ; preds = %193
br i1 %215, label %217, label %297
217: ; preds = %216
%218 = zext i32 %2 to i64
br label %279
219: ; preds = %193
br i1 %215, label %220, label %269
220: ; preds = %219
%221 = zext i32 %2 to i64
%222 = and i64 %221, 1
%223 = icmp eq i32 %2, 1
br i1 %223, label %254, label %224
224: ; preds = %220
%225 = and i64 %221, 4294967294
br label %226
226: ; preds = %226, %224
%227 = phi float [ 0.000000e+00, %224 ], [ %250, %226 ]
%228 = phi i64 [ 0, %224 ], [ %251, %226 ]
%229 = phi i64 [ 0, %224 ], [ %252, %226 ]
%230 = getelementptr inbounds float, float* %1, i64 %228
%231 = load float, float* %230
%232 = load float, float* %5
%233 = fsub float %231, %232
%234 = load float, float* %6
%235 = getelementptr inbounds float, float* %0, i64 %228
%236 = load float, float* %235
%237 = fneg float %234
%238 = tail call float @llvm.fmuladd.f32(float %237, float %236, float %233)
%239 = tail call float @llvm.fmuladd.f32(float %238, float %238, float %227)
store float %239, float* %9
%240 = or i64 %228, 1
%241 = getelementptr inbounds float, float* %1, i64 %240
%242 = load float, float* %241
%243 = load float, float* %5
%244 = fsub float %242, %243
%245 = load float, float* %6
%246 = getelementptr inbounds float, float* %0, i64 %240
%247 = load float, float* %246
%248 = fneg float %245
%249 = tail call float @llvm.fmuladd.f32(float %248, float %247, float %244)
%250 = tail call float @llvm.fmuladd.f32(float %249, float %249, float %239)
store float %250, float* %9
%251 = add nuw nsw i64 %228, 2
%252 = add i64 %229, 2
%253 = icmp eq i64 %252, %225
br i1 %253, label %254, label %226
254: ; preds = %226, %220
%255 = phi float [ 0.000000e+00, %220 ], [ %250, %226 ]
%256 = phi i64 [ 0, %220 ], [ %251, %226 ]
%257 = icmp eq i64 %222, 0
br i1 %257, label %269, label %258
258: ; preds = %254
%259 = getelementptr inbounds float, float* %1, i64 %256
%260 = load float, float* %259
%261 = load float, float* %5
%262 = fsub float %260, %261
%263 = load float, float* %6
%264 = getelementptr inbounds float, float* %0, i64 %256
%265 = load float, float* %264
%266 = fneg float %263
%267 = tail call float @llvm.fmuladd.f32(float %266, float %265, float %262)
%268 = tail call float @llvm.fmuladd.f32(float %267, float %267, float %255)
store float %268, float* %9
br label %269
269: ; preds = %258, %254, %219
store float 1.000000e+00, float* %10
%270 = load float, float* %9
%271 = add nsw i32 %2, -2
%272 = sitofp i32 %271 to float
%273 = fdiv float %270, %272
%274 = tail call float @sqrt(float %273)
%275 = load float, float* %7
%276 = fmul float %274, %275
store float %276, float* %7
%277 = load float, float* %8
%278 = fmul float %274, %277
store float %278, float* %8
br label %298
279: ; preds = %217, %279
%280 = phi float [ 0.000000e+00, %217 ], [ %294, %279 ]
%281 = phi i64 [ 0, %217 ], [ %295, %279 ]
%282 = getelementptr inbounds float, float* %1, i64 %281
%283 = load float, float* %282
%284 = load float, float* %5
%285 = fsub float %283, %284
%286 = load float, float* %6
%287 = getelementptr inbounds float, float* %0, i64 %281
%288 = load float, float* %287
%289 = fneg float %286
%290 = tail call float @llvm.fmuladd.f32(float %289, float %288, float %285)
%291 = getelementptr inbounds float, float* %3, i64 %281
%292 = load float, float* %291
%293 = fdiv float %290, %292
%294 = tail call float @llvm.fmuladd.f32(float %293, float %293, float %280)
store float %294, float* %9
%295 = add nuw nsw i64 %281, 1
%296 = icmp eq i64 %295, %218
br i1 %296, label %297, label %279
297: ; preds = %279, %216
store float 1.000000e+00, float* %10
br label %298
298: ; preds = %297, %269
ret void
}