#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*/ free(float *);
int /*<<< orphan*/ gauss3d_(float *, int *, int *, int *, float *, float *);
int /*<<< orphan*/ imgdap_(float *, int *, int *, int *, float *, int *, int *,
                           int *);
int /*<<< orphan*/ imgpad_(float *, int *, int *, int *, float *, int *, int *,
                           int *);
scalar_t__ malloc(int);
int npad_(int *, int *);
void imag2mask(int *pnx, int *pny, int *pnz, float *imag, short *mask,
               float *mmppix, float *pfhalf, float *pcrit) {
  float vmax;
  float *imgm;
  float *imgmp;
  int i;
  int dimension;
  int margin;
  int nxp, nyp, nzp;
  dimension = (*pnx) * (*pny) * (*pnz);
  margin = (int)(2.0 * 0.1874 / (mmppix[0] * (*pfhalf)));
  nxp = npad_(pnx, &margin);
  margin = (int)(2.0 * 0.1874 / (mmppix[1] * (*pfhalf)));
  nyp = npad_(pny, &margin);
  margin = (int)(2.0 * 0.1874 / (mmppix[2] * (*pfhalf)));
  nzp = npad_(pnz, &margin);
  imgmp = (float *)malloc(nxp * nyp * nzp * sizeof(float));
  imgpad_(imag, pnx, pny, pnz, imgmp, &nxp, &nyp, &nzp);
  gauss3d_(imgmp, &nxp, &nyp, &nzp, mmppix, pfhalf);
  imgm = (float *)malloc(dimension * sizeof(float));
  imgdap_(imgm, pnx, pny, pnz, imgmp, &nxp, &nyp, &nzp);
  free(imgmp);
  vmax = 0.;
  for (i = 0; i < dimension; i++) {
    if (imgm[i] > vmax)
      vmax = imgm[i];
  }
  vmax *= *pcrit;
  for (i = 0; i < dimension; i++) {
    if (imgm[i] >= vmax)
      mask[i] = 1;
    else
      mask[i] = 0;
  }
  free(imgm);
}


