#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 */
scalar_t__ FP_NORMAL;
double NAN;
scalar_t__ fpclassify(double);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *);
double sqrt(double);
int /*<<< orphan*/ stderr;
void standardise_data(double **x, double *y, int N, int p, double *xm,
                      double *xs, double *ym, double *ys) {
  int i, j;
  double delta;
  *ym = *ys = 0;
  for (j = 0; j < p; ++j)
    xm[j] = xs[j] = 0.0;
  for (i = 0; i < N; ++i) {
    delta = y[i] - *ym;
    (*ym) += delta / (i + 1);
    (*ys) += delta * (y[i] - *ym);
    for (j = 0; j < p; ++j) {
      delta = x[i][j] - xm[j];
      xm[j] += delta / (i + 1);
      xs[j] += delta * (x[i][j] - xm[j]);
    }
  }
  if (N < 2) {
    fprintf(stderr, "Warning: fewer than two instances, cannot compute "
                    "standard deviation of values.\n");
    *ys = NAN;
    for (j = 0; j < p; ++j)
      xs[j] = NAN;
  } else {
    delta = 1.0 / (double)(N);
    *ys = sqrt(delta * *ys);
    for (j = 0; j < p; ++j)
      xs[j] = sqrt(delta * xs[j]);
  }
  for (i = 0; i < N; ++i) {
    for (j = 0; j < p; ++j) {
      x[i][j] = (x[i][j] - xm[j]) / xs[j];
      if (fpclassify(x[i][j]) != FP_NORMAL)
        x[i][j] = 0;
    }
    y[i] = (y[i] - *ym) / *ys;
  }
}


