#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*/ bgbmid_free(float *);
scalar_t__ bgbmid_malloc(int);
int floor(int);
int BGBMID_ScaleSampleTrilinear(short *d, int dl, short *s, int sl) {
  float *tb, *tb1, *tb2;
  float s1, s2, s3;
  float r, r1, p1, p2, f, g;
  int v, tl;
  int i, j, k;
  tl = 1;
  while (sl > tl)
    tl <<= 1;
  tb = (float *)bgbmid_malloc(tl * sizeof(float));
  tb1 = (float *)bgbmid_malloc(tl * sizeof(float));
  r = (((float)(sl - 1)) / ((float)tl));
  for (i = 0; i < tl; i++) {
    p1 = i * r;
    j = floor(p1);
    p2 = p1 - j;
    s1 = s[j];
    s2 = s[j + 1];
    f = (s1 * (1 - p2)) + (s2 * p2);
    tb[i] = f;
  }
  while (tl > dl) {
    tl >>= 1;
    for (j = 0; j < tl; j++) {
      f = (tb[j * 2 + 0] + tb[j * 2 + 1]) / 2;
      tb1[j] = f;
    }
    tb2 = tb;
    tb = tb1;
    tb1 = tb2;
  }
  if (tl == dl) {
    for (i = 0; i < dl; i++) {
      v = tb[i];
      if (v < (-32767))
        v = -32767;
      if (v > 32767)
        v = 32767;
      d[i] = v;
    }
    return (0);
  }
  r = (((float)tl) / ((float)dl));
  r1 = (((float)sl) / ((float)dl));
  for (i = 0; i < dl; i++) {
    p1 = i * r;
    j = floor(p1);
    p2 = p1 - j;
    k = floor(i * r1 + 0.5);
    s1 = tb[j];
    s2 = tb[j + 1];
    s3 = s[k];
    v = ((s1 * (1 - p2)) + (s2 * p2)) * 0.5 + (s3 * 0.5);
    if (v < (-32767))
      v = -32767;
    if (v > 32767)
      v = 32767;
    d[i] = v;
  }
  bgbmid_free(tb);
  bgbmid_free(tb1);
  return (0);
}


