void lubksb(double carray[5][5], int *indx, double *barray, int n) {
  int i, j, ii, ll;
  double sum;
  ii = 0;
  for (i = 0; i < n; i++) {
    ll = indx[i];
    sum = barray[ll];
    barray[ll] = barray[i];
    if (ii != 0) {
      for (j = ii; j < i; j++) {
        sum -= carray[i][j] * barray[j];
      }
    } else if (sum != 0.0) {
      ii = i;
    }
    barray[i] = sum;
  }
  for (i = n - 1; i >= 0; i--) {
    sum = barray[i];
    if (i < n - 1) {
      for (j = i + 1; j < n; j++) {
        sum -= carray[i][j] * barray[j];
      }
    }
    barray[i] = sum / carray[i][i];
  }
}


x86 assembly (clang, -O0):
<lubksb>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    122c <lubksb+0x12c>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
cmpl   $0x0,-0x28(%rbp)
je     11e4 <lubksb+0xe4>
mov    -0x28(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11df <lubksb+0xdf>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
imul   $0x28,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1177 <lubksb+0x77>
jmp    120c <lubksb+0x10c>
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1201 <lubksb+0x101>
jp     1201 <lubksb+0x101>
jmp    1207 <lubksb+0x107>
mov    -0x20(%rbp),%eax
mov    %eax,-0x28(%rbp)
jmp    120c <lubksb+0x10c>
movsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1121 <lubksb+0x21>
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x0,-0x20(%rbp)
jl     1310 <lubksb+0x210>
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12d8 <lubksb+0x1d8>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12d3 <lubksb+0x1d3>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
imul   $0x28,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    126b <lubksb+0x16b>
jmp    12d8 <lubksb+0x1d8>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
imul   $0x28,%rcx,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
jmp    1235 <lubksb+0x135>
pop    %rbp
ret

x86 assembly (clang, -O2):
<lubksb>:
push   %r14
push   %rbx
test   %ecx,%ecx
jle    12ac <lubksb+0x1ac>
mov    %ecx,%r8d
lea    0x8(%rdi),%r9
xor    %r11d,%r11d
xorpd  %xmm0,%xmm0
xor    %r10d,%r10d
jmp    1146 <lubksb+0x46>
xchg   %ax,%ax
ucomisd %xmm0,%xmm1
jne    112c <lubksb+0x2c>
jnp    11ce <lubksb+0xce>
mov    %r11d,%r10d
movsd  %xmm1,(%rdx,%r11,8)
add    $0x1,%r11
add    $0x28,%r9
cmp    %r8,%r11
je     11d6 <lubksb+0xd6>
movslq (%rsi,%r11,4),%rax
movsd  (%rdx,%rax,8),%xmm1
movsd  (%rdx,%r11,8),%xmm2
movsd  %xmm2,(%rdx,%rax,8)
test   %r10d,%r10d
je     1120 <lubksb+0x20>
movslq %r10d,%rax
cmp    %rax,%r11
jle    112f <lubksb+0x2f>
mov    %r11d,%ebx
sub    %r10d,%ebx
lea    0x1(%rax),%r14
test   $0x1,%bl
je     118f <lubksb+0x8f>
lea    (%r11,%r11,4),%rbx
lea    (%rdi,%rbx,8),%rbx
movsd  (%rbx,%rax,8),%xmm2
mulsd  (%rdx,%rax,8),%xmm2
subsd  %xmm2,%xmm1
mov    %r14,%rax
cmp    %r14,%r11
je     112f <lubksb+0x2f>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  -0x8(%r9,%rax,8),%xmm2
movsd  (%r9,%rax,8),%xmm3
mulsd  (%rdx,%rax,8),%xmm2
mulsd  0x8(%rdx,%rax,8),%xmm3
subsd  %xmm2,%xmm1
subsd  %xmm3,%xmm1
add    $0x2,%rax
cmp    %rax,%r11
jne    11a0 <lubksb+0xa0>
jmp    112f <lubksb+0x2f>
xor    %r10d,%r10d
jmp    112f <lubksb+0x2f>
test   %ecx,%ecx
jle    12ac <lubksb+0x1ac>
lea    -0x1(%rcx),%eax
mov    %eax,%r9d
lea    (%r9,%r9,4),%rax
lea    (%rdi,%rax,8),%rsi
add    $0x8,%rsi
xor    %r10d,%r10d
mov    %r9,%rax
mov    %r8,%r11
jmp    122e <lubksb+0x12e>
nopl   0x0(%rax,%rax,1)
lea    (%rax,%rax,4),%rbx
lea    (%rdi,%rbx,8),%rbx
divsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,(%rdx,%rax,8)
lea    -0x1(%rax),%rbx
add    $0xffffffffffffffff,%r11
add    $0x1,%r10d
add    $0xffffffffffffffd8,%rsi
test   %rax,%rax
mov    %rbx,%rax
jle    12ac <lubksb+0x1ac>
movsd  (%rdx,%rax,8),%xmm0
cmp    %r9,%rax
jge    1200 <lubksb+0x100>
lea    0x1(%rax),%ebx
cmp    %ecx,%ebx
jge    1200 <lubksb+0x100>
test   $0x1,%r10b
jne    1250 <lubksb+0x150>
mov    %r11,%rbx
cmp    $0x1,%r10d
jne    1280 <lubksb+0x180>
jmp    1200 <lubksb+0x100>
lea    (%rax,%rax,4),%rbx
lea    (%rdi,%rbx,8),%rbx
movsd  (%rbx,%r11,8),%xmm1
mulsd  (%rdx,%r11,8),%xmm1
subsd  %xmm1,%xmm0
lea    0x1(%r11),%rbx
cmp    $0x1,%r10d
je     1200 <lubksb+0x100>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  -0x8(%rsi,%rbx,8),%xmm1
movsd  (%rsi,%rbx,8),%xmm2
mulsd  (%rdx,%rbx,8),%xmm1
mulsd  0x8(%rdx,%rbx,8),%xmm2
subsd  %xmm1,%xmm0
subsd  %xmm2,%xmm0
add    $0x2,%rbx
cmp    %ebx,%r8d
jne    1280 <lubksb+0x180>
jmp    1200 <lubksb+0x100>
pop    %rbx
pop    %r14
ret

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

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