#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 *dvector(int);
int /*<<< orphan*/ free(double *);
int /*<<< orphan*/ gaussj(double **, int, double *, int);
void bayesreg(double **xpx, double *xpy, double *bp, double **priormat,
              double *bpost, double **vpost, int p) {
  int j, k;
  double *bpb;
  bpb = dvector(p);
  for (j = 0; j < p; j++) {
    bpost[j] = 0.0;
    for (k = 0; k < p; k++) {
      vpost[j][k] = xpx[j][k] + priormat[j][k];
    }
  }
  for (j = 0; j < p; j++) {
    bpb[j] = 0.0;
    for (k = 0; k < p; k++) {
      bpb[j] += priormat[j][k] * bp[k];
    }
    bpost[j] = xpy[j] + bpb[j];
  }
  gaussj(vpost, p, bpost, 1);
  free(bpb);
  return;
}


x86 assembly (clang, -O0):
<bayesreg>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    0x10(%rbp),%edi
call   1050 <dvector@plt>
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    11f5 <bayesreg+0xc5>
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    11e2 <bayesreg+0xb2>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1189 <bayesreg+0x59>
jmp    11e7 <bayesreg+0xb7>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1166 <bayesreg+0x36>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    12aa <bayesreg+0x17a>
mov    -0x40(%rbp),%rax
movslq -0x34(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    1275 <bayesreg+0x145>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x40(%rbp),%rax
movslq -0x34(%rbp),%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    121f <bayesreg+0xef>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x34(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11fc <bayesreg+0xcc>
mov    -0x30(%rbp),%rdi
mov    0x10(%rbp),%esi
mov    -0x28(%rbp),%rdx
mov    $0x1,%ecx
call   1040 <gaussj@plt>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<bayesreg>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r9,%r14
mov    %r8,%r15
mov    %rcx,%r12
mov    %rdx,%rbp
mov    %rsi,%rbx
mov    %rdi,0x40(%rsp)
mov    0x80(%rsp),%r13d
mov    %r13d,%edi
call   1050 <dvector@plt>
mov    %rbx,%rcx
mov    %rax,%rbx
test   %r13d,%r13d
jle    149c <bayesreg+0x36c>
mov    %rcx,0x18(%rsp)
mov    %r13d,%eax
mov    %eax,%ecx
and    $0xfffffffc,%ecx
mov    %rcx,0x38(%rsp)
add    $0xfffffffffffffffc,%rcx
mov    %rcx,0x30(%rsp)
mov    %rcx,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
lea    -0x1(%rax),%r10
mov    %rdx,%rcx
mov    %rdx,0x28(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,0x20(%rsp)
mov    %eax,%esi
and    $0x3,%esi
xor    %edi,%edi
mov    %rsi,(%rsp)
mov    %r12,0x10(%rsp)
mov    %r10,0x8(%rsp)
jmp    11dd <bayesreg+0xad>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rdi
cmp    %rax,%rdi
je     13ff <bayesreg+0x2cf>
movq   $0x0,(%r15,%rdi,8)
mov    0x40(%rsp),%rcx
mov    (%rcx,%rdi,8),%rdx
mov    (%r12,%rdi,8),%r8
mov    (%r14,%rdi,8),%r9
cmp    $0x4,%r13d
jae    1210 <bayesreg+0xe0>
xor    %r11d,%r11d
jmp    1360 <bayesreg+0x230>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %r15,%r12
mov    %r13d,%r15d
mov    %r14,%r13
lea    (%r9,%rax,8),%rsi
lea    (%rdx,%rax,8),%rcx
lea    (%r8,%rax,8),%r11
cmp    %rcx,%r9
setb   %r14b
cmp    %rsi,%rdx
setb   %r10b
cmp    %r11,%r9
setb   %cl
cmp    %rsi,%r8
setb   %sil
test   %r10b,%r14b
jne    12f9 <bayesreg+0x1c9>
mov    $0x0,%r11d
and    %sil,%cl
mov    %r13,%r14
mov    %r15d,%r13d
mov    %r12,%r15
mov    0x10(%rsp),%r12
mov    0x8(%rsp),%r10
mov    (%rsp),%rsi
jne    1360 <bayesreg+0x230>
cmpq   $0x0,0x30(%rsp)
je     1315 <bayesreg+0x1e5>
mov    0x20(%rsp),%rsi
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd (%rdx,%rcx,8),%xmm0
movupd 0x10(%rdx,%rcx,8),%xmm1
movupd (%r8,%rcx,8),%xmm2
addpd  %xmm0,%xmm2
movupd 0x10(%r8,%rcx,8),%xmm0
addpd  %xmm1,%xmm0
movupd %xmm2,(%r9,%rcx,8)
movupd %xmm0,0x10(%r9,%rcx,8)
movupd 0x20(%rdx,%rcx,8),%xmm0
movupd 0x30(%rdx,%rcx,8),%xmm1
movupd 0x20(%r8,%rcx,8),%xmm2
addpd  %xmm0,%xmm2
movupd 0x30(%r8,%rcx,8),%xmm0
addpd  %xmm1,%xmm0
movupd %xmm2,0x20(%r9,%rcx,8)
movupd %xmm0,0x30(%r9,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rsi
jne    1290 <bayesreg+0x160>
jmp    1317 <bayesreg+0x1e7>
xor    %r11d,%r11d
mov    %r13,%r14
mov    %r15d,%r13d
mov    %r12,%r15
mov    0x10(%rsp),%r12
mov    0x8(%rsp),%r10
mov    (%rsp),%rsi
jmp    1360 <bayesreg+0x230>
xor    %ecx,%ecx
testb  $0x1,0x28(%rsp)
mov    (%rsp),%rsi
je     134f <bayesreg+0x21f>
movupd (%rdx,%rcx,8),%xmm0
movupd 0x10(%rdx,%rcx,8),%xmm1
movupd (%r8,%rcx,8),%xmm2
addpd  %xmm0,%xmm2
movupd 0x10(%r8,%rcx,8),%xmm0
addpd  %xmm1,%xmm0
movupd %xmm2,(%r9,%rcx,8)
movupd %xmm0,0x10(%r9,%rcx,8)
mov    0x38(%rsp),%rcx
mov    %rcx,%r11
cmp    %rax,%rcx
je     11d0 <bayesreg+0xa0>
mov    %r10,%rcx
sub    %r11,%rcx
test   %rsi,%rsi
je     138c <bayesreg+0x25c>
nopl   0x0(%rax,%rax,1)
movsd  (%rdx,%r11,8),%xmm0
addsd  (%r8,%r11,8),%xmm0
movsd  %xmm0,(%r9,%r11,8)
add    $0x1,%r11
add    $0xffffffffffffffff,%rsi
jne    1370 <bayesreg+0x240>
cmp    $0x3,%rcx
mov    (%rsp),%rsi
jb     11d0 <bayesreg+0xa0>
nopw   0x0(%rax,%rax,1)
movsd  (%rdx,%r11,8),%xmm0
addsd  (%r8,%r11,8),%xmm0
movsd  %xmm0,(%r9,%r11,8)
movsd  0x8(%rdx,%r11,8),%xmm0
addsd  0x8(%r8,%r11,8),%xmm0
movsd  %xmm0,0x8(%r9,%r11,8)
movsd  0x10(%rdx,%r11,8),%xmm0
addsd  0x10(%r8,%r11,8),%xmm0
movsd  %xmm0,0x10(%r9,%r11,8)
movsd  0x18(%rdx,%r11,8),%xmm0
addsd  0x18(%r8,%r11,8),%xmm0
movsd  %xmm0,0x18(%r9,%r11,8)
add    $0x4,%r11
cmp    %r11,%rax
jne    13a0 <bayesreg+0x270>
jmp    11d0 <bayesreg+0xa0>
test   %r13d,%r13d
mov    0x18(%rsp),%r8
jle    149c <bayesreg+0x36c>
mov    %eax,%edx
and    $0xfffffffe,%edx
xor    %esi,%esi
jmp    1435 <bayesreg+0x305>
cs nopw 0x0(%rax,%rax,1)
addsd  (%r8,%rsi,8),%xmm0
movsd  %xmm0,(%r15,%rsi,8)
add    $0x1,%rsi
cmp    %rax,%rsi
je     149c <bayesreg+0x36c>
movq   $0x0,(%rbx,%rsi,8)
mov    (%r12,%rsi,8),%rcx
xorpd  %xmm0,%xmm0
xor    %edi,%edi
test   %r10,%r10
je     1482 <bayesreg+0x352>
nopl   0x0(%rax)
movsd  (%rcx,%rdi,8),%xmm1
mulsd  0x0(%rbp,%rdi,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rsi,8)
movsd  0x8(%rcx,%rdi,8),%xmm0
mulsd  0x8(%rbp,%rdi,8),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rbx,%rsi,8)
add    $0x2,%rdi
cmp    %rdi,%rdx
jne    1450 <bayesreg+0x320>
test   $0x1,%al
je     1420 <bayesreg+0x2f0>
movsd  (%rcx,%rdi,8),%xmm1
mulsd  0x0(%rbp,%rdi,8),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rbx,%rsi,8)
jmp    1420 <bayesreg+0x2f0>
mov    %r14,%rdi
mov    %r13d,%esi
mov    %r15,%rdx
mov    $0x1,%ecx
call   1040 <gaussj@plt>
mov    %rbx,%rdi
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

LLVM IR (-O0):
define void @bayesreg(double** %0, double* %1, double* %2, double** %3, double* %4, double** %5, i32 %6) {
%8 = alloca double**
%9 = alloca double*
%10 = alloca double*
%11 = alloca double**
%12 = alloca double*
%13 = alloca double**
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca double*
store double** %0, double*** %8
store double* %1, double** %9
store double* %2, double** %10
store double** %3, double*** %11
store double* %4, double** %12
store double** %5, double*** %13
store i32 %6, i32* %14
%18 = load i32, i32* %14
%19 = call double* @dvector(i32 %18)
store double* %19, double** %17
store i32 0, i32* %15
br label %20
20: ; preds = %65, %7
%21 = load i32, i32* %15
%22 = load i32, i32* %14
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %68
24: ; preds = %20
%25 = load double*, double** %12
%26 = load i32, i32* %15
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
store double 0.000000e+00, double* %28
store i32 0, i32* %16
br label %29
29: ; preds = %61, %24
%30 = load i32, i32* %16
%31 = load i32, i32* %14
%32 = icmp slt i32 %30, %31
br i1 %32, label %33, label %64
33: ; preds = %29
%34 = load double**, double*** %8
%35 = load i32, i32* %15
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double*, double** %34, i64 %36
%38 = load double*, double** %37
%39 = load i32, i32* %16
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %38, i64 %40
%42 = load double, double* %41
%43 = load double**, double*** %11
%44 = load i32, i32* %15
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds double*, double** %43, i64 %45
%47 = load double*, double** %46
%48 = load i32, i32* %16
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %47, i64 %49
%51 = load double, double* %50
%52 = fadd double %42, %51
%53 = load double**, double*** %13
%54 = load i32, i32* %15
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double*, double** %53, i64 %55
%57 = load double*, double** %56
%58 = load i32, i32* %16
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double, double* %57, i64 %59
store double %52, double* %60
br label %61
61: ; preds = %33
%62 = load i32, i32* %16
%63 = add nsw i32 %62, 1
store i32 %63, i32* %16
br label %29
64: ; preds = %29
br label %65
65: ; preds = %64
%66 = load i32, i32* %15
%67 = add nsw i32 %66, 1
store i32 %67, i32* %15
br label %20
68: ; preds = %20
store i32 0, i32* %15
br label %69
69: ; preds = %122, %68
%70 = load i32, i32* %15
%71 = load i32, i32* %14
%72 = icmp slt i32 %70, %71
br i1 %72, label %73, label %125
73: ; preds = %69
%74 = load double*, double** %17
%75 = load i32, i32* %15
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %74, i64 %76
store double 0.000000e+00, double* %77
store i32 0, i32* %16
br label %78
78: ; preds = %103, %73
%79 = load i32, i32* %16
%80 = load i32, i32* %14
%81 = icmp slt i32 %79, %80
br i1 %81, label %82, label %106
82: ; preds = %78
%83 = load double**, double*** %11
%84 = load i32, i32* %15
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double*, double** %83, i64 %85
%87 = load double*, double** %86
%88 = load i32, i32* %16
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds double, double* %87, i64 %89
%91 = load double, double* %90
%92 = load double*, double** %10
%93 = load i32, i32* %16
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds double, double* %92, i64 %94
%96 = load double, double* %95
%97 = load double*, double** %17
%98 = load i32, i32* %15
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds double, double* %97, i64 %99
%101 = load double, double* %100
%102 = call double @llvm.fmuladd.f64(double %91, double %96, double %101)
store double %102, double* %100
br label %103
103: ; preds = %82
%104 = load i32, i32* %16
%105 = add nsw i32 %104, 1
store i32 %105, i32* %16
br label %78
106: ; preds = %78
%107 = load double*, double** %9
%108 = load i32, i32* %15
%109 = sext i32 %108 to i64
%110 = getelementptr inbounds double, double* %107, i64 %109
%111 = load double, double* %110
%112 = load double*, double** %17
%113 = load i32, i32* %15
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds double, double* %112, i64 %114
%116 = load double, double* %115
%117 = fadd double %111, %116
%118 = load double*, double** %12
%119 = load i32, i32* %15
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds double, double* %118, i64 %120
store double %117, double* %121
br label %122
122: ; preds = %106
%123 = load i32, i32* %15
%124 = add nsw i32 %123, 1
store i32 %124, i32* %15
br label %69
125: ; preds = %69
%126 = load double**, double*** %13
%127 = load i32, i32* %14
%128 = load double*, double** %12
%129 = call i32 @gaussj(double** %126, i32 %127, double* %128, i32 1)
%130 = load double*, double** %17
%131 = call i32 @free(double* %130)
ret void
}

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