void bandet(double *e, int m, int n) {
  int i, k, l, mi, lm, km, m2p1;
  double di, du, dl;
  m2p1 = 2 * m + 1;
  if (m <= 0) {
    return;
  }
  for (i = 1; i <= n; i++) {
    di = e[(i - 1) * m2p1 + m];
    mi = (m < (i - 1)) ? m : i - 1;
    if (mi >= 1) {
      for (k = 1; k <= mi; k++) {
        di -= e[(i - 1) * m2p1 - k + m] * e[(i - k - 1) * m2p1 + k + m];
      }
      e[(i - 1) * m2p1 + m] = di;
    }
    lm = (m < (n - i)) ? m : n - i;
    if (lm >= 1) {
      for (l = 1; l <= lm; l++) {
        dl = e[(i + l - 1) * m2p1 - l + m];
        km = ((m - l) < (i - 1)) ? m - l : i - 1;
        if (km >= 1) {
          du = e[(i - 1) * m2p1 + l + m];
          for (k = 1; k <= km; k++) {
            du -= e[(i - 1) * m2p1 - k + m] * e[(i - k - 1) * m2p1 + l + k + m];
            dl -= e[(l + i - 1) * m2p1 - l - k + m] *
                  e[(i - k - 1) * m2p1 + k + m];
          }
          e[(i - 1) * m2p1 + l + m] = du;
        }
        e[(i + l - 1) * m2p1 - l + m] = dl / di;
      }
    }
  }
}


