#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 pow(float, int);
float sqrt(float);
void slfit(float *x, float *y, float *z, int ndata, float *sig, int weight,
           float *a, float *b, float *siga, float *sigb) {
  float sx = .0f;
  float sy = .0f;
  float st2 = .0f;
  float ss = .0f;
  *b = .0f;
  int q = 0;
  float sigdat = 0.0f;
  float t = 0.0f;
  float wt = 0.0f;
  float chi2 = 0.0f;
  float sxoss = 0.0f;
  int useweight = weight;
  int i = 0;
  for (i = 0; i < ndata; i++) {
    *(z + i) = *(y + i);
    if (*(sig + i) == -999.0) {
      useweight = 0;
      *(z + i) = 0.75 * (*(y + i));
    }
  }
  if (useweight) {
    for (i = 0; i < ndata; i++) {
      wt = 1.0 / (*(sig + i) * (*(sig + i)));
      ss = ss + wt;
      sx = sx + *(x + i) * wt;
      sy = sy + *(z + i) * wt;
    }
  } else {
    for (i = 0; i < ndata; i++) {
      sx = sx + *(x + i);
      sy = sy + *(z + i);
    }
    ss = (float)ndata;
  }
  sxoss = sx / ss;
  if (useweight) {
    for (i = 0; i < ndata; i++) {
      t = (*(x + i) - sxoss) / *(sig + i);
      *b = *b + t * (*(z + i) / *(sig + i));
      st2 = st2 + t * t;
    }
  } else {
    for (i = 0; i < ndata; i++) {
      t = (*(x + i) - sxoss);
      st2 = st2 + t * t;
      *b = *b + t * (*(z + i));
    }
  }
  *b = *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 (!useweight) {
    for (i = 0; i < ndata; i++) {
      chi2 = pow((chi2 + ((*(z + i)) - *a - *b * (*(x + i)))), 2);
    }
    q = 1;
    sigdat = 0.0f;
    if (ndata > 2)
      sigdat = sqrt(chi2 / (ndata - 2));
    *siga = *siga * sigdat;
    *sigb = *sigb * sigdat;
  }
}


