int bandfac(double *a, int ml, int mr, int n) {
  int i, j, k;
  int n1 = n - 1, mt = ml + mr + 1, row, rowi, m, r0, ri0;
  double al;
  for (row = 0; row < n; row++) {
    r0 = row * mt + ml;
    if ((al = a[r0]) == 0.0)
      return (-1 - row);
    al = 1.0 / al;
    m = ((mr) < (n1 - row) ? (mr) : (n1 - row));
    for (j = 1; j <= m; j++)
      a[r0 + j] = a[r0 + j] * al;
    a[r0] = al;
    for (i = 1; i <= ml; i++) {
      rowi = row + i;
      if (rowi > n1)
        break;
      ri0 = rowi * mt + ml;
      al = a[ri0 - i];
      if (al == 0.0)
        continue;
      for (k = 1; k <= m; k++)
        a[ri0 - i + k] = a[ri0 - i + k] - (al * a[r0 + k]);
      a[ri0 - i] = -al;
    }
  }
  return (0);
}


x86 assembly (clang, -O0):
<bandfac>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x14(%rbp),%eax
add    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1352 <bandfac+0x252>
mov    -0x34(%rbp),%eax
imul   -0x30(%rbp),%eax
add    -0x14(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    117b <bandfac+0x7b>
jp     117b <bandfac+0x7b>
mov    $0xffffffff,%eax
sub    -0x34(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1359 <bandfac+0x259>
movsd  0xe7d(%rip),%xmm0
divsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x2c(%rbp),%ecx
sub    -0x34(%rbp),%ecx
cmp    %ecx,%eax
jge    11a9 <bandfac+0xa9>
mov    -0x18(%rbp),%eax
mov    %eax,-0x54(%rbp)
jmp    11b2 <bandfac+0xb2>
mov    -0x2c(%rbp),%eax
sub    -0x34(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    %eax,-0x3c(%rbp)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jg     1202 <bandfac+0x102>
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11bf <bandfac+0xbf>
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     133f <bandfac+0x23f>
mov    -0x34(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jle    1241 <bandfac+0x141>
jmp    133f <bandfac+0x23f>
mov    -0x38(%rbp),%eax
imul   -0x30(%rbp),%eax
add    -0x14(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1282 <bandfac+0x182>
jp     1282 <bandfac+0x182>
jmp    1331 <bandfac+0x231>
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jg     1303 <bandfac+0x203>
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    -0x20(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    -0x20(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1289 <bandfac+0x189>
movsd  -0x50(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    121b <bandfac+0x11b>
jmp    1344 <bandfac+0x244>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    112d <bandfac+0x2d>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<bandfac>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
test   %ecx,%ecx
jle    1535 <bandfac+0x435>
mov    %ecx,%r12d
lea    0x1(%rsi),%eax
lea    (%rax,%rdx,1),%r8d
movslq %r8d,%r9
mov    %ecx,%ebx
movslq %esi,%rbp
mov    %eax,%eax
mov    %rax,0x20(%rsp)
lea    0x8(%rdi),%rax
mov    %rax,-0x60(%rsp)
lea    (%rdx,%rbp,2),%r10d
lea    (%rdx,%rsi,1),%eax
mov    %eax,-0x78(%rsp)
lea    -0x1(%r12),%ecx
lea    0x1(%rbp),%rax
mov    %rax,-0x40(%rsp)
lea    0x38(%rdi,%rbp,8),%r13
mov    %r9,-0x18(%rsp)
lea    0x0(,%r9,8),%rax
mov    %rax,-0x28(%rsp)
mov    %rbp,-0x20(%rsp)
lea    (%rdi,%rbp,8),%r9
lea    0x18(%rdi),%rax
mov    %rax,-0x48(%rsp)
xor    %eax,%eax
xorpd  %xmm0,%xmm0
movsd  0xe77(%rip),%xmm8
movapd 0xe7f(%rip),%xmm2
mov    %r10d,-0x74(%rsp)
mov    %r10d,-0x7c(%rsp)
mov    %ecx,-0x70(%rsp)
mov    %ecx,%ebp
mov    %rbx,-0x50(%rsp)
mov    %r8d,-0x64(%rsp)
mov    %r12,-0x30(%rsp)
mov    %rdx,-0x38(%rsp)
mov    %rsi,0x18(%rsp)
jmp    11ee <bandfac+0xee>
nopl   0x0(%rax)
add    $0x1,%rax
mov    -0x28(%rsp),%rcx
add    %rcx,%r13
mov    -0x68(%rsp),%ebp
add    $0xffffffff,%ebp
add    %rcx,%r9
add    %r8d,-0x7c(%rsp)
cmp    %rbx,%rax
mov    -0x30(%rsp),%r12
mov    -0x38(%rsp),%rdx
je     1535 <bandfac+0x435>
cmp    %ebp,%edx
mov    %ebp,%r15d
cmovl  %edx,%r15d
mov    -0x70(%rsp),%ecx
mov    %ecx,%r14d
sub    %eax,%r14d
cmp    %r14d,%edx
cmovl  %edx,%r14d
mov    %rax,%r11
imul   -0x18(%rsp),%r11
mov    -0x20(%rsp),%rcx
lea    (%r11,%rcx,1),%r10
movsd  (%rdi,%r10,8),%xmm4
ucomisd %xmm0,%xmm4
jne    122c <bandfac+0x12c>
jnp    1539 <bandfac+0x439>
mov    %ebp,-0x68(%rsp)
movapd %xmm8,%xmm3
divsd  %xmm4,%xmm3
mov    %eax,%ecx
not    %ecx
add    %r12d,%ecx
cmp    %edx,%ecx
cmovg  %edx,%ecx
mov    %ecx,-0x80(%rsp)
test   %ecx,%ecx
jle    135b <bandfac+0x25b>
mov    $0x1,%ebx
cmp    $0x4,%r14d
jb     132f <bandfac+0x22f>
mov    %r14d,%ebp
and    $0xfffffffc,%ebp
add    $0xfffffffffffffffc,%rbp
mov    %rbp,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
movapd %xmm3,%xmm4
unpcklpd %xmm3,%xmm4
test   %rbp,%rbp
je     12f4 <bandfac+0x1f4>
mov    %r15d,%ebp
and    $0xfffffffc,%ebp
add    $0xfffffffffffffffc,%rbp
shr    $0x2,%rbp
add    $0x1,%rbp
and    $0xfffffffffffffffe,%rbp
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movupd -0x30(%r13,%rcx,8),%xmm5
movupd -0x20(%r13,%rcx,8),%xmm6
movupd -0x10(%r13,%rcx,8),%xmm7
movupd 0x0(%r13,%rcx,8),%xmm1
mulpd  %xmm4,%xmm5
mulpd  %xmm4,%xmm6
movupd %xmm5,-0x30(%r13,%rcx,8)
movupd %xmm6,-0x20(%r13,%rcx,8)
mulpd  %xmm4,%xmm7
mulpd  %xmm4,%xmm1
movupd %xmm7,-0x10(%r13,%rcx,8)
movupd %xmm1,0x0(%r13,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rbp
jne    12a0 <bandfac+0x1a0>
jmp    12f6 <bandfac+0x1f6>
xor    %ecx,%ecx
mov    %r14d,%ebx
and    $0xfffffffc,%ebx
test   $0x1,%dl
je     1326 <bandfac+0x226>
or     $0x1,%rcx
add    %r10,%rcx
movupd (%rdi,%rcx,8),%xmm1
movupd 0x10(%rdi,%rcx,8),%xmm5
mulpd  %xmm4,%xmm1
mulpd  %xmm4,%xmm5
movupd %xmm1,(%rdi,%rcx,8)
movupd %xmm5,0x10(%rdi,%rcx,8)
cmp    %r14,%rbx
je     135b <bandfac+0x25b>
or     $0x1,%rbx
mov    -0x80(%rsp),%ecx
add    $0xffffffffffffffff,%rbx
nopw   0x0(%rax,%rax,1)
movsd  0x8(%r9,%rbx,8),%xmm1
mulsd  %xmm3,%xmm1
movsd  %xmm1,0x8(%r9,%rbx,8)
add    $0x1,%rbx
cmp    %rcx,%rbx
jb     1340 <bandfac+0x240>
movsd  %xmm3,(%rdi,%r10,8)
test   %esi,%esi
mov    -0x50(%rsp),%rbx
jle    11c0 <bandfac+0xc0>
mov    %r8d,%ecx
imul   %eax,%ecx
add    -0x40(%rsp),%r11
movslq %r15d,%r12
and    $0xfffffffffffffffc,%r12
movslq %r14d,%rbp
add    -0x74(%rsp),%ecx
mov    %ecx,-0x6c(%rsp)
lea    (%rdi,%r11,8),%rcx
mov    %rcx,0x8(%rsp)
lea    (%rcx,%rbp,8),%rcx
mov    %rcx,0x10(%rsp)
movslq -0x80(%rsp),%rdx
mov    -0x60(%rsp),%rcx
lea    (%rcx,%rbp,8),%rcx
mov    %rcx,(%rsp)
mov    %rbp,-0x58(%rsp)
and    $0xfffffffffffffffc,%rbp
mov    %rbp,-0x8(%rsp)
lea    0x1(%rbp),%rcx
mov    %rcx,-0x10(%rsp)
mov    $0x1,%r15d
xor    %r14d,%r14d
mov    -0x7c(%rsp),%r11d
jmp    140d <bandfac+0x30d>
cs nopw 0x0(%rax,%rax,1)
movsd  %xmm3,(%rdi,%r10,8)
mov    0x18(%rsp),%rsi
mov    -0x64(%rsp),%r8d
mov    -0x50(%rsp),%rbx
add    $0x1,%r15
add    $0x1,%r14d
add    -0x78(%rsp),%r11d
cmp    0x20(%rsp),%r15
je     11c0 <bandfac+0xc0>
lea    (%r15,%rax,1),%rcx
cmp    %rbx,%rcx
jae    11c0 <bandfac+0xc0>
mov    %esi,%ebp
sub    %r15d,%ebp
imul   %r8d,%ecx
add    %ebp,%ecx
movslq %ecx,%r10
movsd  (%rdi,%r10,8),%xmm3
ucomisd %xmm0,%xmm3
jne    1436 <bandfac+0x336>
jnp    13f5 <bandfac+0x2f5>
xorpd  %xmm2,%xmm3
cmpl   $0x0,-0x80(%rsp)
jle    13e0 <bandfac+0x2e0>
movslq %r11d,%rcx
mov    $0x1,%ebp
cmpl   $0x4,-0x58(%rsp)
jb     14fb <bandfac+0x3fb>
mov    -0x78(%rsp),%ebx
imul   %r14d,%ebx
add    -0x6c(%rsp),%ebx
movslq %ebx,%rbx
mov    -0x60(%rsp),%rsi
lea    (%rsi,%rbx,8),%rsi
cmp    0x10(%rsp),%rsi
jae    1482 <bandfac+0x382>
mov    (%rsp),%rsi
lea    (%rsi,%rbx,8),%rsi
cmp    %rsi,0x8(%rsp)
jb     14fb <bandfac+0x3fb>
mov    -0x48(%rsp),%rsi
mov    %rcx,%rbx
lea    (%rsi,%rcx,8),%rbp
movapd %xmm3,%xmm4
unpcklpd %xmm3,%xmm4
xor    %r8d,%r8d
nopl   0x0(%rax)
movupd -0x10(%rbp,%r8,8),%xmm1
movupd 0x0(%rbp,%r8,8),%xmm5
movupd 0x8(%r9,%r8,8),%xmm6
movupd 0x18(%r9,%r8,8),%xmm7
mulpd  %xmm4,%xmm6
addpd  %xmm1,%xmm6
mulpd  %xmm4,%xmm7
addpd  %xmm5,%xmm7
movupd %xmm6,-0x10(%rbp,%r8,8)
movupd %xmm7,0x0(%rbp,%r8,8)
add    $0x4,%r8
cmp    %r8,%r12
jne    14a0 <bandfac+0x3a0>
mov    -0x10(%rsp),%rbp
mov    -0x8(%rsp),%rsi
cmp    -0x58(%rsp),%rsi
mov    %rbx,%rcx
je     13e0 <bandfac+0x2e0>
lea    (%rdi,%rcx,8),%rsi
add    $0xffffffffffffffff,%rbp
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  0x8(%r9,%rbp,8),%xmm1
mulsd  %xmm3,%xmm1
addsd  0x8(%rsi,%rbp,8),%xmm1
movsd  %xmm1,0x8(%rsi,%rbp,8)
add    $0x1,%rbp
cmp    %rdx,%rbp
jl     1510 <bandfac+0x410>
jmp    13e0 <bandfac+0x2e0>
xor    %eax,%eax
jmp    153b <bandfac+0x43b>
not    %eax
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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