int matLUbksb(double *aa, int *indx, double *b) {
  double a[4][4];
  int i, j;
  int ii = 0;
  int ll;
  double sum;
  for (i = 0; i < 4; i++)
    for (j = 0; j < 4; j++)
      a[i][j] = aa[i * 4 + j];
  for (i = 0; i < 4; i++) {
    ll = indx[i];
    sum = b[ll];
    b[ll] = b[i];
    if (ii != 0) {
      for (j = ii - 1; j < i - 1; j++)
        sum -= a[i][j] * b[j];
    } else if (sum != 0.0) {
      ii = i;
    }
    b[i] = sum;
  }
  for (i = 3; i >= 0; i--) {
    sum = b[i];
    for (j = i + 1; j < 4; j++)
      sum -= a[i][j] * b[j];
    b[i] = sum / a[i][i];
  }
  return 0;
}


x86 assembly (clang, -O0):
<matLUbksb>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0xac(%rbp)
movl   $0x0,-0xa4(%rbp)
cmpl   $0x4,-0xa4(%rbp)
jge    11b5 <matLUbksb+0xb5>
movl   $0x0,-0xa8(%rbp)
cmpl   $0x4,-0xa8(%rbp)
jge    119c <matLUbksb+0x9c>
mov    -0x8(%rbp),%rax
mov    -0xa4(%rbp),%ecx
shl    $0x2,%ecx
add    -0xa8(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movslq -0xa4(%rbp),%rcx
lea    -0xa0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xa8(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xa8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa8(%rbp)
jmp    113f <matLUbksb+0x3f>
jmp    11a1 <matLUbksb+0xa1>
mov    -0xa4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa4(%rbp)
jmp    1128 <matLUbksb+0x28>
movl   $0x0,-0xa4(%rbp)
cmpl   $0x4,-0xa4(%rbp)
jge    1321 <matLUbksb+0x221>
mov    -0x10(%rbp),%rax
movslq -0xa4(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0xb0(%rbp)
mov    -0x18(%rbp),%rax
movslq -0xb0(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0x18(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0xb0(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
cmpl   $0x0,-0xac(%rbp)
je     12c4 <matLUbksb+0x1c4>
mov    -0xac(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0xa8(%rbp)
mov    -0xa8(%rbp),%eax
mov    -0xa4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12bf <matLUbksb+0x1bf>
movslq -0xa4(%rbp),%rcx
lea    -0xa0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xa8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0xa8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0xb8(%rbp),%xmm1
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0xa8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa8(%rbp)
jmp    1234 <matLUbksb+0x134>
jmp    12f5 <matLUbksb+0x1f5>
movsd  -0xb8(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    12e4 <matLUbksb+0x1e4>
jp     12e4 <matLUbksb+0x1e4>
jmp    12f0 <matLUbksb+0x1f0>
mov    -0xa4(%rbp),%eax
mov    %eax,-0xac(%rbp)
jmp    12f5 <matLUbksb+0x1f5>
movsd  -0xb8(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xa4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa4(%rbp)
jmp    11bf <matLUbksb+0xbf>
movl   $0x3,-0xa4(%rbp)
cmpl   $0x0,-0xa4(%rbp)
jl     142d <matLUbksb+0x32d>
mov    -0x18(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0xa4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa8(%rbp)
cmpl   $0x4,-0xa8(%rbp)
jge    13e0 <matLUbksb+0x2e0>
movslq -0xa4(%rbp),%rcx
lea    -0xa0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xa8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0xa8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0xb8(%rbp),%xmm1
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0xa8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa8(%rbp)
jmp    135f <matLUbksb+0x25f>
movsd  -0xb8(%rbp),%xmm0
ff
movslq -0xa4(%rbp),%rcx
lea    -0xa0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xa4(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xa4(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0xa4(%rbp)
jmp    132b <matLUbksb+0x22b>
xor    %eax,%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<matLUbksb>:
push   %rax
movups 0x70(%rdi),%xmm0
movaps %xmm0,-0x10(%rsp)
movups 0x60(%rdi),%xmm0
movaps %xmm0,-0x20(%rsp)
movups 0x50(%rdi),%xmm0
movaps %xmm0,-0x30(%rsp)
movups 0x40(%rdi),%xmm0
movaps %xmm0,-0x40(%rsp)
movups (%rdi),%xmm0
movups 0x10(%rdi),%xmm1
movups 0x20(%rdi),%xmm2
movups 0x30(%rdi),%xmm3
movaps %xmm3,-0x50(%rsp)
movaps %xmm2,-0x60(%rsp)
movaps %xmm1,-0x70(%rsp)
movaps %xmm0,-0x80(%rsp)
movslq (%rsi),%rax
movsd  (%rdx,%rax,8),%xmm0
movsd  (%rdx),%xmm1
movsd  %xmm1,(%rdx,%rax,8)
movsd  %xmm0,(%rdx)
movslq 0x4(%rsi),%rax
movsd  (%rdx,%rax,8),%xmm0
movsd  0x8(%rdx),%xmm1
movsd  %xmm1,(%rdx,%rax,8)
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
movsd  %xmm0,0x8(%rdx)
movslq 0x8(%rsi),%rax
movsd  (%rdx,%rax,8),%xmm0
movsd  0x10(%rdx),%xmm2
movsd  %xmm2,(%rdx,%rax,8)
jne    1194 <matLUbksb+0x94>
jnp    11a6 <matLUbksb+0xa6>
movsd  -0x40(%rsp),%xmm1
mulsd  (%rdx),%xmm1
subsd  %xmm1,%xmm0
xor    %eax,%eax
jmp    11b3 <matLUbksb+0xb3>
mov    $0x1,%eax
ucomisd %xmm1,%xmm0
jne    11b3 <matLUbksb+0xb3>
jnp    11f1 <matLUbksb+0xf1>
movsd  %xmm0,0x10(%rdx)
movslq 0xc(%rsi),%rcx
movsd  (%rdx,%rcx,8),%xmm0
movsd  0x18(%rdx),%xmm1
movsd  %xmm1,(%rdx,%rcx,8)
movsd  -0x20(%rsp,%rax,8),%xmm1
mulsd  (%rdx,%rax,8),%xmm1
subsd  %xmm1,%xmm0
test   %rax,%rax
jne    1209 <matLUbksb+0x109>
movsd  -0x18(%rsp,%rax,8),%xmm1
mulsd  0x8(%rdx,%rax,8),%xmm1
subsd  %xmm1,%xmm0
jmp    1209 <matLUbksb+0x109>
movsd  %xmm0,0x10(%rdx)
movslq 0xc(%rsi),%rax
movsd  (%rdx,%rax,8),%xmm0
movsd  0x18(%rdx),%xmm1
movsd  %xmm1,(%rdx,%rax,8)
divsd  -0x8(%rsp),%xmm0
movsd  %xmm0,0x18(%rdx)
movsd  0x10(%rdx),%xmm1
movsd  -0x28(%rsp),%xmm2
mulsd  %xmm0,%xmm2
subsd  %xmm2,%xmm1
divsd  -0x30(%rsp),%xmm1
movsd  %xmm1,0x10(%rdx)
movsd  0x8(%rdx),%xmm2
movsd  -0x50(%rsp),%xmm3
mulsd  %xmm1,%xmm3
subsd  %xmm3,%xmm2
movsd  -0x48(%rsp),%xmm3
mulsd  %xmm0,%xmm3
subsd  %xmm3,%xmm2
divsd  -0x58(%rsp),%xmm2
movsd  (%rdx),%xmm3
movsd  %xmm2,0x8(%rdx)
mulsd  -0x78(%rsp),%xmm2
subsd  %xmm2,%xmm3
mulsd  -0x70(%rsp),%xmm1
subsd  %xmm1,%xmm3
mulsd  -0x68(%rsp),%xmm0
subsd  %xmm0,%xmm3
divsd  -0x80(%rsp),%xmm3
movsd  %xmm3,(%rdx)
xor    %eax,%eax
pop    %rcx
ret

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

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