#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 */
unsigned char UCHAR_MAX;
int /*<<< orphan*/ p3dGetElapsedTime_min();
int /*<<< orphan*/ p3dGetElapsedTime_sec();
int /*<<< orphan*/ p3dResetStartTime();
int p3dRidlerThresholding_16(unsigned short *in_im, unsigned char *out_im,
                             const int dimx, const int dimy, const int dimz,
                             unsigned short *thresh,
                             int (*wr_log)(const char *, ...),
                             int (*wr_progress)(const int, ...)) {
  double tt, tb, to, t2;
  int t;
  long N, no, nb;
  int ct;
  if (wr_log != NULL) {
    p3dResetStartTime();
    wr_log("Pore3D - Thresholding image according to Ridler's iterative "
           "method...");
  }
  N = dimx * dimy * dimz;
  tb = 0.0;
  to = 0.0;
  no = 0;
  for (ct = 0; ct < dimx * dimy * dimz; ct++)
    to = to + (in_im[ct]);
  tt = (to / (double)N);
  while (N) {
    no = 0;
    nb = 0;
    tb = 0.0;
    to = 0.0;
    for (ct = 0; ct < dimx * dimy * dimz; ct++) {
      if ((double)(in_im[ct]) >= tt) {
        to = to + (double)(in_im[ct]);
        no++;
      } else {
        tb = tb + (double)(in_im[ct]);
        nb++;
      }
    }
    if (no == 0)
      no = 1;
    if (nb == 0)
      nb = 1;
    t2 = (tb / (double)nb + to / (double)no) / 2.0;
    if (t2 == tt)
      N = 0;
    tt = t2;
  }
  t = (int)tt;
  *thresh = (unsigned short)t;
#pragma omp parallel for
  for (ct = 0; ct < (dimx * dimy * dimz); ct++)
    out_im[ct] = (in_im[ct] > (*thresh)) ? UCHAR_MAX : 0;
  if (wr_log != NULL) {
    wr_log("\tDetermined threshold: %d.", *thresh);
    wr_log("Pore3D - Image thresholded successfully in %dm%0.3fs.",
           p3dGetElapsedTime_min(), p3dGetElapsedTime_sec());
  }
  return 2;
}


