#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ inkernel_quicksort(float *, int /*<<< orphan*/, int);
scalar_t__ malloc(int);
void filter_c_3d_median(float *im, int nz, int ny, int nx, float *res, int nzr,
                        int nyr, int nxr, int w) {
  int nwin = w * w * w;
  float *win = (float *)malloc(nwin * sizeof(float));
  int edgex = w / 2;
  int edgey = w / 2;
  int edgez = w / 2;
  int mpos = nwin / 2;
  int step = ny * nx;
  int x, y, z, wx, wy, wz, ind, indy, indz, indw;
  int nwa;
  for (z = edgez; z < (nz - edgez); ++z) {
    indz = z * step;
    for (y = edgey; y < (ny - edgey); ++y) {
      ind = indz + y * ny;
      for (x = edgex; x < (nx - edgex); ++x) {
        nwa = 0;
        for (wz = 0; wz < w; ++wz) {
          indw = step * (z + wz - edgez);
          for (wy = 0; wy < w; ++wy) {
            indy = indw + ny * (y + wy - edgey);
            for (wx = 0; wx < w; ++wx) {
              win[nwa] = im[indy + x + wx - edgex];
              ++nwa;
            }
          }
        }
        inkernel_quicksort(win, 0, nwin - 1);
        res[ind + x] = win[mpos];
      }
    }
  }
}


