#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 */
float abs(float);
int /*<<< orphan*/ printf(char *, int, int);
double sin(double);
void volume_c_resampling_lanczos2(float *org, int noz, int noy, int nox,
                                  float *trg, int nz, int ny, int nx) {
  float scalez = noz / (float)nz;
  float scaley = noy / (float)ny;
  float scalex = nox / (float)nx;
  int stepo = nox * noy;
  int stept = nx * ny;
  int x, y, z;
  int xi, yi, zi;
  float xt, yt, zt;
  int u, v, w;
  int wz, wy, wx;
  float p, q, r;
  float dx, dy, dz;
  for (z = 0; z < nz; ++z) {
    printf("slice z = %i / %i\n", z + 1, nz);
    zt = z * scalez;
    zi = (int)zt;
    for (y = 0; y < ny; ++y) {
      yt = y * scaley;
      yi = (int)yt;
      for (x = 0; x < nx; ++x) {
        xt = x * scalex;
        xi = (int)xt;
        r = 0;
        for (wz = -1; wz < 3; ++wz) {
          w = zi + wz;
          if (w >= noz) {
            continue;
          }
          if (w < 0) {
            continue;
          }
          dz = zt - w;
          if (abs(dz) > 2.0f) {
            dz = 2.0f;
          }
          q = 0;
          for (wy = -1; wy < 3; ++wy) {
            v = yi + wy;
            if (v >= noy) {
              continue;
            }
            if (v < 0) {
              continue;
            }
            dy = yt - v;
            if (abs(dy) > 2.0f) {
              dy = 2.0f;
            }
            p = 0;
            for (wx = -1; wx < 3; ++wx) {
              u = xi + wx;
              if (u >= nox) {
                continue;
              }
              if (u < 0) {
                continue;
              }
              dx = xt - u;
              if (abs(dx) > 2.0f) {
                dx = 2.0f;
              }
              p = p +
                  org[w * stepo + v * nox + u] *
                      ((dx) == (0)
                           ? 1
                           : sin(3.141592653589793238462643383279 * (dx)) /
                                 (3.141592653589793238462643383279 * (dx))) *
                      ((dx * 0.5f) == (0)
                           ? 1
                           : sin(3.141592653589793238462643383279 *
                                 (dx * 0.5f)) /
                                 (3.141592653589793238462643383279 *
                                  (dx * 0.5f)));
            }
            q = q +
                p *
                    ((dy) == (0)
                         ? 1
                         : sin(3.141592653589793238462643383279 * (dy)) /
                               (3.141592653589793238462643383279 * (dy))) *
                    ((dy * 0.5f) == (0)
                         ? 1
                         : sin(3.141592653589793238462643383279 * (dy * 0.5f)) /
                               (3.141592653589793238462643383279 *
                                (dy * 0.5f)));
          }
          r = r +
              q *
                  ((dz) == (0)
                       ? 1
                       : sin(3.141592653589793238462643383279 * (dz)) /
                             (3.141592653589793238462643383279 * (dz))) *
                  ((dz * 0.5f) == (0)
                       ? 1
                       : sin(3.141592653589793238462643383279 * (dz * 0.5f)) /
                             (3.141592653589793238462643383279 * (dz * 0.5f)));
        }
        trg[z * stept + y * nx + x] = r;
      }
    }
  }
}


