#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double sqrt(double);
void linreg(double *x, double *y, int nitems, double *a, double *b,
            double *minp, double *maxp, double *meanp, double *r) {
  double sx = 0.0;
  double sy = 0.0;
  double st2 = 0.0;
  double min = y[0];
  double max = 0.0;
  double ss, meanx, meany, t;
  double errx, erry, prodsum, sqerrx, sqerry;
  int i;
  *b = 0.0;
  for (i = 0; i < nitems; i++) {
    sx += x[i];
    sy += y[i];
    if (y[i] < min)
      min = y[i];
    if (y[i] > max)
      max = y[i];
  }
  ss = nitems;
  meanx = sx / ss;
  meany = *meanp = sy / ss;
  *minp = min;
  *maxp = max;
  for (i = 0; i < nitems; i++) {
    t = x[i] - meanx;
    st2 += t * t;
    *b += t * y[i];
  }
  *b /= st2;
  *a = (sy - sx * (*b)) / ss;
  prodsum = 0.0;
  sqerrx = 0.0;
  sqerry = 0.0;
  for (i = 0; i < nitems; i++) {
    errx = x[i] - meanx;
    erry = y[i] - meany;
    prodsum += errx * erry;
    sqerrx += errx * errx;
    sqerry += erry * erry;
  }
  *r = prodsum / (sqrt(sqerrx) * sqrt(sqerry));
}