x86 assembly (clang, -O0):
<slfit>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
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    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
movl   $0x0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x48(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x58(%rbp)
movl   $0x0,-0x5c(%rbp)
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1254 <slfit+0x134>
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe0c(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1241 <slfit+0x121>
jp     1241 <slfit+0x121>
movl   $0x0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xde4(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1246 <slfit+0x126>
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    11bd <slfit+0x9d>
cmpl   $0x0,-0x58(%rbp)
je     130e <slfit+0x1ee>
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1309 <slfit+0x1e9>
mov    -0x28(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd69(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x3c(%rbp),%xmm0
addss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x30(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x4c(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x34(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x4c(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    1265 <slfit+0x145>
jmp    1367 <slfit+0x247>
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    135d <slfit+0x23d>
movss  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x34(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    1315 <slfit+0x1f5>
cvtsi2ssl -0x1c(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x30(%rbp),%xmm0
divss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x54(%rbp)
cmpl   $0x0,-0x58(%rbp)
je     141d <slfit+0x2fd>
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1418 <slfit+0x2f8>
mov    -0x8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x54(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x5c(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
mov    0x18(%rbp),%rax
movss  (%rax),%xmm1
movss  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x5c(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    0x18(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x38(%rbp),%xmm1
movss  -0x48(%rbp),%xmm0
movss  -0x48(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    1387 <slfit+0x267>
jmp    14a0 <slfit+0x380>
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    149b <slfit+0x37b>
mov    -0x8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x54(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x38(%rbp),%xmm1
movss  -0x48(%rbp),%xmm0
movss  -0x48(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
mov    0x18(%rbp),%rax
movss  (%rax),%xmm1
movss  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    0x18(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    1424 <slfit+0x304>
jmp    14a0 <slfit+0x380>
mov    0x18(%rbp),%rax
movss  (%rax),%xmm0
divss  -0x38(%rbp),%xmm0
mov    0x18(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x34(%rbp),%xmm1
movss  -0x30(%rbp),%xmm0
mov    0x18(%rbp),%rax
movss  (%rax),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
divss  -0x3c(%rbp),%xmm0
mov    0x10(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x30(%rbp),%xmm0
mulss  -0x30(%rbp),%xmm0
movss  -0x3c(%rbp),%xmm1
mulss  -0x38(%rbp),%xmm1
divss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xaf3(%rip),%xmm0
addsd  %xmm1,%xmm0
movss  -0x3c(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1040 <sqrt@plt>
mov    0x20(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x38(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xac0(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1040 <sqrt@plt>
mov    0x28(%rbp),%rax
movss  %xmm0,(%rax)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
cmpl   $0x0,-0x58(%rbp)
jne    163f <slfit+0x51f>
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    15df <slfit+0x4bf>
movss  -0x50(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    0x10(%rbp),%rax
subss  (%rax),%xmm2
mov    0x18(%rbp),%rax
movss  (%rax),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm3
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
addss  %xmm1,%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movss  %xmm0,-0x50(%rbp)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    156e <slfit+0x44e>
movl   $0x1,-0x40(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
cmpl   $0x2,-0x1c(%rbp)
jle    1615 <slfit+0x4f5>
movss  -0x50(%rbp),%xmm0
mov    -0x1c(%rbp),%eax
sub    $0x2,%eax
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
call   1040 <sqrt@plt>
movss  %xmm0,-0x44(%rbp)
mov    0x20(%rbp),%rax
movss  (%rax),%xmm0
mulss  -0x44(%rbp),%xmm0
mov    0x20(%rbp),%rax
movss  %xmm0,(%rax)
mov    0x28(%rbp),%rax
movss  (%rax),%xmm0
mulss  -0x44(%rbp),%xmm0
mov    0x28(%rbp),%rax
movss  %xmm0,(%rax)
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<slfit>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rdx,%rbp
mov    %rdi,%rbx
mov    0x68(%rsp),%r15
mov    0x60(%rsp),%r12
movl   $0x0,(%r15)
mov    %ecx,0xc(%rsp)
test   %ecx,%ecx
jle    12b6 <slfit+0x196>
mov    0xc(%rsp),%eax
mov    %rax,%rcx
neg    %rcx
xor    %edx,%edx
movss  0xe9b(%rip),%xmm0
movss  0xe9b(%rip),%xmm1
add    $0x1,%rdx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  -0x4(%rsi,%rdx,4),%xmm2
movss  %xmm2,-0x4(%rbp,%rdx,4)
movss  -0x4(%r8,%rdx,4),%xmm2
ucomiss %xmm0,%xmm2
jne    119a <slfit+0x7a>
jnp    11b0 <slfit+0x90>
lea    (%rcx,%rdx,1),%rdi
add    $0x1,%rdi
add    $0x1,%rdx
cmp    $0x1,%rdi
jne    1180 <slfit+0x60>
jmp    11d0 <slfit+0xb0>
xchg   %ax,%ax
movss  -0x4(%rsi,%rdx,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,-0x4(%rbp,%rdx,4)
xor    %r9d,%r9d
mov    %rcx,%rdi
add    %rdx,%rdi
jne    116d <slfit+0x4d>
jmp    12c3 <slfit+0x1a3>
test   %r9d,%r9d
je     12c3 <slfit+0x1a3>
cmpl   $0x0,0xc(%rsp)
jle    12bb <slfit+0x19b>
xorps  %xmm1,%xmm1
xorps  %xmm2,%xmm2
xor    %ecx,%ecx
movss  0xe10(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%r8,%rcx,4),%xmm3
mulss  %xmm3,%xmm3
movaps %xmm0,%xmm4
divss  %xmm3,%xmm4
addss  %xmm4,%xmm2
movss  (%rbx,%rcx,4),%xmm3
movss  0x0(%rbp,%rcx,4),%xmm5
unpcklps %xmm3,%xmm5
shufps $0x0,%xmm4,%xmm4
mulps  %xmm5,%xmm4
addps  %xmm4,%xmm1
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1200 <slfit+0xe0>
movaps %xmm1,%xmm3
shufps $0x55,%xmm1,%xmm3
movaps %xmm3,%xmm0
cmpl   $0x0,0xc(%rsp)
jle    14bd <slfit+0x39d>
divss  %xmm2,%xmm3
movd   (%r15),%xmm4
movq   %xmm4,%xmm7
pxor   %xmm4,%xmm4
shufps $0xe2,%xmm4,%xmm7
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%rbx,%rcx,4),%xmm4
subss  %xmm3,%xmm4
movss  0x0(%rbp,%rcx,4),%xmm5
unpcklps %xmm5,%xmm4
movss  (%r8,%rcx,4),%xmm5
shufps $0x0,%xmm5,%xmm5
divps  %xmm5,%xmm4
movaps %xmm4,%xmm5
shufps $0x0,%xmm4,%xmm5
mulps  %xmm4,%xmm5
addps  %xmm5,%xmm7
movaps %xmm7,%xmm4
shufps $0x55,%xmm7,%xmm4
movss  %xmm4,(%r15)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1270 <slfit+0x150>
xor    %r13d,%r13d
jmp    12ff <slfit+0x1df>
test   %r9d,%r9d
je     12ea <slfit+0x1ca>
xorps  %xmm2,%xmm2
xor    %r13d,%r13d
jmp    12f6 <slfit+0x1d6>
cmpl   $0x0,0xc(%rsp)
jle    12ea <slfit+0x1ca>
lea    -0x1(%rax),%r8
mov    %eax,%edx
and    $0x3,%edx
cmp    $0x3,%r8
jae    14c5 <slfit+0x3a5>
xorps  %xmm1,%xmm1
xor    %esi,%esi
xorps  %xmm0,%xmm0
jmp    1518 <slfit+0x3f8>
xorps  %xmm2,%xmm2
cvtsi2ssl 0xc(%rsp),%xmm2
mov    $0x1,%r13b
xorps  %xmm1,%xmm1
xorps  %xmm0,%xmm0
xorps  %xmm7,%xmm7
movss  (%r15),%xmm3
divss  %xmm7,%xmm3
movss  %xmm3,(%r15)
mulss  %xmm0,%xmm3
subss  %xmm3,%xmm1
divss  %xmm2,%xmm1
mulss  %xmm0,%xmm0
xorps  %xmm3,%xmm3
cvtss2sd %xmm2,%xmm3
mulss  %xmm7,%xmm2
divss  %xmm2,%xmm0
cvtss2sd %xmm0,%xmm0
addsd  0xcd8(%rip),%xmm0
mov    0x70(%rsp),%rax
divsd  %xmm3,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm1,(%r12)
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jb     1359 <slfit+0x239>
sqrtss %xmm0,%xmm0
jmp    137c <slfit+0x25c>
movaps %xmm7,0x10(%rsp)
mov    %r12,%r15
mov    %rax,%r12
call   1040 <sqrt@plt>
xorps  %xmm1,%xmm1
mov    %r12,%rax
mov    %r15,%r12
mov    0x68(%rsp),%r15
movaps 0x10(%rsp),%xmm7
mov    0x78(%rsp),%r14
movss  %xmm0,(%rax)
movss  0xc77(%rip),%xmm0
divss  %xmm7,%xmm0
ucomiss %xmm1,%xmm0
jb     139c <slfit+0x27c>
sqrtss %xmm0,%xmm0
jmp    13a1 <slfit+0x281>
call   1040 <sqrt@plt>
movss  %xmm0,(%r14)
test   %r13b,%r13b
je     15fd <slfit+0x4dd>
cmpl   $0x0,0xc(%rsp)
jle    13d2 <slfit+0x2b2>
mov    0xc(%rsp),%eax
mov    %eax,%r14d
cmp    $0x1,%eax
mov    %r14,0x10(%rsp)
jne    13da <slfit+0x2ba>
xorps  %xmm0,%xmm0
xor    %r13d,%r13d
jmp    1451 <slfit+0x331>
xorps  %xmm1,%xmm1
jmp    15e2 <slfit+0x4c2>
and    $0xfffffffe,%r14d
xorps  %xmm0,%xmm0
xor    %r13d,%r13d
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  0x0(%rbp,%r13,4),%xmm1
subss  (%r12),%xmm1
movss  (%r15),%xmm2
mulss  (%rbx,%r13,4),%xmm2
subss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movaps %xmm1,%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movss  0x4(%rbp,%r13,4),%xmm1
subss  (%r12),%xmm1
movss  (%r15),%xmm2
mulss  0x4(%rbx,%r13,4),%xmm2
subss  %xmm2,%xmm1
addss  %xmm1,%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
add    $0x2,%r13
cmp    %r13,%r14
jne    13f0 <slfit+0x2d0>
testb  $0x1,0x10(%rsp)
je     1482 <slfit+0x362>
movss  0x0(%rbp,%r13,4),%xmm1
subss  (%r12),%xmm1
movss  (%r15),%xmm2
mulss  (%rbx,%r13,4),%xmm2
subss  %xmm2,%xmm1
addss  %xmm1,%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
xorps  %xmm1,%xmm1
mov    0xc(%rsp),%eax
cmp    $0x3,%eax
mov    0x78(%rsp),%r14
jl     15e2 <slfit+0x4c2>
add    $0xfffffffe,%eax
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jb     15da <slfit+0x4ba>
xorps  %xmm1,%xmm1
sqrtss %xmm0,%xmm1
jmp    15e2 <slfit+0x4c2>
xor    %r13d,%r13d
jmp    12fc <slfit+0x1dc>
mov    %eax,%edi
and    $0xfffffffc,%edi
xorps  %xmm1,%xmm1
xor    %esi,%esi
xorps  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
addss  (%rbx,%rsi,4),%xmm0
addss  0x0(%rbp,%rsi,4),%xmm1
addss  0x4(%rbx,%rsi,4),%xmm0
addss  0x4(%rbp,%rsi,4),%xmm1
addss  0x8(%rbx,%rsi,4),%xmm0
addss  0x8(%rbp,%rsi,4),%xmm1
addss  0xc(%rbx,%rsi,4),%xmm0
addss  0xc(%rbp,%rsi,4),%xmm1
add    $0x4,%rsi
cmp    %rsi,%rdi
jne    14e0 <slfit+0x3c0>
test   %rdx,%rdx
je     1543 <slfit+0x423>
lea    0x0(,%rsi,4),%rdi
add    %rbp,%rdi
lea    (%rbx,%rsi,4),%rsi
xor    %ecx,%ecx
xchg   %ax,%ax
addss  (%rsi,%rcx,4),%xmm0
addss  (%rdi,%rcx,4),%xmm1
add    $0x1,%rcx
cmp    %rcx,%rdx
jne    1530 <slfit+0x410>
mov    0xc(%rsp),%ecx
xorps  %xmm2,%xmm2
cvtsi2ss %ecx,%xmm2
mov    $0x1,%r13b
test   %ecx,%ecx
jle    12fc <slfit+0x1dc>
movaps %xmm0,%xmm4
divss  %xmm2,%xmm4
movss  (%r15),%xmm3
test   %r8,%r8
je     160c <slfit+0x4ec>
mov    %eax,%edx
and    $0xfffffffe,%edx
xorps  %xmm7,%xmm7
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movss  (%rbx,%rcx,4),%xmm5
subss  %xmm4,%xmm5
movaps %xmm5,%xmm6
mulss  %xmm5,%xmm6
addss  %xmm7,%xmm6
mulss  0x0(%rbp,%rcx,4),%xmm5
addss  %xmm3,%xmm5
movss  %xmm5,(%r15)
movss  0x4(%rbx,%rcx,4),%xmm3
subss  %xmm4,%xmm3
movaps %xmm3,%xmm7
mulss  %xmm3,%xmm7
addss  %xmm6,%xmm7
mulss  0x4(%rbp,%rcx,4),%xmm3
addss  %xmm5,%xmm3
movss  %xmm3,(%r15)
add    $0x2,%rcx
cmp    %rcx,%rdx
jne    1580 <slfit+0x460>
test   $0x1,%al
je     12ff <slfit+0x1df>
jmp    1619 <slfit+0x4f9>
call   1040 <sqrt@plt>
movaps %xmm0,%xmm1
mov    0x70(%rsp),%rax
movss  (%rax),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,(%rax)
mulss  (%r14),%xmm1
movss  %xmm1,(%r14)
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xorps  %xmm7,%xmm7
xor    %ecx,%ecx
test   $0x1,%al
je     12ff <slfit+0x1df>
movss  (%rbx,%rcx,4),%xmm5
subss  %xmm4,%xmm5
movaps %xmm5,%xmm4
mulss  %xmm5,%xmm4
addss  %xmm4,%xmm7
mulss  0x0(%rbp,%rcx,4),%xmm5
addss  %xmm3,%xmm5
movss  %xmm5,(%r15)
jmp    12ff <slfit+0x1df>

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

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