x86 assembly (clang, -O0):
<volume_c_resampling_lanczos2>:
push   %rbp
mov    %rsp,%rbp
sub    $0x110,%rsp
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)
cvtsi2ssl -0xc(%rbp),%xmm0
cvtsi2ssl -0x24(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
cvtsi2ssl -0x10(%rbp),%xmm0
cvtsi2ssl 0x10(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
cvtsi2ssl -0x14(%rbp),%xmm0
cvtsi2ssl 0x18(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x14(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    0x18(%rbp),%eax
imul   0x10(%rbp),%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1881 <volume_c_resampling_lanczos2+0x751>
mov    -0x44(%rbp),%esi
add    $0x1,%esi
mov    -0x24(%rbp),%edx
lea    0xe52(%rip),%rdi
call   1040 <printf@plt>
cvtsi2ssl -0x44(%rbp),%xmm0
mulss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x5c(%rbp)
cvttss2si -0x5c(%rbp),%eax
mov    %eax,-0x50(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    186e <volume_c_resampling_lanczos2+0x73e>
cvtsi2ssl -0x40(%rbp),%xmm0
mulss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x58(%rbp)
cvttss2si -0x58(%rbp),%eax
mov    %eax,-0x4c(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    185b <volume_c_resampling_lanczos2+0x72b>
cvtsi2ssl -0x3c(%rbp),%xmm0
mulss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x54(%rbp)
cvttss2si -0x54(%rbp),%eax
mov    %eax,-0x48(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x80(%rbp)
movl   $0xffffffff,-0x6c(%rbp)
cmpl   $0x3,-0x6c(%rbp)
jge    1829 <volume_c_resampling_lanczos2+0x6f9>
mov    -0x50(%rbp),%eax
add    -0x6c(%rbp),%eax
mov    %eax,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cmp    -0xc(%rbp),%eax
jl     1269 <volume_c_resampling_lanczos2+0x139>
jmp    181b <volume_c_resampling_lanczos2+0x6eb>
cmpl   $0x0,-0x68(%rbp)
jge    1278 <volume_c_resampling_lanczos2+0x148>
jmp    181b <volume_c_resampling_lanczos2+0x6eb>
movss  -0x5c(%rbp),%xmm0
cvtsi2ssl -0x68(%rbp),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x8c(%rbp)
ff
movss  -0x8c(%rbp),%xmm0
ff
call   1030 <abs@plt>
movss  0xd5d(%rip),%xmm1
ucomiss %xmm1,%xmm0
jbe    12bc <volume_c_resampling_lanczos2+0x18c>
movss  0xd4c(%rip),%xmm0
movss  %xmm0,-0x8c(%rbp)
ff
xorps  %xmm0,%xmm0
movss  %xmm0,-0x7c(%rbp)
movl   $0xffffffff,-0x70(%rbp)
cmpl   $0x3,-0x70(%rbp)
jge    16bc <volume_c_resampling_lanczos2+0x58c>
mov    -0x4c(%rbp),%eax
add    -0x70(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x10(%rbp),%eax
jl     12ef <volume_c_resampling_lanczos2+0x1bf>
jmp    16ae <volume_c_resampling_lanczos2+0x57e>
cmpl   $0x0,-0x64(%rbp)
jge    12fe <volume_c_resampling_lanczos2+0x1ce>
jmp    16ae <volume_c_resampling_lanczos2+0x57e>
movss  -0x58(%rbp),%xmm0
cvtsi2ssl -0x64(%rbp),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x88(%rbp)
ff
movss  -0x88(%rbp),%xmm0
ff
call   1030 <abs@plt>
movss  0xcd7(%rip),%xmm1
ucomiss %xmm1,%xmm0
jbe    1342 <volume_c_resampling_lanczos2+0x212>
movss  0xcc6(%rip),%xmm0
movss  %xmm0,-0x88(%rbp)
ff
xorps  %xmm0,%xmm0
movss  %xmm0,-0x78(%rbp)
movl   $0xffffffff,-0x74(%rbp)
cmpl   $0x3,-0x74(%rbp)
jge    154f <volume_c_resampling_lanczos2+0x41f>
mov    -0x48(%rbp),%eax
add    -0x74(%rbp),%eax
mov    %eax,-0x60(%rbp)
mov    -0x60(%rbp),%eax
cmp    -0x14(%rbp),%eax
jl     1375 <volume_c_resampling_lanczos2+0x245>
jmp    1541 <volume_c_resampling_lanczos2+0x411>
cmpl   $0x0,-0x60(%rbp)
jge    1384 <volume_c_resampling_lanczos2+0x254>
jmp    1541 <volume_c_resampling_lanczos2+0x411>
movss  -0x54(%rbp),%xmm0
cvtsi2ssl -0x60(%rbp),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x84(%rbp)
ff
movss  -0x84(%rbp),%xmm0
ff
call   1030 <abs@plt>
movss  0xc51(%rip),%xmm1
ucomiss %xmm1,%xmm0
jbe    13c8 <volume_c_resampling_lanczos2+0x298>
movss  0xc40(%rip),%xmm0
movss  %xmm0,-0x84(%rbp)
ff
movss  -0x78(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0x8(%rbp),%rax
mov    -0x68(%rbp),%ecx
imul   -0x34(%rbp),%ecx
mov    -0x64(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
add    -0x60(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movss  -0x84(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    1433 <volume_c_resampling_lanczos2+0x303>
jp     1433 <volume_c_resampling_lanczos2+0x303>
movsd  0xbea(%rip),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
jmp    1474 <volume_c_resampling_lanczos2+0x344>
movss  -0x84(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm1
movsd  0xbc1(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <sin@plt>
movss  -0x84(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm2
movsd  0xba4(%rip),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0xa8(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movss  0xb6c(%rip),%xmm0
mulss  -0x84(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    14c7 <volume_c_resampling_lanczos2+0x397>
jp     14c7 <volume_c_resampling_lanczos2+0x397>
movsd  0xb56(%rip),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
jmp    1518 <volume_c_resampling_lanczos2+0x3e8>
movss  0xb35(%rip),%xmm0
mulss  -0x84(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm1
movsd  0xb25(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <sin@plt>
movss  0xb10(%rip),%xmm1
mulss  -0x84(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm2
movsd  0xb00(%rip),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0xa0(%rbp),%xmm1
ff
movsd  -0xb0(%rbp),%xmm0
ff
movsd  -0xb8(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x78(%rbp)
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    1351 <volume_c_resampling_lanczos2+0x221>
movss  -0x7c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movss  -0x78(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movss  -0x88(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    15a0 <volume_c_resampling_lanczos2+0x470>
jp     15a0 <volume_c_resampling_lanczos2+0x470>
movsd  0xa7d(%rip),%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
jmp    15e1 <volume_c_resampling_lanczos2+0x4b1>
movss  -0x88(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm1
movsd  0xa54(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <sin@plt>
movss  -0x88(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm2
movsd  0xa37(%rip),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
movsd  -0xd0(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movss  0x9ff(%rip),%xmm0
mulss  -0x88(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    1634 <volume_c_resampling_lanczos2+0x504>
jp     1634 <volume_c_resampling_lanczos2+0x504>
movsd  0x9e9(%rip),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
jmp    1685 <volume_c_resampling_lanczos2+0x555>
movss  0x9c8(%rip),%xmm0
mulss  -0x88(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm1
movsd  0x9b8(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <sin@plt>
movss  0x9a3(%rip),%xmm1
mulss  -0x88(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm2
movsd  0x993(%rip),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  -0xc8(%rbp),%xmm1
ff
movsd  -0xd8(%rbp),%xmm0
ff
movsd  -0xe0(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x7c(%rbp)
mov    -0x70(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x70(%rbp)
jmp    12cb <volume_c_resampling_lanczos2+0x19b>
movss  -0x80(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movss  -0x7c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
movss  -0x8c(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    170d <volume_c_resampling_lanczos2+0x5dd>
jp     170d <volume_c_resampling_lanczos2+0x5dd>
movsd  0x910(%rip),%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
jmp    174e <volume_c_resampling_lanczos2+0x61e>
movss  -0x8c(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm1
movsd  0x8e7(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <sin@plt>
movss  -0x8c(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm2
movsd  0x8ca(%rip),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
movsd  -0xe8(%rbp),%xmm0
ff
movsd  -0xf8(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
movss  0x892(%rip),%xmm0
mulss  -0x8c(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    17a1 <volume_c_resampling_lanczos2+0x671>
jp     17a1 <volume_c_resampling_lanczos2+0x671>
movsd  0x87c(%rip),%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
jmp    17f2 <volume_c_resampling_lanczos2+0x6c2>
movss  0x85b(%rip),%xmm0
mulss  -0x8c(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm1
movsd  0x84b(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1050 <sin@plt>
movss  0x836(%rip),%xmm1
mulss  -0x8c(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm2
movsd  0x826(%rip),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
movsd  -0xf0(%rbp),%xmm1
ff
movsd  -0x100(%rbp),%xmm0
ff
movsd  -0x108(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x80(%rbp)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    1245 <volume_c_resampling_lanczos2+0x115>
movss  -0x80(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x44(%rbp),%ecx
imul   -0x38(%rbp),%ecx
mov    -0x40(%rbp),%edx
imul   0x18(%rbp),%edx
add    %edx,%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1213 <volume_c_resampling_lanczos2+0xe3>
jmp    1860 <volume_c_resampling_lanczos2+0x730>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11e9 <volume_c_resampling_lanczos2+0xb9>
jmp    1873 <volume_c_resampling_lanczos2+0x743>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    11aa <volume_c_resampling_lanczos2+0x7a>
add    $0x110,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<volume_c_resampling_lanczos2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x138,%rsp
mov    %r8,0xc8(%rsp)
mov    %ecx,0x2c(%rsp)
mov    %rdi,0x108(%rsp)
test   %r9d,%r9d
jle    17c8 <volume_c_resampling_lanczos2+0x698>
mov    %edx,%ebx
mov    %esi,%r14d
cvtsi2ss %esi,%xmm2
cvtsi2ss %r9d,%xmm0
mov    0x178(%rsp),%ecx
cvtsi2ss %edx,%xmm3
cvtsi2ssl 0x170(%rsp),%xmm1
divss  %xmm0,%xmm2
movss  %xmm2,0x70(%rsp)
mov    0x2c(%rsp),%eax
xorps  %xmm2,%xmm2
cvtsi2ss %eax,%xmm2
divss  %xmm1,%xmm3
movss  %xmm3,0x78(%rsp)
xorps  %xmm0,%xmm0
cvtsi2ss %ecx,%xmm0
divss  %xmm0,%xmm2
movss  %xmm2,0x7c(%rsp)
movslq %eax,%r15
mov    %ecx,%eax
mov    %rax,0xc0(%rsp)
mov    %r15d,%eax
mov    %rax,0xf8(%rsp)
mov    %r15d,%r12d
imul   %edx,%r12d
xor    %eax,%eax
mov    %r9d,0x74(%rsp)
mov    %edx,0x90(%rsp)
mov    %esi,0x80(%rsp)
mov    %r12,0xe8(%rsp)
jmp    1206 <volume_c_resampling_lanczos2+0xd6>
xchg   %ax,%ax
mov    0x74(%rsp),%r9d
mov    0xa8(%rsp),%rax
cmp    %r9d,%eax
je     17c8 <volume_c_resampling_lanczos2+0x698>
mov    %rax,%rbp
lea    0x1(%rax),%esi
lea    0xe05(%rip),%rdi
mov    %rsi,0xa8(%rsp)
mov    %r9d,%edx
call   1040 <printf@plt>
cmpl   $0x0,0x170(%rsp)
jle    11f0 <volume_c_resampling_lanczos2+0xc0>
xorps  %xmm0,%xmm0
cvtsi2ss %ebp,%xmm0
mulss  0x70(%rsp),%xmm0
movss  %xmm0,0x84(%rsp)
cvttss2si %xmm0,%ecx
imul   0x170(%rsp),%ebp
mov    %rbp,0xb8(%rsp)
mov    %rcx,%rax
mov    %rcx,0xf0(%rsp)
lea    -0x1(%rcx),%eax
imul   %ebx,%eax
add    $0xffffffff,%eax
mov    %rax,0xb0(%rsp)
xor    %eax,%eax
mov    %rax,0x58(%rsp)
jmp    129a <volume_c_resampling_lanczos2+0x16a>
nopl   0x0(%rax)
mov    0x58(%rsp),%rax
add    $0x1,%eax
mov    %rax,0x58(%rsp)
cmp    0x170(%rsp),%eax
je     11f0 <volume_c_resampling_lanczos2+0xc0>
cmpl   $0x0,0x178(%rsp)
jle    1280 <volume_c_resampling_lanczos2+0x150>
mov    0x58(%rsp),%rcx
xorps  %xmm0,%xmm0
cvtsi2ss %ecx,%xmm0
mulss  0x78(%rsp),%xmm0
movss  %xmm0,0x94(%rsp)
cvttss2si %xmm0,%edx
mov    0xb8(%rsp),%rax
add    %ecx,%eax
imul   0x178(%rsp),%eax
mov    %rax,0xd8(%rsp)
mov    0xb0(%rsp),%rax
mov    %rdx,0x100(%rsp)
add    %edx,%eax
imul   0x2c(%rsp),%eax
mov    %rax,0xd0(%rsp)
xor    %edx,%edx
jmp    1333 <volume_c_resampling_lanczos2+0x203>
mov    0xd8(%rsp),%rax
mov    0xe0(%rsp),%rdx
add    %edx,%eax
cltq
mov    0xc8(%rsp),%rcx
movss  %xmm0,(%rcx,%rax,4)
add    $0x1,%rdx
cmp    0xc0(%rsp),%rdx
je     1280 <volume_c_resampling_lanczos2+0x150>
mov    %rdx,0xe0(%rsp)
xorps  %xmm0,%xmm0
cvtsi2ss %edx,%xmm0
mulss  0x7c(%rsp),%xmm0
movss  %xmm0,0x9c(%rsp)
cvttss2si %xmm0,%r13d
movslq %r13d,%rbp
mov    0xd0(%rsp),%rax
add    %ebp,%eax
mov    %rax,0xa0(%rsp)
mov    $0xffffffff,%ecx
xorps  %xmm0,%xmm0
jmp    13be <volume_c_resampling_lanczos2+0x28e>
cs nopw 0x0(%rax,%rax,1)
nop
movss  0x8c(%rsp),%xmm1
cvtss2sd %xmm1,%xmm1
movss  0xc(%rsp),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm1,%xmm3
xorps  %xmm0,%xmm0
cvtsd2ss %xmm3,%xmm0
add    $0x1,%ecx
add    %r12,0xa0(%rsp)
cmp    $0x3,%ecx
je     1300 <volume_c_resampling_lanczos2+0x1d0>
mov    0xf0(%rsp),%rax
add    %ecx,%eax
cmp    %r14d,%eax
jge    13aa <volume_c_resampling_lanczos2+0x27a>
test   %eax,%eax
js     13aa <volume_c_resampling_lanczos2+0x27a>
movss  %xmm0,0x8c(%rsp)
mov    %rcx,0x60(%rsp)
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
movss  0x84(%rsp),%xmm0
subss  %xmm1,%xmm0
movaps %xmm0,0x110(%rsp)
call   1030 <abs@plt>
movss  %xmm0,0x88(%rsp)
mov    $0xffffffff,%r12d
xorps  %xmm0,%xmm0
movss  %xmm0,0xc(%rsp)
mov    0xa0(%rsp),%r14
jmp    146d <volume_c_resampling_lanczos2+0x33d>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  0xc(%rsp),%xmm1
cvtss2sd %xmm1,%xmm1
cvtss2sd %xmm2,%xmm2
mulsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm1,%xmm3
xorps  %xmm0,%xmm0
cvtsd2ss %xmm3,%xmm0
movss  %xmm0,0xc(%rsp)
add    $0x1,%r12d
add    0xf8(%rsp),%r14
cmp    $0x3,%r12d
je     16f0 <volume_c_resampling_lanczos2+0x5c0>
mov    0x100(%rsp),%rax
add    %r12d,%eax
cmp    %ebx,%eax
jge    1457 <volume_c_resampling_lanczos2+0x327>
test   %eax,%eax
js     1457 <volume_c_resampling_lanczos2+0x327>
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
movss  0x94(%rsp),%xmm0
subss  %xmm1,%xmm0
movaps %xmm0,0x120(%rsp)
call   1030 <abs@plt>
movss  %xmm0,0x98(%rsp)
xorps  %xmm2,%xmm2
mov    $0xffffffffffffffff,%rbx
jmp    14f2 <volume_c_resampling_lanczos2+0x3c2>
cs nopw 0x0(%rax,%rax,1)
movss  0x30(%rsp),%xmm1
cvtss2sd %xmm1,%xmm1
xorps  %xmm2,%xmm2
cvtss2sd %xmm3,%xmm2
mulsd  %xmm2,%xmm4
mulsd  %xmm0,%xmm4
addsd  %xmm1,%xmm4
xorps  %xmm2,%xmm2
cvtsd2ss %xmm4,%xmm2
add    $0x1,%rbx
cmp    $0x3,%rbx
je     1610 <volume_c_resampling_lanczos2+0x4e0>
lea    (%rbx,%rbp,1),%rax
cmp    %r15,%rax
jge    14e4 <volume_c_resampling_lanczos2+0x3b4>
test   %rax,%rax
js     14e4 <volume_c_resampling_lanczos2+0x3b4>
movss  %xmm2,0x30(%rsp)
lea    (%rbx,%r13,1),%eax
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
movss  0x9c(%rsp),%xmm0
subss  %xmm1,%xmm0
movaps %xmm0,0x10(%rsp)
call   1030 <abs@plt>
movss  0xad0(%rip),%xmm1
movaps %xmm1,%xmm2
cmpltss %xmm0,%xmm2
movaps %xmm2,%xmm0
andnps 0x10(%rsp),%xmm0
andps  %xmm1,%xmm2
orps   %xmm0,%xmm2
lea    (%r14,%rbx,1),%eax
cltq
mov    0x108(%rsp),%rcx
movss  (%rcx,%rax,4),%xmm3
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm2
movsd  0xaa1(%rip),%xmm0
movapd %xmm0,%xmm4
jne    156f <volume_c_resampling_lanczos2+0x43f>
jnp    15b4 <volume_c_resampling_lanczos2+0x484>
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
mulsd  0xa92(%rip),%xmm0
movsd  %xmm0,0x68(%rsp)
movaps %xmm2,0x40(%rsp)
movss  %xmm3,0x10(%rsp)
call   1050 <sin@plt>
movss  0x10(%rsp),%xmm3
movaps 0x40(%rsp),%xmm2
xorps  %xmm1,%xmm1
movapd %xmm0,%xmm4
movsd  0xa5a(%rip),%xmm0
divsd  0x68(%rsp),%xmm4
mulss  0xa48(%rip),%xmm2
ucomiss %xmm1,%xmm2
jne    15c7 <volume_c_resampling_lanczos2+0x497>
jnp    14c0 <volume_c_resampling_lanczos2+0x390>
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
mulsd  0xa3a(%rip),%xmm0
movsd  %xmm0,0x68(%rsp)
movss  %xmm3,0x10(%rsp)
movsd  %xmm4,0x40(%rsp)
call   1050 <sin@plt>
movsd  0x40(%rsp),%xmm4
movss  0x10(%rsp),%xmm3
divsd  0x68(%rsp),%xmm0
jmp    14c0 <volume_c_resampling_lanczos2+0x390>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  0x9e8(%rip),%xmm1
movaps %xmm1,%xmm4
cmpltss 0x98(%rsp),%xmm4
movaps %xmm4,%xmm0
andps  %xmm1,%xmm0
andnps 0x120(%rsp),%xmm4
orps   %xmm0,%xmm4
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm4
movsd  0x9c4(%rip),%xmm0
movapd %xmm0,%xmm3
jne    164c <volume_c_resampling_lanczos2+0x51c>
jnp    1691 <volume_c_resampling_lanczos2+0x561>
xorps  %xmm0,%xmm0
cvtss2sd %xmm4,%xmm0
mulsd  0x9b5(%rip),%xmm0
movsd  %xmm0,0x40(%rsp)
movss  %xmm2,0x30(%rsp)
movaps %xmm4,0x10(%rsp)
call   1050 <sin@plt>
movaps 0x10(%rsp),%xmm4
movss  0x30(%rsp),%xmm2
xorps  %xmm1,%xmm1
movapd %xmm0,%xmm3
movsd  0x97d(%rip),%xmm0
divsd  0x40(%rsp),%xmm3
mulss  0x96b(%rip),%xmm4
ucomiss %xmm1,%xmm4
mov    0x90(%rsp),%ebx
jne    16ab <volume_c_resampling_lanczos2+0x57b>
jnp    1430 <volume_c_resampling_lanczos2+0x300>
xorps  %xmm0,%xmm0
cvtss2sd %xmm4,%xmm0
mulsd  0x956(%rip),%xmm0
movsd  %xmm0,0x40(%rsp)
movss  %xmm2,0x30(%rsp)
movsd  %xmm3,0x10(%rsp)
call   1050 <sin@plt>
movsd  0x10(%rsp),%xmm3
movss  0x30(%rsp),%xmm2
divsd  0x40(%rsp),%xmm0
jmp    1430 <volume_c_resampling_lanczos2+0x300>
nopl   0x0(%rax,%rax,1)
movss  0x908(%rip),%xmm1
movaps %xmm1,%xmm2
cmpltss 0x88(%rsp),%xmm2
movaps %xmm2,%xmm0
andps  %xmm1,%xmm0
andnps 0x110(%rsp),%xmm2
orps   %xmm0,%xmm2
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm2
movsd  0x8e4(%rip),%xmm0
movapd %xmm0,%xmm3
mov    0x60(%rsp),%rcx
jne    1731 <volume_c_resampling_lanczos2+0x601>
jnp    176f <volume_c_resampling_lanczos2+0x63f>
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
mulsd  0x8d0(%rip),%xmm0
movsd  %xmm0,0x10(%rsp)
movaps %xmm2,0x30(%rsp)
call   1050 <sin@plt>
movaps 0x30(%rsp),%xmm2
mov    0x60(%rsp),%rcx
xorps  %xmm1,%xmm1
movapd %xmm0,%xmm3
movsd  0x89f(%rip),%xmm0
divsd  0x10(%rsp),%xmm3
mulss  0x88d(%rip),%xmm2
ucomiss %xmm1,%xmm2
mov    0x80(%rsp),%r14d
mov    0xe8(%rsp),%r12
jne    1792 <volume_c_resampling_lanczos2+0x662>
jnp    1380 <volume_c_resampling_lanczos2+0x250>
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
mulsd  0x86f(%rip),%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  %xmm3,0x30(%rsp)
call   1050 <sin@plt>
movsd  0x30(%rsp),%xmm3
mov    0x60(%rsp),%rcx
divsd  0x10(%rsp),%xmm0
jmp    1380 <volume_c_resampling_lanczos2+0x250>
add    $0x138,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [19 x i8] c"slice z = %i / %i\0A\00"
define void @volume_c_resampling_lanczos2(float* %0, i32 %1, i32 %2, i32 %3, float* %4, i32 %5, i32 %6, i32 %7) {
%9 = alloca float*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca float*
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%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
%26 = alloca i32
%27 = alloca i32
%28 = alloca float
%29 = alloca float
%30 = alloca float
%31 = alloca i32
%32 = alloca i32
%33 = alloca i32
%34 = alloca i32
%35 = alloca i32
%36 = alloca i32
%37 = alloca float
%38 = alloca float
%39 = alloca float
%40 = alloca float
%41 = alloca float
%42 = alloca float
store float* %0, float** %9
store i32 %1, i32* %10
store i32 %2, i32* %11
store i32 %3, i32* %12
store float* %4, float** %13
store i32 %5, i32* %14
store i32 %6, i32* %15
store i32 %7, i32* %16
%43 = load i32, i32* %10
%44 = sitofp i32 %43 to float
%45 = load i32, i32* %14
%46 = sitofp i32 %45 to float
%47 = fdiv float %44, %46
store float %47, float* %17
%48 = load i32, i32* %11
%49 = sitofp i32 %48 to float
%50 = load i32, i32* %15
%51 = sitofp i32 %50 to float
%52 = fdiv float %49, %51
store float %52, float* %18
%53 = load i32, i32* %12
%54 = sitofp i32 %53 to float
%55 = load i32, i32* %16
%56 = sitofp i32 %55 to float
%57 = fdiv float %54, %56
store float %57, float* %19
%58 = load i32, i32* %12
%59 = load i32, i32* %11
%60 = mul nsw i32 %58, %59
store i32 %60, i32* %20
%61 = load i32, i32* %16
%62 = load i32, i32* %15
%63 = mul nsw i32 %61, %62
store i32 %63, i32* %21
store i32 0, i32* %24
br label %64
64: ; preds = %335, %8
%65 = load i32, i32* %24
%66 = load i32, i32* %14
%67 = icmp slt i32 %65, %66
br i1 %67, label %68, label %338
68: ; preds = %64
%69 = load i32, i32* %24
%70 = add nsw i32 %69, 1
%71 = load i32, i32* %14
%72 = call i32 @printf(i8* getelementptr inbounds ([19 x i8], [19 x i8]* @.str, i64 0, i64 0), i32 %70, i32 %71)
%73 = load i32, i32* %24
%74 = sitofp i32 %73 to float
%75 = load float, float* %17
%76 = fmul float %74, %75
store float %76, float* %30
%77 = load float, float* %30
%78 = fptosi float %77 to i32
store i32 %78, i32* %27
store i32 0, i32* %23
br label %79
79: ; preds = %331, %68
%80 = load i32, i32* %23
%81 = load i32, i32* %15
%82 = icmp slt i32 %80, %81
br i1 %82, label %83, label %334
83: ; preds = %79
%84 = load i32, i32* %23
%85 = sitofp i32 %84 to float
%86 = load float, float* %18
%87 = fmul float %85, %86
store float %87, float* %29
%88 = load float, float* %29
%89 = fptosi float %88 to i32
store i32 %89, i32* %26
store i32 0, i32* %22
br label %90
90: ; preds = %327, %83
%91 = load i32, i32* %22
%92 = load i32, i32* %16
%93 = icmp slt i32 %91, %92
br i1 %93, label %94, label %330
94: ; preds = %90
%95 = load i32, i32* %22
%96 = sitofp i32 %95 to float
%97 = load float, float* %19
%98 = fmul float %96, %97
store float %98, float* %28
%99 = load float, float* %28
%100 = fptosi float %99 to i32
store i32 %100, i32* %25
store float 0.000000e+00, float* %39
store i32 -1, i32* %34
br label %101
101: ; preds = %310, %94
%102 = load i32, i32* %34
%103 = icmp slt i32 %102, 3
br i1 %103, label %104, label %313
104: ; preds = %101
%105 = load i32, i32* %27
%106 = load i32, i32* %34
%107 = add nsw i32 %105, %106
store i32 %107, i32* %33
%108 = load i32, i32* %33
%109 = load i32, i32* %10
%110 = icmp sge i32 %108, %109
br i1 %110, label %111, label %112
111: ; preds = %104
br label %310
112: ; preds = %104
%113 = load i32, i32* %33
%114 = icmp slt i32 %113, 0
br i1 %114, label %115, label %116
115: ; preds = %112
br label %310
116: ; preds = %112
%117 = load float, float* %30
%118 = load i32, i32* %33
%119 = sitofp i32 %118 to float
%120 = fsub float %117, %119
store float %120, float* %42
%121 = load float, float* %42
%122 = call float @abs(float %121)
%123 = fcmp ogt float %122, 2.000000e+00
br i1 %123, label %124, label %125
124: ; preds = %116
store float 2.000000e+00, float* %42
br label %125
125: ; preds = %124, %116
store float 0.000000e+00, float* %38
store i32 -1, i32* %35
br label %126
126: ; preds = %268, %125
%127 = load i32, i32* %35
%128 = icmp slt i32 %127, 3
br i1 %128, label %129, label %271
129: ; preds = %126
%130 = load i32, i32* %26
%131 = load i32, i32* %35
%132 = add nsw i32 %130, %131
store i32 %132, i32* %32
%133 = load i32, i32* %32
%134 = load i32, i32* %11
%135 = icmp sge i32 %133, %134
br i1 %135, label %136, label %137
136: ; preds = %129
br label %268
137: ; preds = %129
%138 = load i32, i32* %32
%139 = icmp slt i32 %138, 0
br i1 %139, label %140, label %141
140: ; preds = %137
br label %268
141: ; preds = %137
%142 = load float, float* %29
%143 = load i32, i32* %32
%144 = sitofp i32 %143 to float
%145 = fsub float %142, %144
store float %145, float* %41
%146 = load float, float* %41
%147 = call float @abs(float %146)
%148 = fcmp ogt float %147, 2.000000e+00
br i1 %148, label %149, label %150
149: ; preds = %141
store float 2.000000e+00, float* %41
br label %150
150: ; preds = %149, %141
store float 0.000000e+00, float* %37
store i32 -1, i32* %36
br label %151
151: ; preds = %226, %150
%152 = load i32, i32* %36
%153 = icmp slt i32 %152, 3
br i1 %153, label %154, label %229
154: ; preds = %151
%155 = load i32, i32* %25
%156 = load i32, i32* %36
%157 = add nsw i32 %155, %156
store i32 %157, i32* %31
%158 = load i32, i32* %31
%159 = load i32, i32* %12
%160 = icmp sge i32 %158, %159
br i1 %160, label %161, label %162
161: ; preds = %154
br label %226
162: ; preds = %154
%163 = load i32, i32* %31
%164 = icmp slt i32 %163, 0
br i1 %164, label %165, label %166
165: ; preds = %162
br label %226
166: ; preds = %162
%167 = load float, float* %28
%168 = load i32, i32* %31
%169 = sitofp i32 %168 to float
%170 = fsub float %167, %169
store float %170, float* %40
%171 = load float, float* %40
%172 = call float @abs(float %171)
%173 = fcmp ogt float %172, 2.000000e+00
br i1 %173, label %174, label %175
174: ; preds = %166
store float 2.000000e+00, float* %40
br label %175
175: ; preds = %174, %166
%176 = load float, float* %37
%177 = fpext float %176 to double
%178 = load float*, float** %9
%179 = load i32, i32* %33
%180 = load i32, i32* %20
%181 = mul nsw i32 %179, %180
%182 = load i32, i32* %32
%183 = load i32, i32* %12
%184 = mul nsw i32 %182, %183
%185 = add nsw i32 %181, %184
%186 = load i32, i32* %31
%187 = add nsw i32 %185, %186
%188 = sext i32 %187 to i64
%189 = getelementptr inbounds float, float* %178, i64 %188
%190 = load float, float* %189
%191 = fpext float %190 to double
%192 = load float, float* %40
%193 = fcmp oeq float %192, 0.000000e+00
br i1 %193, label %194, label %195
194: ; preds = %175
br label %204
195: ; preds = %175
%196 = load float, float* %40
%197 = fpext float %196 to double
%198 = fmul double 0x400921FB54442D18, %197
%199 = call double @sin(double %198)
%200 = load float, float* %40
%201 = fpext float %200 to double
%202 = fmul double 0x400921FB54442D18, %201
%203 = fdiv double %199, %202
br label %204
204: ; preds = %195, %194
%205 = phi double [ 1.000000e+00, %194 ], [ %203, %195 ]
%206 = fmul double %191, %205
%207 = load float, float* %40
%208 = fmul float %207, 5.000000e-01
%209 = fcmp oeq float %208, 0.000000e+00
br i1 %209, label %210, label %211
210: ; preds = %204
br label %222
211: ; preds = %204
%212 = load float, float* %40
%213 = fmul float %212, 5.000000e-01
%214 = fpext float %213 to double
%215 = fmul double 0x400921FB54442D18, %214
%216 = call double @sin(double %215)
%217 = load float, float* %40
%218 = fmul float %217, 5.000000e-01
%219 = fpext float %218 to double
%220 = fmul double 0x400921FB54442D18, %219
%221 = fdiv double %216, %220
br label %222
222: ; preds = %211, %210
%223 = phi double [ 1.000000e+00, %210 ], [ %221, %211 ]
%224 = call double @llvm.fmuladd.f64(double %206, double %223, double %177)
%225 = fptrunc double %224 to float
store float %225, float* %37
br label %226
226: ; preds = %222, %165, %161
%227 = load i32, i32* %36
%228 = add nsw i32 %227, 1
store i32 %228, i32* %36
br label %151
229: ; preds = %151
%230 = load float, float* %38
%231 = fpext float %230 to double
%232 = load float, float* %37
%233 = fpext float %232 to double
%234 = load float, float* %41
%235 = fcmp oeq float %234, 0.000000e+00
br i1 %235, label %236, label %237
236: ; preds = %229
br label %246
237: ; preds = %229
%238 = load float, float* %41
%239 = fpext float %238 to double
%240 = fmul double 0x400921FB54442D18, %239
%241 = call double @sin(double %240)
%242 = load float, float* %41
%243 = fpext float %242 to double
%244 = fmul double 0x400921FB54442D18, %243
%245 = fdiv double %241, %244
br label %246
246: ; preds = %237, %236
%247 = phi double [ 1.000000e+00, %236 ], [ %245, %237 ]
%248 = fmul double %233, %247
%249 = load float, float* %41
%250 = fmul float %249, 5.000000e-01
%251 = fcmp oeq float %250, 0.000000e+00
br i1 %251, label %252, label %253
252: ; preds = %246
br label %264
253: ; preds = %246
%254 = load float, float* %41
%255 = fmul float %254, 5.000000e-01
%256 = fpext float %255 to double
%257 = fmul double 0x400921FB54442D18, %256
%258 = call double @sin(double %257)
%259 = load float, float* %41
%260 = fmul float %259, 5.000000e-01
%261 = fpext float %260 to double
%262 = fmul double 0x400921FB54442D18, %261
%263 = fdiv double %258, %262
br label %264
264: ; preds = %253, %252
%265 = phi double [ 1.000000e+00, %252 ], [ %263, %253 ]
%266 = call double @llvm.fmuladd.f64(double %248, double %265, double %231)
%267 = fptrunc double %266 to float
store float %267, float* %38
br label %268
268: ; preds = %264, %140, %136
%269 = load i32, i32* %35
%270 = add nsw i32 %269, 1
store i32 %270, i32* %35
br label %126
271: ; preds = %126
%272 = load float, float* %39
%273 = fpext float %272 to double
%274 = load float, float* %38
%275 = fpext float %274 to double
%276 = load float, float* %42
%277 = fcmp oeq float %276, 0.000000e+00
br i1 %277, label %278, label %279
278: ; preds = %271
br label %288
279: ; preds = %271
%280 = load float, float* %42
%281 = fpext float %280 to double
%282 = fmul double 0x400921FB54442D18, %281
%283 = call double @sin(double %282)
%284 = load float, float* %42
%285 = fpext float %284 to double
%286 = fmul double 0x400921FB54442D18, %285
%287 = fdiv double %283, %286
br label %288
288: ; preds = %279, %278
%289 = phi double [ 1.000000e+00, %278 ], [ %287, %279 ]
%290 = fmul double %275, %289
%291 = load float, float* %42
%292 = fmul float %291, 5.000000e-01
%293 = fcmp oeq float %292, 0.000000e+00
br i1 %293, label %294, label %295
294: ; preds = %288
br label %306
295: ; preds = %288
%296 = load float, float* %42
%297 = fmul float %296, 5.000000e-01
%298 = fpext float %297 to double
%299 = fmul double 0x400921FB54442D18, %298
%300 = call double @sin(double %299)
%301 = load float, float* %42
%302 = fmul float %301, 5.000000e-01
%303 = fpext float %302 to double
%304 = fmul double 0x400921FB54442D18, %303
%305 = fdiv double %300, %304
br label %306
306: ; preds = %295, %294
%307 = phi double [ 1.000000e+00, %294 ], [ %305, %295 ]
%308 = call double @llvm.fmuladd.f64(double %290, double %307, double %273)
%309 = fptrunc double %308 to float
store float %309, float* %39
br label %310
310: ; preds = %306, %115, %111
%311 = load i32, i32* %34
%312 = add nsw i32 %311, 1
store i32 %312, i32* %34
br label %101
313: ; preds = %101
%314 = load float, float* %39
%315 = load float*, float** %13
%316 = load i32, i32* %24
%317 = load i32, i32* %21
%318 = mul nsw i32 %316, %317
%319 = load i32, i32* %23
%320 = load i32, i32* %16
%321 = mul nsw i32 %319, %320
%322 = add nsw i32 %318, %321
%323 = load i32, i32* %22
%324 = add nsw i32 %322, %323
%325 = sext i32 %324 to i64
%326 = getelementptr inbounds float, float* %315, i64 %325
store float %314, float* %326
br label %327
327: ; preds = %313
%328 = load i32, i32* %22
%329 = add nsw i32 %328, 1
store i32 %329, i32* %22
br label %90
330: ; preds = %90
br label %331
331: ; preds = %330
%332 = load i32, i32* %23
%333 = add nsw i32 %332, 1
store i32 %333, i32* %23
br label %79
334: ; preds = %79
br label %335
335: ; preds = %334
%336 = load i32, i32* %24
%337 = add nsw i32 %336, 1
store i32 %337, i32* %24
br label %64
338: ; preds = %64
ret void
}

LLVM IR (-O2):
@.str = [19 x i8] c"slice z = %i / %i\0A\00"
define void @volume_c_resampling_lanczos2(float* nocapture %0, i32 %1, i32 %2, i32 %3, float* nocapture writeonly %4, i32 %5, i32 %6, i32 %7) {
%9 = sitofp i32 %1 to float
%10 = sitofp i32 %5 to float
%11 = fdiv float %9, %10
%12 = sitofp i32 %2 to float
%13 = sitofp i32 %6 to float
%14 = fdiv float %12, %13
%15 = sitofp i32 %3 to float
%16 = sitofp i32 %7 to float
%17 = fdiv float %15, %16
%18 = icmp sgt i32 %5, 0
br i1 %18, label %19, label %186
19: ; preds = %8
%20 = icmp sgt i32 %6, 0
%21 = icmp sgt i32 %7, 0
%22 = sext i32 %3 to i64
%23 = zext i32 %7 to i64
br label %26
24: ; preds = %183, %26
%25 = icmp eq i32 %28, %5
br i1 %25, label %186, label %26
26: ; preds = %19, %24
%27 = phi i32 [ 0, %19 ], [ %28, %24 ]
%28 = add nuw nsw i32 %27, 1
%29 = tail call i32 @printf(i8* getelementptr inbounds ([19 x i8], [19 x i8]* @.str, i64 0, i64 0), i32 %28, i32 %5)
%30 = sitofp i32 %27 to float
%31 = fmul float %11, %30
%32 = fptosi float %31 to i32
br i1 %20, label %33, label %24
33: ; preds = %26
%34 = mul i32 %27, %6
br label %35
35: ; preds = %33, %183
%36 = phi i32 [ 0, %33 ], [ %184, %183 ]
%37 = sitofp i32 %36 to float
%38 = fmul float %14, %37
%39 = fptosi float %38 to i32
br i1 %21, label %40, label %183
40: ; preds = %35
%41 = add i32 %36, %34
%42 = mul i32 %41, %7
br label %43
43: ; preds = %40, %177
%44 = phi i64 [ 0, %40 ], [ %181, %177 ]
%45 = trunc i64 %44 to i32
%46 = sitofp i32 %45 to float
%47 = fmul float %17, %46
%48 = fptosi float %47 to i32
%49 = sext i32 %48 to i64
br label %50
50: ; preds = %43, %173
%51 = phi float [ 0.000000e+00, %43 ], [ %174, %173 ]
%52 = phi i32 [ -1, %43 ], [ %175, %173 ]
%53 = add nsw i32 %52, %32
%54 = icmp sge i32 %53, %1
%55 = icmp slt i32 %53, 0
%56 = or i1 %54, %55
br i1 %56, label %173, label %57
57: ; preds = %50
%58 = sitofp i32 %53 to float
%59 = fsub float %31, %58
%60 = tail call float @abs(float %59)
%61 = fcmp ogt float %60, 2.000000e+00
%62 = mul i32 %53, %2
br label %63
63: ; preds = %57, %145
%64 = phi float [ 0.000000e+00, %57 ], [ %146, %145 ]
%65 = phi i32 [ -1, %57 ], [ %147, %145 ]
%66 = add nsw i32 %65, %39
%67 = icmp sge i32 %66, %2
%68 = icmp slt i32 %66, 0
%69 = or i1 %67, %68
br i1 %69, label %145, label %70
70: ; preds = %63
%71 = sitofp i32 %66 to float
%72 = fsub float %38, %71
%73 = tail call float @abs(float %72)
%74 = fcmp ogt float %73, 2.000000e+00
%75 = add i32 %66, %62
%76 = mul i32 %75, %3
br label %77
77: ; preds = %70, %117
%78 = phi i64 [ -1, %70 ], [ %119, %117 ]
%79 = phi float [ 0.000000e+00, %70 ], [ %118, %117 ]
%80 = add nsw i64 %78, %49
%81 = icmp sge i64 %80, %22
%82 = icmp slt i64 %80, 0
%83 = or i1 %81, %82
br i1 %83, label %117, label %84
84: ; preds = %77
%85 = trunc i64 %80 to i32
%86 = sitofp i32 %85 to float
%87 = fsub float %47, %86
%88 = tail call float @abs(float %87)
%89 = fcmp ogt float %88, 2.000000e+00
%90 = select i1 %89, float 2.000000e+00, float %87
%91 = fpext float %79 to double
%92 = add i32 %76, %85
%93 = sext i32 %92 to i64
%94 = getelementptr inbounds float, float* %0, i64 %93
%95 = load float, float* %94
%96 = fpext float %95 to double
%97 = fcmp oeq float %90, 0.000000e+00
br i1 %97, label %103, label %98
98: ; preds = %84
%99 = fpext float %90 to double
%100 = fmul double %99, 0x400921FB54442D18
%101 = tail call double @sin(double %100)
%102 = fdiv double %101, %100
br label %103
103: ; preds = %84, %98
%104 = phi double [ %102, %98 ], [ 1.000000e+00, %84 ]
%105 = fmul double %104, %96
%106 = fmul float %90, 5.000000e-01
%107 = fcmp oeq float %106, 0.000000e+00
br i1 %107, label %113, label %108
108: ; preds = %103
%109 = fpext float %106 to double
%110 = fmul double %109, 0x400921FB54442D18
%111 = tail call double @sin(double %110)
%112 = fdiv double %111, %110
br label %113
113: ; preds = %103, %108
%114 = phi double [ %112, %108 ], [ 1.000000e+00, %103 ]
%115 = tail call double @llvm.fmuladd.f64(double %105, double %114, double %91)
%116 = fptrunc double %115 to float
br label %117
117: ; preds = %77, %113
%118 = phi float [ %79, %77 ], [ %116, %113 ]
%119 = add nsw i64 %78, 1
%120 = icmp eq i64 %119, 3
br i1 %120, label %121, label %77
121: ; preds = %117
%122 = select i1 %74, float 2.000000e+00, float %72
%123 = fpext float %64 to double
%124 = fpext float %118 to double
%125 = fcmp oeq float %122, 0.000000e+00
br i1 %125, label %131, label %126
126: ; preds = %121
%127 = fpext float %122 to double
%128 = fmul double %127, 0x400921FB54442D18
%129 = tail call double @sin(double %128)
%130 = fdiv double %129, %128
br label %131
131: ; preds = %121, %126
%132 = phi double [ %130, %126 ], [ 1.000000e+00, %121 ]
%133 = fmul double %132, %124
%134 = fmul float %122, 5.000000e-01
%135 = fcmp oeq float %134, 0.000000e+00
br i1 %135, label %141, label %136
136: ; preds = %131
%137 = fpext float %134 to double
%138 = fmul double %137, 0x400921FB54442D18
%139 = tail call double @sin(double %138)
%140 = fdiv double %139, %138
br label %141
141: ; preds = %131, %136
%142 = phi double [ %140, %136 ], [ 1.000000e+00, %131 ]
%143 = tail call double @llvm.fmuladd.f64(double %133, double %142, double %123)
%144 = fptrunc double %143 to float
br label %145
145: ; preds = %63, %141
%146 = phi float [ %64, %63 ], [ %144, %141 ]
%147 = add nsw i32 %65, 1
%148 = icmp eq i32 %147, 3
br i1 %148, label %149, label %63
149: ; preds = %145
%150 = select i1 %61, float 2.000000e+00, float %59
%151 = fpext float %51 to double
%152 = fpext float %146 to double
%153 = fcmp oeq float %150, 0.000000e+00
br i1 %153, label %159, label %154
154: ; preds = %149
%155 = fpext float %150 to double
%156 = fmul double %155, 0x400921FB54442D18
%157 = tail call double @sin(double %156)
%158 = fdiv double %157, %156
br label %159
159: ; preds = %149, %154
%160 = phi double [ %158, %154 ], [ 1.000000e+00, %149 ]
%161 = fmul double %160, %152
%162 = fmul float %150, 5.000000e-01
%163 = fcmp oeq float %162, 0.000000e+00
br i1 %163, label %169, label %164
164: ; preds = %159
%165 = fpext float %162 to double
%166 = fmul double %165, 0x400921FB54442D18
%167 = tail call double @sin(double %166)
%168 = fdiv double %167, %166
br label %169
169: ; preds = %159, %164
%170 = phi double [ %168, %164 ], [ 1.000000e+00, %159 ]
%171 = tail call double @llvm.fmuladd.f64(double %161, double %170, double %151)
%172 = fptrunc double %171 to float
br label %173
173: ; preds = %50, %169
%174 = phi float [ %51, %50 ], [ %172, %169 ]
%175 = add nsw i32 %52, 1
%176 = icmp eq i32 %175, 3
br i1 %176, label %177, label %50
177: ; preds = %173
%178 = add i32 %42, %45
%179 = sext i32 %178 to i64
%180 = getelementptr inbounds float, float* %4, i64 %179
store float %174, float* %180
%181 = add nuw nsw i64 %44, 1
%182 = icmp eq i64 %181, %23
br i1 %182, label %183, label %43
183: ; preds = %177, %35
%184 = add nuw nsw i32 %36, 1
%185 = icmp eq i32 %184, %6
br i1 %185, label %24, label %35
186: ; preds = %24, %8
ret void
}