x86 assembly (clang, -O0):
<BGBMID_ScaleSampleTrilinear>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %ecx,-0x24(%rbp)
movl   $0x1,-0x6c(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x6c(%rbp),%eax
jle    116a <BGBMID_ScaleSampleTrilinear+0x3a>
mov    -0x6c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    1150 <BGBMID_ScaleSampleTrilinear+0x20>
movslq -0x6c(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1030 <bgbmid_malloc@plt>
mov    %rax,-0x30(%rbp)
movslq -0x6c(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1030 <bgbmid_malloc@plt>
mov    %rax,-0x38(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
cvtsi2ss %eax,%xmm0
cvtsi2ssl -0x6c(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x50(%rbp)
movl   $0x0,-0x70(%rbp)
mov    -0x70(%rbp),%eax
cmp    -0x6c(%rbp),%eax
jge    1262 <BGBMID_ScaleSampleTrilinear+0x132>
cvtsi2ssl -0x70(%rbp),%xmm0
movss  -0x50(%rbp),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x58(%rbp)
cvttss2si -0x58(%rbp),%edi
call   1050 <floor@plt>
mov    %eax,-0x74(%rbp)
movss  -0x58(%rbp),%xmm0
cvtsi2ssl -0x74(%rbp),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x5c(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x74(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x74(%rbp),%rcx
movswl 0x2(%rax,%rcx,2),%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x44(%rbp),%xmm0
movss  0xde2(%rip),%xmm2
subss  -0x5c(%rbp),%xmm2
movss  -0x48(%rbp),%xmm1
mulss  -0x5c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x60(%rbp)
movss  -0x60(%rbp),%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)
jmp    11af <BGBMID_ScaleSampleTrilinear+0x7f>
jmp    1267 <BGBMID_ScaleSampleTrilinear+0x137>
mov    -0x6c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jle    1307 <BGBMID_ScaleSampleTrilinear+0x1d7>
mov    -0x6c(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x6c(%rbp)
movl   $0x0,-0x74(%rbp)
mov    -0x74(%rbp),%eax
cmp    -0x6c(%rbp),%eax
jge    12ea <BGBMID_ScaleSampleTrilinear+0x1ba>
mov    -0x30(%rbp),%rax
mov    -0x74(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x74(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  0xd4b(%rip),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x60(%rbp)
movss  -0x60(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x74(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    1283 <BGBMID_ScaleSampleTrilinear+0x153>
mov    -0x30(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x40(%rbp),%rax
mov    %rax,-0x38(%rbp)
jmp    1267 <BGBMID_ScaleSampleTrilinear+0x137>
mov    -0x6c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jne    138a <BGBMID_ScaleSampleTrilinear+0x25a>
movl   $0x0,-0x70(%rbp)
mov    -0x70(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    137e <BGBMID_ScaleSampleTrilinear+0x24e>
mov    -0x30(%rbp),%rax
movslq -0x70(%rbp),%rcx
cvttss2si (%rax,%rcx,4),%eax
mov    %eax,-0x68(%rbp)
cmpl   $0xffff8001,-0x68(%rbp)
jge    134a <BGBMID_ScaleSampleTrilinear+0x21a>
movl   $0xffff8001,-0x68(%rbp)
cmpl   $0x7fff,-0x68(%rbp)
jle    135e <BGBMID_ScaleSampleTrilinear+0x22e>
movl   $0x7fff,-0x68(%rbp)
mov    -0x68(%rbp),%eax
mov    %ax,%dx
mov    -0x10(%rbp),%rax
movslq -0x70(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x70(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x70(%rbp)
jmp    131a <BGBMID_ScaleSampleTrilinear+0x1ea>
movl   $0x0,-0x4(%rbp)
jmp    150f <BGBMID_ScaleSampleTrilinear+0x3df>
cvtsi2ssl -0x6c(%rbp),%xmm0
cvtsi2ssl -0x14(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x50(%rbp)
cvtsi2ssl -0x24(%rbp),%xmm0
cvtsi2ssl -0x14(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x54(%rbp)
movl   $0x0,-0x70(%rbp)
mov    -0x70(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    14f6 <BGBMID_ScaleSampleTrilinear+0x3c6>
cvtsi2ssl -0x70(%rbp),%xmm0
movss  -0x50(%rbp),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x58(%rbp)
cvttss2si -0x58(%rbp),%edi
call   1050 <floor@plt>
mov    %eax,-0x74(%rbp)
movss  -0x58(%rbp),%xmm0
cvtsi2ssl -0x74(%rbp),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x5c(%rbp)
cvtsi2ssl -0x70(%rbp),%xmm0
movss  -0x54(%rbp),%xmm1
mulss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xbf0(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1050 <floor@plt>
mov    %eax,-0x78(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x74(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x74(%rbp),%rcx
movss  0x4(%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x78(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x44(%rbp),%xmm0
movss  0xba1(%rip),%xmm2
subss  -0x5c(%rbp),%xmm2
movss  -0x48(%rbp),%xmm1
mulss  -0x5c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x4c(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xb6d(%rip),%xmm2
mulsd  %xmm2,%xmm1
movsd  0xb61(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x68(%rbp)
cmpl   $0xffff8001,-0x68(%rbp)
jge    14c2 <BGBMID_ScaleSampleTrilinear+0x392>
movl   $0xffff8001,-0x68(%rbp)
cmpl   $0x7fff,-0x68(%rbp)
jle    14d6 <BGBMID_ScaleSampleTrilinear+0x3a6>
movl   $0x7fff,-0x68(%rbp)
mov    -0x68(%rbp),%eax
mov    %ax,%dx
mov    -0x10(%rbp),%rax
movslq -0x70(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x70(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x70(%rbp)
jmp    13b7 <BGBMID_ScaleSampleTrilinear+0x287>
mov    -0x30(%rbp),%rdi
call   1040 <bgbmid_free@plt>
mov    -0x38(%rbp),%rdi
call   1040 <bgbmid_free@plt>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<BGBMID_ScaleSampleTrilinear>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rdx,%r15
mov    %esi,0x10(%rsp)
mov    %rdi,0x20(%rsp)
mov    $0x1,%eax
nop
mov    %eax,%r14d
lea    (%r14,%r14,1),%eax
cmp    %ecx,%r14d
jl     1150 <BGBMID_ScaleSampleTrilinear+0x20>
mov    %rcx,0x8(%rsp)
lea    0x0(,%r14,4),%ebx
mov    %ebx,%edi
call   1030 <bgbmid_malloc@plt>
mov    %rax,%r12
mov    %ebx,%edi
call   1030 <bgbmid_malloc@plt>
mov    %rax,%rbx
test   %r14d,%r14d
jle    121d <BGBMID_ScaleSampleTrilinear+0xed>
mov    0x8(%rsp),%rax
add    $0xffffffff,%eax
cvtsi2ss %eax,%xmm1
cvtsi2ss %r14d,%xmm0
divss  %xmm0,%xmm1
movss  %xmm1,0x4(%rsp)
mov    %r14d,%r13d
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
xorps  %xmm0,%xmm0
cvtsi2ss %ebp,%xmm0
mulss  0x4(%rsp),%xmm0
movss  %xmm0,(%rsp)
cvttss2si %xmm0,%edi
call   1050 <floor@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
cltq
movswl (%r15,%rax,2),%ecx
xorps  %xmm1,%xmm1
cvtsi2ss %ecx,%xmm1
movswl 0x2(%r15,%rax,2),%eax
xorps  %xmm2,%xmm2
cvtsi2ss %eax,%xmm2
movss  (%rsp),%xmm3
subss  %xmm0,%xmm3
movss  0xe02(%rip),%xmm0
subss  %xmm3,%xmm0
mulss  %xmm1,%xmm0
mulss  %xmm3,%xmm2
addss  %xmm0,%xmm2
movss  %xmm2,(%r12,%rbp,4)
add    $0x1,%rbp
cmp    %rbp,%r13
jne    11b0 <BGBMID_ScaleSampleTrilinear+0x80>
mov    0x10(%rsp),%esi
cmp    %esi,%r14d
jle    12e7 <BGBMID_ScaleSampleTrilinear+0x1b7>
movss  0xdd2(%rip),%xmm0
mov    %r12,%rax
jmp    124c <BGBMID_ScaleSampleTrilinear+0x11c>
nopw   0x0(%rax,%rax,1)
mov    %r12,%rax
cmp    %esi,%r14d
jle    12e7 <BGBMID_ScaleSampleTrilinear+0x1b7>
mov    %rbx,%r12
mov    %rax,%rbx
mov    %r14d,%ecx
sar    %r14d
cmp    $0x2,%ecx
jl     1240 <BGBMID_ScaleSampleTrilinear+0x110>
cmp    $0x2,%r14d
mov    $0x1,%eax
cmovge %r14d,%eax
cmp    $0x4,%ecx
jge    1280 <BGBMID_ScaleSampleTrilinear+0x150>
xor    %ecx,%ecx
jmp    12c5 <BGBMID_ScaleSampleTrilinear+0x195>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %eax,%edx
and    $0x7ffffffe,%edx
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movss  (%rbx,%rcx,8),%xmm1
addss  0x4(%rbx,%rcx,8),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%r12,%rcx,4)
movss  0x8(%rbx,%rcx,8),%xmm1
addss  0xc(%rbx,%rcx,8),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,0x4(%r12,%rcx,4)
add    $0x2,%rcx
cmp    %rcx,%rdx
jne    1290 <BGBMID_ScaleSampleTrilinear+0x160>
test   $0x1,%al
je     1240 <BGBMID_ScaleSampleTrilinear+0x110>
movss  (%rbx,%rcx,8),%xmm1
addss  0x4(%rbx,%rcx,8),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%r12,%rcx,4)
jmp    1240 <BGBMID_ScaleSampleTrilinear+0x110>
cmp    %esi,%r14d
jne    130b <BGBMID_ScaleSampleTrilinear+0x1db>
test   %esi,%esi
mov    0x20(%rsp),%rdi
jle    1445 <BGBMID_ScaleSampleTrilinear+0x315>
mov    %esi,%eax
cmp    $0x8,%esi
jae    1456 <BGBMID_ScaleSampleTrilinear+0x326>
xor    %ecx,%ecx
jmp    14ee <BGBMID_ScaleSampleTrilinear+0x3be>
test   %esi,%esi
mov    0x20(%rsp),%rax
jle    1435 <BGBMID_ScaleSampleTrilinear+0x305>
xorps  %xmm1,%xmm1
cvtsi2ss %r14d,%xmm1
xorps  %xmm0,%xmm0
cvtsi2ss %esi,%xmm0
xorps  %xmm2,%xmm2
cvtsi2ssl 0x8(%rsp),%xmm2
divss  %xmm0,%xmm1
movss  %xmm1,0x8(%rsp)
divss  %xmm0,%xmm2
movss  %xmm2,0x1c(%rsp)
mov    %esi,%ecx
mov    %rcx,0x10(%rsp)
xor    %ebp,%ebp
nopl   (%rax)
xorps  %xmm1,%xmm1
cvtsi2ss %ebp,%xmm1
movss  %xmm1,0x4(%rsp)
movss  0x8(%rsp),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,(%rsp)
cvttss2si %xmm0,%edi
mov    %rax,%r14
call   1050 <floor@plt>
mov    %eax,%r13d
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
movss  (%rsp),%xmm1
subss  %xmm0,%xmm1
movss  %xmm1,(%rsp)
movss  0x4(%rsp),%xmm0
mulss  0x1c(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xc60(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1050 <floor@plt>
movslq %r13d,%rcx
cltq
movss  0xc3e(%rip),%xmm0
movss  (%rsp),%xmm1
subss  %xmm1,%xmm0
mulss  0x4(%r12,%rcx,4),%xmm1
movswl (%r15,%rax,2),%eax
mulss  (%r12,%rcx,4),%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0xc14(%rip),%xmm2
mulsd  %xmm2,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
cmp    $0xffff8002,%eax
mov    $0xffff8001,%ecx
cmovl  %ecx,%eax
cmp    $0x7fff,%eax
mov    $0x7fff,%ecx
cmovge %ecx,%eax
mov    %ax,(%r14,%rbp,2)
mov    %r14,%rax
add    $0x1,%rbp
cmp    %rbp,0x10(%rsp)
jne    1350 <BGBMID_ScaleSampleTrilinear+0x220>
mov    %r12,%rdi
call   1040 <bgbmid_free@plt>
mov    %rbx,%rdi
call   1040 <bgbmid_free@plt>
xor    %eax,%eax
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %eax,%ecx
and    $0xfffffff8,%ecx
xor    %edx,%edx
movdqa 0xbab(%rip),%xmm0
movdqa 0xbb3(%rip),%xmm1
nopl   (%rax)
movups (%r12,%rdx,4),%xmm2
movups 0x10(%r12,%rdx,4),%xmm3
cvttps2dq %xmm2,%xmm2
cvttps2dq %xmm3,%xmm3
movdqa %xmm2,%xmm4
pcmpgtd %xmm0,%xmm4
pand   %xmm4,%xmm2
pandn  %xmm0,%xmm4
por    %xmm2,%xmm4
movdqa %xmm3,%xmm2
pcmpgtd %xmm0,%xmm2
pand   %xmm2,%xmm3
pandn  %xmm0,%xmm2
por    %xmm3,%xmm2
movdqa %xmm1,%xmm3
pcmpgtd %xmm4,%xmm3
pand   %xmm3,%xmm4
pandn  %xmm1,%xmm3
por    %xmm4,%xmm3
movdqa %xmm1,%xmm4
pcmpgtd %xmm2,%xmm4
pand   %xmm4,%xmm2
pandn  %xmm1,%xmm4
por    %xmm2,%xmm4
packssdw %xmm4,%xmm3
movdqu %xmm3,(%rdi,%rdx,2)
add    $0x8,%rdx
cmp    %rdx,%rcx
jne    1470 <BGBMID_ScaleSampleTrilinear+0x340>
cmp    %rax,%rcx
je     1445 <BGBMID_ScaleSampleTrilinear+0x315>
mov    $0xffff8001,%edx
jmp    1511 <BGBMID_ScaleSampleTrilinear+0x3e1>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %si,(%rdi,%rcx,2)
add    $0x1,%rcx
cmp    %rcx,%rax
je     1445 <BGBMID_ScaleSampleTrilinear+0x315>
cvttss2si (%r12,%rcx,4),%esi
cmp    $0xffff8002,%esi
cmovl  %edx,%esi
cmp    $0x7fff,%esi
jl     1500 <BGBMID_ScaleSampleTrilinear+0x3d0>
mov    $0x7fff,%esi
jmp    1500 <BGBMID_ScaleSampleTrilinear+0x3d0>

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

LLVM IR (-O2):
define i32 @BGBMID_ScaleSampleTrilinear(i16* nocapture writeonly %0, i32 %1, i16* nocapture %2, i32 %3) {
br label %5
5: ; preds = %5, %4
%6 = phi i32 [ 1, %4 ], [ %8, %5 ]
%7 = icmp slt i32 %6, %3
%8 = shl i32 %6, 1
br i1 %7, label %5, label %9
9: ; preds = %5
%10 = shl i32 %6, 2
%11 = tail call i64 @bgbmid_malloc(i32 %10)
%12 = inttoptr i64 %11 to float*
%13 = tail call i64 @bgbmid_malloc(i32 %10)
%14 = inttoptr i64 %13 to float*
%15 = add nsw i32 %3, -1
%16 = sitofp i32 %15 to float
%17 = sitofp i32 %6 to float
%18 = fdiv float %16, %17
%19 = icmp sgt i32 %6, 0
br i1 %19, label %20, label %22
20: ; preds = %9
%21 = zext i32 %6 to i64
br label %24
22: ; preds = %24, %9
%23 = icmp sgt i32 %6, %1
br i1 %23, label %63, label %101
24: ; preds = %20, %24
%25 = phi i64 [ 0, %20 ], [ %46, %24 ]
%26 = trunc i64 %25 to i32
%27 = sitofp i32 %26 to float
%28 = fmul float %18, %27
%29 = fptosi float %28 to i32
%30 = tail call i32 @floor(i32 %29)
%31 = sitofp i32 %30 to float
%32 = fsub float %28, %31
%33 = sext i32 %30 to i64
%34 = getelementptr inbounds i16, i16* %2, i64 %33
%35 = load i16, i16* %34
%36 = sitofp i16 %35 to float
%37 = add nsw i32 %30, 1
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds i16, i16* %2, i64 %38
%40 = load i16, i16* %39
%41 = sitofp i16 %40 to float
%42 = fsub float 1.000000e+00, %32
%43 = fmul float %32, %41
%44 = tail call float @llvm.fmuladd.f32(float %36, float %42, float %43)
%45 = getelementptr inbounds float, float* %12, i64 %25
store float %44, float* %45
%46 = add nuw nsw i64 %25, 1
%47 = icmp eq i64 %46, %21
br i1 %47, label %22, label %24
48: ; preds = %76, %69
%49 = phi i64 [ 0, %69 ], [ %98, %76 ]
%50 = icmp eq i64 %72, 0
br i1 %50, label %61, label %51
51: ; preds = %48
%52 = shl nuw nsw i64 %49, 1
%53 = getelementptr inbounds float, float* %65, i64 %52
%54 = load float, float* %53
%55 = or i64 %52, 1
%56 = getelementptr inbounds float, float* %65, i64 %55
%57 = load float, float* %56
%58 = fadd float %54, %57
%59 = fmul float %58, 5.000000e-01
%60 = getelementptr inbounds float, float* %66, i64 %49
store float %59, float* %60
br label %61
61: ; preds = %51, %48, %63
%62 = icmp sgt i32 %67, %1
br i1 %62, label %63, label %101
63: ; preds = %22, %61
%64 = phi i32 [ %67, %61 ], [ %6, %22 ]
%65 = phi float* [ %66, %61 ], [ %12, %22 ]
%66 = phi float* [ %65, %61 ], [ %14, %22 ]
%67 = ashr i32 %64, 1
%68 = icmp sgt i32 %64, 1
br i1 %68, label %69, label %61
69: ; preds = %63
%70 = call i32 @llvm.smax.i32(i32 %67, i32 1)
%71 = zext i32 %70 to i64
%72 = and i64 %71, 1
%73 = icmp slt i32 %64, 4
br i1 %73, label %48, label %74
74: ; preds = %69
%75 = and i64 %71, 2147483646
br label %76
76: ; preds = %76, %74
%77 = phi i64 [ 0, %74 ], [ %98, %76 ]
%78 = phi i64 [ 0, %74 ], [ %99, %76 ]
%79 = shl nuw nsw i64 %77, 1
%80 = getelementptr inbounds float, float* %65, i64 %79
%81 = load float, float* %80
%82 = or i64 %79, 1
%83 = getelementptr inbounds float, float* %65, i64 %82
%84 = load float, float* %83
%85 = fadd float %81, %84
%86 = fmul float %85, 5.000000e-01
%87 = getelementptr inbounds float, float* %66, i64 %77
store float %86, float* %87
%88 = or i64 %77, 1
%89 = shl nuw nsw i64 %88, 1
%90 = getelementptr inbounds float, float* %65, i64 %89
%91 = load float, float* %90
%92 = or i64 %89, 1
%93 = getelementptr inbounds float, float* %65, i64 %92
%94 = load float, float* %93
%95 = fadd float %91, %94
%96 = fmul float %95, 5.000000e-01
%97 = getelementptr inbounds float, float* %66, i64 %88
store float %96, float* %97
%98 = add nuw nsw i64 %77, 2
%99 = add nuw i64 %78, 2
%100 = icmp eq i64 %99, %75
br i1 %100, label %48, label %76
101: ; preds = %61, %22
%102 = phi float* [ %14, %22 ], [ %65, %61 ]
%103 = phi float* [ %12, %22 ], [ %66, %61 ]
%104 = phi i32 [ %6, %22 ], [ %67, %61 ]
%105 = icmp eq i32 %104, %1
br i1 %105, label %106, label %156
106: ; preds = %101
%107 = icmp sgt i32 %1, 0
br i1 %107, label %108, label %208
108: ; preds = %106
%109 = zext i32 %1 to i64
%110 = icmp ult i32 %1, 8
br i1 %110, label %141, label %111
111: ; preds = %108
%112 = and i64 %109, 4294967288
br label %113
113: ; preds = %113, %111
%114 = phi i64 [ 0, %111 ], [ %137, %113 ]
%115 = getelementptr inbounds float, float* %103, i64 %114
%116 = bitcast float* %115 to <4 x float>*
%117 = load <4 x float>, <4 x float>* %116
%118 = getelementptr inbounds float, float* %115, i64 4
%119 = bitcast float* %118 to <4 x float>*
%120 = load <4 x float>, <4 x float>* %119
%121 = fptosi <4 x float> %117 to <4 x i32>
%122 = fptosi <4 x float> %120 to <4 x i32>
%123 = icmp sgt <4 x i32> %121, <i32 -32767, i32 -32767, i32 -32767, i32 -32767>
%124 = icmp sgt <4 x i32> %122, <i32 -32767, i32 -32767, i32 -32767, i32 -32767>
%125 = select <4 x i1> %123, <4 x i32> %121, <4 x i32> <i32 -32767, i32 -32767, i32 -32767, i32 -32767>
%126 = select <4 x i1> %124, <4 x i32> %122, <4 x i32> <i32 -32767, i32 -32767, i32 -32767, i32 -32767>
%127 = icmp slt <4 x i32> %125, <i32 32767, i32 32767, i32 32767, i32 32767>
%128 = icmp slt <4 x i32> %126, <i32 32767, i32 32767, i32 32767, i32 32767>
%129 = select <4 x i1> %127, <4 x i32> %125, <4 x i32> <i32 32767, i32 32767, i32 32767, i32 32767>
%130 = select <4 x i1> %128, <4 x i32> %126, <4 x i32> <i32 32767, i32 32767, i32 32767, i32 32767>
%131 = trunc <4 x i32> %129 to <4 x i16>
%132 = trunc <4 x i32> %130 to <4 x i16>
%133 = getelementptr inbounds i16, i16* %0, i64 %114
%134 = bitcast i16* %133 to <4 x i16>*
store <4 x i16> %131, <4 x i16>* %134
%135 = getelementptr inbounds i16, i16* %133, i64 4
%136 = bitcast i16* %135 to <4 x i16>*
store <4 x i16> %132, <4 x i16>* %136
%137 = add nuw i64 %114, 8
%138 = icmp eq i64 %137, %112
br i1 %138, label %139, label %113
139: ; preds = %113
%140 = icmp eq i64 %112, %109
br i1 %140, label %208, label %141
141: ; preds = %108, %139
%142 = phi i64 [ 0, %108 ], [ %112, %139 ]
br label %143
143: ; preds = %141, %143
%144 = phi i64 [ %154, %143 ], [ %142, %141 ]
%145 = getelementptr inbounds float, float* %103, i64 %144
%146 = load float, float* %145
%147 = fptosi float %146 to i32
%148 = icmp sgt i32 %147, -32767
%149 = select i1 %148, i32 %147, i32 -32767
%150 = icmp slt i32 %149, 32767
%151 = select i1 %150, i32 %149, i32 32767
%152 = trunc i32 %151 to i16
%153 = getelementptr inbounds i16, i16* %0, i64 %144
store i16 %152, i16* %153
%154 = add nuw nsw i64 %144, 1
%155 = icmp eq i64 %154, %109
br i1 %155, label %208, label %143
156: ; preds = %101
%157 = sitofp i32 %104 to float
%158 = sitofp i32 %1 to float
%159 = fdiv float %157, %158
%160 = sitofp i32 %3 to float
%161 = fdiv float %160, %158
%162 = icmp sgt i32 %1, 0
br i1 %162, label %163, label %205
163: ; preds = %156
%164 = zext i32 %1 to i64
br label %165
165: ; preds = %163, %165
%166 = phi i64 [ 0, %163 ], [ %203, %165 ]
%167 = trunc i64 %166 to i32
%168 = sitofp i32 %167 to float
%169 = fmul float %159, %168
%170 = fptosi float %169 to i32
%171 = tail call i32 @floor(i32 %170)
%172 = sitofp i32 %171 to float
%173 = fsub float %169, %172
%174 = fmul float %161, %168
%175 = fpext float %174 to double
%176 = fadd double %175, 5.000000e-01
%177 = fptosi double %176 to i32
%178 = tail call i32 @floor(i32 %177)
%179 = sext i32 %171 to i64
%180 = getelementptr inbounds float, float* %103, i64 %179
%181 = load float, float* %180
%182 = add nsw i32 %171, 1
%183 = sext i32 %182 to i64
%184 = getelementptr inbounds float, float* %103, i64 %183
%185 = load float, float* %184
%186 = sext i32 %178 to i64
%187 = getelementptr inbounds i16, i16* %2, i64 %186
%188 = load i16, i16* %187
%189 = fsub float 1.000000e+00, %173
%190 = fmul float %173, %185
%191 = tail call float @llvm.fmuladd.f32(float %181, float %189, float %190)
%192 = fpext float %191 to double
%193 = sitofp i16 %188 to double
%194 = fmul double %193, 5.000000e-01
%195 = tail call double @llvm.fmuladd.f64(double %192, double 5.000000e-01, double %194)
%196 = fptosi double %195 to i32
%197 = icmp sgt i32 %196, -32767
%198 = select i1 %197, i32 %196, i32 -32767
%199 = icmp slt i32 %198, 32767
%200 = select i1 %199, i32 %198, i32 32767
%201 = trunc i32 %200 to i16
%202 = getelementptr inbounds i16, i16* %0, i64 %166
store i16 %201, i16* %202
%203 = add nuw nsw i64 %166, 1
%204 = icmp eq i64 %203, %164
br i1 %204, label %205, label %165
205: ; preds = %165, %156
%206 = tail call i32 @bgbmid_free(float* %103)
%207 = tail call i32 @bgbmid_free(float* %102)
br label %208
208: ; preds = %143, %139, %106, %205
ret i32 0
}