#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 KineticEnergy(int, double *, double *, double *, double *);
double *calloc(int, int);
int /*<<< orphan*/ free(double *);
double sqrt(double);
void VerletUpdate(int nAtom, double dt, double *mass, double temp, double *x,
                  double *y, double *z, double *px, double *py, double *pz,
                  double *vx, double *vy, double *vz, double *fx, double *fy,
                  double *fz) {
  int n;
  double *cx, *cy, *cz;
  cx = calloc(nAtom, sizeof(double));
  cy = calloc(nAtom, sizeof(double));
  cz = calloc(nAtom, sizeof(double));
  for (n = 0; n < nAtom; n++) {
    cx[n] = x[n];
    cy[n] = y[n];
    cz[n] = z[n];
  }
  for (n = 0; n < nAtom; n++) {
    vx[n] = (cx[n] - px[n]) / dt + 0.5 * dt * fx[n] / mass[n];
    vy[n] = (cy[n] - py[n]) / dt + 0.5 * dt * fy[n] / mass[n];
    vz[n] = (cz[n] - pz[n]) / dt + 0.5 * dt * fz[n] / mass[n];
  }
  if (temp > 0.0) {
    double alpha = sqrt(1.5 * 0.0083144621 * temp * nAtom /
                        KineticEnergy(nAtom, mass, vx, vy, vz));
    for (n = 0; n < nAtom; n++) {
      vx[n] *= alpha;
      vy[n] *= alpha;
      vz[n] *= alpha;
    }
  }
  for (n = 0; n < nAtom; n++) {
    x[n] = cx[n] + dt * vx[n] + 0.5 * dt * dt * fx[n] / mass[n];
    y[n] = cy[n] + dt * vy[n] + 0.5 * dt * dt * fy[n] / mass[n];
    z[n] = cz[n] + dt * vz[n] + 0.5 * dt * dt * fz[n] / mass[n];
  }
  for (n = 0; n < nAtom; n++) {
    px[n] = cx[n];
    py[n] = cy[n];
    pz[n] = cz[n];
  }
  free(cx);
  free(cy);
  free(cz);
}


