int linSystemSolve(int dim, int order, double *a, double *b, int *pivot) {
  int i, k, m;
  double temp;
  if (order == 1) {
    if (*(a + (0) * dim + (0)) == 0.0) {
      return 1;
    } else {
      b[0] /= *(a + (0) * dim + (0));
      return 0;
    }
  }
  for (k = 0; k < order - 1; ++k) {
    m = pivot[k];
    temp = b[m];
    b[m] = b[k];
    b[k] = temp;
    for (i = k + 1; i < order; ++i) {
      b[i] += *(a + (i)*dim + (k)) * temp;
    }
  }
  for (k = order - 1; k >= 0; --k) {
    if (*(a + (k)*dim + (k)) == 0) {
      return 1;
    }
    b[k] /= *(a + (k)*dim + (k));
    temp = -b[k];
    for (i = 0; i < k; ++i) {
      b[i] += *(a + (i)*dim + (k)) * temp;
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<linSystemSolve>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
cmpl   $0x1,-0xc(%rbp)
jne    117b <linSystemSolve+0x7b>
mov    -0x18(%rbp),%rax
imul   $0x0,-0x8(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    114f <linSystemSolve+0x4f>
jp     114f <linSystemSolve+0x4f>
movl   $0x1,-0x4(%rbp)
jmp    1377 <linSystemSolve+0x277>
mov    -0x18(%rbp),%rax
imul   $0x0,-0x8(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x4(%rbp)
jmp    1377 <linSystemSolve+0x277>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1252 <linSystemSolve+0x152>
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x34(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    123f <linSystemSolve+0x13f>
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x40(%rbp),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11e8 <linSystemSolve+0xe8>
jmp    1244 <linSystemSolve+0x144>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1182 <linSystemSolve+0x82>
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x0,-0x30(%rbp)
jl     1370 <linSystemSolve+0x270>
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x8(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    12a2 <linSystemSolve+0x1a2>
jp     12a2 <linSystemSolve+0x1a2>
movl   $0x1,-0x4(%rbp)
jmp    1377 <linSystemSolve+0x277>
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x8(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jge    135d <linSystemSolve+0x25d>
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x40(%rbp),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1306 <linSystemSolve+0x206>
jmp    1362 <linSystemSolve+0x262>
mov    -0x30(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x30(%rbp)
jmp    125b <linSystemSolve+0x15b>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<linSystemSolve>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
cmp    $0x1,%esi
jne    113b <linSystemSolve+0x3b>
movsd  (%rdx),%xmm0
mov    $0x1,%eax
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1128 <linSystemSolve+0x28>
jnp    1327 <linSystemSolve+0x227>
movsd  (%rcx),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rcx)
xor    %eax,%eax
jmp    1327 <linSystemSolve+0x227>
xor    %eax,%eax
cmp    $0x2,%esi
jl     1327 <linSystemSolve+0x227>
lea    -0x1(%rsi),%r10d
movslq %edi,%r11
mov    %esi,%r9d
mov    %rdx,-0x10(%rsp)
lea    (%rdx,%r11,8),%r15
lea    0x0(,%r11,8),%r13
lea    0x8(,%r11,8),%r14
xor    %r12d,%r12d
mov    %r9,%rdi
mov    %rcx,%rax
jmp    1190 <linSystemSolve+0x90>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    %r14,%r15
add    $0x8,%rax
add    $0xffffffffffffffff,%rdi
cmp    %r10,%r12
je     11e0 <linSystemSolve+0xe0>
movslq (%r8,%r12,4),%rbx
movsd  (%rcx,%rbx,8),%xmm0
movsd  (%rcx,%r12,8),%xmm1
movsd  %xmm1,(%rcx,%rbx,8)
movsd  %xmm0,(%rcx,%r12,8)
add    $0x1,%r12
cmp    %r9,%r12
jae    1180 <linSystemSolve+0x80>
mov    $0x1,%ebp
mov    %r15,%rbx
nopl   0x0(%rax,%rax,1)
movsd  (%rbx),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rax,%rbp,8),%xmm1
movsd  %xmm1,(%rax,%rbp,8)
add    %r13,%rbx
add    $0x1,%rbp
cmp    %rbp,%rdi
jne    11c0 <linSystemSolve+0xc0>
jmp    1180 <linSystemSolve+0x80>
test   %esi,%esi
jle    131e <linSystemSolve+0x21e>
lea    -0x2(%r9),%rsi
lea    -0x1(%r9),%r10
mov    -0x10(%rsp),%rdx
lea    (%rdx,%r9,8),%r8
add    $0xfffffffffffffff8,%r8
mov    %r11,%rdi
shl    $0x4,%rdi
xor    %r12d,%r12d
xorpd  %xmm0,%xmm0
movapd 0xded(%rip),%xmm1
mov    %r9,%r15
mov    %r9,-0x8(%rsp)
jmp    123e <linSystemSolve+0x13e>
nopl   (%rax)
add    $0x1,%r12
add    $0xfffffffffffffff8,%r8
add    $0xffffffffffffffff,%r10
cmp    $0x1,%r15
mov    %r13,%r15
mov    -0x8(%rsp),%r9
jbe    1327 <linSystemSolve+0x227>
lea    -0x1(%r15),%r13
mov    %r13,%rax
imul   %r11,%rax
lea    (%rdx,%rax,8),%rax
movsd  -0x8(%rax,%r15,8),%xmm3
ucomisd %xmm0,%xmm3
jne    1260 <linSystemSolve+0x160>
jnp    1322 <linSystemSolve+0x222>
movsd  (%rcx,%r13,8),%xmm2
divsd  %xmm3,%xmm2
movsd  %xmm2,(%rcx,%r13,8)
xor    %eax,%eax
cmp    $0x2,%r15
jb     1327 <linSystemSolve+0x227>
mov    %r12,%rbx
not    %rbx
add    %r9,%rbx
xorpd  %xmm1,%xmm2
cmp    %r12,%rsi
jne    12a0 <linSystemSolve+0x1a0>
xor    %r9d,%r9d
jmp    12e9 <linSystemSolve+0x1e9>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r10,%rbp
and    $0xfffffffffffffffe,%rbp
mov    %r8,%r14
xor    %r9d,%r9d
nopl   (%rax)
movsd  (%r14),%xmm3
mulsd  %xmm2,%xmm3
addsd  (%rcx,%r9,8),%xmm3
movsd  %xmm3,(%rcx,%r9,8)
movsd  (%r14,%r11,8),%xmm3
mulsd  %xmm2,%xmm3
addsd  0x8(%rcx,%r9,8),%xmm3
movsd  %xmm3,0x8(%rcx,%r9,8)
add    $0x2,%r9
add    %rdi,%r14
cmp    %r9,%rbp
jne    12b0 <linSystemSolve+0x1b0>
test   $0x1,%bl
mov    -0x10(%rsp),%rdx
je     1220 <linSystemSolve+0x120>
lea    (%rdx,%r13,8),%rbx
mov    %rsi,%rbp
mov    %r9,%rsi
imul   %r11,%rsi
mulsd  (%rbx,%rsi,8),%xmm2
mov    %rbp,%rsi
addsd  (%rcx,%r9,8),%xmm2
movsd  %xmm2,(%rcx,%r9,8)
jmp    1220 <linSystemSolve+0x120>
xor    %eax,%eax
jmp    1327 <linSystemSolve+0x227>
mov    $0x1,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @linSystemSolve(i32 %0, i32 %1, double* nocapture %2, double* nocapture %3, i32* nocapture %4) {
%6 = icmp eq i32 %1, 1
br i1 %6, label %15, label %7
7: ; preds = %5
%8 = icmp sgt i32 %1, 1
br i1 %8, label %9, label %108
9: ; preds = %7
%10 = add nsw i32 %1, -1
%11 = sext i32 %0 to i64
%12 = zext i32 %1 to i64
%13 = zext i32 %10 to i64
%14 = zext i32 %1 to i64
br label %30
15: ; preds = %5
%16 = load double, double* %2
%17 = fcmp oeq double %16, 0.000000e+00
br i1 %17, label %108, label %18
18: ; preds = %15
%19 = load double, double* %3
%20 = fdiv double %19, %16
store double %20, double* %3
br label %108
21: ; preds = %44, %30
%22 = add nuw nsw i64 %32, 1
%23 = icmp eq i64 %40, %13
br i1 %23, label %24, label %30
24: ; preds = %21
%25 = icmp sgt i32 %1, 0
br i1 %25, label %26, label %108
26: ; preds = %24
%27 = sext i32 %0 to i64
%28 = zext i32 %1 to i64
%29 = add nsw i64 %14, -2
br label %54
30: ; preds = %9, %21
%31 = phi i64 [ 0, %9 ], [ %40, %21 ]
%32 = phi i64 [ 1, %9 ], [ %22, %21 ]
%33 = getelementptr inbounds i32, i32* %4, i64 %31
%34 = load i32, i32* %33
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double, double* %3, i64 %35
%37 = load double, double* %36
%38 = getelementptr inbounds double, double* %3, i64 %31
%39 = load double, double* %38
store double %39, double* %36
store double %37, double* %38
%40 = add nuw nsw i64 %31, 1
%41 = icmp ult i64 %40, %12
br i1 %41, label %42, label %21
42: ; preds = %30
%43 = getelementptr inbounds double, double* %2, i64 %31
br label %44
44: ; preds = %42, %44
%45 = phi i64 [ %32, %42 ], [ %52, %44 ]
%46 = mul nsw i64 %45, %11
%47 = getelementptr inbounds double, double* %43, i64 %46
%48 = load double, double* %47
%49 = getelementptr inbounds double, double* %3, i64 %45
%50 = load double, double* %49
%51 = tail call double @llvm.fmuladd.f64(double %48, double %37, double %50)
store double %51, double* %49
%52 = add nuw nsw i64 %45, 1
%53 = icmp eq i64 %52, %14
br i1 %53, label %21, label %44
54: ; preds = %26, %106
%55 = phi i64 [ 0, %26 ], [ %107, %106 ]
%56 = phi i64 [ %28, %26 ], [ %59, %106 ]
%57 = xor i64 %55, -1
%58 = add i64 %57, %14
%59 = add nsw i64 %56, -1
%60 = mul nsw i64 %59, %27
%61 = getelementptr inbounds double, double* %2, i64 %60
%62 = getelementptr inbounds double, double* %61, i64 %59
%63 = load double, double* %62
%64 = fcmp oeq double %63, 0.000000e+00
br i1 %64, label %108, label %65
65: ; preds = %54
%66 = getelementptr inbounds double, double* %3, i64 %59
%67 = load double, double* %66
%68 = fdiv double %67, %63
store double %68, double* %66
%69 = fneg double %68
%70 = icmp ugt i64 %56, 1
br i1 %70, label %71, label %108
71: ; preds = %65
%72 = getelementptr inbounds double, double* %2, i64 %59
%73 = and i64 %58, 1
%74 = icmp eq i64 %29, %55
br i1 %74, label %96, label %75
75: ; preds = %71
%76 = and i64 %58, -2
br label %77
77: ; preds = %77, %75
%78 = phi i64 [ 0, %75 ], [ %93, %77 ]
%79 = phi i64 [ 0, %75 ], [ %94, %77 ]
%80 = mul nsw i64 %78, %27
%81 = getelementptr inbounds double, double* %72, i64 %80
%82 = load double, double* %81
%83 = getelementptr inbounds double, double* %3, i64 %78
%84 = load double, double* %83
%85 = tail call double @llvm.fmuladd.f64(double %82, double %69, double %84)
store double %85, double* %83
%86 = or i64 %78, 1
%87 = mul nsw i64 %86, %27
%88 = getelementptr inbounds double, double* %72, i64 %87
%89 = load double, double* %88
%90 = getelementptr inbounds double, double* %3, i64 %86
%91 = load double, double* %90
%92 = tail call double @llvm.fmuladd.f64(double %89, double %69, double %91)
store double %92, double* %90
%93 = add nuw nsw i64 %78, 2
%94 = add i64 %79, 2
%95 = icmp eq i64 %94, %76
br i1 %95, label %96, label %77
96: ; preds = %77, %71
%97 = phi i64 [ 0, %71 ], [ %93, %77 ]
%98 = icmp eq i64 %73, 0
br i1 %98, label %106, label %99
99: ; preds = %96
%100 = mul nsw i64 %97, %27
%101 = getelementptr inbounds double, double* %72, i64 %100
%102 = load double, double* %101
%103 = getelementptr inbounds double, double* %3, i64 %97
%104 = load double, double* %103
%105 = tail call double @llvm.fmuladd.f64(double %102, double %69, double %104)
store double %105, double* %103
br label %106
106: ; preds = %96, %99
%107 = add i64 %55, 1
br i1 %70, label %54, label %108
108: ; preds = %65, %54, %106, %7, %24, %15, %18
%109 = phi i32 [ 0, %18 ], [ 1, %15 ], [ 0, %24 ], [ 0, %7 ], [ 0, %65 ], [ 1, %54 ], [ 0, %106 ]
ret i32 %109
}