x86 assembly (clang, -O0):
<filter_c_3d_median>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
mov    0x20(%rbp),%eax
imul   0x20(%rbp),%eax
imul   0x20(%rbp),%eax
mov    %eax,-0x28(%rbp)
movslq -0x28(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x30(%rbp)
mov    0x20(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x34(%rbp)
mov    0x20(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x38(%rbp)
mov    0x20(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x3c(%rbp)
mov    -0x28(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x40(%rbp)
mov    -0x10(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x50(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    -0x3c(%rbp),%ecx
cmp    %ecx,%eax
jge    1336 <filter_c_3d_median+0x216>
mov    -0x50(%rbp),%eax
imul   -0x44(%rbp),%eax
mov    %eax,-0x68(%rbp)
mov    -0x38(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    -0x38(%rbp),%ecx
cmp    %ecx,%eax
jge    1323 <filter_c_3d_median+0x203>
mov    -0x68(%rbp),%eax
mov    -0x4c(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x60(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x34(%rbp),%ecx
cmp    %ecx,%eax
jge    1310 <filter_c_3d_median+0x1f0>
movl   $0x0,-0x70(%rbp)
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    12d2 <filter_c_3d_median+0x1b2>
mov    -0x44(%rbp),%eax
mov    -0x50(%rbp),%ecx
add    -0x5c(%rbp),%ecx
sub    -0x3c(%rbp),%ecx
imul   %ecx,%eax
mov    %eax,-0x6c(%rbp)
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    12bf <filter_c_3d_median+0x19f>
mov    -0x6c(%rbp),%eax
mov    -0x10(%rbp),%ecx
mov    -0x4c(%rbp),%edx
add    -0x58(%rbp),%edx
sub    -0x38(%rbp),%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x64(%rbp)
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    12ac <filter_c_3d_median+0x18c>
mov    -0x8(%rbp),%rax
mov    -0x64(%rbp),%ecx
add    -0x48(%rbp),%ecx
add    -0x54(%rbp),%ecx
sub    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x70(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x70(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x70(%rbp)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1264 <filter_c_3d_median+0x144>
jmp    12b1 <filter_c_3d_median+0x191>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    123a <filter_c_3d_median+0x11a>
jmp    12c4 <filter_c_3d_median+0x1a4>
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    1215 <filter_c_3d_median+0xf5>
mov    -0x30(%rbp),%rdi
mov    -0x28(%rbp),%edx
sub    $0x1,%edx
xor    %esi,%esi
call   1030 <inkernel_quicksort@plt>
mov    -0x30(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x60(%rbp),%ecx
add    -0x48(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    11f6 <filter_c_3d_median+0xd6>
jmp    1315 <filter_c_3d_median+0x1f5>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11d0 <filter_c_3d_median+0xb0>
jmp    1328 <filter_c_3d_median+0x208>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    11af <filter_c_3d_median+0x8f>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<filter_c_3d_median>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x178,%rsp
mov    %r8,0xf0(%rsp)
mov    %ecx,0xc(%rsp)
mov    %edx,%r12d
mov    %esi,%ebx
mov    %rdi,%rbp
mov    0x1c0(%rsp),%r13d
mov    %r13d,%eax
imul   %eax,%eax
imul   %r13d,%eax
mov    %rax,%r14
lea    0x0(,%rax,4),%edi
call   1040 <malloc@plt>
mov    %rax,%r15
mov    %r13d,%eax
shr    $0x1f,%eax
add    %r13d,%eax
sar    %eax
sub    %eax,%ebx
mov    %ebx,0x4c(%rsp)
mov    %eax,0x8(%rsp)
cmp    %ebx,%eax
jge    17b5 <filter_c_3d_median+0x695>
mov    0xc(%rsp),%ebx
mov    %ebx,%edx
imul   %r12d,%edx
mov    %r14d,%eax
shr    $0x1f,%eax
add    %r14d,%eax
sar    %eax
mov    %r12d,%esi
mov    0x8(%rsp),%ecx
sub    %ecx,%esi
mov    %esi,0x20(%rsp)
sub    %ecx,%ebx
mov    %ebx,0xc(%rsp)
mov    0x1c0(%rsp),%ebx
add    $0xffffffff,%r14d
mov    %r14,0xe8(%rsp)
cltq
mov    %rax,0xe0(%rsp)
movslq %ecx,%rax
movslq %r12d,%rcx
mov    %edx,0x54(%rsp)
movslq %edx,%rdx
mov    %ebx,%r13d
lea    -0x1(%r13),%rsi
mov    %rsi,0x160(%rsp)
lea    (%r15,%r13,4),%rsi
mov    %rsi,0x150(%rsp)
lea    0x0(,%r13,4),%rsi
add    %rbp,%rsi
mov    %rsi,0x158(%rsp)
mov    %r13d,%esi
and    $0xfffffff8,%esi
mov    %rsi,0x148(%rsp)
add    $0xfffffffffffffff8,%rsi
mov    %rsi,0x140(%rsp)
mov    %rsi,%rdi
shr    $0x3,%rdi
add    $0x1,%rdi
mov    %rdi,%rsi
mov    %rdi,0x138(%rsp)
and    $0xfffffffffffffffe,%rdi
mov    %rdi,0x130(%rsp)
mov    %r13,%rsi
neg    %rsi
mov    %rsi,0x168(%rsp)
lea    0x30(%r15),%rsi
mov    %rsi,0x128(%rsp)
mov    %rdx,0xc0(%rsp)
mov    %edx,%edx
mov    %rdx,0xb8(%rsp)
mov    %rcx,0xc8(%rsp)
mov    %ecx,%r14d
xor    %ecx,%ecx
mov    %rcx,0x78(%rsp)
mov    $0x8,%ecx
mov    %rcx,0x70(%rsp)
mov    $0x1,%ecx
mov    %rcx,0x68(%rsp)
movl   $0x0,0x1c(%rsp)
mov    %rax,0x80(%rsp)
mov    %rax,0x28(%rsp)
mov    %r12d,0x60(%rsp)
mov    %r14,0x10(%rsp)
jmp    12f4 <filter_c_3d_median+0x1d4>
cs nopw 0x0(%rax,%rax,1)
nop
mov    0x28(%rsp),%rcx
add    $0x1,%rcx
addl   $0x1,0x1c(%rsp)
mov    0xb8(%rsp),%rax
add    %rax,0x70(%rsp)
add    %rax,0x78(%rsp)
add    %rax,0x68(%rsp)
mov    %rcx,0x28(%rsp)
cmp    %ecx,0x4c(%rsp)
je     17b5 <filter_c_3d_median+0x695>
mov    0x20(%rsp),%eax
cmp    %eax,0x8(%rsp)
jge    12c0 <filter_c_3d_median+0x1a0>
mov    0xc0(%rsp),%rax
mov    %eax,%ecx
imul   0x1c(%rsp),%ecx
mov    %ecx,0x50(%rsp)
mov    0x28(%rsp),%rcx
mov    %rcx,%rdx
imul   %rax,%rdx
mov    %rdx,0xd0(%rsp)
sub    0x8(%rsp),%ecx
mov    %rcx,0x100(%rsp)
movl   $0x0,0x24(%rsp)
mov    0x68(%rsp),%rax
mov    %rax,0x98(%rsp)
mov    0x78(%rsp),%rax
mov    %rax,0x90(%rsp)
mov    0x70(%rsp),%rax
mov    %rax,0x88(%rsp)
mov    0x80(%rsp),%rax
mov    %rax,0x30(%rsp)
jmp    13a5 <filter_c_3d_median+0x285>
nop
mov    0x30(%rsp),%rax
add    $0x1,%rax
addl   $0x1,0x24(%rsp)
add    %r14,0x88(%rsp)
add    %r14,0x90(%rsp)
add    %r14,0x98(%rsp)
mov    %rax,0x30(%rsp)
cmp    %eax,0x20(%rsp)
je     12c0 <filter_c_3d_median+0x1a0>
mov    0x8(%rsp),%eax
cmp    0xc(%rsp),%eax
jge    1370 <filter_c_3d_median+0x250>
mov    0x24(%rsp),%eax
mov    %eax,%ecx
mov    0xc8(%rsp),%rax
imul   %eax,%ecx
add    0x50(%rsp),%ecx
mov    %rcx,0xd8(%rsp)
mov    0x30(%rsp),%rcx
mov    %rcx,%rdx
imul   %rax,%rdx
add    0xd0(%rsp),%rdx
mov    %rdx,0xf8(%rsp)
sub    0x8(%rsp),%ecx
mov    %rcx,0x170(%rsp)
xor    %eax,%eax
mov    %rax,0x40(%rsp)
mov    0x98(%rsp),%rax
mov    %rax,0xb0(%rsp)
mov    0x90(%rsp),%rax
mov    %rax,0xa8(%rsp)
mov    0x88(%rsp),%rax
mov    %rax,0xa0(%rsp)
mov    0x80(%rsp),%rax
mov    %rax,0x38(%rsp)
jmp    14b8 <filter_c_3d_median+0x398>
nopw   0x0(%rax,%rax,1)
mov    %r15,%rdi
xor    %esi,%esi
mov    0xe8(%rsp),%rdx
call   1030 <inkernel_quicksort@plt>
mov    0xe0(%rsp),%rax
movss  (%r15,%rax,4),%xmm0
mov    0xf8(%rsp),%rax
mov    0x38(%rsp),%rdx
add    %rdx,%rax
mov    0xf0(%rsp),%rcx
movss  %xmm0,(%rcx,%rax,4)
add    $0x1,%rdx
mov    0x40(%rsp),%rax
add    $0x1,%eax
mov    %rax,0x40(%rsp)
addq   $0x1,0xa0(%rsp)
addq   $0x1,0xa8(%rsp)
addq   $0x1,0xb0(%rsp)
mov    %rdx,0x38(%rsp)
cmp    %edx,0xc(%rsp)
je     1370 <filter_c_3d_median+0x250>
test   %ebx,%ebx
jle    1440 <filter_c_3d_median+0x320>
mov    0xd8(%rsp),%rax
mov    0x40(%rsp),%rcx
add    %ecx,%eax
mov    %eax,0x5c(%rsp)
mov    0x38(%rsp),%rax
sub    0x8(%rsp),%eax
mov    %eax,0x58(%rsp)
xor    %edx,%edx
mov    0xb0(%rsp),%r10
mov    0xa8(%rsp),%r11
mov    0xa0(%rsp),%rsi
xor    %edi,%edi
jmp    153c <filter_c_3d_median+0x41c>
nopw   0x0(%rax,%rax,1)
mov    0x108(%rsp),%rdi
add    $0x1,%edi
mov    0xb8(%rsp),%rax
mov    0x110(%rsp),%rsi
add    %rax,%rsi
mov    0x118(%rsp),%r11
add    %rax,%r11
mov    0x120(%rsp),%r10
add    %rax,%r10
cmp    %ebx,%edi
je     1440 <filter_c_3d_median+0x320>
mov    0x54(%rsp),%eax
mov    %eax,%ecx
imul   %edi,%ecx
add    0x5c(%rsp),%ecx
mov    %ecx,0x64(%rsp)
mov    0x100(%rsp),%rcx
mov    %rdi,0x108(%rsp)
add    %ecx,%edi
imul   %eax,%edi
add    0x58(%rsp),%edi
mov    %r10,0x120(%rsp)
mov    %r11,0x118(%rsp)
mov    %rsi,0x110(%rsp)
mov    %rsi,%rax
xor    %r9d,%r9d
jmp    15a6 <filter_c_3d_median+0x486>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r9d
add    %r14,%rax
add    %r14,%r11
add    %r14,%r10
cmp    %ebx,%r9d
je     1500 <filter_c_3d_median+0x3e0>
mov    0x170(%rsp),%rcx
lea    (%rcx,%r9,1),%r8d
imul   %r12d,%r8d
add    %edi,%r8d
movslq %edx,%rdx
cmp    $0x8,%ebx
jae    15d0 <filter_c_3d_median+0x4b0>
xor    %esi,%esi
jmp    165d <filter_c_3d_median+0x53d>
nopl   0x0(%rax,%rax,1)
mov    %edi,%r14d
mov    %r9d,%ecx
mov    %r12d,%edi
imul   %r12d,%ecx
add    0x64(%rsp),%ecx
mov    0x160(%rsp),%r12
lea    (%rcx,%r12,1),%esi
cmp    %ecx,%esi
jl     1650 <filter_c_3d_median+0x530>
shr    $0x20,%r12
jne    1650 <filter_c_3d_median+0x530>
movslq %ecx,%rcx
mov    0x158(%rsp),%rsi
lea    (%rsi,%rcx,4),%rsi
lea    (%r15,%rdx,4),%rbx
cmp    %rsi,%rbx
jae    16d9 <filter_c_3d_median+0x5b9>
lea    0x0(,%rcx,4),%rcx
add    %rbp,%rcx
mov    0x150(%rsp),%rsi
lea    (%rsi,%rdx,4),%rsi
cmp    %rsi,%rcx
jae    16d9 <filter_c_3d_median+0x5b9>
mov    %edi,%r12d
xor    %esi,%esi
mov    0x1c0(%rsp),%ebx
jmp    1655 <filter_c_3d_median+0x535>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xor    %esi,%esi
mov    %edi,%r12d
mov    %r14d,%edi
mov    0x10(%rsp),%r14
test   $0x1,%r13b
jne    1670 <filter_c_3d_median+0x550>
mov    %rsi,%rcx
jmp    168d <filter_c_3d_median+0x56d>
nopl   0x0(%rax,%rax,1)
add    %esi,%r8d
movslq %r8d,%rcx
movss  0x0(%rbp,%rcx,4),%xmm0
movss  %xmm0,(%r15,%rdx,4)
add    $0x1,%rdx
mov    %rsi,%rcx
or     $0x1,%rcx
not    %rsi
cmp    0x168(%rsp),%rsi
je     1590 <filter_c_3d_median+0x470>
xchg   %ax,%ax
lea    (%r11,%rcx,1),%esi
movslq %esi,%rsi
movss  0x0(%rbp,%rsi,4),%xmm0
movss  %xmm0,(%r15,%rdx,4)
lea    (%r10,%rcx,1),%esi
movslq %esi,%rsi
movss  0x0(%rbp,%rsi,4),%xmm0
movss  %xmm0,0x4(%r15,%rdx,4)
add    $0x2,%rdx
add    $0x2,%rcx
cmp    %rcx,%r13
jne    16a0 <filter_c_3d_median+0x580>
jmp    1590 <filter_c_3d_median+0x470>
cmpq   $0x0,0x140(%rsp)
je     179f <filter_c_3d_median+0x67f>
mov    0x128(%rsp),%rcx
lea    (%rcx,%rdx,4),%rsi
mov    0x130(%rsp),%r12
xor    %ecx,%ecx
mov    %r14d,%edi
mov    0x10(%rsp),%r14
cs nopw 0x0(%rax,%rax,1)
lea    (%r11,%rcx,1),%ebx
movslq %ebx,%rbx
movups 0x0(%rbp,%rbx,4),%xmm0
movups 0x10(%rbp,%rbx,4),%xmm1
movups %xmm0,-0x30(%rsi,%rcx,4)
movups %xmm1,-0x20(%rsi,%rcx,4)
lea    (%rax,%rcx,1),%ebx
movslq %ebx,%rbx
movups 0x0(%rbp,%rbx,4),%xmm0
movups 0x10(%rbp,%rbx,4),%xmm1
movups %xmm0,-0x10(%rsi,%rcx,4)
movups %xmm1,(%rsi,%rcx,4)
add    $0x10,%rcx
add    $0xfffffffffffffffe,%r12
jne    1710 <filter_c_3d_median+0x5f0>
testb  $0x1,0x138(%rsp)
01
je     1777 <filter_c_3d_median+0x657>
lea    (%rdx,%rcx,1),%rsi
add    %r8d,%ecx
movslq %ecx,%rcx
movups 0x0(%rbp,%rcx,4),%xmm0
movups 0x10(%rbp,%rcx,4),%xmm1
movups %xmm0,(%r15,%rsi,4)
movups %xmm1,0x10(%r15,%rsi,4)
mov    0x148(%rsp),%rcx
add    %rcx,%rdx
mov    %rcx,%rsi
cmp    %r13,%rcx
mov    0x60(%rsp),%r12d
mov    0x1c0(%rsp),%ebx
je     1590 <filter_c_3d_median+0x470>
jmp    165d <filter_c_3d_median+0x53d>
xor    %ecx,%ecx
mov    %r14d,%edi
mov    0x10(%rsp),%r14
testb  $0x1,0x138(%rsp)
01
jne    1758 <filter_c_3d_median+0x638>
jmp    1777 <filter_c_3d_median+0x657>
add    $0x178,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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