x86 assembly (clang, -O0):
<imag2mask>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    -0x10(%rbp),%rcx
imul   (%rcx),%eax
mov    -0x18(%rbp),%rcx
imul   (%rcx),%eax
mov    %eax,-0x50(%rbp)
mov    -0x30(%rbp),%rax
movss  (%rax),%xmm0
mov    0x10(%rbp),%rax
mulss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe45(%rip),%xmm0
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x54(%rbp)
mov    -0x8(%rbp),%rdi
lea    -0x54(%rbp),%rsi
call   1080 <npad_@plt>
mov    %eax,-0x58(%rbp)
mov    -0x30(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    0x10(%rbp),%rax
mulss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe0d(%rip),%xmm0
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x54(%rbp)
mov    -0x10(%rbp),%rdi
lea    -0x54(%rbp),%rsi
call   1080 <npad_@plt>
mov    %eax,-0x5c(%rbp)
mov    -0x30(%rbp),%rax
movss  0x8(%rax),%xmm0
mov    0x10(%rbp),%rax
mulss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdd5(%rip),%xmm0
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x54(%rbp)
mov    -0x18(%rbp),%rdi
lea    -0x54(%rbp),%rsi
call   1080 <npad_@plt>
mov    %eax,-0x60(%rbp)
mov    -0x58(%rbp),%eax
imul   -0x5c(%rbp),%eax
imul   -0x60(%rbp),%eax
cltq
shl    $0x2,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x48(%rbp)
mov    -0x20(%rbp),%rdi
mov    -0x8(%rbp),%rsi
mov    -0x10(%rbp),%rdx
mov    -0x18(%rbp),%rcx
mov    -0x48(%rbp),%r8
lea    -0x58(%rbp),%r9
lea    -0x5c(%rbp),%r10
lea    -0x60(%rbp),%rax
mov    %r10,(%rsp)
mov    %rax,0x8(%rsp)
call   1070 <imgpad_@plt>
mov    -0x48(%rbp),%rdi
mov    -0x30(%rbp),%r8
mov    0x10(%rbp),%r9
lea    -0x58(%rbp),%rsi
lea    -0x5c(%rbp),%rdx
lea    -0x60(%rbp),%rcx
call   1030 <gauss3d_@plt>
movslq -0x50(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rdi
mov    -0x8(%rbp),%rsi
mov    -0x10(%rbp),%rdx
mov    -0x18(%rbp),%rcx
mov    -0x48(%rbp),%r8
lea    -0x58(%rbp),%r9
lea    -0x5c(%rbp),%r10
lea    -0x60(%rbp),%rax
mov    %r10,(%rsp)
mov    %rax,0x8(%rsp)
call   1060 <imgdap_@plt>
mov    -0x48(%rbp),%rdi
call   1040 <free@plt>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x50(%rbp),%eax
jge    134e <imag2mask+0x1ee>
mov    -0x40(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x34(%rbp),%xmm0
jbe    133b <imag2mask+0x1db>
mov    -0x40(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x34(%rbp)
jmp    1340 <imag2mask+0x1e0>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1306 <imag2mask+0x1a6>
mov    0x18(%rbp),%rax
movss  (%rax),%xmm0
mulss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x50(%rbp),%eax
jge    13be <imag2mask+0x25e>
mov    -0x40(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x34(%rbp),%xmm0
jb     139d <imag2mask+0x23d>
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movw   $0x1,(%rax,%rcx,2)
jmp    13ab <imag2mask+0x24b>
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movw   $0x0,(%rax,%rcx,2)
jmp    13b0 <imag2mask+0x250>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1367 <imag2mask+0x207>
mov    -0x40(%rbp),%rdi
call   1040 <free@plt>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<imag2mask>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r8,%r14
mov    %rcx,0x30(%rsp)
mov    %rdx,0x10(%rsp)
mov    %rsi,%r12
mov    (%rsi),%eax
mov    %rax,%rcx
imul   (%rdi),%ecx
mov    %rdi,0x18(%rsp)
mov    0x70(%rsp),%rax
imul   (%rdx),%ecx
mov    %rcx,0x28(%rsp)
movss  (%r9),%xmm0
mov    %r9,%rbx
mulss  (%rax),%xmm0
mov    %rax,%r13
cvtss2sd %xmm0,%xmm0
movsd  0xe4d(%rip),%xmm1
divsd  %xmm0,%xmm1
cvttsd2si %xmm1,%eax
mov    %eax,(%rsp)
mov    %rsp,%r15
mov    %r15,%rsi
call   1080 <npad_@plt>
mov    %eax,%ebp
mov    %eax,0xc(%rsp)
movss  0x4(%rbx),%xmm0
mov    %rbx,0x20(%rsp)
mulss  0x0(%r13),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe15(%rip),%xmm1
divsd  %xmm0,%xmm1
cvttsd2si %xmm1,%eax
mov    %eax,(%rsp)
mov    %r12,%rdi
mov    %r15,%rsi
call   1080 <npad_@plt>
mov    %eax,%r15d
mov    %eax,0x8(%rsp)
movss  0x8(%rbx),%xmm0
mulss  0x0(%r13),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xde1(%rip),%xmm1
divsd  %xmm0,%xmm1
cvttsd2si %xmm1,%eax
mov    %eax,(%rsp)
mov    0x10(%rsp),%rbx
mov    %rbx,%rdi
mov    %rsp,%rsi
call   1080 <npad_@plt>
mov    %eax,0x4(%rsp)
imul   %eax,%ebp
imul   %r15d,%ebp
shl    $0x2,%ebp
mov    %ebp,%edi
call   1050 <malloc@plt>
mov    %rax,%rbp
lea    0x4(%rsp),%rax
lea    0x8(%rsp),%r10
lea    0xc(%rsp),%r13
mov    0x30(%rsp),%rdi
mov    0x18(%rsp),%rsi
mov    %r12,%r15
mov    %r12,%rdx
mov    %rbx,%rcx
mov    %rbp,%r8
mov    %r13,%r9
push   %rax
mov    %rax,%r12
push   %r10
mov    %r10,%rbx
call   1070 <imgpad_@plt>
add    $0x10,%rsp
mov    %rbp,%rdi
mov    %r13,%rsi
mov    %rbx,%rdx
mov    %r12,%rcx
mov    0x20(%rsp),%r8
mov    0x70(%rsp),%r9
call   1030 <gauss3d_@plt>
mov    0x28(%rsp),%r12
lea    0x0(,%r12,4),%edi
call   1050 <malloc@plt>
mov    %rax,%rbx
mov    %rax,%rdi
mov    0x18(%rsp),%rsi
mov    %r15,%rdx
mov    0x10(%rsp),%rcx
mov    %rbp,%r8
mov    %r13,%r9
lea    0x4(%rsp),%rax
push   %rax
lea    0x10(%rsp),%rax
push   %rax
call   1060 <imgdap_@plt>
add    $0x10,%rsp
mov    %rbp,%rdi
call   1040 <free@plt>
test   %r12d,%r12d
jle    1543 <imag2mask+0x3e3>
mov    %r12d,%eax
lea    -0x1(%rax),%rdx
mov    %eax,%ecx
and    $0x3,%ecx
cmp    $0x3,%rdx
jae    1312 <imag2mask+0x1b2>
xorpd  %xmm1,%xmm1
xor    %edx,%edx
jmp    1350 <imag2mask+0x1f0>
mov    %eax,%esi
and    $0xfffffffc,%esi
xorpd  %xmm1,%xmm1
xor    %edx,%edx
nopl   (%rax)
movss  (%rbx,%rdx,4),%xmm0
maxss  %xmm1,%xmm0
movss  0x4(%rbx,%rdx,4),%xmm1
maxss  %xmm0,%xmm1
movss  0x8(%rbx,%rdx,4),%xmm0
maxss  %xmm1,%xmm0
movss  0xc(%rbx,%rdx,4),%xmm1
maxss  %xmm0,%xmm1
add    $0x4,%rdx
cmp    %rdx,%rsi
jne    1320 <imag2mask+0x1c0>
movaps %xmm1,%xmm0
test   %rcx,%rcx
je     1375 <imag2mask+0x215>
lea    (%rbx,%rdx,4),%rdx
xor    %esi,%esi
xchg   %ax,%ax
movss  (%rdx,%rsi,4),%xmm0
maxss  %xmm1,%xmm0
add    $0x1,%rsi
movaps %xmm0,%xmm1
cmp    %rsi,%rcx
jne    1360 <imag2mask+0x200>
test   %r12d,%r12d
jle    1543 <imag2mask+0x3e3>
mov    0x78(%rsp),%rcx
mulss  (%rcx),%xmm0
cmp    $0x8,%r12d
jb     13a7 <imag2mask+0x247>
lea    (%rbx,%rax,4),%rcx
cmp    %r14,%rcx
jbe    144c <imag2mask+0x2ec>
lea    (%r14,%rax,2),%rcx
cmp    %rcx,%rbx
jae    144c <imag2mask+0x2ec>
xor    %ecx,%ecx
mov    %rcx,%rdx
not    %rdx
add    %rax,%rdx
mov    %rax,%rsi
and    $0x3,%rsi
je     13dd <imag2mask+0x27d>
nopl   0x0(%rax,%rax,1)
movss  (%rbx,%rcx,4),%xmm1
xor    %edi,%edi
ucomiss %xmm0,%xmm1
setae  %dil
mov    %di,(%r14,%rcx,2)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    13c0 <imag2mask+0x260>
cmp    $0x3,%rdx
jb     1543 <imag2mask+0x3e3>
nopw   0x0(%rax,%rax,1)
movss  (%rbx,%rcx,4),%xmm1
xor    %edx,%edx
ucomiss %xmm0,%xmm1
setae  %dl
mov    %dx,(%r14,%rcx,2)
movss  0x4(%rbx,%rcx,4),%xmm1
xor    %edx,%edx
ucomiss %xmm0,%xmm1
setae  %dl
mov    %dx,0x2(%r14,%rcx,2)
movss  0x8(%rbx,%rcx,4),%xmm1
xor    %edx,%edx
ucomiss %xmm0,%xmm1
setae  %dl
mov    %dx,0x4(%r14,%rcx,2)
movss  0xc(%rbx,%rcx,4),%xmm1
xor    %edx,%edx
ucomiss %xmm0,%xmm1
setae  %dl
mov    %dx,0x6(%r14,%rcx,2)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    13f0 <imag2mask+0x290>
jmp    1543 <imag2mask+0x3e3>
mov    %eax,%ecx
and    $0xfffffff8,%ecx
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
lea    -0x8(%rcx),%rdx
mov    %rdx,%rsi
shr    $0x3,%rsi
add    $0x1,%rsi
test   %rdx,%rdx
je     155a <imag2mask+0x3fa>
mov    %rsi,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %edx,%edx
movdqa 0xb8f(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movups (%rbx,%rdx,4),%xmm3
movups 0x10(%rbx,%rdx,4),%xmm4
movaps %xmm1,%xmm5
cmpleps %xmm3,%xmm5
packssdw %xmm5,%xmm5
pand   %xmm2,%xmm5
movaps %xmm1,%xmm3
cmpleps %xmm4,%xmm3
packssdw %xmm3,%xmm3
pand   %xmm2,%xmm3
punpcklqdq %xmm3,%xmm5
movdqu %xmm5,(%r14,%rdx,2)
movups 0x20(%rbx,%rdx,4),%xmm3
movups 0x30(%rbx,%rdx,4),%xmm4
movaps %xmm1,%xmm5
cmpleps %xmm3,%xmm5
packssdw %xmm5,%xmm5
pand   %xmm2,%xmm5
movaps %xmm1,%xmm3
cmpleps %xmm4,%xmm3
packssdw %xmm3,%xmm3
pand   %xmm2,%xmm3
punpcklqdq %xmm3,%xmm5
movdqu %xmm5,0x10(%r14,%rdx,2)
add    $0x10,%rdx
add    $0xfffffffffffffffe,%rdi
jne    1490 <imag2mask+0x330>
test   $0x1,%sil
je     153a <imag2mask+0x3da>
movups (%rbx,%rdx,4),%xmm2
movups 0x10(%rbx,%rdx,4),%xmm3
movaps %xmm1,%xmm4
cmpleps %xmm2,%xmm4
packssdw %xmm4,%xmm4
movdqa 0xaf0(%rip),%xmm2
pand   %xmm2,%xmm4
cmpleps %xmm3,%xmm1
packssdw %xmm1,%xmm1
pand   %xmm2,%xmm1
punpcklqdq %xmm1,%xmm4
movdqu %xmm4,(%r14,%rdx,2)
cmp    %rax,%rcx
jne    13a9 <imag2mask+0x249>
mov    %rbx,%rdi
call   1040 <free@plt>
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %edx,%edx
test   $0x1,%sil
jne    1504 <imag2mask+0x3a4>
jmp    153a <imag2mask+0x3da>

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

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