#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 */
double sqrt(double);
int pcapack_scale(const bool centerx, const bool scalex, const int m,
                  const int n, double *restrict x) {
  int i, j;
  double colmean, colvar;
  double dt, tmp;
  if (m == 0 || n == 0)
    return 0;
  ;
  if (centerx && scalex) {
    tmp = 1. / ((double)m - 1);
#pragma omp parallel for private(i, j, colmean, colvar, dt)                    \
    shared(x) if (m * n > OMP_MIN_SIZE)
    for (j = 0; j < n; j++) {
      colmean = 0;
      colvar = 0;

      for (i = 0; i < m; i++) {
        dt = x[i + m * j] - colmean;
        colmean += dt / ((double)i + 1);
        colvar += dt * (x[i + m * j] - colmean);
      }
      colvar = sqrt(colvar * tmp);

      for (i = 0; i < m; i++)
        x[i + m * j] = (x[i + m * j] - colmean) / colvar;
    }
  } else if (centerx) {
    const double div = 1. / ((double)m);
#pragma omp parallel for private(i, j, colmean)                                \
    shared(x) if (m * n > OMP_MIN_SIZE)
    for (j = 0; j < n; j++) {
      colmean = 0;

      for (i = 0; i < m; i++)
        colmean += x[i + m * j] * div;

      for (i = 0; i < m; i++)
        x[i + m * j] -= colmean;
    }
  } else if (scalex) {
    const double div = 1. / ((double)m - 1);
#pragma omp parallel for private(i, j, colvar, tmp)                            \
    shared(x) if (m * n > OMP_MIN_SIZE)
    for (j = 0; j < n; j++) {
      colvar = 0;

      for (i = 0; i < m; i++) {
        tmp = x[i + m * j];
        colvar += tmp * tmp * div;
      }
      colvar = sqrt(colvar);

      for (i = 0; i < m; i++)
        x[i + m * j] /= colvar;
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<pcapack_scale>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8,-0x20(%rbp)
cmpl   $0x0,-0x10(%rbp)
je     113c <pcapack_scale+0x2c>
cmpl   $0x0,-0x14(%rbp)
jne    1148 <pcapack_scale+0x38>
movl   $0x0,-0x4(%rbp)
jmp    14d2 <pcapack_scale+0x3c2>
cmpl   $0x0,-0x8(%rbp)
je     12c3 <pcapack_scale+0x1b3>
cmpl   $0x0,-0xc(%rbp)
je     12c3 <pcapack_scale+0x1b3>
cvtsi2sdl -0x10(%rbp),%xmm1
movsd  0xe97(%rip),%xmm0
subsd  %xmm0,%xmm1
movsd  0xe8b(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    12be <pcapack_scale+0x1ae>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    123c <pcapack_scale+0x12c>
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x28(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
cvtsi2sdl -0x24(%rbp),%xmm1
movsd  0xe18(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x28(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
subsd  -0x30(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11a8 <pcapack_scale+0x98>
movsd  -0x38(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    12ab <pcapack_scale+0x19b>
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x28(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x30(%rbp),%xmm0
divsd  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x28(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1257 <pcapack_scale+0x147>
jmp    12b0 <pcapack_scale+0x1a0>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1185 <pcapack_scale+0x75>
jmp    14cb <pcapack_scale+0x3bb>
cmpl   $0x0,-0x8(%rbp)
je     13ad <pcapack_scale+0x29d>
cvtsi2sdl -0x10(%rbp),%xmm1
movsd  0xd26(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    13a8 <pcapack_scale+0x298>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    134e <pcapack_scale+0x23e>
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x28(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1305 <pcapack_scale+0x1f5>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1395 <pcapack_scale+0x285>
movsd  -0x30(%rbp),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x28(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1355 <pcapack_scale+0x245>
jmp    139a <pcapack_scale+0x28a>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    12ea <pcapack_scale+0x1da>
jmp    14c6 <pcapack_scale+0x3b6>
cmpl   $0x0,-0xc(%rbp)
je     14c1 <pcapack_scale+0x3b1>
cvtsi2sdl -0x10(%rbp),%xmm1
movsd  0xc3c(%rip),%xmm0
subsd  %xmm0,%xmm1
movsd  0xc30(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    14bc <pcapack_scale+0x3ac>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1453 <pcapack_scale+0x343>
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x28(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    13fb <pcapack_scale+0x2eb>
movsd  -0x38(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    14a9 <pcapack_scale+0x399>
movsd  -0x38(%rbp),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x28(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1469 <pcapack_scale+0x359>
jmp    14ae <pcapack_scale+0x39e>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    13e0 <pcapack_scale+0x2d0>
jmp    14c1 <pcapack_scale+0x3b1>
jmp    14c6 <pcapack_scale+0x3b6>
jmp    14cb <pcapack_scale+0x3bb>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<pcapack_scale>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %r8,0x8(%rsp)
mov    %edx,0x4(%rsp)
test   %edx,%edx
je     17f9 <pcapack_scale+0x6e9>
test   %ecx,%ecx
je     17f9 <pcapack_scale+0x6e9>
test   %edi,%edi
je     138f <pcapack_scale+0x27f>
test   %esi,%esi
je     138f <pcapack_scale+0x27f>
test   %ecx,%ecx
jle    17f9 <pcapack_scale+0x6e9>
mov    0x4(%rsp),%eax
cvtsi2sd %eax,%xmm0
addsd  0xea1(%rip),%xmm0
movsd  0xea1(%rip),%xmm5
movapd %xmm5,%xmm1
divsd  %xmm0,%xmm1
movslq %eax,%rdx
mov    %ecx,%ebp
xorpd  %xmm6,%xmm6
movsd  %xmm1,0x38(%rsp)
mulsd  %xmm6,%xmm1
movsd  %xmm1,0x10(%rsp)
mov    %eax,%ebx
mov    %ebx,%r15d
and    $0xfffffffe,%r15d
lea    -0x2(%r15),%rax
mov    %rax,0x30(%rsp)
mov    %rax,%rcx
shr    %rcx
add    $0x1,%rcx
mov    %rcx,%rax
mov    %rcx,0x28(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,0x20(%rsp)
mov    %rdx,0x18(%rsp)
lea    0x0(,%rdx,8),%r14
xor    %r12d,%r12d
mov    0x8(%rsp),%r13
jmp    11e0 <pcapack_scale+0xd0>
nopl   0x0(%rax)
add    $0x1,%r12
add    %r14,%r13
cmp    %rbp,%r12
je     17f9 <pcapack_scale+0x6e9>
cmpl   $0x0,0x4(%rsp)
jle    1260 <pcapack_scale+0x150>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
xorpd  %xmm7,%xmm7
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movapd %xmm7,%xmm1
movsd  0x0(%r13,%rax,8),%xmm2
movapd %xmm2,%xmm3
subsd  %xmm7,%xmm3
xorps  %xmm4,%xmm4
cvtsi2sd %eax,%xmm4
addsd  %xmm5,%xmm4
movapd %xmm3,%xmm7
divsd  %xmm4,%xmm7
addsd  %xmm1,%xmm7
subsd  %xmm7,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm2,%xmm0
add    $0x1,%rax
cmp    %rax,%rbx
jne    1200 <pcapack_scale+0xf0>
mulsd  0x38(%rsp),%xmm0
ucomisd %xmm6,%xmm0
jb     1290 <pcapack_scale+0x180>
sqrtsd %xmm0,%xmm0
cmpl   $0x0,0x4(%rsp)
jg     12b8 <pcapack_scale+0x1a8>
jmp    11d0 <pcapack_scale+0xc0>
nopl   0x0(%rax,%rax,1)
movsd  0x10(%rsp),%xmm0
ucomisd %xmm6,%xmm0
jae    11d0 <pcapack_scale+0xc0>
movsd  0x10(%rsp),%xmm0
call   1030 <sqrt@plt>
xorpd  %xmm6,%xmm6
movsd  0xd81(%rip),%xmm5
jmp    11d0 <pcapack_scale+0xc0>
nopl   0x0(%rax)
movapd %xmm7,0x40(%rsp)
call   1030 <sqrt@plt>
movapd 0x40(%rsp),%xmm7
xorpd  %xmm6,%xmm6
movsd  0xd5b(%rip),%xmm5
cmpl   $0x0,0x4(%rsp)
jle    11d0 <pcapack_scale+0xc0>
cmpl   $0x1,0x4(%rsp)
jne    12c6 <pcapack_scale+0x1b6>
xor    %eax,%eax
jmp    1360 <pcapack_scale+0x250>
movapd %xmm7,%xmm1
unpcklpd %xmm7,%xmm1
movapd %xmm0,%xmm2
unpcklpd %xmm0,%xmm2
cmpq   $0x0,0x30(%rsp)
je     1384 <pcapack_scale+0x274>
mov    0x20(%rsp),%rcx
xor    %eax,%eax
nopl   0x0(%rax)
movupd 0x0(%r13,%rax,8),%xmm3
movupd 0x10(%r13,%rax,8),%xmm4
subpd  %xmm1,%xmm3
divpd  %xmm2,%xmm3
movupd %xmm3,0x0(%r13,%rax,8)
subpd  %xmm1,%xmm4
divpd  %xmm2,%xmm4
movupd %xmm4,0x10(%r13,%rax,8)
add    $0x4,%rax
add    $0xfffffffffffffffe,%rcx
jne    12f0 <pcapack_scale+0x1e0>
testb  $0x1,0x28(%rsp)
je     1350 <pcapack_scale+0x240>
mov    %r12,%rcx
imul   0x18(%rsp),%rcx
add    %rcx,%rax
mov    0x8(%rsp),%rcx
movupd (%rcx,%rax,8),%xmm3
subpd  %xmm1,%xmm3
divpd  %xmm2,%xmm3
movupd %xmm3,(%rcx,%rax,8)
mov    %r15,%rax
cmp    %rbx,%r15
je     11d0 <pcapack_scale+0xc0>
nopl   0x0(%rax)
movsd  0x0(%r13,%rax,8),%xmm1
subsd  %xmm7,%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%r13,%rax,8)
add    $0x1,%rax
cmp    %rax,%rbx
jne    1360 <pcapack_scale+0x250>
jmp    11d0 <pcapack_scale+0xc0>
xor    %eax,%eax
testb  $0x1,0x28(%rsp)
jne    132d <pcapack_scale+0x21d>
jmp    1350 <pcapack_scale+0x240>
test   %edi,%edi
je     1586 <pcapack_scale+0x476>
test   %ecx,%ecx
jle    17f9 <pcapack_scale+0x6e9>
mov    0x4(%rsp),%eax
cvtsi2sd %eax,%xmm1
movsd  0xc59(%rip),%xmm0
divsd  %xmm1,%xmm0
movslq %eax,%r8
mov    %ecx,%r15d
mov    %eax,%ebp
lea    -0x1(%rbp),%r14
mov    %ebp,%edi
and    $0xfffffffc,%edi
lea    -0x4(%rdi),%r10
mov    %r10,%r11
shr    $0x2,%r11
add    $0x1,%r11
mov    %ebp,%edx
and    $0x3,%edx
mov    %r11,%r9
and    $0xfffffffffffffffe,%r9
mov    0x8(%rsp),%rcx
lea    0x18(%rcx),%rax
lea    0x0(,%r8,8),%r12
xor    %r13d,%r13d
jmp    1413 <pcapack_scale+0x303>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r13
add    %r12,%rax
add    %r12,%rcx
cmp    %r15,%r13
je     17f9 <pcapack_scale+0x6e9>
cmpl   $0x0,0x4(%rsp)
jle    1400 <pcapack_scale+0x2f0>
xorpd  %xmm1,%xmm1
xor    %ebx,%ebx
cmp    $0x3,%r14
jb     1470 <pcapack_scale+0x360>
cs nopw 0x0(%rax,%rax,1)
movsd  -0x18(%rax,%rbx,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
movsd  -0x10(%rax,%rbx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
movsd  -0x8(%rax,%rbx,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
movsd  (%rax,%rbx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
add    $0x4,%rbx
cmp    %rbx,%rdi
jne    1430 <pcapack_scale+0x320>
test   %rdx,%rdx
je     1496 <pcapack_scale+0x386>
lea    (%rcx,%rbx,8),%rsi
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
movsd  (%rsi,%rbx,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm2,%xmm1
add    $0x1,%rbx
cmp    %rbx,%rdx
jne    1480 <pcapack_scale+0x370>
cmpl   $0x0,0x4(%rsp)
jle    1400 <pcapack_scale+0x2f0>
cmpl   $0x4,0x4(%rsp)
jae    14af <pcapack_scale+0x39f>
xor    %esi,%esi
jmp    1560 <pcapack_scale+0x450>
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
test   %r10,%r10
je     157c <pcapack_scale+0x46c>
mov    %r9,%rsi
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
movupd (%rcx,%rbx,8),%xmm3
movupd 0x10(%rcx,%rbx,8),%xmm4
movupd 0x20(%rcx,%rbx,8),%xmm5
movupd 0x30(%rcx,%rbx,8),%xmm6
subpd  %xmm2,%xmm3
subpd  %xmm2,%xmm4
movupd %xmm3,(%rcx,%rbx,8)
movupd %xmm4,0x10(%rcx,%rbx,8)
subpd  %xmm2,%xmm5
subpd  %xmm2,%xmm6
movupd %xmm5,0x20(%rcx,%rbx,8)
movupd %xmm6,0x30(%rcx,%rbx,8)
add    $0x8,%rbx
add    $0xfffffffffffffffe,%rsi
jne    14d0 <pcapack_scale+0x3c0>
test   $0x1,%r11b
je     154b <pcapack_scale+0x43b>
mov    %r13,%rsi
imul   %r8,%rsi
add    %rsi,%rbx
mov    0x8(%rsp),%rsi
movupd (%rsi,%rbx,8),%xmm3
movupd 0x10(%rsi,%rbx,8),%xmm4
subpd  %xmm2,%xmm3
subpd  %xmm2,%xmm4
movupd %xmm3,(%rsi,%rbx,8)
movupd %xmm4,0x10(%rsi,%rbx,8)
mov    %rdi,%rsi
cmp    %rbp,%rdi
je     1400 <pcapack_scale+0x2f0>
nopw   0x0(%rax,%rax,1)
movsd  (%rcx,%rsi,8),%xmm2
subsd  %xmm1,%xmm2
movsd  %xmm2,(%rcx,%rsi,8)
add    $0x1,%rsi
cmp    %rsi,%rbp
jne    1560 <pcapack_scale+0x450>
jmp    1400 <pcapack_scale+0x2f0>
xor    %ebx,%ebx
test   $0x1,%r11b
jne    151e <pcapack_scale+0x40e>
jmp    154b <pcapack_scale+0x43b>
test   %esi,%esi
je     17f9 <pcapack_scale+0x6e9>
test   %ecx,%ecx
jle    17f9 <pcapack_scale+0x6e9>
mov    0x4(%rsp),%eax
cvtsi2sd %eax,%xmm0
addsd  0xa5a(%rip),%xmm0
movsd  0xa5a(%rip),%xmm4
divsd  %xmm0,%xmm4
movslq %eax,%rdx
mov    %ecx,%esi
mov    %eax,%r13d
lea    -0x1(%r13),%rdi
mov    %r13d,%eax
and    $0xfffffffe,%eax
mov    %rax,0x28(%rsp)
add    $0xfffffffffffffffe,%rax
mov    %rax,0x20(%rsp)
mov    %rax,%rcx
shr    %rcx
add    $0x1,%rcx
mov    %r13d,%ebx
and    $0x3,%ebx
mov    %r13d,%r14d
and    $0xfffffffc,%r14d
mov    %rcx,%rax
mov    %rcx,0x18(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,0x60(%rsp)
mov    0x8(%rsp),%rbp
lea    0x18(%rbp),%r15
mov    %rdx,0x58(%rsp)
lea    0x0(,%rdx,8),%rdx
xor    %r12d,%r12d
xorpd  %xmm5,%xmm5
movsd  %xmm4,0x38(%rsp)
mov    %rsi,0x40(%rsp)
mov    %rdi,0x10(%rsp)
mov    %rdx,0x30(%rsp)
jmp    1643 <pcapack_scale+0x533>
xchg   %ax,%ax
add    $0x1,%r12
add    %rdx,%r15
add    %rdx,%rbp
cmp    %rsi,%r12
je     17f9 <pcapack_scale+0x6e9>
cmpl   $0x0,0x4(%rsp)
jle    1630 <pcapack_scale+0x520>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
cmp    $0x3,%rdi
jb     16b4 <pcapack_scale+0x5a4>
cs nopw 0x0(%rax,%rax,1)
movsd  -0x18(%r15,%rax,8),%xmm1
movsd  -0x10(%r15,%rax,8),%xmm2
mulsd  %xmm1,%xmm1
mulsd  %xmm4,%xmm1
addsd  %xmm0,%xmm1
mulsd  %xmm2,%xmm2
mulsd  %xmm4,%xmm2
addsd  %xmm1,%xmm2
movsd  -0x8(%r15,%rax,8),%xmm1
mulsd  %xmm1,%xmm1
mulsd  %xmm4,%xmm1
addsd  %xmm2,%xmm1
movsd  (%r15,%rax,8),%xmm0
mulsd  %xmm0,%xmm0
mulsd  %xmm4,%xmm0
addsd  %xmm1,%xmm0
add    $0x4,%rax
cmp    %rax,%r14
jne    1660 <pcapack_scale+0x550>
test   %rbx,%rbx
je     16ea <pcapack_scale+0x5da>
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm1
mulsd  %xmm4,%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rcx
cmp    %rcx,%rbx
jne    16d0 <pcapack_scale+0x5c0>
ucomisd %xmm5,%xmm0
jb     1700 <pcapack_scale+0x5f0>
sqrtsd %xmm0,%xmm0
cmpl   $0x0,0x4(%rsp)
jg     1729 <pcapack_scale+0x619>
jmp    1630 <pcapack_scale+0x520>
call   1030 <sqrt@plt>
xorpd  %xmm5,%xmm5
mov    0x30(%rsp),%rdx
mov    0x10(%rsp),%rdi
mov    0x40(%rsp),%rsi
movsd  0x38(%rsp),%xmm4
cmpl   $0x0,0x4(%rsp)
jle    1630 <pcapack_scale+0x520>
cmpl   $0x1,0x4(%rsp)
jne    1737 <pcapack_scale+0x627>
xor    %eax,%eax
jmp    17d0 <pcapack_scale+0x6c0>
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
cmpq   $0x0,0x20(%rsp)
je     17ee <pcapack_scale+0x6de>
mov    0x60(%rsp),%rcx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd 0x0(%rbp,%rax,8),%xmm2
movupd 0x10(%rbp,%rax,8),%xmm3
divpd  %xmm1,%xmm2
movupd %xmm2,0x0(%rbp,%rax,8)
divpd  %xmm1,%xmm3
movupd %xmm3,0x10(%rbp,%rax,8)
add    $0x4,%rax
add    $0xfffffffffffffffe,%rcx
jne    1760 <pcapack_scale+0x650>
testb  $0x1,0x18(%rsp)
je     17b0 <pcapack_scale+0x6a0>
mov    %r12,%rcx
imul   0x58(%rsp),%rcx
add    %rcx,%rax
mov    0x8(%rsp),%rcx
movupd (%rcx,%rax,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rcx,%rax,8)
mov    0x28(%rsp),%rcx
mov    %rcx,%rax
cmp    %r13,%rcx
je     1630 <pcapack_scale+0x520>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%rbp,%rax,8)
add    $0x1,%rax
cmp    %rax,%r13
jne    17d0 <pcapack_scale+0x6c0>
jmp    1630 <pcapack_scale+0x520>
xor    %eax,%eax
testb  $0x1,0x18(%rsp)
jne    1791 <pcapack_scale+0x681>
jmp    17b0 <pcapack_scale+0x6a0>
xor    %eax,%eax
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @pcapack_scale(i32 %0, i32 %1, i32 %2, i32 %3, double* noalias nocapture %4) {
%6 = icmp eq i32 %2, 0
%7 = icmp eq i32 %3, 0
%8 = or i1 %6, %7
br i1 %8, label %381, label %9
9: ; preds = %5
%10 = icmp ne i32 %0, 0
%11 = icmp ne i32 %1, 0
%12 = and i1 %10, %11
br i1 %12, label %13, label %115
13: ; preds = %9
%14 = sitofp i32 %2 to double
%15 = fadd double %14, -1.000000e+00
%16 = fdiv double 1.000000e+00, %15
%17 = icmp sgt i32 %3, 0
br i1 %17, label %18, label %381
18: ; preds = %13
%19 = icmp sgt i32 %2, 0
%20 = sext i32 %2 to i64
%21 = zext i32 %3 to i64
%22 = fmul double %16, 0.000000e+00
%23 = zext i32 %2 to i64
%24 = zext i32 %2 to i64
%25 = and i64 %23, 4294967294
%26 = add nsw i64 %25, -2
%27 = lshr exact i64 %26, 1
%28 = add nuw i64 %27, 1
%29 = icmp eq i32 %2, 1
%30 = and i64 %23, 4294967294
%31 = and i64 %28, 1
%32 = icmp eq i64 %26, 0
%33 = and i64 %28, -2
%34 = icmp eq i64 %31, 0
%35 = icmp eq i64 %30, %23
br label %36
36: ; preds = %18, %112
%37 = phi i64 [ 0, %18 ], [ %113, %112 ]
br i1 %19, label %40, label %38
38: ; preds = %36
%39 = tail call double @sqrt(double %22)
br label %112
40: ; preds = %36
%41 = mul nsw i64 %37, %20
br label %42
42: ; preds = %40, %42
%43 = phi i64 [ 0, %40 ], [ %57, %42 ]
%44 = phi double [ 0.000000e+00, %40 ], [ %56, %42 ]
%45 = phi double [ 0.000000e+00, %40 ], [ %54, %42 ]
%46 = add nsw i64 %43, %41
%47 = getelementptr inbounds double, double* %4, i64 %46
%48 = load double, double* %47
%49 = fsub double %48, %45
%50 = trunc i64 %43 to i32
%51 = sitofp i32 %50 to double
%52 = fadd double %51, 1.000000e+00
%53 = fdiv double %49, %52
%54 = fadd double %45, %53
%55 = fsub double %48, %54
%56 = tail call double @llvm.fmuladd.f64(double %49, double %55, double %44)
%57 = add nuw nsw i64 %43, 1
%58 = icmp eq i64 %57, %23
br i1 %58, label %59, label %42
59: ; preds = %42
%60 = fmul double %16, %56
%61 = tail call double @sqrt(double %60)
br i1 %19, label %62, label %112
62: ; preds = %59
%63 = mul nsw i64 %37, %20
br i1 %29, label %101, label %64
64: ; preds = %62
%65 = insertelement <2 x double> poison, double %54, i64 0
%66 = shufflevector <2 x double> %65, <2 x double> poison, <2 x i32> zeroinitializer
%67 = insertelement <2 x double> poison, double %61, i64 0
%68 = shufflevector <2 x double> %67, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %32, label %90, label %69
69: ; preds = %64, %69
%70 = phi i64 [ %87, %69 ], [ 0, %64 ]
%71 = phi i64 [ %88, %69 ], [ 0, %64 ]
%72 = add nsw i64 %70, %63
%73 = getelementptr inbounds double, double* %4, i64 %72
%74 = bitcast double* %73 to <2 x double>*
%75 = load <2 x double>, <2 x double>* %74
%76 = fsub <2 x double> %75, %66
%77 = fdiv <2 x double> %76, %68
%78 = bitcast double* %73 to <2 x double>*
store <2 x double> %77, <2 x double>* %78
%79 = or i64 %70, 2
%80 = add nsw i64 %79, %63
%81 = getelementptr inbounds double, double* %4, i64 %80
%82 = bitcast double* %81 to <2 x double>*
%83 = load <2 x double>, <2 x double>* %82
%84 = fsub <2 x double> %83, %66
%85 = fdiv <2 x double> %84, %68
%86 = bitcast double* %81 to <2 x double>*
store <2 x double> %85, <2 x double>* %86
%87 = add nuw i64 %70, 4
%88 = add nuw i64 %71, 2
%89 = icmp eq i64 %88, %33
br i1 %89, label %90, label %69
90: ; preds = %69, %64
%91 = phi i64 [ 0, %64 ], [ %87, %69 ]
br i1 %34, label %100, label %92
92: ; preds = %90
%93 = add nsw i64 %91, %63
%94 = getelementptr inbounds double, double* %4, i64 %93
%95 = bitcast double* %94 to <2 x double>*
%96 = load <2 x double>, <2 x double>* %95
%97 = fsub <2 x double> %96, %66
%98 = fdiv <2 x double> %97, %68
%99 = bitcast double* %94 to <2 x double>*
store <2 x double> %98, <2 x double>* %99
br label %100
100: ; preds = %90, %92
br i1 %35, label %112, label %101
101: ; preds = %62, %100
%102 = phi i64 [ 0, %62 ], [ %30, %100 ]
br label %103
103: ; preds = %101, %103
%104 = phi i64 [ %110, %103 ], [ %102, %101 ]
%105 = add nsw i64 %104, %63
%106 = getelementptr inbounds double, double* %4, i64 %105
%107 = load double, double* %106
%108 = fsub double %107, %54
%109 = fdiv double %108, %61
store double %109, double* %106
%110 = add nuw nsw i64 %104, 1
%111 = icmp eq i64 %110, %24
br i1 %111, label %112, label %103
112: ; preds = %103, %100, %38, %59
%113 = add nuw nsw i64 %37, 1
%114 = icmp eq i64 %113, %21
br i1 %114, label %381, label %36
115: ; preds = %9
br i1 %10, label %116, label %253
116: ; preds = %115
%117 = sitofp i32 %2 to double
%118 = fdiv double 1.000000e+00, %117
%119 = icmp sgt i32 %3, 0
br i1 %119, label %120, label %381
120: ; preds = %116
%121 = icmp sgt i32 %2, 0
%122 = sext i32 %2 to i64
%123 = zext i32 %3 to i64
%124 = zext i32 %2 to i64
%125 = zext i32 %2 to i64
%126 = add nsw i64 %124, -1
%127 = and i64 %124, 4294967292
%128 = add nsw i64 %127, -4
%129 = lshr exact i64 %128, 2
%130 = add nuw nsw i64 %129, 1
%131 = and i64 %124, 3
%132 = icmp ult i64 %126, 3
%133 = and i64 %124, 4294967292
%134 = icmp eq i64 %131, 0
%135 = icmp ult i32 %2, 4
%136 = and i64 %124, 4294967292
%137 = and i64 %130, 1
%138 = icmp eq i64 %128, 0
%139 = and i64 %130, 9223372036854775806
%140 = icmp eq i64 %137, 0
%141 = icmp eq i64 %136, %124
br label %142
142: ; preds = %120, %250
%143 = phi i64 [ 0, %120 ], [ %251, %250 ]
br i1 %121, label %144, label %250
144: ; preds = %142
%145 = mul nsw i64 %143, %122
br i1 %132, label %146, label %216
146: ; preds = %216, %144
%147 = phi double [ undef, %144 ], [ %238, %216 ]
%148 = phi i64 [ 0, %144 ], [ %239, %216 ]
%149 = phi double [ 0.000000e+00, %144 ], [ %238, %216 ]
br i1 %134, label %161, label %150
150: ; preds = %146, %150
%151 = phi i64 [ %158, %150 ], [ %148, %146 ]
%152 = phi double [ %157, %150 ], [ %149, %146 ]
%153 = phi i64 [ %159, %150 ], [ 0, %146 ]
%154 = add nsw i64 %151, %145
%155 = getelementptr inbounds double, double* %4, i64 %154
%156 = load double, double* %155
%157 = tail call double @llvm.fmuladd.f64(double %156, double %118, double %152)
%158 = add nuw nsw i64 %151, 1
%159 = add i64 %153, 1
%160 = icmp eq i64 %159, %131
br i1 %160, label %161, label %150
161: ; preds = %150, %146
%162 = phi double [ %147, %146 ], [ %157, %150 ]
br i1 %121, label %163, label %250
163: ; preds = %161
%164 = mul nsw i64 %143, %122
br i1 %135, label %214, label %165
165: ; preds = %163
%166 = insertelement <2 x double> poison, double %162, i64 0
%167 = shufflevector <2 x double> %166, <2 x double> poison, <2 x i32> zeroinitializer
%168 = insertelement <2 x double> poison, double %162, i64 0
%169 = shufflevector <2 x double> %168, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %138, label %199, label %170
170: ; preds = %165, %170
%171 = phi i64 [ %196, %170 ], [ 0, %165 ]
%172 = phi i64 [ %197, %170 ], [ 0, %165 ]
%173 = add nsw i64 %171, %164
%174 = getelementptr inbounds double, double* %4, i64 %173
%175 = bitcast double* %174 to <2 x double>*
%176 = load <2 x double>, <2 x double>* %175
%177 = getelementptr inbounds double, double* %174, i64 2
%178 = bitcast double* %177 to <2 x double>*
%179 = load <2 x double>, <2 x double>* %178
%180 = fsub <2 x double> %176, %167
%181 = fsub <2 x double> %179, %169
%182 = bitcast double* %174 to <2 x double>*
store <2 x double> %180, <2 x double>* %182
%183 = bitcast double* %177 to <2 x double>*
store <2 x double> %181, <2 x double>* %183
%184 = or i64 %171, 4
%185 = add nsw i64 %184, %164
%186 = getelementptr inbounds double, double* %4, i64 %185
%187 = bitcast double* %186 to <2 x double>*
%188 = load <2 x double>, <2 x double>* %187
%189 = getelementptr inbounds double, double* %186, i64 2
%190 = bitcast double* %189 to <2 x double>*
%191 = load <2 x double>, <2 x double>* %190
%192 = fsub <2 x double> %188, %167
%193 = fsub <2 x double> %191, %169
%194 = bitcast double* %186 to <2 x double>*
store <2 x double> %192, <2 x double>* %194
%195 = bitcast double* %189 to <2 x double>*
store <2 x double> %193, <2 x double>* %195
%196 = add nuw i64 %171, 8
%197 = add i64 %172, 2
%198 = icmp eq i64 %197, %139
br i1 %198, label %199, label %170
199: ; preds = %170, %165
%200 = phi i64 [ 0, %165 ], [ %196, %170 ]
br i1 %140, label %213, label %201
201: ; preds = %199
%202 = add nsw i64 %200, %164
%203 = getelementptr inbounds double, double* %4, i64 %202
%204 = bitcast double* %203 to <2 x double>*
%205 = load <2 x double>, <2 x double>* %204
%206 = getelementptr inbounds double, double* %203, i64 2
%207 = bitcast double* %206 to <2 x double>*
%208 = load <2 x double>, <2 x double>* %207
%209 = fsub <2 x double> %205, %167
%210 = fsub <2 x double> %208, %169
%211 = bitcast double* %203 to <2 x double>*
store <2 x double> %209, <2 x double>* %211
%212 = bitcast double* %206 to <2 x double>*
store <2 x double> %210, <2 x double>* %212
br label %213
213: ; preds = %199, %201
br i1 %141, label %250, label %214
214: ; preds = %163, %213
%215 = phi i64 [ 0, %163 ], [ %136, %213 ]
br label %242
216: ; preds = %144, %216
%217 = phi i64 [ %239, %216 ], [ 0, %144 ]
%218 = phi double [ %238, %216 ], [ 0.000000e+00, %144 ]
%219 = phi i64 [ %240, %216 ], [ 0, %144 ]
%220 = add nsw i64 %217, %145
%221 = getelementptr inbounds double, double* %4, i64 %220
%222 = load double, double* %221
%223 = tail call double @llvm.fmuladd.f64(double %222, double %118, double %218)
%224 = or i64 %217, 1
%225 = add nsw i64 %224, %145
%226 = getelementptr inbounds double, double* %4, i64 %225
%227 = load double, double* %226
%228 = tail call double @llvm.fmuladd.f64(double %227, double %118, double %223)
%229 = or i64 %217, 2
%230 = add nsw i64 %229, %145
%231 = getelementptr inbounds double, double* %4, i64 %230
%232 = load double, double* %231
%233 = tail call double @llvm.fmuladd.f64(double %232, double %118, double %228)
%234 = or i64 %217, 3
%235 = add nsw i64 %234, %145
%236 = getelementptr inbounds double, double* %4, i64 %235
%237 = load double, double* %236
%238 = tail call double @llvm.fmuladd.f64(double %237, double %118, double %233)
%239 = add nuw nsw i64 %217, 4
%240 = add i64 %219, 4
%241 = icmp eq i64 %240, %133
br i1 %241, label %146, label %216
242: ; preds = %214, %242
%243 = phi i64 [ %248, %242 ], [ %215, %214 ]
%244 = add nsw i64 %243, %164
%245 = getelementptr inbounds double, double* %4, i64 %244
%246 = load double, double* %245
%247 = fsub double %246, %162
store double %247, double* %245
%248 = add nuw nsw i64 %243, 1
%249 = icmp eq i64 %248, %125
br i1 %249, label %250, label %242
250: ; preds = %242, %213, %142, %161
%251 = add nuw nsw i64 %143, 1
%252 = icmp eq i64 %251, %123
br i1 %252, label %381, label %142
253: ; preds = %115
br i1 %11, label %254, label %381
254: ; preds = %253
%255 = sitofp i32 %2 to double
%256 = fadd double %255, -1.000000e+00
%257 = fdiv double 1.000000e+00, %256
%258 = icmp sgt i32 %3, 0
br i1 %258, label %259, label %381
259: ; preds = %254
%260 = icmp sgt i32 %2, 0
%261 = sext i32 %2 to i64
%262 = zext i32 %3 to i64
%263 = zext i32 %2 to i64
%264 = zext i32 %2 to i64
%265 = add nsw i64 %263, -1
%266 = and i64 %263, 4294967294
%267 = add nsw i64 %266, -2
%268 = lshr exact i64 %267, 1
%269 = add nuw i64 %268, 1
%270 = and i64 %263, 3
%271 = icmp ult i64 %265, 3
%272 = and i64 %263, 4294967292
%273 = icmp eq i64 %270, 0
%274 = icmp eq i32 %2, 1
%275 = and i64 %263, 4294967294
%276 = and i64 %269, 1
%277 = icmp eq i64 %267, 0
%278 = and i64 %269, -2
%279 = icmp eq i64 %276, 0
%280 = icmp eq i64 %275, %263
br label %281
281: ; preds = %259, %378
%282 = phi i64 [ 0, %259 ], [ %379, %378 ]
br i1 %260, label %283, label %378
283: ; preds = %281
%284 = mul nsw i64 %282, %261
br i1 %271, label %315, label %285
285: ; preds = %283, %285
%286 = phi i64 [ %312, %285 ], [ 0, %283 ]
%287 = phi double [ %311, %285 ], [ 0.000000e+00, %283 ]
%288 = phi i64 [ %313, %285 ], [ 0, %283 ]
%289 = add nsw i64 %286, %284
%290 = getelementptr inbounds double, double* %4, i64 %289
%291 = load double, double* %290
%292 = fmul double %291, %291
%293 = tail call double @llvm.fmuladd.f64(double %292, double %257, double %287)
%294 = or i64 %286, 1
%295 = add nsw i64 %294, %284
%296 = getelementptr inbounds double, double* %4, i64 %295
%297 = load double, double* %296
%298 = fmul double %297, %297
%299 = tail call double @llvm.fmuladd.f64(double %298, double %257, double %293)
%300 = or i64 %286, 2
%301 = add nsw i64 %300, %284
%302 = getelementptr inbounds double, double* %4, i64 %301
%303 = load double, double* %302
%304 = fmul double %303, %303
%305 = tail call double @llvm.fmuladd.f64(double %304, double %257, double %299)
%306 = or i64 %286, 3
%307 = add nsw i64 %306, %284
%308 = getelementptr inbounds double, double* %4, i64 %307
%309 = load double, double* %308
%310 = fmul double %309, %309
%311 = tail call double @llvm.fmuladd.f64(double %310, double %257, double %305)
%312 = add nuw nsw i64 %286, 4
%313 = add i64 %288, 4
%314 = icmp eq i64 %313, %272
br i1 %314, label %315, label %285
315: ; preds = %285, %283
%316 = phi double [ undef, %283 ], [ %311, %285 ]
%317 = phi i64 [ 0, %283 ], [ %312, %285 ]
%318 = phi double [ 0.000000e+00, %283 ], [ %311, %285 ]
br i1 %273, label %331, label %319
319: ; preds = %315, %319
%320 = phi i64 [ %328, %319 ], [ %317, %315 ]
%321 = phi double [ %327, %319 ], [ %318, %315 ]
%322 = phi i64 [ %329, %319 ], [ 0, %315 ]
%323 = add nsw i64 %320, %284
%324 = getelementptr inbounds double, double* %4, i64 %323
%325 = load double, double* %324
%326 = fmul double %325, %325
%327 = tail call double @llvm.fmuladd.f64(double %326, double %257, double %321)
%328 = add nuw nsw i64 %320, 1
%329 = add i64 %322, 1
%330 = icmp eq i64 %329, %270
br i1 %330, label %331, label %319
331: ; preds = %319, %315
%332 = phi double [ %316, %315 ], [ %327, %319 ]
%333 = tail call double @sqrt(double %332)
br i1 %260, label %334, label %378
334: ; preds = %331
%335 = mul nsw i64 %282, %261
br i1 %274, label %368, label %336
336: ; preds = %334
%337 = insertelement <2 x double> poison, double %333, i64 0
%338 = shufflevector <2 x double> %337, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %277, label %358, label %339
339: ; preds = %336, %339
%340 = phi i64 [ %355, %339 ], [ 0, %336 ]
%341 = phi i64 [ %356, %339 ], [ 0, %336 ]
%342 = add nsw i64 %340, %335
%343 = getelementptr inbounds double, double* %4, i64 %342
%344 = bitcast double* %343 to <2 x double>*
%345 = load <2 x double>, <2 x double>* %344
%346 = fdiv <2 x double> %345, %338
%347 = bitcast double* %343 to <2 x double>*
store <2 x double> %346, <2 x double>* %347
%348 = or i64 %340, 2
%349 = add nsw i64 %348, %335
%350 = getelementptr inbounds double, double* %4, i64 %349
%351 = bitcast double* %350 to <2 x double>*
%352 = load <2 x double>, <2 x double>* %351
%353 = fdiv <2 x double> %352, %338
%354 = bitcast double* %350 to <2 x double>*
store <2 x double> %353, <2 x double>* %354
%355 = add nuw i64 %340, 4
%356 = add nuw i64 %341, 2
%357 = icmp eq i64 %356, %278
br i1 %357, label %358, label %339
358: ; preds = %339, %336
%359 = phi i64 [ 0, %336 ], [ %355, %339 ]
br i1 %279, label %367, label %360
360: ; preds = %358
%361 = add nsw i64 %359, %335
%362 = getelementptr inbounds double, double* %4, i64 %361
%363 = bitcast double* %362 to <2 x double>*
%364 = load <2 x double>, <2 x double>* %363
%365 = fdiv <2 x double> %364, %338
%366 = bitcast double* %362 to <2 x double>*
store <2 x double> %365, <2 x double>* %366
br label %367
367: ; preds = %358, %360
br i1 %280, label %378, label %368
368: ; preds = %334, %367
%369 = phi i64 [ 0, %334 ], [ %275, %367 ]
br label %370
370: ; preds = %368, %370
%371 = phi i64 [ %376, %370 ], [ %369, %368 ]
%372 = add nsw i64 %371, %335
%373 = getelementptr inbounds double, double* %4, i64 %372
%374 = load double, double* %373
%375 = fdiv double %374, %333
store double %375, double* %373
%376 = add nuw nsw i64 %371, 1
%377 = icmp eq i64 %376, %264
br i1 %377, label %378, label %370
378: ; preds = %370, %367, %281, %331
%379 = add nuw nsw i64 %282, 1
%380 = icmp eq i64 %379, %262
br i1 %380, label %381, label %281
381: ; preds = %378, %250, %112, %254, %116, %13, %253, %5
ret i32 0
}