x86 assembly (clang, -O0):
<bandet>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    -0xc(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0xc(%rbp)
jg     1129 <bandet+0x29>
jmp    145c <bandet+0x35c>
movl   $0x1,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     145c <bandet+0x35c>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0xc(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1176 <bandet+0x76>
mov    -0xc(%rbp),%eax
mov    %eax,-0x4c(%rbp)
jmp    117f <bandet+0x7f>
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x1,-0x20(%rbp)
jl     1232 <bandet+0x132>
movl   $0x1,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     1214 <bandet+0x114>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
sub    -0x18(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    -0x18(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
add    -0x18(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%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    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1196 <bandet+0x96>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xc(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    -0x14(%rbp),%ecx
cmp    %ecx,%eax
jge    124e <bandet+0x14e>
mov    -0xc(%rbp),%eax
mov    %eax,-0x50(%rbp)
jmp    1257 <bandet+0x157>
mov    -0x10(%rbp),%eax
sub    -0x14(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x50(%rbp),%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x1,-0x24(%rbp)
jl     1449 <bandet+0x349>
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     1444 <bandet+0x344>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0xc(%rbp),%eax
sub    -0x1c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12c0 <bandet+0x1c0>
mov    -0xc(%rbp),%eax
sub    -0x1c(%rbp),%eax
mov    %eax,-0x54(%rbp)
jmp    12c9 <bandet+0x1c9>
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    %eax,-0x28(%rbp)
cmpl   $0x1,-0x28(%rbp)
jl     140d <bandet+0x30d>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
add    -0x1c(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x1,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x28(%rbp),%eax
jg     13ec <bandet+0x2ec>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
sub    -0x18(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    -0x18(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
add    -0x1c(%rbp),%ecx
add    -0x18(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    -0x14(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
sub    -0x18(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    -0x18(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
add    -0x18(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1301 <bandet+0x201>
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
add    -0x1c(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
divsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    126e <bandet+0x16e>
jmp    1449 <bandet+0x349>
jmp    144e <bandet+0x34e>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1130 <bandet+0x30>
pop    %rbp
ret

x86 assembly (clang, -O2):
<bandet>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %edx,-0x54(%rsp)
mov    %rsi,-0x48(%rsp)
test   %esi,%esi
jle    134a <bandet+0x24a>
cmpl   $0x0,-0x54(%rsp)
jle    134a <bandet+0x24a>
mov    -0x48(%rsp),%rcx
lea    (%rcx,%rcx,1),%eax
lea    (%rcx,%rcx,1),%edx
add    $0x1,%edx
mov    %edx,-0x34(%rsp)
movslq %edx,%rdx
mov    %rdx,-0x18(%rsp)
mov    %ecx,%edx
mov    %rdx,-0x28(%rsp)
mov    -0x54(%rsp),%edx
mov    %rdx,-0x20(%rsp)
mov    %ecx,%edx
neg    %edx
mov    %edx,-0x4c(%rsp)
lea    -0x1(%rcx),%edx
mov    %edx,-0x50(%rsp)
mov    $0x1,%edx
sub    %ecx,%edx
mov    %edx,-0x38(%rsp)
lea    (%rcx,%rcx,2),%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x3c(%rsp)
xor    %ebx,%ebx
mov    $0x1,%r9d
jmp    11af <bandet+0xaf>
xchg   %ax,%ax
mov    -0x10(%rsp),%r9
add    $0x1,%r9
mov    -0x8(%rsp),%rbx
add    $0x1,%ebx
mov    -0x34(%rsp),%ecx
add    %ecx,-0x4c(%rsp)
add    %ecx,-0x50(%rsp)
add    %ecx,-0x38(%rsp)
add    %ecx,-0x3c(%rsp)
cmp    -0x54(%rsp),%ebx
je     134a <bandet+0x24a>
mov    -0x48(%rsp),%r10
cmp    %ebx,%r10d
mov    %ebx,%edx
cmovl  %r10d,%edx
lea    -0x1(%r9),%rbp
mov    %rbp,%rcx
imul   -0x18(%rsp),%rcx
mov    -0x28(%rsp),%r11
mov    %rcx,%rsi
mov    %rcx,-0x30(%rsp)
lea    (%rcx,%r11,1),%r8
movsd  (%rdi,%r8,8),%xmm0
cmp    %r11,%rbp
mov    %ebp,%esi
cmovg  %r10d,%esi
test   %esi,%esi
mov    %rbx,-0x8(%rsp)
jle    1235 <bandet+0x135>
add    $0x1,%edx
add    $0xffffffffffffffff,%rdx
mov    -0x50(%rsp),%esi
mov    -0x4c(%rsp),%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movslq %esi,%rsi
movsd  (%rdi,%rsi,8),%xmm1
movslq %ebx,%rbx
mulsd  (%rdi,%rbx,8),%xmm1
subsd  %xmm1,%xmm0
sub    %eax,%ebx
add    $0xffffffff,%esi
add    $0xffffffffffffffff,%rdx
jne    1210 <bandet+0x110>
movsd  %xmm0,(%rdi,%r8,8)
mov    -0x20(%rsp),%rcx
mov    %r9,-0x10(%rsp)
sub    %r9,%rcx
cmp    -0x28(%rsp),%rcx
cmovg  -0x48(%rsp),%ecx
test   %ecx,%ecx
jle    1180 <bandet+0x80>
mov    -0x30(%rsp),%rdx
add    -0x48(%rsp),%edx
mov    %rdx,-0x30(%rsp)
mov    %ecx,%r10d
mov    $0x1,%r14d
mov    -0x3c(%rsp),%r15d
mov    -0x38(%rsp),%edx
jmp    12a0 <bandet+0x1a0>
cs nopw 0x0(%rax,%rax,1)
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi,%r12,8)
lea    0x1(%r14),%rcx
add    $0x1,%edx
add    %eax,%r15d
cmp    %r10,%r14
mov    %rcx,%r14
jae    1180 <bandet+0x80>
lea    (%r14,%rbp,1),%ecx
imul   -0x34(%rsp),%ecx
sub    %r14d,%ecx
add    -0x48(%rsp),%ecx
movslq %ecx,%r12
movsd  (%rdi,%r12,8),%xmm1
mov    -0x28(%rsp),%rcx
sub    %r14,%rcx
cmp    %rbp,%rcx
cmovge %rbp,%rcx
test   %ecx,%ecx
jle    1280 <bandet+0x180>
mov    -0x30(%rsp),%rsi
add    %r14d,%esi
movslq %esi,%r13
movslq %ecx,%r11
movhpd (%rdi,%r13,8),%xmm1
mov    -0x50(%rsp),%ecx
mov    %ecx,%ebx
mov    %r15d,%r8d
mov    %edx,%r9d
mov    -0x4c(%rsp),%ecx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movslq %ebx,%rbx
movslq %r9d,%r9
movslq %r8d,%r8
movsd  (%rdi,%r8,8),%xmm2
movhpd (%rdi,%rbx,8),%xmm2
movslq %ecx,%rcx
movsd  (%rdi,%rcx,8),%xmm3
movhpd (%rdi,%r9,8),%xmm3
mulpd  %xmm2,%xmm3
subpd  %xmm3,%xmm1
add    $0x1,%rsi
sub    %eax,%ecx
sub    %eax,%r9d
add    $0xffffffff,%r8d
add    $0xffffffff,%ebx
cmp    %r11,%rsi
jl     1300 <bandet+0x200>
movhpd %xmm1,(%rdi,%r13,8)
jmp    1280 <bandet+0x180>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @bandet(double* %0, i32 %1, i32 %2) {
%4 = alloca double*
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca double
%15 = alloca double
%16 = alloca double
store double* %0, double** %4
store i32 %1, i32* %5
store i32 %2, i32* %6
%17 = load i32, i32* %5
%18 = mul nsw i32 2, %17
%19 = add nsw i32 %18, 1
store i32 %19, i32* %13
%20 = load i32, i32* %5
%21 = icmp sle i32 %20, 0
br i1 %21, label %22, label %23
22: ; preds = %3
br label %277
23: ; preds = %3
store i32 1, i32* %7
br label %24
24: ; preds = %274, %23
%25 = load i32, i32* %7
%26 = load i32, i32* %6
%27 = icmp sle i32 %25, %26
br i1 %27, label %28, label %277
28: ; preds = %24
%29 = load double*, double** %4
%30 = load i32, i32* %7
%31 = sub nsw i32 %30, 1
%32 = load i32, i32* %13
%33 = mul nsw i32 %31, %32
%34 = load i32, i32* %5
%35 = add nsw i32 %33, %34
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %29, i64 %36
%38 = load double, double* %37
store double %38, double* %14
%39 = load i32, i32* %5
%40 = load i32, i32* %7
%41 = sub nsw i32 %40, 1
%42 = icmp slt i32 %39, %41
br i1 %42, label %43, label %45
43: ; preds = %28
%44 = load i32, i32* %5
br label %48
45: ; preds = %28
%46 = load i32, i32* %7
%47 = sub nsw i32 %46, 1
br label %48
48: ; preds = %45, %43
%49 = phi i32 [ %44, %43 ], [ %47, %45 ]
store i32 %49, i32* %10
%50 = load i32, i32* %10
%51 = icmp sge i32 %50, 1
br i1 %51, label %52, label %101
52: ; preds = %48
store i32 1, i32* %8
br label %53
53: ; preds = %87, %52
%54 = load i32, i32* %8
%55 = load i32, i32* %10
%56 = icmp sle i32 %54, %55
br i1 %56, label %57, label %90
57: ; preds = %53
%58 = load double*, double** %4
%59 = load i32, i32* %7
%60 = sub nsw i32 %59, 1
%61 = load i32, i32* %13
%62 = mul nsw i32 %60, %61
%63 = load i32, i32* %8
%64 = sub nsw i32 %62, %63
%65 = load i32, i32* %5
%66 = add nsw i32 %64, %65
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds double, double* %58, i64 %67
%69 = load double, double* %68
%70 = load double*, double** %4
%71 = load i32, i32* %7
%72 = load i32, i32* %8
%73 = sub nsw i32 %71, %72
%74 = sub nsw i32 %73, 1
%75 = load i32, i32* %13
%76 = mul nsw i32 %74, %75
%77 = load i32, i32* %8
%78 = add nsw i32 %76, %77
%79 = load i32, i32* %5
%80 = add nsw i32 %78, %79
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %70, i64 %81
%83 = load double, double* %82
%84 = load double, double* %14
%85 = fneg double %69
%86 = call double @llvm.fmuladd.f64(double %85, double %83, double %84)
store double %86, double* %14
br label %87
87: ; preds = %57
%88 = load i32, i32* %8
%89 = add nsw i32 %88, 1
store i32 %89, i32* %8
br label %53
90: ; preds = %53
%91 = load double, double* %14
%92 = load double*, double** %4
%93 = load i32, i32* %7
%94 = sub nsw i32 %93, 1
%95 = load i32, i32* %13
%96 = mul nsw i32 %94, %95
%97 = load i32, i32* %5
%98 = add nsw i32 %96, %97
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds double, double* %92, i64 %99
store double %91, double* %100
br label %101
101: ; preds = %90, %48
%102 = load i32, i32* %5
%103 = load i32, i32* %6
%104 = load i32, i32* %7
%105 = sub nsw i32 %103, %104
%106 = icmp slt i32 %102, %105
br i1 %106, label %107, label %109
107: ; preds = %101
%108 = load i32, i32* %5
br label %113
109: ; preds = %101
%110 = load i32, i32* %6
%111 = load i32, i32* %7
%112 = sub nsw i32 %110, %111
br label %113
113: ; preds = %109, %107
%114 = phi i32 [ %108, %107 ], [ %112, %109 ]
store i32 %114, i32* %11
%115 = load i32, i32* %11
%116 = icmp sge i32 %115, 1
br i1 %116, label %117, label %273
117: ; preds = %113
store i32 1, i32* %9
br label %118
118: ; preds = %269, %117
%119 = load i32, i32* %9
%120 = load i32, i32* %11
%121 = icmp sle i32 %119, %120
br i1 %121, label %122, label %272
122: ; preds = %118
%123 = load double*, double** %4
%124 = load i32, i32* %7
%125 = load i32, i32* %9
%126 = add nsw i32 %124, %125
%127 = sub nsw i32 %126, 1
%128 = load i32, i32* %13
%129 = mul nsw i32 %127, %128
%130 = load i32, i32* %9
%131 = sub nsw i32 %129, %130
%132 = load i32, i32* %5
%133 = add nsw i32 %131, %132
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds double, double* %123, i64 %134
%136 = load double, double* %135
store double %136, double* %16
%137 = load i32, i32* %5
%138 = load i32, i32* %9
%139 = sub nsw i32 %137, %138
%140 = load i32, i32* %7
%141 = sub nsw i32 %140, 1
%142 = icmp slt i32 %139, %141
br i1 %142, label %143, label %147
143: ; preds = %122
%144 = load i32, i32* %5
%145 = load i32, i32* %9
%146 = sub nsw i32 %144, %145
br label %150
147: ; preds = %122
%148 = load i32, i32* %7
%149 = sub nsw i32 %148, 1
br label %150
150: ; preds = %147, %143
%151 = phi i32 [ %146, %143 ], [ %149, %147 ]
store i32 %151, i32* %12
%152 = load i32, i32* %12
%153 = icmp sge i32 %152, 1
br i1 %153, label %154, label %252
154: ; preds = %150
%155 = load double*, double** %4
%156 = load i32, i32* %7
%157 = sub nsw i32 %156, 1
%158 = load i32, i32* %13
%159 = mul nsw i32 %157, %158
%160 = load i32, i32* %9
%161 = add nsw i32 %159, %160
%162 = load i32, i32* %5
%163 = add nsw i32 %161, %162
%164 = sext i32 %163 to i64
%165 = getelementptr inbounds double, double* %155, i64 %164
%166 = load double, double* %165
store double %166, double* %15
store i32 1, i32* %8
br label %167
167: ; preds = %236, %154
%168 = load i32, i32* %8
%169 = load i32, i32* %12
%170 = icmp sle i32 %168, %169
br i1 %170, label %171, label %239
171: ; preds = %167
%172 = load double*, double** %4
%173 = load i32, i32* %7
%174 = sub nsw i32 %173, 1
%175 = load i32, i32* %13
%176 = mul nsw i32 %174, %175
%177 = load i32, i32* %8
%178 = sub nsw i32 %176, %177
%179 = load i32, i32* %5
%180 = add nsw i32 %178, %179
%181 = sext i32 %180 to i64
%182 = getelementptr inbounds double, double* %172, i64 %181
%183 = load double, double* %182
%184 = load double*, double** %4
%185 = load i32, i32* %7
%186 = load i32, i32* %8
%187 = sub nsw i32 %185, %186
%188 = sub nsw i32 %187, 1
%189 = load i32, i32* %13
%190 = mul nsw i32 %188, %189
%191 = load i32, i32* %9
%192 = add nsw i32 %190, %191
%193 = load i32, i32* %8
%194 = add nsw i32 %192, %193
%195 = load i32, i32* %5
%196 = add nsw i32 %194, %195
%197 = sext i32 %196 to i64
%198 = getelementptr inbounds double, double* %184, i64 %197
%199 = load double, double* %198
%200 = load double, double* %15
%201 = fneg double %183
%202 = call double @llvm.fmuladd.f64(double %201, double %199, double %200)
store double %202, double* %15
%203 = load double*, double** %4
%204 = load i32, i32* %9
%205 = load i32, i32* %7
%206 = add nsw i32 %204, %205
%207 = sub nsw i32 %206, 1
%208 = load i32, i32* %13
%209 = mul nsw i32 %207, %208
%210 = load i32, i32* %9
%211 = sub nsw i32 %209, %210
%212 = load i32, i32* %8
%213 = sub nsw i32 %211, %212
%214 = load i32, i32* %5
%215 = add nsw i32 %213, %214
%216 = sext i32 %215 to i64
%217 = getelementptr inbounds double, double* %203, i64 %216
%218 = load double, double* %217
%219 = load double*, double** %4
%220 = load i32, i32* %7
%221 = load i32, i32* %8
%222 = sub nsw i32 %220, %221
%223 = sub nsw i32 %222, 1
%224 = load i32, i32* %13
%225 = mul nsw i32 %223, %224
%226 = load i32, i32* %8
%227 = add nsw i32 %225, %226
%228 = load i32, i32* %5
%229 = add nsw i32 %227, %228
%230 = sext i32 %229 to i64
%231 = getelementptr inbounds double, double* %219, i64 %230
%232 = load double, double* %231
%233 = load double, double* %16
%234 = fneg double %218
%235 = call double @llvm.fmuladd.f64(double %234, double %232, double %233)
store double %235, double* %16
br label %236
236: ; preds = %171
%237 = load i32, i32* %8
%238 = add nsw i32 %237, 1
store i32 %238, i32* %8
br label %167
239: ; preds = %167
%240 = load double, double* %15
%241 = load double*, double** %4
%242 = load i32, i32* %7
%243 = sub nsw i32 %242, 1
%244 = load i32, i32* %13
%245 = mul nsw i32 %243, %244
%246 = load i32, i32* %9
%247 = add nsw i32 %245, %246
%248 = load i32, i32* %5
%249 = add nsw i32 %247, %248
%250 = sext i32 %249 to i64
%251 = getelementptr inbounds double, double* %241, i64 %250
store double %240, double* %251
br label %252
252: ; preds = %239, %150
%253 = load double, double* %16
%254 = load double, double* %14
%255 = fdiv double %253, %254
%256 = load double*, double** %4
%257 = load i32, i32* %7
%258 = load i32, i32* %9
%259 = add nsw i32 %257, %258
%260 = sub nsw i32 %259, 1
%261 = load i32, i32* %13
%262 = mul nsw i32 %260, %261
%263 = load i32, i32* %9
%264 = sub nsw i32 %262, %263
%265 = load i32, i32* %5
%266 = add nsw i32 %264, %265
%267 = sext i32 %266 to i64
%268 = getelementptr inbounds double, double* %256, i64 %267
store double %255, double* %268
br label %269
269: ; preds = %252
%270 = load i32, i32* %9
%271 = add nsw i32 %270, 1
store i32 %271, i32* %9
br label %118
272: ; preds = %118
br label %273
273: ; preds = %272, %113
br label %274
274: ; preds = %273
%275 = load i32, i32* %7
%276 = add nsw i32 %275, 1
store i32 %276, i32* %7
br label %24
277: ; preds = %22, %24
ret void
}

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