void dwt2d(void (*wavelet)(float *, float *, int, int, int), float *out,
           float *in, int N0, int W, int H, int SO, int SI, int SW) {
  for (int j = 0; j < H; ++j) {
    wavelet(out + SO * SW * j, in + SI * SW * j, W, SO, SI);
    for (int i = 0; i < W; ++i)
      in[(SW * j + i) * SI] = out[(SW * j + i) * SO];
  }
  for (int i = 0; i < W; ++i)
    wavelet(out + SO * i, in + SI * i, H, SO * SW, SI * SW);
  int W2 = (W + 1) / 2, H2 = (H + 1) / 2;
  for (int j = 0; j < H2; ++j)
    for (int i = 0; i < W2; ++i)
      in[(SW * j + i) * SI] = out[(SW * j + i) * SO];
  if (W2 >= N0 && H2 >= N0)
    dwt2d(wavelet, out, in, N0, W2, H2, SO, SI, SW);
}


x86 assembly (clang, -O0):
<dwt2d>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11f5 <dwt2d+0xe5>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rdi
mov    0x10(%rbp),%ecx
imul   0x20(%rbp),%ecx
imul   -0x28(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rdi
mov    -0x18(%rbp),%rsi
mov    0x18(%rbp),%ecx
imul   0x20(%rbp),%ecx
imul   -0x28(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rsi
mov    -0x20(%rbp),%edx
mov    0x10(%rbp),%ecx
mov    0x18(%rbp),%r8d
call   *%rax
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11e2 <dwt2d+0xd2>
mov    -0x10(%rbp),%rax
mov    0x20(%rbp),%ecx
imul   -0x28(%rbp),%ecx
add    -0x2c(%rbp),%ecx
imul   0x10(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    0x20(%rbp),%ecx
imul   -0x28(%rbp),%ecx
add    -0x2c(%rbp),%ecx
imul   0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1194 <dwt2d+0x84>
jmp    11e7 <dwt2d+0xd7>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    113f <dwt2d+0x2f>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1259 <dwt2d+0x149>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rdi
mov    0x10(%rbp),%ecx
imul   -0x30(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rdi
mov    -0x18(%rbp),%rsi
mov    0x18(%rbp),%ecx
imul   -0x30(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rsi
mov    -0x24(%rbp),%edx
mov    0x10(%rbp),%ecx
imul   0x20(%rbp),%ecx
mov    0x18(%rbp),%r8d
imul   0x20(%rbp),%r8d
call   *%rax
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11fc <dwt2d+0xec>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x34(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    12f6 <dwt2d+0x1e6>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    12e3 <dwt2d+0x1d3>
mov    -0x10(%rbp),%rax
mov    0x20(%rbp),%ecx
imul   -0x3c(%rbp),%ecx
add    -0x40(%rbp),%ecx
imul   0x10(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    0x20(%rbp),%ecx
imul   -0x3c(%rbp),%ecx
add    -0x40(%rbp),%ecx
imul   0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1295 <dwt2d+0x185>
jmp    12e8 <dwt2d+0x1d8>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1282 <dwt2d+0x172>
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jl     1342 <dwt2d+0x232>
mov    -0x38(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jl     1342 <dwt2d+0x232>
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    -0x1c(%rbp),%ecx
mov    -0x34(%rbp),%r8d
mov    -0x38(%rbp),%r9d
mov    0x10(%rbp),%r11d
mov    0x18(%rbp),%r10d
mov    0x20(%rbp),%eax
mov    %r11d,(%rsp)
mov    %r10d,0x8(%rsp)
mov    %eax,0x10(%rsp)
call   1030 <dwt2d@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dwt2d>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x108,%rsp
mov    %r9d,%ebx
mov    %r8d,%r10d
mov    %rdi,%r13
mov    0x150(%rsp),%ebp
mov    0x148(%rsp),%r8d
mov    0x140(%rsp),%edi
mov    %ebp,%eax
imul   %edi,%eax
mov    %eax,0xa0(%rsp)
mov    %ebp,%eax
imul   %r8d,%eax
mov    %eax,0x3c(%rsp)
movslq %edi,%r11
movslq %r8d,%rax
mov    %eax,%r8d
xor    $0x1,%r8d
mov    %r11d,%edi
xor    $0x1,%edi
or     %r8d,%edi
setne  0x7(%rsp)
mov    %rax,%rdi
shl    $0x5,%rdi
mov    %rdi,0xc0(%rsp)
movslq %ebp,%rdi
mov    %rdi,%rbp
imul   %rax,%rbp
shl    $0x2,%rbp
mov    %rbp,0xe8(%rsp)
mov    %rax,%rbp
shl    $0x6,%rbp
mov    %rbp,0xc8(%rsp)
mov    %r11,%rbp
shl    $0x5,%rbp
mov    %rbp,0xd8(%rsp)
mov    %rdi,0x78(%rsp)
imul   %r11,%rdi
shl    $0x2,%rdi
mov    %rdi,0xe0(%rsp)
mov    %r11,%rdi
shl    $0x6,%rdi
mov    %rdi,0xb8(%rsp)
mov    %rdx,0x8(%rsp)
add    $0x10,%rdx
mov    %rdx,0xb0(%rsp)
lea    0x10(%rsi),%rdx
mov    %rdx,0xa8(%rsp)
lea    0x0(,%rax,4),%r12
mov    %rax,0x50(%rsp)
lea    0x0(,%rax,8),%r14
lea    0x0(,%r11,4),%r9
mov    %r11,0x58(%rsp)
lea    0x0(,%r11,8),%r15
mov    %r13,%r11
mov    %r13,0xd0(%rsp)
mov    %r9,0x68(%rsp)
mov    %ecx,0xa4(%rsp)
mov    %rsi,0x10(%rsp)
mov    %r12,0x70(%rsp)
mov    %r14,0xf8(%rsp)
mov    %r15,0xf0(%rsp)
mov    %r10d,%eax
mov    %rax,0x18(%rsp)
mov    %rbx,0x30(%rsp)
test   %ebx,%ebx
mov    %r10,0x98(%rsp)
jle    1580 <dwt2d+0x480>
mov    0x30(%rsp),%eax
mov    %rax,0x90(%rsp)
mov    0x18(%rsp),%rcx
mov    %ecx,%eax
and    $0xfffffff8,%eax
mov    %rax,0x48(%rsp)
add    $0xfffffffffffffff8,%rax
mov    %rax,0x40(%rsp)
mov    %rax,%rdx
shr    $0x3,%rdx
add    $0x1,%rdx
mov    %rdx,%rax
mov    %rdx,0x80(%rsp)
movabs $0x3ffffffffffffffe,%rax
ff ff 3f
and    %rax,%rdx
mov    %rdx,0x100(%rsp)
neg    %rcx
mov    %rcx,0x88(%rsp)
mov    $0x1,%eax
mov    %rax,0x60(%rsp)
xor    %r13d,%r13d
mov    0xa8(%rsp),%rax
mov    %rax,0x20(%rsp)
mov    0xb0(%rsp),%rax
mov    %rax,0x28(%rsp)
xor    %r12d,%r12d
jmp    133b <dwt2d+0x23b>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r12
mov    0x28(%rsp),%rax
add    0xe8(%rsp),%rax
mov    %rax,0x28(%rsp)
mov    0x20(%rsp),%rax
add    0xe0(%rsp),%rax
mov    %rax,0x20(%rsp)
mov    0x78(%rsp),%rax
add    %rax,%r13
add    %rax,0x60(%rsp)
cmp    0x90(%rsp),%r12
mov    0xd0(%rsp),%r11
je     1580 <dwt2d+0x480>
mov    %r12d,%eax
imul   0x150(%rsp),%eax
mov    %eax,%ecx
mov    0x140(%rsp),%ebp
imul   %ebp,%ecx
movslq %ecx,%rcx
mov    0x10(%rsp),%rdx
lea    (%rdx,%rcx,4),%rdi
mov    0x148(%rsp),%r8d
imul   %r8d,%eax
cltq
mov    0x8(%rsp),%rcx
lea    (%rcx,%rax,4),%rsi
mov    %r10d,%edx
mov    %ebp,%ecx
call   *%r11
mov    0x68(%rsp),%r9
mov    0x98(%rsp),%r10
test   %r10d,%r10d
mov    0xd8(%rsp),%rax
jle    12f0 <dwt2d+0x1f0>
mov    %r12,%r11
imul   0x78(%rsp),%r11
cmp    $0x8,%r10d
setb   %cl
or     0x7(%rsp),%cl
jne    13f3 <dwt2d+0x2f3>
mov    0x18(%rsp),%rcx
add    %r11,%rcx
mov    0x8(%rsp),%rdi
lea    (%rdi,%r11,4),%rdx
mov    0x10(%rsp),%rbp
lea    0x0(,%rcx,4),%rsi
add    %rbp,%rsi
cmp    %rsi,%rdx
jae    14a7 <dwt2d+0x3a7>
lea    (%rdi,%rcx,4),%rcx
lea    0x0(,%r11,4),%rdx
add    %rbp,%rdx
cmp    %rcx,%rdx
jae    14a7 <dwt2d+0x3a7>
xor    %edx,%edx
mov    %rdx,%rcx
testb  $0x1,0x18(%rsp)
je     142d <dwt2d+0x32d>
add    %rdx,%r11
mov    %r11,%rcx
imul   0x58(%rsp),%rcx
mov    0x10(%rsp),%rsi
movss  (%rsi,%rcx,4),%xmm0
imul   0x50(%rsp),%r11
mov    0x8(%rsp),%rcx
movss  %xmm0,(%rcx,%r11,4)
mov    %rdx,%rcx
or     $0x1,%rcx
not    %rdx
cmp    0x88(%rsp),%rdx
je     12f0 <dwt2d+0x1f0>
mov    0x18(%rsp),%rax
sub    %rcx,%rax
lea    (%rcx,%r13,1),%rdx
mov    0x70(%rsp),%rdi
mov    %rdi,%rsi
imul   %rdx,%rsi
imul   %r9,%rdx
add    0x60(%rsp),%rcx
imul   %rcx,%rdi
imul   %r9,%rcx
mov    0x10(%rsp),%rbp
mov    0x8(%rsp),%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  0x0(%rbp,%rdx,1),%xmm0
movss  %xmm0,(%rbx,%rsi,1)
movss  0x0(%rbp,%rcx,1),%xmm0
movss  %xmm0,(%rbx,%rdi,1)
add    %r14,%rbx
add    %r15,%rbp
add    $0xfffffffffffffffe,%rax
jne    1480 <dwt2d+0x380>
jmp    12f0 <dwt2d+0x1f0>
mov    %rbp,%r9
cmpq   $0x0,0x40(%rsp)
mov    0xc8(%rsp),%rbx
je     1572 <dwt2d+0x472>
mov    0x100(%rsp),%rdx
mov    0x20(%rsp),%rsi
mov    0x28(%rsp),%rdi
xor    %ecx,%ecx
mov    0xc0(%rsp),%rbp
mov    0xb8(%rsp),%r8
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movups -0x10(%rsi),%xmm0
movups (%rsi),%xmm1
movups %xmm0,-0x10(%rdi)
movups %xmm1,(%rdi)
movups -0x10(%rsi,%rax,1),%xmm0
movups (%rsi,%rax,1),%xmm1
movups %xmm0,-0x10(%rdi,%rbp,1)
movups %xmm1,(%rdi,%rbp,1)
add    $0x10,%rcx
add    %rbx,%rdi
add    %r8,%rsi
add    $0xfffffffffffffffe,%rdx
jne    14f0 <dwt2d+0x3f0>
testb  $0x1,0x80(%rsp)
01
je     1555 <dwt2d+0x455>
add    %r11,%rcx
mov    %rcx,%rdx
imul   0x58(%rsp),%rdx
movups (%r9,%rdx,4),%xmm0
movups 0x10(%r9,%rdx,4),%xmm1
imul   0x50(%rsp),%rcx
mov    0x8(%rsp),%rax
movups %xmm0,(%rax,%rcx,4)
movups %xmm1,0x10(%rax,%rcx,4)
mov    0x48(%rsp),%rax
mov    %rax,%rdx
cmp    0x18(%rsp),%rax
mov    0x68(%rsp),%r9
je     12f0 <dwt2d+0x1f0>
jmp    13f5 <dwt2d+0x2f5>
xor    %ecx,%ecx
testb  $0x1,0x80(%rsp)
01
jne    152a <dwt2d+0x42a>
jmp    1555 <dwt2d+0x455>
test   %r10d,%r10d
mov    0xa0(%rsp),%r15d
mov    0x3c(%rsp),%r8d
mov    0x70(%rsp),%rsi
mov    0x30(%rsp),%rdx
mov    0x18(%rsp),%r12
jle    15dd <dwt2d+0x4dd>
mov    %r9,%rbx
mov    0x10(%rsp),%rbp
mov    0x8(%rsp),%r13
xchg   %ax,%ax
mov    %rbp,%rdi
mov    %r13,%rsi
mov    %r15d,%ecx
mov    %r11,%r14
call   *%r11
mov    0x3c(%rsp),%r8d
mov    0x30(%rsp),%rdx
mov    %r14,%r11
mov    0x70(%rsp),%rsi
add    %rsi,%r13
add    %rbx,%rbp
add    $0xffffffffffffffff,%r12
jne    15b0 <dwt2d+0x4b0>
mov    0x98(%rsp),%rdi
lea    0x1(%rdi),%eax
shr    $0x1f,%eax
add    %edi,%eax
add    $0x1,%eax
sar    %eax
mov    %eax,0x48(%rsp)
lea    0x1(%rdx),%eax
shr    $0x1f,%eax
add    %edx,%eax
add    $0x1,%eax
sar    %eax
mov    %eax,0x40(%rsp)
test   %edx,%edx
mov    0xf8(%rsp),%r14
mov    0xf0(%rsp),%r15
jle    18e0 <dwt2d+0x7e0>
mov    %rsi,%r12
mov    0x48(%rsp),%eax
cmp    $0x2,%eax
mov    $0x1,%ecx
cmovge %eax,%ecx
mov    0x40(%rsp),%eax
cmp    $0x2,%eax
mov    $0x1,%esi
cmovge %eax,%esi
mov    %rsi,0x18(%rsp)
mov    %ecx,%eax
and    $0x7ffffff8,%eax
mov    %rax,0x30(%rsp)
add    $0xfffffffffffffff8,%rax
mov    %rax,0x90(%rsp)
mov    %rax,%rsi
shr    $0x3,%rsi
add    $0x1,%rsi
mov    %rsi,%rax
mov    %rsi,0x88(%rsp)
movabs $0x3ffffffffffffffe,%rax
ff ff 3f
and    %rax,%rsi
mov    %rsi,0x80(%rsp)
mov    %rcx,0x20(%rsp)
neg    %rcx
mov    %rcx,0x60(%rsp)
mov    $0x1,%eax
mov    %rax,0x28(%rsp)
xor    %ecx,%ecx
mov    0xa8(%rsp),%r13
mov    0xb0(%rsp),%r11
xor    %r8d,%r8d
jmp    16f4 <dwt2d+0x5f4>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r8
add    0xe8(%rsp),%r11
add    0xe0(%rsp),%r13
mov    0x78(%rsp),%rax
add    %rax,%rcx
add    %rax,0x28(%rsp)
cmp    0x18(%rsp),%r8
mov    0x98(%rsp),%rdi
je     18e0 <dwt2d+0x7e0>
test   %edi,%edi
jle    16c0 <dwt2d+0x5c0>
mov    %r8,%r9
imul   0x78(%rsp),%r9
mov    0x20(%rsp),%rdx
cmp    $0x8,%edx
setb   %al
or     0x7(%rsp),%al
jne    1750 <dwt2d+0x650>
lea    (%r9,%rdx,1),%rax
mov    0x8(%rsp),%rdi
lea    (%rdi,%r9,4),%rdx
mov    0x10(%rsp),%rbp
lea    0x0(,%rax,4),%rsi
add    %rbp,%rsi
cmp    %rsi,%rdx
jae    17f8 <dwt2d+0x6f8>
lea    (%rdi,%rax,4),%rax
lea    0x0(,%r9,4),%rdx
add    %rbp,%rdx
cmp    %rax,%rdx
jae    17f8 <dwt2d+0x6f8>
xor    %edx,%edx
mov    %rdx,%r10
testb  $0x1,0x20(%rsp)
je     178a <dwt2d+0x68a>
add    %rdx,%r9
mov    %r9,%rax
imul   0x58(%rsp),%rax
mov    0x10(%rsp),%rsi
movss  (%rsi,%rax,4),%xmm0
imul   0x50(%rsp),%r9
mov    0x8(%rsp),%rax
movss  %xmm0,(%rax,%r9,4)
mov    %rdx,%r10
or     $0x1,%r10
not    %rdx
cmp    0x60(%rsp),%rdx
je     16c0 <dwt2d+0x5c0>
mov    0x20(%rsp),%rdx
sub    %r10,%rdx
lea    (%r10,%rcx,1),%rbx
mov    %r12,%rsi
imul   %rbx,%rsi
mov    0x68(%rsp),%rax
imul   %rax,%rbx
add    0x28(%rsp),%r10
mov    %r12,%rdi
imul   %r10,%rdi
imul   %rax,%r10
mov    0x10(%rsp),%rax
mov    0x8(%rsp),%rbp
xchg   %ax,%ax
movss  (%rax,%rbx,1),%xmm0
movss  %xmm0,0x0(%rbp,%rsi,1)
movss  (%rax,%r10,1),%xmm0
movss  %xmm0,0x0(%rbp,%rdi,1)
add    %r14,%rbp
add    %r15,%rax
add    $0xfffffffffffffffe,%rdx
jne    17d0 <dwt2d+0x6d0>
jmp    16c0 <dwt2d+0x5c0>
cmpq   $0x0,0x90(%rsp)
je     18c5 <dwt2d+0x7c5>
mov    0x80(%rsp),%rdx
mov    %r13,%rsi
mov    %r11,%rdi
xor    %ebx,%ebx
mov    0xc0(%rsp),%rax
mov    0xc8(%rsp),%r10
mov    0xd8(%rsp),%rbp
mov    0xb8(%rsp),%r12
nopw   0x0(%rax,%rax,1)
movups -0x10(%rsi),%xmm0
movups (%rsi),%xmm1
movups %xmm0,-0x10(%rdi)
movups %xmm1,(%rdi)
movups -0x10(%rsi,%rbp,1),%xmm0
movups (%rsi,%rbp,1),%xmm1
movups %xmm0,-0x10(%rdi,%rax,1)
movups %xmm1,(%rdi,%rax,1)
add    $0x10,%rbx
add    %r10,%rdi
add    %r12,%rsi
add    $0xfffffffffffffffe,%rdx
jne    1840 <dwt2d+0x740>
testb  $0x1,0x88(%rsp)
01
je     18a8 <dwt2d+0x7a8>
add    %r9,%rbx
mov    %rbx,%rax
imul   0x58(%rsp),%rax
mov    0x10(%rsp),%rdx
movups (%rdx,%rax,4),%xmm0
movups 0x10(%rdx,%rax,4),%xmm1
imul   0x50(%rsp),%rbx
mov    0x8(%rsp),%rax
movups %xmm0,(%rax,%rbx,4)
movups %xmm1,0x10(%rax,%rbx,4)
mov    0x30(%rsp),%rax
mov    %rax,%rdx
cmp    0x20(%rsp),%rax
mov    0x70(%rsp),%r12
je     16c0 <dwt2d+0x5c0>
jmp    1752 <dwt2d+0x652>
xor    %ebx,%ebx
testb  $0x1,0x88(%rsp)
01
jne    187a <dwt2d+0x77a>
jmp    18a8 <dwt2d+0x7a8>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    0xa4(%rsp),%ecx
mov    0x48(%rsp),%eax
cmp    %ecx,%eax
mov    0xd0(%rsp),%r11
mov    0x68(%rsp),%r9
jl     190d <dwt2d+0x80d>
mov    %eax,%r10d
mov    0x40(%rsp),%eax
mov    %eax,%ebx
cmp    %ecx,%eax
jge    1240 <dwt2d+0x140>
add    $0x108,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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