void banbks(float **a, unsigned long n, int m1, int m2, float **al,
            unsigned long indx[], float b[]) {
  unsigned long i, k, l;
  int mm;
  float dum;
  mm = m1 + m2 + 1;
  l = m1;
  for (k = 1; k <= n; k++) {
    i = indx[k];
    if (i != k) {
      dum = (b[k]);
      (b[k]) = (b[i]);
      (b[i]) = dum;
    }
    if (l < n)
      l++;
    for (i = k + 1; i <= l; i++)
      b[i] -= al[k][i - k] * b[k];
  }
  l = 1;
  for (i = n; i >= 1; i--) {
    dum = b[i];
    for (k = 2; k <= l; k++)
      dum -= a[i][k] * b[k + i - 1];
    b[i] = dum / a[i][1];
    if (l < mm)
      l++;
  }
}


x86 assembly (clang, -O0):
<banbks>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    -0x14(%rbp),%eax
add    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
movslq -0x14(%rbp),%rax
mov    %rax,-0x40(%rbp)
movq   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
ja     124d <banbks+0x14d>
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x38(%rbp),%rax
je     11a4 <banbks+0xa4>
mov    0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
mov    0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x48(%rbp),%xmm0
mov    0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    11be <banbks+0xbe>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x40(%rbp),%rax
ja     1237 <banbks+0x137>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x30(%rbp),%rcx
sub    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movd   %xmm0,%edx
xor    $0x80000000,%edx
movd   %edx,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    11ca <banbks+0xca>
jmp    123c <banbks+0x13c>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    113a <banbks+0x3a>
movq   $0x1,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x30(%rbp)
cmpq   $0x1,-0x30(%rbp)
jb     1340 <banbks+0x240>
mov    0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
movq   $0x2,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x40(%rbp),%rax
ja     12ea <banbks+0x1ea>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
add    -0x30(%rbp),%rcx
sub    $0x1,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x48(%rbp),%xmm1
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1282 <banbks+0x182>
movss  -0x48(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
divss  0x4(%rax),%xmm0
mov    0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
cmp    %rcx,%rax
jae    132a <banbks+0x22a>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    132f <banbks+0x22f>
mov    -0x30(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x30(%rbp)
jmp    125d <banbks+0x15d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<banbks>:
push   %r15
push   %r14
push   %rbx
test   %rsi,%rsi
je     126f <banbks+0x16f>
mov    0x20(%rsp),%r15
lea    (%rdx,%rcx,1),%r10d
add    $0x1,%r10d
movslq %edx,%rdx
mov    $0x1,%r11d
jmp    1135 <banbks+0x35>
cs nopw 0x0(%rax,%rax,1)
cmp    %rsi,%rbx
je     11a9 <banbks+0xa9>
mov    %r11,%rbx
mov    (%r9,%r11,8),%rax
cmp    %r11,%rax
je     1159 <banbks+0x59>
movss  (%r15,%rbx,4),%xmm0
movss  (%r15,%rax,4),%xmm1
movss  %xmm1,(%r15,%rbx,4)
movss  %xmm0,(%r15,%rax,4)
cmp    %rsi,%rdx
adc    $0x0,%rdx
lea    0x1(%rbx),%r11
cmp    %rdx,%rbx
jae    1130 <banbks+0x30>
mov    (%r8,%rbx,8),%rcx
add    $0x4,%rcx
mov    %r11,%rax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rcx),%xmm0
movss  (%r15,%rax,4),%xmm1
mulss  (%r15,%rbx,4),%xmm0
subss  %xmm0,%xmm1
movss  %xmm1,(%r15,%rax,4)
add    $0x1,%rax
add    $0x4,%rcx
cmp    %rdx,%rax
jbe    1180 <banbks+0x80>
jmp    1130 <banbks+0x30>
test   %rsi,%rsi
je     126f <banbks+0x16f>
movslq %r10d,%r8
lea    (%r15,%rsi,4),%r11
add    $0x8,%r11
mov    $0x1,%r10d
jmp    11ef <banbks+0xef>
cs nopw 0x0(%rax,%rax,1)
nop
divss  0x4(%rcx),%xmm0
movss  %xmm0,(%r15,%r9,4)
cmp    %r8,%r10
adc    $0x0,%r10
add    $0xfffffffffffffffc,%r11
test   %rsi,%rsi
je     126f <banbks+0x16f>
mov    %rsi,%r9
movss  (%r15,%rsi,4),%xmm0
add    $0xffffffffffffffff,%rsi
mov    (%rdi,%r9,8),%rcx
cmp    $0x2,%r10
jb     11d0 <banbks+0xd0>
lea    -0x1(%r10),%r14
mov    $0x2,%edx
cmp    $0x2,%r10
je     124e <banbks+0x14e>
mov    %r14,%rbx
and    $0xfffffffffffffffe,%rbx
xor    %edx,%edx
xchg   %ax,%ax
movss  0x8(%rcx,%rdx,4),%xmm1
movss  0xc(%rcx,%rdx,4),%xmm2
mulss  -0x4(%r11,%rdx,4),%xmm1
mulss  (%r11,%rdx,4),%xmm2
subss  %xmm1,%xmm0
subss  %xmm2,%xmm0
add    $0x2,%rdx
cmp    %rdx,%rbx
jne    1220 <banbks+0x120>
add    $0x2,%rdx
test   $0x1,%r14b
je     11d0 <banbks+0xd0>
movss  (%rcx,%rdx,4),%xmm1
add    %rsi,%rdx
mulss  (%r15,%rdx,4),%xmm1
subss  %xmm1,%xmm0
jmp    11d0 <banbks+0xd0>
pop    %rbx
pop    %r14
pop    %r15
ret

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

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