x86 assembly (clang, -O0):
<linreg>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
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    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x28(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0xa4(%rbp)
mov    -0xa4(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1240 <linreg+0x130>
mov    -0x8(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x10(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x50(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11f7 <linreg+0xe7>
mov    -0x10(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x58(%rbp),%xmm0
jbe    1227 <linreg+0x117>
mov    -0x10(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    122c <linreg+0x11c>
mov    -0xa4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa4(%rbp)
jmp    1180 <linreg+0x70>
cvtsi2sdl -0x14(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x38(%rbp),%xmm0
divsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x40(%rbp),%xmm0
divsd  -0x60(%rbp),%xmm0
mov    0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  %xmm0,-0x70(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x58(%rbp),%xmm0
mov    0x10(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x0,-0xa4(%rbp)
mov    -0xa4(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1316 <linreg+0x206>
mov    -0x8(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x78(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0xa4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa4(%rbp)
jmp    1294 <linreg+0x184>
movsd  -0x48(%rbp),%xmm1
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
movsd  -0x40(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x60(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movl   $0x0,-0xa4(%rbp)
mov    -0xa4(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    145d <linreg+0x34d>
mov    -0x8(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  -0x88(%rbp),%xmm2
ff
movsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
movsd  -0x98(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0x88(%rbp),%xmm2
ff
movsd  -0xa0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0xa4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa4(%rbp)
jmp    1394 <linreg+0x284>
movsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
call   1030 <sqrt@plt>
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0xa0(%rbp),%xmm0
ff
call   1030 <sqrt@plt>
movsd  -0xb8(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xb0(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    0x20(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<linreg>:
push   %r15
push   %r14
push   %rbx
sub    $0x20,%rsp
mov    0x48(%rsp),%r11
mov    0x40(%rsp),%r10
movsd  (%rsi),%xmm3
movq   $0x0,(%r8)
mov    %edx,%r14d
test   %edx,%edx
jle    1289 <linreg+0x179>
mov    %r14,%r15
add    $0xffffffffffffffff,%r15
je     12b9 <linreg+0x1a9>
mov    %r14d,%ebx
and    $0xfffffffe,%ebx
xorpd  %xmm1,%xmm1
xorpd  %xmm0,%xmm0
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
movsd  (%rsi,%rax,8),%xmm2
movsd  0x8(%rsi,%rax,8),%xmm4
movapd %xmm2,%xmm5
movhpd (%rdi,%rax,8),%xmm5
addpd  %xmm1,%xmm5
movapd %xmm2,%xmm6
minsd  %xmm3,%xmm6
maxsd  %xmm0,%xmm2
movapd %xmm4,%xmm1
movhpd 0x8(%rdi,%rax,8),%xmm1
addpd  %xmm5,%xmm1
movapd %xmm4,%xmm3
minsd  %xmm6,%xmm3
maxsd  %xmm2,%xmm4
add    $0x2,%rax
movapd %xmm4,%xmm0
cmp    %rax,%rbx
jne    1160 <linreg+0x50>
test   $0x1,%r14b
je     11d7 <linreg+0xc7>
movsd  (%rsi,%rax,8),%xmm0
movapd %xmm0,%xmm2
movhpd (%rdi,%rax,8),%xmm2
addpd  %xmm2,%xmm1
movapd %xmm0,%xmm2
minsd  %xmm3,%xmm2
maxsd  %xmm4,%xmm0
movapd %xmm2,%xmm3
movapd %xmm0,%xmm4
xorps  %xmm2,%xmm2
cvtsi2sd %edx,%xmm2
movapd %xmm2,%xmm5
unpcklpd %xmm2,%xmm5
movapd %xmm1,%xmm0
divpd  %xmm5,%xmm0
movlpd %xmm0,(%r11)
movsd  %xmm3,(%r9)
movsd  %xmm4,(%r10)
test   %edx,%edx
jle    12ae <linreg+0x19e>
movsd  (%r8),%xmm3
movapd %xmm0,%xmm5
unpckhpd %xmm0,%xmm5
test   %r15,%r15
je     12d2 <linreg+0x1c2>
mov    %r14d,%ebx
and    $0xfffffffe,%ebx
xorpd  %xmm4,%xmm4
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rax,8),%xmm6
subsd  %xmm5,%xmm6
movapd %xmm6,%xmm7
mulsd  %xmm6,%xmm7
addsd  %xmm4,%xmm7
mulsd  (%rsi,%rax,8),%xmm6
addsd  %xmm3,%xmm6
movsd  %xmm6,(%r8)
movsd  0x8(%rdi,%rax,8),%xmm3
subsd  %xmm5,%xmm3
movapd %xmm3,%xmm4
mulsd  %xmm3,%xmm4
addsd  %xmm7,%xmm4
mulsd  0x8(%rsi,%rax,8),%xmm3
addsd  %xmm6,%xmm3
movsd  %xmm3,(%r8)
add    $0x2,%rax
cmp    %rax,%rbx
jne    1230 <linreg+0x120>
test   $0x1,%r14b
jne    12de <linreg+0x1ce>
jmp    1301 <linreg+0x1f1>
cvtsi2sd %edx,%xmm2
xorpd  %xmm0,%xmm0
divsd  %xmm2,%xmm0
movsd  %xmm0,(%r11)
unpcklpd %xmm0,%xmm0
movsd  %xmm3,(%r9)
movq   $0x0,(%r10)
xorpd  %xmm1,%xmm1
movsd  (%r8),%xmm3
xorpd  %xmm4,%xmm4
jmp    1301 <linreg+0x1f1>
xor    %eax,%eax
xorpd  %xmm1,%xmm1
xorpd  %xmm4,%xmm4
test   $0x1,%r14b
jne    11b1 <linreg+0xa1>
jmp    11d7 <linreg+0xc7>
xorpd  %xmm4,%xmm4
xor    %eax,%eax
test   $0x1,%r14b
je     1301 <linreg+0x1f1>
movsd  (%rdi,%rax,8),%xmm6
subsd  %xmm5,%xmm6
movapd %xmm6,%xmm5
mulsd  %xmm6,%xmm5
addsd  %xmm5,%xmm4
mulsd  (%rsi,%rax,8),%xmm6
addsd  %xmm6,%xmm3
movsd  %xmm3,(%r8)
divsd  %xmm4,%xmm3
movsd  %xmm3,(%r8)
movapd %xmm1,%xmm4
unpckhpd %xmm1,%xmm4
mulsd  %xmm3,%xmm4
subsd  %xmm4,%xmm1
divsd  %xmm2,%xmm1
movsd  %xmm1,(%rcx)
xorpd  %xmm1,%xmm1
test   %edx,%edx
jle    1348 <linreg+0x238>
cmp    $0x1,%edx
jne    1355 <linreg+0x245>
xorpd  %xmm2,%xmm2
xorpd  %xmm5,%xmm5
xor    %eax,%eax
test   $0x1,%r14b
jne    13d5 <linreg+0x2c5>
jmp    13fb <linreg+0x2eb>
xorpd  %xmm6,%xmm6
xorpd  %xmm5,%xmm5
jmp    13fb <linreg+0x2eb>
mov    %r14d,%ecx
and    $0xfffffffe,%ecx
xorpd  %xmm2,%xmm2
xorpd  %xmm5,%xmm5
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rsi,%rax,8),%xmm4
movhpd (%rdi,%rax,8),%xmm4
subpd  %xmm0,%xmm4
movapd %xmm4,%xmm3
unpckhpd %xmm4,%xmm3
mulpd  %xmm4,%xmm3
mulsd  %xmm4,%xmm4
addsd  %xmm5,%xmm4
movsd  0x8(%rsi,%rax,8),%xmm5
addpd  %xmm2,%xmm3
movhpd 0x8(%rdi,%rax,8),%xmm5
subpd  %xmm0,%xmm5
movapd %xmm5,%xmm6
unpckhpd %xmm5,%xmm6
mulpd  %xmm5,%xmm6
movapd %xmm6,%xmm2
addpd  %xmm3,%xmm2
mulsd  %xmm5,%xmm5
addsd  %xmm4,%xmm5
add    $0x2,%rax
cmp    %rax,%rcx
jne    1370 <linreg+0x260>
addpd  %xmm3,%xmm6
test   $0x1,%r14b
je     13fb <linreg+0x2eb>
movsd  (%rsi,%rax,8),%xmm3
movhpd (%rdi,%rax,8),%xmm3
subpd  %xmm0,%xmm3
movapd %xmm3,%xmm6
unpckhpd %xmm3,%xmm6
mulpd  %xmm3,%xmm6
addpd  %xmm2,%xmm6
mulsd  %xmm3,%xmm3
addsd  %xmm3,%xmm5
movapd %xmm6,%xmm0
unpckhpd %xmm6,%xmm0
ucomisd %xmm1,%xmm0
movapd %xmm6,0x10(%rsp)
jb     1418 <linreg+0x308>
xorps  %xmm2,%xmm2
sqrtsd %xmm0,%xmm2
jmp    142b <linreg+0x31b>
movapd %xmm5,(%rsp)
call   1030 <sqrt@plt>
movapd (%rsp),%xmm5
movapd %xmm0,%xmm2
mov    0x50(%rsp),%rbx
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm5
jb     1443 <linreg+0x333>
xorps  %xmm0,%xmm0
sqrtsd %xmm5,%xmm0
jmp    1456 <linreg+0x346>
movapd %xmm5,%xmm0
movsd  %xmm2,(%rsp)
call   1030 <sqrt@plt>
movsd  (%rsp),%xmm2
movapd 0x10(%rsp),%xmm1
mulsd  %xmm0,%xmm2
divsd  %xmm2,%xmm1
movsd  %xmm1,(%rbx)
add    $0x20,%rsp
pop    %rbx
pop    %r14
pop    %r15
ret

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

LLVM IR (-O2):
define void @linreg(double* nocapture %0, double* nocapture %1, i32 %2, double* nocapture writeonly %3, double* nocapture %4, double* nocapture writeonly %5, double* nocapture writeonly %6, double* nocapture writeonly %7, double* nocapture writeonly %8) {
%10 = load double, double* %1
store double 0.000000e+00, double* %4
%11 = icmp sgt i32 %2, 0
br i1 %11, label %18, label %12
12: ; preds = %9
%13 = sitofp i32 %2 to double
%14 = insertelement <2 x double> poison, double %13, i64 0
%15 = fdiv <2 x double> <double 0.000000e+00, double poison>, %14
%16 = shufflevector <2 x double> %15, <2 x double> poison, <2 x i32> zeroinitializer
%17 = extractelement <2 x double> %15, i64 0
store double %17, double* %7
store double %10, double* %5
store double 0.000000e+00, double* %6
br label %87
18: ; preds = %9
%19 = zext i32 %2 to i64
%20 = add nsw i64 %19, -1
%21 = and i64 %19, 1
%22 = icmp eq i64 %20, 0
br i1 %22, label %57, label %23
23: ; preds = %18
%24 = and i64 %19, 4294967294
br label %25
25: ; preds = %25, %23
%26 = phi i64 [ 0, %23 ], [ %54, %25 ]
%27 = phi double [ %10, %23 ], [ %51, %25 ]
%28 = phi double [ 0.000000e+00, %23 ], [ %53, %25 ]
%29 = phi <2 x double> [ zeroinitializer, %23 ], [ %49, %25 ]
%30 = phi i64 [ 0, %23 ], [ %55, %25 ]
%31 = getelementptr inbounds double, double* %0, i64 %26
%32 = load double, double* %31
%33 = getelementptr inbounds double, double* %1, i64 %26
%34 = load double, double* %33
%35 = insertelement <2 x double> poison, double %34, i64 0
%36 = insertelement <2 x double> %35, double %32, i64 1
%37 = fadd <2 x double> %29, %36
%38 = fcmp olt double %34, %27
%39 = select i1 %38, double %34, double %27
%40 = fcmp ogt double %34, %28
%41 = select i1 %40, double %34, double %28
%42 = or i64 %26, 1
%43 = getelementptr inbounds double, double* %0, i64 %42
%44 = load double, double* %43
%45 = getelementptr inbounds double, double* %1, i64 %42
%46 = load double, double* %45
%47 = insertelement <2 x double> poison, double %46, i64 0
%48 = insertelement <2 x double> %47, double %44, i64 1
%49 = fadd <2 x double> %37, %48
%50 = fcmp olt double %46, %39
%51 = select i1 %50, double %46, double %39
%52 = fcmp ogt double %46, %41
%53 = select i1 %52, double %46, double %41
%54 = add nuw nsw i64 %26, 2
%55 = add i64 %30, 2
%56 = icmp eq i64 %55, %24
br i1 %56, label %57, label %25
57: ; preds = %25, %18
%58 = phi <2 x double> [ undef, %18 ], [ %49, %25 ]
%59 = phi double [ undef, %18 ], [ %51, %25 ]
%60 = phi double [ undef, %18 ], [ %53, %25 ]
%61 = phi i64 [ 0, %18 ], [ %54, %25 ]
%62 = phi double [ %10, %18 ], [ %51, %25 ]
%63 = phi double [ 0.000000e+00, %18 ], [ %53, %25 ]
%64 = phi <2 x double> [ zeroinitializer, %18 ], [ %49, %25 ]
%65 = icmp eq i64 %21, 0
br i1 %65, label %78, label %66
66: ; preds = %57
%67 = getelementptr inbounds double, double* %0, i64 %61
%68 = load double, double* %67
%69 = getelementptr inbounds double, double* %1, i64 %61
%70 = load double, double* %69
%71 = insertelement <2 x double> poison, double %70, i64 0
%72 = insertelement <2 x double> %71, double %68, i64 1
%73 = fadd <2 x double> %64, %72
%74 = fcmp olt double %70, %62
%75 = select i1 %74, double %70, double %62
%76 = fcmp ogt double %70, %63
%77 = select i1 %76, double %70, double %63
br label %78
78: ; preds = %57, %66
%79 = phi <2 x double> [ %58, %57 ], [ %73, %66 ]
%80 = phi double [ %59, %57 ], [ %75, %66 ]
%81 = phi double [ %60, %57 ], [ %77, %66 ]
%82 = sitofp i32 %2 to double
%83 = insertelement <2 x double> poison, double %82, i64 0
%84 = shufflevector <2 x double> %83, <2 x double> poison, <2 x i32> zeroinitializer
%85 = fdiv <2 x double> %79, %84
%86 = extractelement <2 x double> %85, i64 0
store double %86, double* %7
store double %80, double* %5
store double %81, double* %6
br i1 %11, label %92, label %87
87: ; preds = %12, %78
%88 = phi double [ %13, %12 ], [ %82, %78 ]
%89 = phi <2 x double> [ %16, %12 ], [ %85, %78 ]
%90 = phi <2 x double> [ zeroinitializer, %12 ], [ %79, %78 ]
%91 = load double, double* %4
br label %139
92: ; preds = %78
%93 = load double, double* %4
%94 = and i64 %19, 1
%95 = icmp eq i64 %20, 0
br i1 %95, label %123, label %96
96: ; preds = %92
%97 = and i64 %19, 4294967294
%98 = extractelement <2 x double> %85, i64 1
%99 = extractelement <2 x double> %85, i64 1
br label %100
100: ; preds = %100, %96
%101 = phi double [ %93, %96 ], [ %119, %100 ]
%102 = phi i64 [ 0, %96 ], [ %120, %100 ]
%103 = phi double [ 0.000000e+00, %96 ], [ %116, %100 ]
%104 = phi i64 [ 0, %96 ], [ %121, %100 ]
%105 = getelementptr inbounds double, double* %0, i64 %102
%106 = load double, double* %105
%107 = fsub double %106, %98
%108 = tail call double @llvm.fmuladd.f64(double %107, double %107, double %103)
%109 = getelementptr inbounds double, double* %1, i64 %102
%110 = load double, double* %109
%111 = tail call double @llvm.fmuladd.f64(double %107, double %110, double %101)
store double %111, double* %4
%112 = or i64 %102, 1
%113 = getelementptr inbounds double, double* %0, i64 %112
%114 = load double, double* %113
%115 = fsub double %114, %99
%116 = tail call double @llvm.fmuladd.f64(double %115, double %115, double %108)
%117 = getelementptr inbounds double, double* %1, i64 %112
%118 = load double, double* %117
%119 = tail call double @llvm.fmuladd.f64(double %115, double %118, double %111)
store double %119, double* %4
%120 = add nuw nsw i64 %102, 2
%121 = add i64 %104, 2
%122 = icmp eq i64 %121, %97
br i1 %122, label %123, label %100
123: ; preds = %100, %92
%124 = phi double [ undef, %92 ], [ %116, %100 ]
%125 = phi double [ undef, %92 ], [ %119, %100 ]
%126 = phi double [ %93, %92 ], [ %119, %100 ]
%127 = phi i64 [ 0, %92 ], [ %120, %100 ]
%128 = phi double [ 0.000000e+00, %92 ], [ %116, %100 ]
%129 = icmp eq i64 %94, 0
br i1 %129, label %139, label %130
130: ; preds = %123
%131 = getelementptr inbounds double, double* %0, i64 %127
%132 = load double, double* %131
%133 = extractelement <2 x double> %85, i64 1
%134 = fsub double %132, %133
%135 = tail call double @llvm.fmuladd.f64(double %134, double %134, double %128)
%136 = getelementptr inbounds double, double* %1, i64 %127
%137 = load double, double* %136
%138 = tail call double @llvm.fmuladd.f64(double %134, double %137, double %126)
store double %138, double* %4
br label %139
139: ; preds = %130, %123, %87
%140 = phi double [ %88, %87 ], [ %82, %123 ], [ %82, %130 ]
%141 = phi double [ %91, %87 ], [ %125, %123 ], [ %138, %130 ]
%142 = phi double [ 0.000000e+00, %87 ], [ %124, %123 ], [ %135, %130 ]
%143 = phi <2 x double> [ %89, %87 ], [ %85, %123 ], [ %85, %130 ]
%144 = phi <2 x double> [ %90, %87 ], [ %79, %123 ], [ %79, %130 ]
%145 = fdiv double %141, %142
store double %145, double* %4
%146 = extractelement <2 x double> %144, i64 1
%147 = fneg double %146
%148 = extractelement <2 x double> %144, i64 0
%149 = tail call double @llvm.fmuladd.f64(double %147, double %145, double %148)
%150 = fdiv double %149, %140
store double %150, double* %3
br i1 %11, label %151, label %209
151: ; preds = %139
%152 = zext i32 %2 to i64
%153 = and i64 %152, 1
%154 = icmp eq i32 %2, 1
br i1 %154, label %190, label %155
155: ; preds = %151
%156 = and i64 %152, 4294967294
br label %157
157: ; preds = %157, %155
%158 = phi i64 [ 0, %155 ], [ %185, %157 ]
%159 = phi double [ 0.000000e+00, %155 ], [ %184, %157 ]
%160 = phi <2 x double> [ zeroinitializer, %155 ], [ %182, %157 ]
%161 = phi i64 [ 0, %155 ], [ %186, %157 ]
%162 = getelementptr inbounds double, double* %0, i64 %158
%163 = load double, double* %162
%164 = getelementptr inbounds double, double* %1, i64 %158
%165 = load double, double* %164
%166 = insertelement <2 x double> poison, double %165, i64 0
%167 = insertelement <2 x double> %166, double %163, i64 1
%168 = fsub <2 x double> %167, %143
%169 = shufflevector <2 x double> %168, <2 x double> undef, <2 x i32> <i32 1, i32 1>
%170 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %169, <2 x double> %168, <2 x double> %160)
%171 = extractelement <2 x double> %168, i64 0
%172 = tail call double @llvm.fmuladd.f64(double %171, double %171, double %159)
%173 = or i64 %158, 1
%174 = getelementptr inbounds double, double* %0, i64 %173
%175 = load double, double* %174
%176 = getelementptr inbounds double, double* %1, i64 %173
%177 = load double, double* %176
%178 = insertelement <2 x double> poison, double %177, i64 0
%179 = insertelement <2 x double> %178, double %175, i64 1
%180 = fsub <2 x double> %179, %143
%181 = shufflevector <2 x double> %180, <2 x double> undef, <2 x i32> <i32 1, i32 1>
%182 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %181, <2 x double> %180, <2 x double> %170)
%183 = extractelement <2 x double> %180, i64 0
%184 = tail call double @llvm.fmuladd.f64(double %183, double %183, double %172)
%185 = add nuw nsw i64 %158, 2
%186 = add i64 %161, 2
%187 = icmp eq i64 %186, %156
br i1 %187, label %188, label %157
188: ; preds = %157
%189 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %181, <2 x double> %180, <2 x double> %170)
br label %190
190: ; preds = %188, %151
%191 = phi <2 x double> [ undef, %151 ], [ %189, %188 ]
%192 = phi double [ undef, %151 ], [ %184, %188 ]
%193 = phi i64 [ 0, %151 ], [ %185, %188 ]
%194 = phi double [ 0.000000e+00, %151 ], [ %184, %188 ]
%195 = phi <2 x double> [ zeroinitializer, %151 ], [ %182, %188 ]
%196 = icmp eq i64 %153, 0
br i1 %196, label %209, label %197
197: ; preds = %190
%198 = getelementptr inbounds double, double* %0, i64 %193
%199 = load double, double* %198
%200 = getelementptr inbounds double, double* %1, i64 %193
%201 = load double, double* %200
%202 = insertelement <2 x double> poison, double %201, i64 0
%203 = insertelement <2 x double> %202, double %199, i64 1
%204 = fsub <2 x double> %203, %143
%205 = shufflevector <2 x double> %204, <2 x double> undef, <2 x i32> <i32 1, i32 1>
%206 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %205, <2 x double> %204, <2 x double> %195)
%207 = extractelement <2 x double> %204, i64 0
%208 = tail call double @llvm.fmuladd.f64(double %207, double %207, double %194)
br label %209
209: ; preds = %197, %190, %139
%210 = phi double [ 0.000000e+00, %139 ], [ %192, %190 ], [ %208, %197 ]
%211 = phi <2 x double> [ zeroinitializer, %139 ], [ %191, %190 ], [ %206, %197 ]
%212 = extractelement <2 x double> %211, i64 1
%213 = tail call double @sqrt(double %212)
%214 = tail call double @sqrt(double %210)
%215 = fmul double %213, %214
%216 = extractelement <2 x double> %211, i64 0
%217 = fdiv double %216, %215
store double %217, double* %8
ret void
}