x86 assembly (clang, -O0):
<VerletUpdate>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    0x48(%rbp),%rax
mov    0x40(%rbp),%rax
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
mov    %rcx,-0x30(%rbp)
mov    %r8,-0x38(%rbp)
mov    %r9,-0x40(%rbp)
mov    -0x4(%rbp),%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x50(%rbp)
mov    -0x4(%rbp),%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x58(%rbp)
mov    -0x4(%rbp),%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x60(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    122b <VerletUpdate+0xeb>
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x58(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x60(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    11c3 <VerletUpdate+0x83>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1351 <VerletUpdate+0x211>
mov    -0x50(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  0xda3(%rip),%xmm1
mulsd  -0x10(%rbp),%xmm1
mov    0x38(%rbp),%rax
movslq -0x44(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm1
addsd  %xmm1,%xmm0
mov    0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  0xd4c(%rip),%xmm1
mulsd  -0x10(%rbp),%xmm1
mov    0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm1
addsd  %xmm1,%xmm0
mov    0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x60(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  0xcf5(%rip),%xmm1
mulsd  -0x10(%rbp),%xmm1
mov    0x48(%rbp),%rax
movslq -0x44(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm1
addsd  %xmm1,%xmm0
mov    0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1232 <VerletUpdate+0xf2>
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1417 <VerletUpdate+0x2d7>
movsd  0xc95(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0x4(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x4(%rbp),%edi
mov    -0x18(%rbp),%rsi
mov    0x20(%rbp),%rdx
mov    0x28(%rbp),%rcx
mov    0x30(%rbp),%r8
call   1030 <KineticEnergy@plt>
movaps %xmm0,%xmm1
movsd  -0x70(%rbp),%xmm0
divsd  %xmm1,%xmm0
call   1060 <sqrt@plt>
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1412 <VerletUpdate+0x2d2>
movsd  -0x68(%rbp),%xmm0
mov    0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x68(%rbp),%xmm0
mov    0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x68(%rbp),%xmm0
mov    0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    13b3 <VerletUpdate+0x273>
jmp    1417 <VerletUpdate+0x2d7>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1564 <VerletUpdate+0x424>
mov    -0x50(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x10(%rbp),%xmm0
mov    0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xbaf(%rip),%xmm1
mulsd  -0x10(%rbp),%xmm1
mulsd  -0x10(%rbp),%xmm1
mov    0x38(%rbp),%rax
movslq -0x44(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm1
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x10(%rbp),%xmm0
mov    0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xb4b(%rip),%xmm1
mulsd  -0x10(%rbp),%xmm1
mulsd  -0x10(%rbp),%xmm1
mov    0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm1
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x60(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x10(%rbp),%xmm0
mov    0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xae7(%rip),%xmm1
mulsd  -0x10(%rbp),%xmm1
mulsd  -0x10(%rbp),%xmm1
mov    0x48(%rbp),%rax
movslq -0x44(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm1
addsd  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    141e <VerletUpdate+0x2de>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    15d3 <VerletUpdate+0x493>
mov    -0x50(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x60(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    156b <VerletUpdate+0x42b>
mov    -0x50(%rbp),%rdi
call   1040 <free@plt>
mov    -0x58(%rbp),%rdi
call   1040 <free@plt>
mov    -0x60(%rbp),%rdi
call   1040 <free@plt>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<VerletUpdate>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %r9,0x40(%rsp)
mov    %r8,0x78(%rsp)
mov    %rcx,0x30(%rsp)
mov    %rdx,0x70(%rsp)
movsd  %xmm1,0x60(%rsp)
mov    %rsi,%rbp
movsd  %xmm0,0x50(%rsp)
mov    %edi,%r12d
mov    0xf8(%rsp),%rax
mov    %rax,0x68(%rsp)
mov    0xf0(%rsp),%rax
mov    %rax,0x58(%rsp)
mov    0xe8(%rsp),%rbx
mov    0xe0(%rsp),%r13
mov    0xd8(%rsp),%r14
mov    0xd0(%rsp),%r15
mov    0xc8(%rsp),%rax
mov    %rax,0x38(%rsp)
mov    0xc0(%rsp),%rax
mov    %rax,0x48(%rsp)
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,0x28(%rsp)
mov    %r12d,%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,0x18(%rsp)
mov    %r12d,%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    0x58(%rsp),%r9
movsd  0x50(%rsp),%xmm3
mov    %r12d,%r10d
mov    0x18(%rsp),%r11
mov    0x28(%rsp),%rdi
mov    %r12d,%ecx
mov    %rcx,0x8(%rsp)
test   %r12d,%r12d
mov    %r12d,0x14(%rsp)
mov    %rax,0x20(%rsp)
jle    1350 <VerletUpdate+0x200>
mov    0x8(%rsp),%rcx
lea    0x0(,%rcx,8),%r12
mov    0x70(%rsp),%rsi
mov    %r12,%rdx
call   1060 <memcpy@plt>
mov    0x18(%rsp),%rdi
mov    0x30(%rsp),%rsi
mov    %r12,%rdx
call   1060 <memcpy@plt>
mov    0x20(%rsp),%rdi
mov    0x78(%rsp),%rsi
mov    %r12,%rdx
call   1060 <memcpy@plt>
mov    0x58(%rsp),%r9
mov    0x68(%rsp),%rsi
movsd  0x50(%rsp),%xmm3
mov    0x14(%rsp),%r10d
mov    0x20(%rsp),%rax
mov    0x18(%rsp),%r11
mov    0x28(%rsp),%rdi
mov    0x48(%rsp),%r8
movsd  0xd55(%rip),%xmm0
mulsd  %xmm3,%xmm0
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%rcx,1),%xmm1
mov    0x40(%rsp),%rdx
subsd  (%rdx,%rcx,1),%xmm1
divsd  %xmm3,%xmm1
movsd  (%rbx,%rcx,1),%xmm2
mulsd  %xmm0,%xmm2
divsd  0x0(%rbp,%rcx,1),%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%r15,%rcx,1)
movsd  (%r11,%rcx,1),%xmm1
subsd  (%r8,%rcx,1),%xmm1
divsd  %xmm3,%xmm1
movsd  (%r9,%rcx,1),%xmm2
mulsd  %xmm0,%xmm2
divsd  0x0(%rbp,%rcx,1),%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%r14,%rcx,1)
movsd  (%rax,%rcx,1),%xmm1
mov    0x38(%rsp),%rdx
subsd  (%rdx,%rcx,1),%xmm1
divsd  %xmm3,%xmm1
movsd  (%rsi,%rcx,1),%xmm2
mulsd  %xmm0,%xmm2
divsd  0x0(%rbp,%rcx,1),%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,0x0(%r13,%rcx,1)
add    $0x8,%rcx
cmp    %rcx,%r12
jne    12c0 <VerletUpdate+0x170>
mov    %rbx,0x80(%rsp)
xorpd  %xmm1,%xmm1
movsd  0x60(%rsp),%xmm0
ucomisd %xmm1,%xmm0
mov    0x30(%rsp),%r12
mov    0x8(%rsp),%rbx
jbe    1642 <VerletUpdate+0x4f2>
mulsd  0xc8a(%rip),%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %r10d,%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x60(%rsp)
mov    %r10d,%edi
mov    %rbp,%rsi
mov    %r15,%rdx
mov    %r14,%rcx
mov    %r13,%r8
call   1030 <KineticEnergy@plt>
movsd  0x60(%rsp),%xmm1
divsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
ucomisd 0xc56(%rip),%xmm1
jb     13c2 <VerletUpdate+0x272>
sqrtsd %xmm0,%xmm0
jmp    13c7 <VerletUpdate+0x277>
call   1070 <sqrt@plt>
mov    0x28(%rsp),%rdi
mov    0x18(%rsp),%r11
mov    0x20(%rsp),%rax
movsd  0x50(%rsp),%xmm3
mov    0x14(%rsp),%ecx
test   %ecx,%ecx
mov    0x8(%rsp),%rbx
mov    0x30(%rsp),%r12
mov    0x58(%rsp),%r9
jle    1978 <VerletUpdate+0x828>
cmp    $0x1,%ecx
jne    14d6 <VerletUpdate+0x386>
xor    %ecx,%ecx
mov    %rcx,%rdx
not    %rdx
test   $0x1,%bl
je     1443 <VerletUpdate+0x2f3>
movsd  (%r15,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%r15,%rcx,8)
movsd  (%r14,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%r14,%rcx,8)
movsd  0x0(%r13,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%r13,%rcx,8)
or     $0x1,%rcx
add    %rbx,%rdx
mov    0x14(%rsp),%r10d
je     1642 <VerletUpdate+0x4f2>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%r15,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%r15,%rcx,8)
movsd  (%r14,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%r14,%rcx,8)
movsd  0x0(%r13,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%r13,%rcx,8)
movsd  0x8(%r15,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%r15,%rcx,8)
movsd  0x8(%r14,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%r14,%rcx,8)
movsd  0x8(%r13,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%r13,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%rbx
jne    1460 <VerletUpdate+0x310>
jmp    1642 <VerletUpdate+0x4f2>
lea    (%r14,%rbx,8),%rax
cmp    %r15,%rax
seta   %r8b
lea    (%r15,%rbx,8),%rdx
cmp    %r14,%rdx
seta   %r12b
lea    0x0(,%rbx,8),%rdi
add    %r13,%rdi
cmp    %r15,%rdi
seta   %bl
cmp    %r13,%rdx
seta   %sil
cmp    %r14,%rdi
seta   %dl
cmp    %r13,%rax
seta   %dil
xor    %ecx,%ecx
test   %r12b,%r8b
jne    19a4 <VerletUpdate+0x854>
and    %sil,%bl
jne    19a4 <VerletUpdate+0x854>
and    %dil,%dl
mov    0x28(%rsp),%rdi
mov    0x18(%rsp),%r11
mov    0x20(%rsp),%rax
mov    0x8(%rsp),%rbx
mov    0x30(%rsp),%r12
jne    1402 <VerletUpdate+0x2b2>
mov    %ebx,%ecx
and    $0xfffffffe,%ecx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x2(%rcx),%rsi
mov    %rsi,%r9
shr    %r9
add    $0x1,%r9
test   %rsi,%rsi
je     19d7 <VerletUpdate+0x887>
mov    %r9,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd (%r15,%rdx,8),%xmm2
mulpd  %xmm1,%xmm2
movupd %xmm2,(%r15,%rdx,8)
movupd (%r14,%rdx,8),%xmm2
mulpd  %xmm1,%xmm2
movupd %xmm2,(%r14,%rdx,8)
movupd 0x0(%r13,%rdx,8),%xmm2
mulpd  %xmm1,%xmm2
movupd %xmm2,0x0(%r13,%rdx,8)
movupd 0x10(%r15,%rdx,8),%xmm2
mulpd  %xmm1,%xmm2
movupd %xmm2,0x10(%r15,%rdx,8)
movupd 0x10(%r14,%rdx,8),%xmm2
mulpd  %xmm1,%xmm2
movupd %xmm2,0x10(%r14,%rdx,8)
movupd 0x10(%r13,%rdx,8),%xmm2
mulpd  %xmm1,%xmm2
movupd %xmm2,0x10(%r13,%rdx,8)
add    $0x4,%rdx
add    $0xfffffffffffffffe,%rsi
jne    1580 <VerletUpdate+0x430>
test   $0x1,%r9b
je     162a <VerletUpdate+0x4da>
movupd (%r15,%rdx,8),%xmm2
mulpd  %xmm1,%xmm2
movupd %xmm2,(%r15,%rdx,8)
movupd (%r14,%rdx,8),%xmm2
mulpd  %xmm1,%xmm2
movupd %xmm2,(%r14,%rdx,8)
movupd 0x0(%r13,%rdx,8),%xmm2
mulpd  %xmm1,%xmm2
movupd %xmm2,0x0(%r13,%rdx,8)
mov    0x8(%rsp),%rbx
cmp    %rbx,%rcx
mov    0x14(%rsp),%r10d
mov    0x58(%rsp),%r9
jne    1402 <VerletUpdate+0x2b2>
test   %r10d,%r10d
mov    0x78(%rsp),%rsi
jle    1978 <VerletUpdate+0x828>
movsd  0x9a8(%rip),%xmm0
mulsd  %xmm3,%xmm0
mulsd  %xmm3,%xmm0
lea    0x0(,%rbx,8),%rdx
xor    %ecx,%ecx
mov    %r9,%r10
mov    0x70(%rsp),%r9
mov    %r12,%r8
mov    0x68(%rsp),%r12
mov    0x80(%rsp),%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%r15,%rcx,1),%xmm1
mulsd  %xmm3,%xmm1
addsd  (%rdi,%rcx,1),%xmm1
movsd  (%rbx,%rcx,1),%xmm2
mulsd  %xmm0,%xmm2
divsd  0x0(%rbp,%rcx,1),%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%r9,%rcx,1)
movsd  (%r14,%rcx,1),%xmm1
mulsd  %xmm3,%xmm1
addsd  (%r11,%rcx,1),%xmm1
movsd  (%r10,%rcx,1),%xmm2
mulsd  %xmm0,%xmm2
divsd  0x0(%rbp,%rcx,1),%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%r8,%rcx,1)
movsd  0x0(%r13,%rcx,1),%xmm1
mulsd  %xmm3,%xmm1
addsd  (%rax,%rcx,1),%xmm1
movsd  (%r12,%rcx,1),%xmm2
mulsd  %xmm0,%xmm2
divsd  0x0(%rbp,%rcx,1),%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%rsi,%rcx,1)
add    $0x8,%rcx
cmp    %rcx,%rdx
jne    1690 <VerletUpdate+0x540>
mov    0x14(%rsp),%ecx
test   %ecx,%ecx
mov    0x8(%rsp),%rsi
mov    0x48(%rsp),%r8
jle    1978 <VerletUpdate+0x828>
cmp    $0x4,%ecx
jae    17e4 <VerletUpdate+0x694>
xor    %ebp,%ebp
mov    %rbp,%rcx
not    %rcx
test   $0x1,%sil
je     1773 <VerletUpdate+0x623>
movsd  (%rdi,%rbp,8),%xmm0
mov    0x40(%rsp),%rdx
movsd  %xmm0,(%rdx,%rbp,8)
movsd  (%r11,%rbp,8),%xmm0
movsd  %xmm0,(%r8,%rbp,8)
movsd  (%rax,%rbp,8),%xmm0
mov    0x38(%rsp),%rdx
movsd  %xmm0,(%rdx,%rbp,8)
or     $0x1,%rbp
add    %rsi,%rcx
je     1978 <VerletUpdate+0x828>
mov    0x40(%rsp),%rcx
mov    0x38(%rsp),%rdx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rdi,%rbp,8),%xmm0
movsd  %xmm0,(%rcx,%rbp,8)
movsd  (%r11,%rbp,8),%xmm0
movsd  %xmm0,(%r8,%rbp,8)
movsd  (%rax,%rbp,8),%xmm0
movsd  %xmm0,(%rdx,%rbp,8)
movsd  0x8(%rdi,%rbp,8),%xmm0
movsd  %xmm0,0x8(%rcx,%rbp,8)
movsd  0x8(%r11,%rbp,8),%xmm0
movsd  %xmm0,0x8(%r8,%rbp,8)
movsd  0x8(%rax,%rbp,8),%xmm0
movsd  %xmm0,0x8(%rdx,%rbp,8)
add    $0x2,%rbp
cmp    %rbp,%rsi
jne    1790 <VerletUpdate+0x640>
jmp    1978 <VerletUpdate+0x828>
mov    %rax,%r12
mov    %rdi,%r14
lea    (%r8,%rsi,8),%rax
mov    0x40(%rsp),%rcx
cmp    %rcx,%rax
seta   %r15b
lea    (%rcx,%rsi,8),%rdx
cmp    %r8,%rdx
seta   %r13b
mov    %r8,%rbp
mov    0x38(%rsp),%rbx
lea    (%rbx,%rsi,8),%rdi
cmp    %rcx,%rdi
seta   %cl
cmp    %rbx,%rdx
seta   %sil
cmp    %r8,%rdi
seta   %dl
cmp    %rbx,%rax
seta   %dil
xor    %ebp,%ebp
test   %r13b,%r15b
jne    19c2 <VerletUpdate+0x872>
and    %sil,%cl
jne    19c2 <VerletUpdate+0x872>
and    %dil,%dl
mov    0x48(%rsp),%r8
mov    %r14,%rdi
mov    %r12,%rax
mov    0x8(%rsp),%rsi
jne    1739 <VerletUpdate+0x5e9>
mov    %esi,%ebp
and    $0xfffffffc,%ebp
lea    -0x4(%rbp),%rcx
mov    %rcx,%r9
shr    $0x2,%r9
add    $0x1,%r9
test   %rcx,%rcx
je     19e8 <VerletUpdate+0x898>
mov    %r9,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %ecx,%ecx
mov    0x40(%rsp),%rdx
mov    0x38(%rsp),%rbx
nopl   0x0(%rax,%rax,1)
movups (%rdi,%rcx,8),%xmm0
movups 0x10(%rdi,%rcx,8),%xmm1
movups %xmm0,(%rdx,%rcx,8)
movups %xmm1,0x10(%rdx,%rcx,8)
movups (%r11,%rcx,8),%xmm0
movups 0x10(%r11,%rcx,8),%xmm1
movups %xmm0,(%r8,%rcx,8)
movups %xmm1,0x10(%r8,%rcx,8)
movups (%rax,%rcx,8),%xmm0
movups 0x10(%rax,%rcx,8),%xmm1
movups %xmm0,(%rbx,%rcx,8)
movups %xmm1,0x10(%rbx,%rcx,8)
movups 0x20(%rdi,%rcx,8),%xmm0
movups 0x30(%rdi,%rcx,8),%xmm1
movups %xmm0,0x20(%rdx,%rcx,8)
movups %xmm1,0x30(%rdx,%rcx,8)
movups 0x20(%r11,%rcx,8),%xmm0
movups 0x30(%r11,%rcx,8),%xmm1
movups %xmm0,0x20(%r8,%rcx,8)
movups %xmm1,0x30(%r8,%rcx,8)
movupd 0x20(%rax,%rcx,8),%xmm0
movupd 0x30(%rax,%rcx,8),%xmm1
movupd %xmm0,0x20(%rbx,%rcx,8)
movupd %xmm1,0x30(%rbx,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rsi
jne    1890 <VerletUpdate+0x740>
test   $0x1,%r9b
je     196a <VerletUpdate+0x81a>
movups (%rdi,%rcx,8),%xmm0
movups 0x10(%rdi,%rcx,8),%xmm1
mov    0x40(%rsp),%rdx
movups %xmm0,(%rdx,%rcx,8)
movups %xmm1,0x10(%rdx,%rcx,8)
movups (%r11,%rcx,8),%xmm0
movups 0x10(%r11,%rcx,8),%xmm1
movups %xmm0,(%r8,%rcx,8)
movups %xmm1,0x10(%r8,%rcx,8)
movupd (%rax,%rcx,8),%xmm0
movupd 0x10(%rax,%rcx,8),%xmm1
mov    0x38(%rsp),%rdx
movupd %xmm0,(%rdx,%rcx,8)
movupd %xmm1,0x10(%rdx,%rcx,8)
mov    0x8(%rsp),%rsi
cmp    %rsi,%rbp
jne    1739 <VerletUpdate+0x5e9>
mov    %r11,%rbp
mov    %rax,%rbx
call   1040 <free@plt>
mov    %rbp,%rdi
call   1040 <free@plt>
mov    %rbx,%rdi
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <free@plt>
mov    0x28(%rsp),%rdi
mov    0x18(%rsp),%r11
mov    0x20(%rsp),%rax
mov    0x8(%rsp),%rbx
mov    0x30(%rsp),%r12
jmp    1402 <VerletUpdate+0x2b2>
mov    0x48(%rsp),%r8
mov    %r14,%rdi
mov    %r12,%rax
mov    0x8(%rsp),%rsi
jmp    1739 <VerletUpdate+0x5e9>
xor    %edx,%edx
test   $0x1,%r9b
jne    15f8 <VerletUpdate+0x4a8>
jmp    162a <VerletUpdate+0x4da>
xor    %ecx,%ecx
test   $0x1,%r9b
jne    1922 <VerletUpdate+0x7d2>
jmp    196a <VerletUpdate+0x81a>

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

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