#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 */
int /*<<< orphan*/ fast_inverse(double *, int, double *);
void solveRegression(double *AtA, double *AtB, int *selectK, int *selectN,
                     int sK, int N, int K, double *Z, double *sAtA,
                     double *inVsAtA) {
  int k1, k2, sk1, sk2, i;
  sk1 = 0;
  for (k1 = 0; k1 < K; k1++) {
    if (selectK[k1]) {
      sk2 = 0;
      for (k2 = 0; k2 < K; k2++) {
        if (selectK[k2]) {
          sAtA[sk1 * sK + sk2] = AtA[k1 * K + k2];
          sk2++;
        }
      }
      sk1++;
    }
  }
  if (sK > 1)
    fast_inverse(sAtA, sK, inVsAtA);
  else
    inVsAtA[0] = 1 / sAtA[0];
  sk1 = 0;
  for (k1 = 0; k1 < K; k1++) {
    if (selectK[k1]) {
      for (i = 0; i < N; i++) {
        if (selectN[i]) {
          Z[k1 * N + i] = 0;
          sk2 = 0;
          for (k2 = 0; k2 < K; k2++) {
            if (selectK[k2]) {
              Z[k1 * N + i] += inVsAtA[sk1 * sK + sk2] * AtB[k2 * N + i];
              sk2++;
            }
          }
        }
      }
      sk1++;
    } else {
      for (i = 0; i < N; i++) {
        if (selectN[i]) {
          Z[k1 * N + i] = 0;
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<solveRegression>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9d,-0x28(%rbp)
movl   $0x0,-0x34(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    11fb <solveRegression+0xeb>
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     11e8 <solveRegression+0xd8>
movl   $0x0,-0x38(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    11df <solveRegression+0xcf>
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     11cc <solveRegression+0xbc>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   0x10(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x24(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11d1 <solveRegression+0xc1>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1179 <solveRegression+0x69>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11ed <solveRegression+0xdd>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    114d <solveRegression+0x3d>
cmpl   $0x1,-0x24(%rbp)
jle    121a <solveRegression+0x10a>
mov    0x20(%rbp),%rdi
mov    -0x24(%rbp),%esi
mov    0x28(%rbp),%rdx
call   1030 <fast_inverse@plt>
jmp    1232 <solveRegression+0x122>
mov    0x20(%rbp),%rax
movsd  0xdda(%rip),%xmm0
divsd  (%rax),%xmm0
mov    0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x0,-0x34(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    13c2 <solveRegression+0x2b2>
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1359 <solveRegression+0x249>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    134b <solveRegression+0x23b>
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1338 <solveRegression+0x228>
mov    0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x28(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x38(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    1333 <solveRegression+0x223>
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1320 <solveRegression+0x210>
mov    0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x24(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x28(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x28(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1325 <solveRegression+0x215>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12aa <solveRegression+0x19a>
jmp    1338 <solveRegression+0x228>
jmp    133d <solveRegression+0x22d>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1265 <solveRegression+0x155>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    13af <solveRegression+0x29f>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    13aa <solveRegression+0x29a>
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1397 <solveRegression+0x287>
mov    0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x28(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    139c <solveRegression+0x28c>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1360 <solveRegression+0x250>
jmp    13af <solveRegression+0x29f>
jmp    13b4 <solveRegression+0x2a4>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1240 <solveRegression+0x130>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<solveRegression>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r9d,%r15d
mov    %rcx,0x10(%rsp)
mov    %rdx,%rbx
mov    %rsi,0x30(%rsp)
mov    %rdi,(%rsp)
mov    0x90(%rsp),%rdi
mov    0x80(%rsp),%eax
mov    %eax,%r12d
mov    %eax,%r13d
test   %eax,%eax
mov    %r8d,0xc(%rsp)
jle    121f <solveRegression+0x10f>
mov    %r12d,%r14d
and    $0xfffffffe,%r14d
mov    (%rsp),%rax
lea    0x8(%rax),%rdx
lea    0x0(,%r12,8),%r9
xor    %r11d,%r11d
xor    %r10d,%r10d
jmp    1199 <solveRegression+0x89>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r10d
mov    0xc(%rsp),%r8d
add    $0x1,%r11
add    %r9,%rdx
cmp    %r12,%r11
je     121f <solveRegression+0x10f>
cmpl   $0x0,(%rbx,%r11,4)
je     1189 <solveRegression+0x79>
mov    %r10d,%esi
imul   %r8d,%esi
xor    %eax,%eax
xor    %ebp,%ebp
cmp    $0x1,%r13d
jne    11e9 <solveRegression+0xd9>
test   $0x1,%r12b
je     1180 <solveRegression+0x70>
cmpl   $0x0,(%rbx,%rax,4)
je     1180 <solveRegression+0x70>
mov    %r11,%rcx
imul   %r12,%rcx
add    %rcx,%rax
mov    (%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm0
add    %esi,%ebp
movslq %ebp,%rax
movsd  %xmm0,(%rdi,%rax,8)
jmp    1180 <solveRegression+0x70>
nopl   0x0(%rax)
add    $0x2,%rax
cmp    %rax,%r14
je     11b1 <solveRegression+0xa1>
cmpl   $0x0,(%rbx,%rax,4)
je     1203 <solveRegression+0xf3>
movsd  -0x8(%rdx,%rax,8),%xmm0
lea    (%rsi,%rbp,1),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rdi,%rcx,8)
add    $0x1,%ebp
cmpl   $0x0,0x4(%rbx,%rax,4)
je     11e0 <solveRegression+0xd0>
movsd  (%rdx,%rax,8),%xmm0
lea    (%rsi,%rbp,1),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rdi,%rcx,8)
add    $0x1,%ebp
jmp    11e0 <solveRegression+0xd0>
mov    0x98(%rsp),%rbp
cmp    $0x2,%r8d
jl     1242 <solveRegression+0x132>
mov    %r8d,%esi
mov    %rbp,%rdx
call   1030 <fast_inverse@plt>
test   %r13d,%r13d
jg     125c <solveRegression+0x14c>
jmp    1413 <solveRegression+0x303>
movsd  0xdb6(%rip),%xmm0
divsd  (%rdi),%xmm0
movsd  %xmm0,0x0(%rbp)
test   %r13d,%r13d
jle    1413 <solveRegression+0x303>
mov    0x88(%rsp),%rax
movslq %r15d,%rcx
mov    %ecx,%edx
mov    %rdx,(%rsp)
mov    %edx,%r8d
and    $0xfffffffe,%r8d
mov    %rcx,0x20(%rsp)
lea    0x0(,%rcx,8),%rdi
lea    0x8(%rax),%rsi
xor    %r9d,%r9d
movl   $0x0,0x8(%rsp)
mov    %r15d,0x1c(%rsp)
mov    %r8,0x38(%rsp)
jmp    12c9 <solveRegression+0x1b9>
nopl   0x0(%rax)
addl   $0x1,0x8(%rsp)
mov    0x1c(%rsp),%r15d
mov    0x38(%rsp),%r8
mov    0x40(%rsp),%rsi
mov    0x28(%rsp),%r9
add    $0x1,%r9
add    %rdi,%rsi
cmp    %r12,%r9
je     1413 <solveRegression+0x303>
cmpl   $0x0,(%rbx,%r9,4)
je     1390 <solveRegression+0x280>
mov    %r9,0x28(%rsp)
mov    %rsi,0x40(%rsp)
test   %r15d,%r15d
jle    12a0 <solveRegression+0x190>
mov    0x28(%rsp),%rsi
imul   0x20(%rsp),%rsi
mov    0x8(%rsp),%ecx
mov    %ecx,%r8d
imul   0xc(%rsp),%r8d
mov    0x30(%rsp),%r11
xor    %r15d,%r15d
jmp    131e <solveRegression+0x20e>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r15
add    $0x8,%r11
cmp    (%rsp),%r15
je     12a0 <solveRegression+0x190>
mov    0x10(%rsp),%rcx
cmpl   $0x0,(%rcx,%r15,4)
je     1310 <solveRegression+0x200>
lea    (%r15,%rsi,1),%r10
movq   $0x0,(%rax,%r10,8)
xor    %r9d,%r9d
xorpd  %xmm0,%xmm0
mov    %r11,%r14
xor    %ecx,%ecx
jmp    135c <solveRegression+0x24c>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%r9
add    %rdi,%r14
cmp    %r9,%r12
je     1310 <solveRegression+0x200>
cmpl   $0x0,(%rbx,%r9,4)
je     1350 <solveRegression+0x240>
lea    (%rcx,%r8,1),%r13d
movslq %r13d,%rdx
movsd  0x0(%rbp,%rdx,8),%xmm1
mulsd  (%r14),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%r10,8)
add    $0x1,%ecx
jmp    1350 <solveRegression+0x240>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
test   %r15d,%r15d
jle    12b9 <solveRegression+0x1a9>
xor    %ecx,%ecx
cmp    $0x1,%r15d
jne    13e9 <solveRegression+0x2d9>
testb  $0x1,(%rsp)
je     12b9 <solveRegression+0x1a9>
mov    0x10(%rsp),%rdx
cmpl   $0x0,(%rdx,%rcx,4)
je     12b9 <solveRegression+0x1a9>
mov    %r9,%rdx
imul   0x20(%rsp),%rdx
add    %rdx,%rcx
movq   $0x0,(%rax,%rcx,8)
jmp    12b9 <solveRegression+0x1a9>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x2,%rcx
cmp    %rcx,%r8
je     13a1 <solveRegression+0x291>
mov    0x10(%rsp),%rdx
cmpl   $0x0,(%rdx,%rcx,4)
je     13fd <solveRegression+0x2ed>
movq   $0x0,-0x8(%rsi,%rcx,8)
mov    0x10(%rsp),%rdx
cmpl   $0x0,0x4(%rdx,%rcx,4)
je     13e0 <solveRegression+0x2d0>
movq   $0x0,(%rsi,%rcx,8)
jmp    13e0 <solveRegression+0x2d0>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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