x86 assembly (clang, -O0):
<p3dRidlerThresholding_16>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9,-0x28(%rbp)
cmpq   $0x0,0x10(%rbp)
je     1177 <p3dRidlerThresholding_16+0x47>
mov    $0x0,%al
call   1050 <p3dResetStartTime@plt>
lea    0xe96(%rip),%rdi
mov    $0x0,%al
call   *0x10(%rbp)
mov    -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
imul   -0x1c(%rbp),%eax
cltq
mov    %rax,-0x58(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movq   $0x0,-0x60(%rbp)
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
mov    -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jge    11e9 <p3dRidlerThresholding_16+0xb9>
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movzwl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    11a7 <p3dRidlerThresholding_16+0x77>
movsd  -0x40(%rbp),%xmm0
cvtsi2sdq -0x58(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
cmpq   $0x0,-0x58(%rbp)
je     1352 <p3dRidlerThresholding_16+0x222>
movq   $0x0,-0x60(%rbp)
movq   $0x0,-0x68(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
mov    -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jge    12cc <p3dRidlerThresholding_16+0x19c>
mov    -0x8(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movzwl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm0
ucomisd -0x30(%rbp),%xmm0
jb     128f <p3dRidlerThresholding_16+0x15f>
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movzwl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x60(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x60(%rbp)
jmp    12b9 <p3dRidlerThresholding_16+0x189>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movzwl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x68(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x68(%rbp)
jmp    12be <p3dRidlerThresholding_16+0x18e>
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    122f <p3dRidlerThresholding_16+0xff>
cmpq   $0x0,-0x60(%rbp)
jne    12df <p3dRidlerThresholding_16+0x1af>
movq   $0x1,-0x60(%rbp)
cmpq   $0x0,-0x68(%rbp)
jne    12f2 <p3dRidlerThresholding_16+0x1c2>
movq   $0x1,-0x68(%rbp)
movsd  -0x38(%rbp),%xmm0
cvtsi2sdq -0x68(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
cvtsi2sdq -0x60(%rbp),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  0xce4(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
ucomisd -0x30(%rbp),%xmm0
jne    1343 <p3dRidlerThresholding_16+0x213>
jp     1343 <p3dRidlerThresholding_16+0x213>
movq   $0x0,-0x58(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    11fd <p3dRidlerThresholding_16+0xcd>
cvttsd2si -0x30(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    %ax,%cx
mov    -0x28(%rbp),%rax
mov    %cx,(%rax)
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
mov    -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jge    13d9 <p3dRidlerThresholding_16+0x2a9>
mov    -0x8(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movzwl (%rax,%rcx,2),%eax
mov    -0x28(%rbp),%rcx
movzwl (%rcx),%ecx
cmp    %ecx,%eax
jle    13b1 <p3dRidlerThresholding_16+0x281>
mov    0x2c4a(%rip),%rax
movzbl (%rax),%eax
mov    %eax,-0x70(%rbp)
jmp    13bb <p3dRidlerThresholding_16+0x28b>
xor    %eax,%eax
mov    %eax,-0x70(%rbp)
jmp    13bb <p3dRidlerThresholding_16+0x28b>
mov    -0x70(%rbp),%eax
mov    %al,%dl
mov    -0x10(%rbp),%rax
movslq -0x6c(%rbp),%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    136e <p3dRidlerThresholding_16+0x23e>
cmpq   $0x0,0x10(%rbp)
je     1427 <p3dRidlerThresholding_16+0x2f7>
mov    0x10(%rbp),%rcx
mov    -0x28(%rbp),%rax
movzwl (%rax),%esi
lea    0xc58(%rip),%rdi
mov    $0x0,%al
call   *%rcx
mov    0x10(%rbp),%rax
mov    %rax,-0x78(%rbp)
mov    $0x0,%al
call   1040 <p3dGetElapsedTime_min@plt>
mov    %eax,-0x7c(%rbp)
mov    $0x0,%al
call   1030 <p3dGetElapsedTime_sec@plt>
mov    -0x7c(%rbp),%esi
mov    -0x78(%rbp),%rcx
mov    %eax,%edx
lea    0xc46(%rip),%rdi
mov    $0x0,%al
call   *%rcx
mov    $0x2,%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<p3dRidlerThresholding_16>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %r8d,%r12d
mov    %ecx,%ebp
mov    %edx,%r13d
mov    %rsi,%r15
mov    %rdi,%rbx
mov    0x40(%rsp),%r14
test   %r14,%r14
je     1174 <p3dRidlerThresholding_16+0x44>
xor    %eax,%eax
mov    %r15,(%rsp)
mov    %r9,%r15
call   1050 <p3dResetStartTime@plt>
lea    0xeb8(%rip),%rdi
xor    %eax,%eax
call   *%r14
mov    %r15,%r9
mov    (%rsp),%r15
imul   %r13d,%ebp
imul   %r12d,%ebp
mov    %ebp,%r8d
test   %ebp,%ebp
jle    119b <p3dRidlerThresholding_16+0x6b>
lea    -0x1(%r8),%rcx
mov    %r8d,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    11a4 <p3dRidlerThresholding_16+0x74>
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
jmp    11f8 <p3dRidlerThresholding_16+0xc8>
xorpd  %xmm0,%xmm0
jmp    1228 <p3dRidlerThresholding_16+0xf8>
mov    %r8d,%edx
and    $0xfffffffc,%edx
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
movzwl (%rbx,%rcx,2),%esi
xorps  %xmm1,%xmm1
cvtsi2sd %esi,%xmm1
movzwl 0x2(%rbx,%rcx,2),%esi
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
addsd  %xmm0,%xmm1
movzwl 0x4(%rbx,%rcx,2),%esi
xorps  %xmm3,%xmm3
cvtsi2sd %esi,%xmm3
addsd  %xmm1,%xmm2
movzwl 0x6(%rbx,%rcx,2),%esi
xorps  %xmm0,%xmm0
cvtsi2sd %esi,%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    11b0 <p3dRidlerThresholding_16+0x80>
test   %rax,%rax
je     1228 <p3dRidlerThresholding_16+0xf8>
lea    (%rbx,%rcx,2),%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movzwl (%rcx,%rdx,2),%esi
xorps  %xmm1,%xmm1
cvtsi2sd %esi,%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1210 <p3dRidlerThresholding_16+0xe0>
xorps  %xmm1,%xmm1
cvtsi2sd %ebp,%xmm1
divsd  %xmm1,%xmm0
test   %ebp,%ebp
sete   %r10b
mov    %r8d,%edx
and    $0xfffffffe,%edx
mov    $0x1,%r11d
movsd  0xdb3(%rip),%xmm1
movapd %xmm0,%xmm2
jmp    12a2 <p3dRidlerThresholding_16+0x172>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xor    %edi,%edi
xorpd  %xmm4,%xmm4
mov    $0x1,%eax
test   %rdi,%rdi
cmove  %r11,%rdi
xorps  %xmm5,%xmm5
cvtsi2sd %rdi,%xmm5
divsd  %xmm5,%xmm3
xorps  %xmm5,%xmm5
cvtsi2sd %rax,%xmm5
divsd  %xmm5,%xmm4
addsd  %xmm3,%xmm4
mulsd  %xmm1,%xmm4
ucomisd %xmm2,%xmm4
movapd %xmm4,%xmm2
jne    12a2 <p3dRidlerThresholding_16+0x172>
jnp    1360 <p3dRidlerThresholding_16+0x230>
test   $0x1,%r10b
jne    136c <p3dRidlerThresholding_16+0x23c>
xorpd  %xmm3,%xmm3
test   %ebp,%ebp
jle    1260 <p3dRidlerThresholding_16+0x130>
xor    %esi,%esi
xor    %edi,%edi
xor    %eax,%eax
xorpd  %xmm4,%xmm4
cmp    $0x1,%ebp
jne    1301 <p3dRidlerThresholding_16+0x1d1>
test   $0x1,%r8b
je     12e2 <p3dRidlerThresholding_16+0x1b2>
movzwl (%rbx,%rsi,2),%ecx
xorps  %xmm5,%xmm5
cvtsi2sd %ecx,%xmm5
ucomisd %xmm2,%xmm5
jae    1344 <p3dRidlerThresholding_16+0x214>
addsd  %xmm5,%xmm3
add    $0x1,%rdi
test   %rax,%rax
jne    126b <p3dRidlerThresholding_16+0x13b>
jmp    1266 <p3dRidlerThresholding_16+0x136>
nopl   0x0(%rax)
addsd  %xmm5,%xmm3
add    $0x1,%rdi
add    $0x2,%rsi
cmp    %rsi,%rdx
je     12c3 <p3dRidlerThresholding_16+0x193>
movzwl (%rbx,%rsi,2),%ecx
xorps  %xmm5,%xmm5
cvtsi2sd %ecx,%xmm5
ucomisd %xmm2,%xmm5
jae    1320 <p3dRidlerThresholding_16+0x1f0>
addsd  %xmm5,%xmm3
add    $0x1,%rdi
jmp    1328 <p3dRidlerThresholding_16+0x1f8>
nopl   0x0(%rax)
addsd  %xmm5,%xmm4
add    $0x1,%rax
movzwl 0x2(%rbx,%rsi,2),%ecx
xorps  %xmm5,%xmm5
cvtsi2sd %ecx,%xmm5
ucomisd %xmm2,%xmm5
jb     12f0 <p3dRidlerThresholding_16+0x1c0>
addsd  %xmm5,%xmm4
add    $0x1,%rax
jmp    12f8 <p3dRidlerThresholding_16+0x1c8>
addsd  %xmm5,%xmm4
add    $0x1,%rax
test   %rax,%rax
jne    126b <p3dRidlerThresholding_16+0x13b>
jmp    1266 <p3dRidlerThresholding_16+0x136>
nopw   0x0(%rax,%rax,1)
mov    $0x1,%r10b
movapd %xmm4,%xmm0
jmp    124d <p3dRidlerThresholding_16+0x11d>
cvttsd2si %xmm0,%eax
mov    %ax,(%r9)
test   %ebp,%ebp
jle    14c3 <p3dRidlerThresholding_16+0x393>
xor    %edx,%edx
cmp    %ax,(%rbx)
mov    0x2c68(%rip),%rcx
movzbl (%rcx),%eax
cmovbe %edx,%eax
mov    %al,(%r15)
cmp    $0x1,%ebp
je     14c3 <p3dRidlerThresholding_16+0x393>
lea    -0x1(%r8),%r10
mov    $0x1,%eax
cmp    $0x10,%r10
jb     149a <p3dRidlerThresholding_16+0x36a>
lea    0x1(%r15),%rdx
lea    (%r15,%r8,1),%rbp
lea    0x2(%r9),%rsi
lea    0x2(%rbx),%rdi
cmp    %rsi,%rdx
setb   %r11b
mov    %r9,%r14
cmp    %r9,%rbp
seta   %r12b
lea    (%rbx,%r8,2),%rsi
cmp    %rsi,%rdx
setb   %dl
cmp    %rdi,%rbp
seta   %dil
cmp    %rcx,%r15
setb   %sil
cmp    %rcx,%rbp
seta   %r9b
test   %r12b,%r11b
jne    1497 <p3dRidlerThresholding_16+0x367>
and    %dil,%dl
jne    1497 <p3dRidlerThresholding_16+0x367>
and    %r9b,%sil
mov    %r14,%r9
jne    149a <p3dRidlerThresholding_16+0x36a>
mov    %r10,%rsi
and    $0xfffffffffffffff0,%rsi
lea    0x1(%rsi),%rax
movzwl (%r9),%edx
movd   %edx,%xmm0
pshuflw $0x0,%xmm0,%xmm0
pshufd $0x0,%xmm0,%xmm0
movzbl (%rcx),%edx
movd   %edx,%xmm1
punpcklbw %xmm1,%xmm1
pshuflw $0x0,%xmm1,%xmm1
xor    %edi,%edi
movdqa 0xbce(%rip),%xmm2
pxor   %xmm2,%xmm0
cs nopw 0x0(%rax,%rax,1)
movdqu 0x2(%rbx,%rdi,2),%xmm3
movdqu 0x12(%rbx,%rdi,2),%xmm4
pxor   %xmm2,%xmm3
pcmpgtw %xmm0,%xmm3
packsswb %xmm3,%xmm3
pxor   %xmm2,%xmm4
pcmpgtw %xmm0,%xmm4
packsswb %xmm4,%xmm4
pand   %xmm1,%xmm3
pand   %xmm1,%xmm4
punpcklqdq %xmm4,%xmm3
movdqu %xmm3,0x1(%r15,%rdi,1)
add    $0x10,%rdi
cmp    %rdi,%rsi
jne    1450 <p3dRidlerThresholding_16+0x320>
cmp    %rsi,%r10
jne    149a <p3dRidlerThresholding_16+0x36a>
jmp    14c3 <p3dRidlerThresholding_16+0x393>
mov    %r14,%r9
mov    %rax,%rdx
not    %rdx
test   $0x1,%r8b
jne    14be <p3dRidlerThresholding_16+0x38e>
movzwl (%rbx,%rax,2),%esi
xor    %edi,%edi
cmp    (%r9),%si
movzbl (%rcx),%esi
cmovbe %edi,%esi
mov    %sil,(%r15,%rax,1)
add    $0x1,%rax
add    %r8,%rdx
jne    151e <p3dRidlerThresholding_16+0x3ee>
mov    0x40(%rsp),%rbp
test   %rbp,%rbp
je     14fb <p3dRidlerThresholding_16+0x3cb>
movzwl (%r9),%esi
lea    0xb8e(%rip),%rdi
xor    %eax,%eax
call   *%rbp
xor    %eax,%eax
call   1040 <p3dGetElapsedTime_min@plt>
mov    %eax,%ebx
xor    %eax,%eax
call   1030 <p3dGetElapsedTime_sec@plt>
lea    0xb8e(%rip),%rdi
mov    %ebx,%esi
mov    %eax,%edx
xor    %eax,%eax
call   *%rbp
mov    $0x2,%eax
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
nop
mov    %dl,0x1(%r15,%rax,1)
add    $0x2,%rax
cmp    %rax,%r8
je     14c3 <p3dRidlerThresholding_16+0x393>
movzwl (%rbx,%rax,2),%esi
xor    %edx,%edx
cmp    (%r9),%si
ja     1530 <p3dRidlerThresholding_16+0x400>
xor    %esi,%esi
jmp    1533 <p3dRidlerThresholding_16+0x403>
xchg   %ax,%ax
movzbl (%rcx),%esi
mov    %sil,(%r15,%rax,1)
movzwl 0x2(%rbx,%rax,2),%esi
cmp    (%r9),%si
jbe    1510 <p3dRidlerThresholding_16+0x3e0>
movzbl (%rcx),%edx
jmp    1510 <p3dRidlerThresholding_16+0x3e0>

LLVM IR (-O0):
@.str = [70 x i8] c"Pore3D - Thresholding image according to Ridler's iterative method...\00"
@UCHAR_MAX = global i8 0
@.str.1 = [27 x i8] c"\09Determined threshold: %d.\00"
@.str.2 = [54 x i8] c"Pore3D - Image thresholded successfully in %dm%0.3fs.\00"
define i32 @p3dRidlerThresholding_16(i16* %0, i8* %1, i32 %2, i32 %3, i32 %4, i16* %5, i32 (i8*, ...)* %6, i32 (i32, ...)* %7) {
%9 = alloca i16*
%10 = alloca i8*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i16*
%15 = alloca i32 (i8*, ...)*
%16 = alloca i32 (i32, ...)*
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca i32
%22 = alloca i64
%23 = alloca i64
%24 = alloca i64
%25 = alloca i32
store i16* %0, i16** %9
store i8* %1, i8** %10
store i32 %2, i32* %11
store i32 %3, i32* %12
store i32 %4, i32* %13
store i16* %5, i16** %14
store i32 (i8*, ...)* %6, i32 (i8*, ...)** %15
store i32 (i32, ...)* %7, i32 (i32, ...)** %16
%26 = load i32 (i8*, ...)*, i32 (i8*, ...)** %15
%27 = icmp ne i32 (i8*, ...)* %26, null
br i1 %27, label %28, label %32
28: ; preds = %8
%29 = call i32 (...) @p3dResetStartTime()
%30 = load i32 (i8*, ...)*, i32 (i8*, ...)** %15
%31 = call i32 (i8*, ...) %30(i8* getelementptr inbounds ([70 x i8], [70 x i8]* @.str, i64 0, i64 0))
br label %32
32: ; preds = %28, %8
%33 = load i32, i32* %11
%34 = load i32, i32* %12
%35 = mul nsw i32 %33, %34
%36 = load i32, i32* %13
%37 = mul nsw i32 %35, %36
%38 = sext i32 %37 to i64
store i64 %38, i64* %22
store double 0.000000e+00, double* %18
store double 0.000000e+00, double* %19
store i64 0, i64* %23
store i32 0, i32* %25
br label %39
39: ; preds = %57, %32
%40 = load i32, i32* %25
%41 = load i32, i32* %11
%42 = load i32, i32* %12
%43 = mul nsw i32 %41, %42
%44 = load i32, i32* %13
%45 = mul nsw i32 %43, %44
%46 = icmp slt i32 %40, %45
br i1 %46, label %47, label %60
47: ; preds = %39
%48 = load double, double* %19
%49 = load i16*, i16** %9
%50 = load i32, i32* %25
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds i16, i16* %49, i64 %51
%53 = load i16, i16* %52
%54 = zext i16 %53 to i32
%55 = sitofp i32 %54 to double
%56 = fadd double %48, %55
store double %56, double* %19
br label %57
57: ; preds = %47
%58 = load i32, i32* %25
%59 = add nsw i32 %58, 1
store i32 %59, i32* %25
br label %39
60: ; preds = %39
%61 = load double, double* %19
%62 = load i64, i64* %22
%63 = sitofp i64 %62 to double
%64 = fdiv double %61, %63
store double %64, double* %17
br label %65
65: ; preds = %135, %60
%66 = load i64, i64* %22
%67 = icmp ne i64 %66, 0
br i1 %67, label %68, label %137
68: ; preds = %65
store i64 0, i64* %23
store i64 0, i64* %24
store double 0.000000e+00, double* %18
store double 0.000000e+00, double* %19
store i32 0, i32* %25
br label %69
69: ; preds = %109, %68
%70 = load i32, i32* %25
%71 = load i32, i32* %11
%72 = load i32, i32* %12
%73 = mul nsw i32 %71, %72
%74 = load i32, i32* %13
%75 = mul nsw i32 %73, %74
%76 = icmp slt i32 %70, %75
br i1 %76, label %77, label %112
77: ; preds = %69
%78 = load i16*, i16** %9
%79 = load i32, i32* %25
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds i16, i16* %78, i64 %80
%82 = load i16, i16* %81
%83 = uitofp i16 %82 to double
%84 = load double, double* %17
%85 = fcmp oge double %83, %84
br i1 %85, label %86, label %97
86: ; preds = %77
%87 = load double, double* %19
%88 = load i16*, i16** %9
%89 = load i32, i32* %25
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds i16, i16* %88, i64 %90
%92 = load i16, i16* %91
%93 = uitofp i16 %92 to double
%94 = fadd double %87, %93
store double %94, double* %19
%95 = load i64, i64* %23
%96 = add nsw i64 %95, 1
store i64 %96, i64* %23
br label %108
97: ; preds = %77
%98 = load double, double* %18
%99 = load i16*, i16** %9
%100 = load i32, i32* %25
%101 = sext i32 %100 to i64
%102 = getelementptr inbounds i16, i16* %99, i64 %101
%103 = load i16, i16* %102
%104 = uitofp i16 %103 to double
%105 = fadd double %98, %104
store double %105, double* %18
%106 = load i64, i64* %24
%107 = add nsw i64 %106, 1
store i64 %107, i64* %24
br label %108
108: ; preds = %97, %86
br label %109
109: ; preds = %108
%110 = load i32, i32* %25
%111 = add nsw i32 %110, 1
store i32 %111, i32* %25
br label %69
112: ; preds = %69
%113 = load i64, i64* %23
%114 = icmp eq i64 %113, 0
br i1 %114, label %115, label %116
115: ; preds = %112
store i64 1, i64* %23
br label %116
116: ; preds = %115, %112
%117 = load i64, i64* %24
%118 = icmp eq i64 %117, 0
br i1 %118, label %119, label %120
119: ; preds = %116
store i64 1, i64* %24
br label %120
120: ; preds = %119, %116
%121 = load double, double* %18
%122 = load i64, i64* %24
%123 = sitofp i64 %122 to double
%124 = fdiv double %121, %123
%125 = load double, double* %19
%126 = load i64, i64* %23
%127 = sitofp i64 %126 to double
%128 = fdiv double %125, %127
%129 = fadd double %124, %128
%130 = fdiv double %129, 2.000000e+00
store double %130, double* %20
%131 = load double, double* %20
%132 = load double, double* %17
%133 = fcmp oeq double %131, %132
br i1 %133, label %134, label %135
134: ; preds = %120
store i64 0, i64* %22
br label %135
135: ; preds = %134, %120
%136 = load double, double* %20
store double %136, double* %17
br label %65
137: ; preds = %65
%138 = load double, double* %17
%139 = fptosi double %138 to i32
store i32 %139, i32* %21
%140 = load i32, i32* %21
%141 = trunc i32 %140 to i16
%142 = load i16*, i16** %14
store i16 %141, i16* %142
store i32 0, i32* %25
br label %143
143: ; preds = %173, %137
%144 = load i32, i32* %25
%145 = load i32, i32* %11
%146 = load i32, i32* %12
%147 = mul nsw i32 %145, %146
%148 = load i32, i32* %13
%149 = mul nsw i32 %147, %148
%150 = icmp slt i32 %144, %149
br i1 %150, label %151, label %176
151: ; preds = %143
%152 = load i16*, i16** %9
%153 = load i32, i32* %25
%154 = sext i32 %153 to i64
%155 = getelementptr inbounds i16, i16* %152, i64 %154
%156 = load i16, i16* %155
%157 = zext i16 %156 to i32
%158 = load i16*, i16** %14
%159 = load i16, i16* %158
%160 = zext i16 %159 to i32
%161 = icmp sgt i32 %157, %160
br i1 %161, label %162, label %165
162: ; preds = %151
%163 = load i8, i8* @UCHAR_MAX
%164 = zext i8 %163 to i32
br label %166
165: ; preds = %151
br label %166
166: ; preds = %165, %162
%167 = phi i32 [ %164, %162 ], [ 0, %165 ]
%168 = trunc i32 %167 to i8
%169 = load i8*, i8** %10
%170 = load i32, i32* %25
%171 = sext i32 %170 to i64
%172 = getelementptr inbounds i8, i8* %169, i64 %171
store i8 %168, i8* %172
br label %173
173: ; preds = %166
%174 = load i32, i32* %25
%175 = add nsw i32 %174, 1
store i32 %175, i32* %25
br label %143
176: ; preds = %143
%177 = load i32 (i8*, ...)*, i32 (i8*, ...)** %15
%178 = icmp ne i32 (i8*, ...)* %177, null
br i1 %178, label %179, label %189
179: ; preds = %176
%180 = load i32 (i8*, ...)*, i32 (i8*, ...)** %15
%181 = load i16*, i16** %14
%182 = load i16, i16* %181
%183 = zext i16 %182 to i32
%184 = call i32 (i8*, ...) %180(i8* getelementptr inbounds ([27 x i8], [27 x i8]* @.str.1, i64 0, i64 0), i32 %183)
%185 = load i32 (i8*, ...)*, i32 (i8*, ...)** %15
%186 = call i32 (...) @p3dGetElapsedTime_min()
%187 = call i32 (...) @p3dGetElapsedTime_sec()
%188 = call i32 (i8*, ...) %185(i8* getelementptr inbounds ([54 x i8], [54 x i8]* @.str.2, i64 0, i64 0), i32 %186, i32 %187)
br label %189
189: ; preds = %179, %176
ret i32 2
}

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