x86 assembly (clang, -O0):
<standardise_data>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11ad <standardise_data+0x7d>
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1173 <standardise_data+0x43>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    12e8 <standardise_data+0x1b8>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
subsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    0x10(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x10(%rbp),%rax
subsd  (%rax),%xmm2
mov    0x18(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    12d5 <standardise_data+0x1a5>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    122e <standardise_data+0xfe>
jmp    12da <standardise_data+0x1aa>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11b4 <standardise_data+0x84>
cmpl   $0x2,-0x14(%rbp)
jge    1358 <standardise_data+0x228>
mov    0x2cd7(%rip),%rax
mov    (%rax),%edi
lea    0xd06(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2cea(%rip),%rax
movsd  (%rax),%xmm0
mov    0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1353 <standardise_data+0x223>
mov    0x2cc4(%rip),%rax
movsd  (%rax),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1321 <standardise_data+0x1f1>
jmp    13d2 <standardise_data+0x2a2>
cvtsi2sdl -0x14(%rbp),%xmm1
movsd  0xc9b(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    0x18(%rbp),%rax
mulsd  (%rax),%xmm0
call   1050 <sqrt@plt>
mov    0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    13cd <standardise_data+0x29d>
movsd  -0x38(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
call   1050 <sqrt@plt>
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    138f <standardise_data+0x25f>
jmp    13d2 <standardise_data+0x2a2>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    14c9 <standardise_data+0x399>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1491 <standardise_data+0x361>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1040 <fpclassify@plt>
mov    0x2b8b(%rip),%rcx
cmp    (%rcx),%rax
je     147e <standardise_data+0x34e>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1483 <standardise_data+0x353>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    13ec <standardise_data+0x2bc>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
subsd  (%rax),%xmm0
mov    0x18(%rbp),%rax
divsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    13d9 <standardise_data+0x2a9>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<standardise_data>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r9,%r15
mov    %r8,%r14
mov    %ecx,%r12d
mov    %rsi,0x30(%rsp)
mov    %rdi,0x10(%rsp)
mov    0x78(%rsp),%rax
mov    0x70(%rsp),%rcx
movq   $0x0,(%rax)
movq   $0x0,(%rcx)
mov    %r12d,%r13d
test   %r12d,%r12d
jle    1339 <standardise_data+0x209>
cmp    $0x4,%r12d
jb     1195 <standardise_data+0x65>
lea    (%r14,%r13,8),%rax
cmp    %r15,%rax
jbe    1234 <standardise_data+0x104>
lea    (%r15,%r13,8),%rax
cmp    %r14,%rax
jbe    1234 <standardise_data+0x104>
xor    %eax,%eax
mov    %rax,%rcx
not    %rcx
add    %r13,%rcx
mov    %r13,%rsi
and    $0x3,%rsi
je     11ca <standardise_data+0x9a>
nopl   0x0(%rax)
movq   $0x0,(%r15,%rax,8)
movq   $0x0,(%r14,%rax,8)
add    $0x1,%rax
add    $0xffffffffffffffff,%rsi
jne    11b0 <standardise_data+0x80>
cmp    $0x3,%rcx
jb     1339 <standardise_data+0x209>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movq   $0x0,(%r15,%rax,8)
movq   $0x0,(%r14,%rax,8)
movq   $0x0,0x8(%r15,%rax,8)
movq   $0x0,0x8(%r14,%rax,8)
movq   $0x0,0x10(%r15,%rax,8)
movq   $0x0,0x10(%r14,%rax,8)
movq   $0x0,0x18(%r15,%rax,8)
movq   $0x0,0x18(%r14,%rax,8)
add    $0x4,%rax
cmp    %rax,%r13
jne    11e0 <standardise_data+0xb0>
jmp    1339 <standardise_data+0x209>
mov    %r13d,%eax
and    $0xfffffffc,%eax
lea    -0x4(%rax),%rsi
mov    %rsi,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
mov    %edi,%ecx
and    $0x3,%ecx
cmp    $0xc,%rsi
jae    125b <standardise_data+0x12b>
xor    %esi,%esi
jmp    12e8 <standardise_data+0x1b8>
and    $0xfffffffffffffffc,%rdi
xor    %esi,%esi
xorpd  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
nop
movupd %xmm0,(%r15,%rsi,8)
movupd %xmm0,0x10(%r15,%rsi,8)
movupd %xmm0,(%r14,%rsi,8)
movupd %xmm0,0x10(%r14,%rsi,8)
movupd %xmm0,0x20(%r15,%rsi,8)
movupd %xmm0,0x30(%r15,%rsi,8)
movupd %xmm0,0x20(%r14,%rsi,8)
movupd %xmm0,0x30(%r14,%rsi,8)
movupd %xmm0,0x40(%r15,%rsi,8)
movupd %xmm0,0x50(%r15,%rsi,8)
movupd %xmm0,0x40(%r14,%rsi,8)
movupd %xmm0,0x50(%r14,%rsi,8)
movupd %xmm0,0x60(%r15,%rsi,8)
movupd %xmm0,0x70(%r15,%rsi,8)
movupd %xmm0,0x60(%r14,%rsi,8)
movupd %xmm0,0x70(%r14,%rsi,8)
add    $0x10,%rsi
add    $0xfffffffffffffffc,%rdi
jne    1270 <standardise_data+0x140>
test   %rcx,%rcx
je     1330 <standardise_data+0x200>
lea    (%r14,%rsi,8),%rbp
add    $0x10,%rbp
lea    (%r15,%rsi,8),%rsi
add    $0x10,%rsi
shl    $0x5,%rcx
xor    %edi,%edi
xorpd  %xmm0,%xmm0
nopw   0x0(%rax,%rax,1)
movupd %xmm0,-0x10(%rsi,%rdi,1)
movupd %xmm0,(%rsi,%rdi,1)
movupd %xmm0,-0x10(%rbp,%rdi,1)
movupd %xmm0,0x0(%rbp,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%rcx
jne    1310 <standardise_data+0x1e0>
cmp    %r13,%rax
jne    1197 <standardise_data+0x67>
mov    %edx,%eax
mov    %rax,0x28(%rsp)
test   %edx,%edx
mov    %edx,0xc(%rsp)
jle    14dc <standardise_data+0x3ac>
lea    (%r15,%r13,8),%r8
cmp    %r14,%r8
seta   %al
lea    (%r14,%r13,8),%r9
cmp    %r15,%r9
seta   %r10b
and    %al,%r10b
mov    %r13d,%ebp
and    $0xfffffffe,%ebp
xor    %edi,%edi
jmp    137b <standardise_data+0x24b>
xchg   %ax,%ax
cmp    0x28(%rsp),%rdi
je     14cf <standardise_data+0x39f>
mov    %rdi,%rax
add    $0x1,%rdi
mov    0x30(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm1
mov    0x70(%rsp),%rdx
movsd  (%rdx),%xmm2
subsd  %xmm2,%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %edi,%xmm0
movapd %xmm1,%xmm3
divsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%rdx)
movsd  (%rcx,%rax,8),%xmm2
subsd  %xmm3,%xmm2
mulsd  %xmm1,%xmm2
mov    0x78(%rsp),%rcx
addsd  (%rcx),%xmm2
movsd  %xmm2,(%rcx)
test   %r12d,%r12d
jle    1370 <standardise_data+0x240>
mov    0x10(%rsp),%rcx
mov    (%rcx,%rax,8),%rax
cmp    $0x1,%r12d
jne    1430 <standardise_data+0x300>
xor    %esi,%esi
movsd  (%rax,%rsi,8),%xmm1
movsd  (%r14,%rsi,8),%xmm2
subsd  %xmm2,%xmm1
movapd %xmm1,%xmm3
divsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%r14,%rsi,8)
movsd  (%rax,%rsi,8),%xmm2
subsd  %xmm3,%xmm2
mulsd  %xmm1,%xmm2
addsd  (%r15,%rsi,8),%xmm2
movsd  %xmm2,(%r15,%rsi,8)
add    $0x1,%rsi
cmp    %rsi,%r13
jne    13e0 <standardise_data+0x2b0>
jmp    1370 <standardise_data+0x240>
nopl   0x0(%rax,%rax,1)
lea    (%rax,%r13,8),%rcx
cmp    %r14,%rcx
seta   %bl
cmp    %r9,%rax
setb   %dl
and    %bl,%dl
or     %r10b,%dl
cmp    %r15,%rcx
seta   %bl
cmp    %r8,%rax
setb   %cl
test   %dl,%dl
jne    14c8 <standardise_data+0x398>
mov    $0x0,%esi
and    %cl,%bl
jne    13e0 <standardise_data+0x2b0>
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
movupd (%rax,%rsi,8),%xmm2
movupd (%r14,%rsi,8),%xmm3
subpd  %xmm3,%xmm2
movapd %xmm2,%xmm4
divpd  %xmm1,%xmm4
addpd  %xmm3,%xmm4
movupd %xmm4,(%r14,%rsi,8)
movupd (%rax,%rsi,8),%xmm3
subpd  %xmm4,%xmm3
movupd (%r15,%rsi,8),%xmm4
mulpd  %xmm2,%xmm3
addpd  %xmm4,%xmm3
movupd %xmm3,(%r15,%rsi,8)
add    $0x2,%rsi
cmp    %rsi,%rbp
jne    1470 <standardise_data+0x340>
mov    %rbp,%rsi
cmp    %r13,%rbp
jne    13e0 <standardise_data+0x2b0>
jmp    1370 <standardise_data+0x240>
xor    %esi,%esi
jmp    13e0 <standardise_data+0x2b0>
mov    0xc(%rsp),%eax
cmp    $0x1,%eax
jg     15a9 <standardise_data+0x479>
mov    0x2aed(%rip),%rax
mov    (%rax),%edi
lea    0xb1c(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2b00(%rip),%rax
movsd  (%rax),%xmm0
mov    0x78(%rsp),%rcx
movsd  %xmm0,(%rcx)
test   %r12d,%r12d
jle    180d <standardise_data+0x6dd>
cmp    $0x4,%r12d
jb     1532 <standardise_data+0x402>
lea    0x8(%rax),%rcx
cmp    %rcx,%r15
jae    1667 <standardise_data+0x537>
lea    (%r15,%r13,8),%rcx
cmp    0x2acc(%rip),%rcx
jbe    1667 <standardise_data+0x537>
xor    %ecx,%ecx
mov    %rcx,%rdx
not    %rdx
add    %r13,%rdx
mov    %r13,%rsi
and    $0x3,%rsi
je     1564 <standardise_data+0x434>
cs nopw 0x0(%rax,%rax,1)
movsd  (%rax),%xmm0
movsd  %xmm0,(%r15,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    1550 <standardise_data+0x420>
cmp    $0x3,%rdx
jb     180d <standardise_data+0x6dd>
xchg   %ax,%ax
movsd  (%rax),%xmm0
movsd  %xmm0,(%r15,%rcx,8)
movsd  (%rax),%xmm0
movsd  %xmm0,0x8(%r15,%rcx,8)
movsd  (%rax),%xmm0
movsd  %xmm0,0x10(%r15,%rcx,8)
movsd  (%rax),%xmm0
movsd  %xmm0,0x18(%r15,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%r13
jne    1570 <standardise_data+0x440>
jmp    180d <standardise_data+0x6dd>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0xa48(%rip),%xmm2
divsd  %xmm0,%xmm2
mov    0x78(%rsp),%rax
movsd  (%rax),%xmm0
mulsd  %xmm2,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
movsd  %xmm2,(%rsp)
jb     15de <standardise_data+0x4ae>
sqrtsd %xmm0,%xmm0
jmp    15e8 <standardise_data+0x4b8>
call   1050 <sqrt@plt>
movsd  (%rsp),%xmm2
mov    0x78(%rsp),%rax
movsd  %xmm0,(%rax)
test   %r12d,%r12d
jle    180d <standardise_data+0x6dd>
lea    -0x1(%r13),%rax
mov    %r13d,%ecx
and    $0x3,%ecx
mov    %rcx,0x20(%rsp)
cmp    $0x3,%rax
jae    1692 <standardise_data+0x562>
xor    %ebx,%ebx
cmpq   $0x0,0x20(%rsp)
je     180d <standardise_data+0x6dd>
lea    (%r15,%rbx,8),%rbx
xor    %ebp,%ebp
xorpd  %xmm1,%xmm1
jmp    1652 <standardise_data+0x522>
nopl   (%rax)
call   1050 <sqrt@plt>
xorpd  %xmm1,%xmm1
movsd  (%rsp),%xmm2
movsd  %xmm0,(%rbx,%rbp,8)
add    $0x1,%rbp
cmp    %rbp,0x20(%rsp)
je     180d <standardise_data+0x6dd>
movsd  (%rbx,%rbp,8),%xmm0
mulsd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jb     1630 <standardise_data+0x500>
sqrtsd %xmm0,%xmm0
jmp    163e <standardise_data+0x50e>
mov    %r13d,%ecx
and    $0xfffffffc,%ecx
lea    -0x4(%rcx),%rdi
mov    %rdi,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %esi,%edx
and    $0x3,%edx
cmp    $0xc,%rdi
jae    1773 <standardise_data+0x643>
xor    %edi,%edi
jmp    17d1 <standardise_data+0x6a1>
mov    %r13d,%ebp
and    $0xfffffffc,%ebp
xor    %ebx,%ebx
xorpd  %xmm1,%xmm1
jmp    16c2 <standardise_data+0x592>
call   1050 <sqrt@plt>
xorpd  %xmm1,%xmm1
movsd  (%rsp),%xmm2
movsd  %xmm0,0x18(%r15,%rbx,8)
add    $0x4,%rbx
cmp    %rbx,%rbp
je     1615 <standardise_data+0x4e5>
movsd  (%r15,%rbx,8),%xmm0
mulsd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jb     16e0 <standardise_data+0x5b0>
sqrtsd %xmm0,%xmm0
jmp    16ee <standardise_data+0x5be>
nopl   0x0(%rax,%rax,1)
call   1050 <sqrt@plt>
xorpd  %xmm1,%xmm1
movsd  (%rsp),%xmm2
movsd  %xmm0,(%r15,%rbx,8)
movsd  0x8(%r15,%rbx,8),%xmm0
mulsd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jb     1710 <standardise_data+0x5e0>
sqrtsd %xmm0,%xmm0
jmp    171e <standardise_data+0x5ee>
nopl   0x0(%rax,%rax,1)
call   1050 <sqrt@plt>
xorpd  %xmm1,%xmm1
movsd  (%rsp),%xmm2
movsd  %xmm0,0x8(%r15,%rbx,8)
movsd  0x10(%r15,%rbx,8),%xmm0
mulsd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jb     1740 <standardise_data+0x610>
sqrtsd %xmm0,%xmm0
jmp    174e <standardise_data+0x61e>
nopl   0x0(%rax)
call   1050 <sqrt@plt>
xorpd  %xmm1,%xmm1
movsd  (%rsp),%xmm2
movsd  %xmm0,0x10(%r15,%rbx,8)
movsd  0x18(%r15,%rbx,8),%xmm0
mulsd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jb     16a0 <standardise_data+0x570>
sqrtsd %xmm0,%xmm0
jmp    16ae <standardise_data+0x57e>
and    $0xfffffffffffffffc,%rsi
movsd  (%rax),%xmm0
unpcklpd %xmm0,%xmm0
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movupd %xmm0,(%r15,%rdi,8)
movupd %xmm0,0x10(%r15,%rdi,8)
movupd %xmm0,0x20(%r15,%rdi,8)
movupd %xmm0,0x30(%r15,%rdi,8)
movupd %xmm0,0x40(%r15,%rdi,8)
movupd %xmm0,0x50(%r15,%rdi,8)
movupd %xmm0,0x60(%r15,%rdi,8)
movupd %xmm0,0x70(%r15,%rdi,8)
add    $0x10,%rdi
add    $0xfffffffffffffffc,%rsi
jne    1790 <standardise_data+0x660>
test   %rdx,%rdx
je     1804 <standardise_data+0x6d4>
movsd  (%rax),%xmm0
unpcklpd %xmm0,%xmm0
lea    (%r15,%rdi,8),%rsi
add    $0x10,%rsi
shl    $0x5,%rdx
xor    %edi,%edi
nopl   0x0(%rax)
movupd %xmm0,-0x10(%rsi,%rdi,1)
movupd %xmm0,(%rsi,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%rdx
jne    17f0 <standardise_data+0x6c0>
cmp    %r13,%rcx
jne    1534 <standardise_data+0x404>
cmpl   $0x0,0xc(%rsp)
jle    18b7 <standardise_data+0x787>
xor    %ebp,%ebp
mov    0x27c7(%rip),%rbx
mov    %r12d,0x1c(%rsp)
jmp    1861 <standardise_data+0x731>
nopl   0x0(%rax,%rax,1)
mov    0x30(%rsp),%rax
movsd  (%rax,%rbp,8),%xmm0
mov    0x70(%rsp),%rcx
subsd  (%rcx),%xmm0
mov    0x78(%rsp),%rcx
divsd  (%rcx),%xmm0
movsd  %xmm0,(%rax,%rbp,8)
add    $0x1,%rbp
cmp    0x28(%rsp),%rbp
mov    0x1c(%rsp),%r12d
je     18b7 <standardise_data+0x787>
test   %r12d,%r12d
jle    1830 <standardise_data+0x700>
xor    %r12d,%r12d
jmp    1879 <standardise_data+0x749>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r12
cmp    %r12,%r13
je     1830 <standardise_data+0x700>
mov    0x10(%rsp),%rax
mov    (%rax,%rbp,8),%rax
movsd  (%rax,%r12,8),%xmm0
subsd  (%r14,%r12,8),%xmm0
divsd  (%r15,%r12,8),%xmm0
movsd  %xmm0,(%rax,%r12,8)
call   1040 <fpclassify@plt>
cmp    (%rbx),%rax
je     1870 <standardise_data+0x740>
mov    0x10(%rsp),%rax
mov    (%rax,%rbp,8),%rax
movq   $0x0,(%rax,%r12,8)
jmp    1870 <standardise_data+0x740>
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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