void hessian_edmiston(const int Norbs, double *eri, double *vector_in,
                      double *vector_out) {
  int rowindex, p, q, r, s;
  const int num_vars = (Norbs * (Norbs - 1)) / 2;
#pragma omp parallel for schedule(static)
  for (rowindex = 0; rowindex < num_vars; rowindex++) {
    s = num_vars - 1 - rowindex;
    r = 2;
    while ((r * (r - 1)) / 2 <= s) {
      r++;
    }
    p = Norbs - r;
    q = rowindex - p * Norbs + ((p + 1) * (p + 2)) / 2;
    double value = 0.0;
    for (s = p + 1; s < Norbs; s++) {
      const int pqps = p + Norbs * (q + Norbs * (p + Norbs * s));
      const int ppqs = p + Norbs * (p + Norbs * (q + Norbs * s));
      const int qqqs = q + Norbs * (q + Norbs * (q + Norbs * s));
      const int sssq = s + Norbs * (s + Norbs * (s + Norbs * q));
      const double prefactor =
          2 * (4 * eri[pqps] + 2 * eri[ppqs] - eri[qqqs] - eri[sssq]);
      const int colindex = rowindex + s - q;
      value += prefactor * vector_in[colindex];
    }
    for (r = 0; r < q; r++) {
      const int qpqr = q + Norbs * (p + Norbs * (q + Norbs * r));
      const int qqpr = q + Norbs * (q + Norbs * (p + Norbs * r));
      const int pppr = p + Norbs * (p + Norbs * (p + Norbs * r));
      const int rrrp = r + Norbs * (r + Norbs * (r + Norbs * p));
      const double prefactor =
          2 * (4 * eri[qpqr] + 2 * eri[qqpr] - eri[pppr] - eri[rrrp]);
      const int colindex = q + r * Norbs - ((r + 1) * (r + 2)) / 2;
      value += prefactor * vector_in[colindex];
    }
    for (s = q + 1; s < Norbs; s++) {
      const int qpqs = q + Norbs * (p + Norbs * (q + Norbs * s));
      const int qqps = q + Norbs * (q + Norbs * (p + Norbs * s));
      const int ppps = p + Norbs * (p + Norbs * (p + Norbs * s));
      const int sssp = s + Norbs * (s + Norbs * (s + Norbs * p));
      const double prefactor =
          2 * (4 * eri[qpqs] + 2 * eri[qqps] - eri[ppps] - eri[sssp]);
      const int colindex = s + q * Norbs - ((q + 1) * (q + 2)) / 2;
      value -= prefactor * vector_in[colindex];
    }
    for (r = 0; r < p; r++) {
      const int pqpr = p + Norbs * (q + Norbs * (p + Norbs * r));
      const int ppqr = p + Norbs * (p + Norbs * (q + Norbs * r));
      const int qqqr = q + Norbs * (q + Norbs * (q + Norbs * r));
      const int rrrq = r + Norbs * (r + Norbs * (r + Norbs * q));
      const double prefactor =
          2 * (4 * eri[pqpr] + 2 * eri[ppqr] - eri[qqqr] - eri[rrrq]);
      const int colindex = p + r * Norbs - ((r + 1) * (r + 2)) / 2;
      value -= prefactor * vector_in[colindex];
    }
    vector_out[rowindex] = -value;
  }
}


x86 assembly (clang, -O0):
<hessian_edmiston>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x4(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    1827 <hessian_edmiston+0x727>
mov    -0x38(%rbp),%eax
sub    $0x1,%eax
sub    -0x24(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x2,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    -0x34(%rbp),%eax
jg     117f <hessian_edmiston+0x7f>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1154 <hessian_edmiston+0x54>
mov    -0x4(%rbp),%eax
sub    -0x30(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x28(%rbp),%ecx
imul   -0x4(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,-0xc0(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    -0x28(%rbp),%ecx
add    $0x2,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0xc0(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x2c(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1316 <hessian_edmiston+0x216>
mov    -0x28(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x2c(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x28(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x34(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x44(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x28(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x2c(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x34(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x2c(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x2c(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x34(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x34(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x34(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x2c(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  0xd68(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movsd  0xd63(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
movsd  0xd31(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x24(%rbp),%eax
add    -0x34(%rbp),%eax
sub    -0x2c(%rbp),%eax
mov    %eax,-0x5c(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11cf <hessian_edmiston+0xcf>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    148e <hessian_edmiston+0x38e>
mov    -0x2c(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x28(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x2c(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x30(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x60(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x2c(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x28(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x30(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x64(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x28(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x28(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x30(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x68(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x30(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x30(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x28(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x6c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  0xc1a(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movsd  0xc15(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movslq -0x68(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x6c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
movsd  0xbe3(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x30(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0xc4(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    -0x30(%rbp),%ecx
add    $0x2,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0xc4(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x7c(%rbp)
movsd  -0x78(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x7c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    131d <hessian_edmiston+0x21d>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    163d <hessian_edmiston+0x53d>
mov    -0x2c(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x28(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x2c(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x34(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x80(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x2c(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x28(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x34(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x84(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x28(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x28(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x34(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x88(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x34(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x34(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x28(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x8c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x80(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  0xa94(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movsd  0xa8f(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movslq -0x88(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x8c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
movsd  0xa57(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x34(%rbp),%eax
mov    -0x2c(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0xc8(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    -0x2c(%rbp),%ecx
add    $0x2,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0xc8(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x9c(%rbp)
movsd  -0x98(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x9c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1497 <hessian_edmiston+0x397>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    17f0 <hessian_edmiston+0x6f0>
mov    -0x28(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x2c(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x28(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x30(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0xa0(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x28(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x2c(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x30(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0xa4(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x2c(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x2c(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x30(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0xa8(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x4(%rbp),%ecx
mov    -0x30(%rbp),%edx
mov    -0x4(%rbp),%esi
mov    -0x30(%rbp),%edi
mov    -0x4(%rbp),%r8d
imul   -0x2c(%rbp),%r8d
add    %r8d,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0xac(%rbp)
mov    -0x10(%rbp),%rax
movslq -0xa0(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  0x8e1(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movsd  0x8dc(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movslq -0xa8(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0xac(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
movsd  0x8a4(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0x28(%rbp),%eax
mov    -0x30(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0xcc(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    -0x30(%rbp),%ecx
add    $0x2,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0xcc(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0xbc(%rbp)
movsd  -0xb8(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0xbc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1644 <hessian_edmiston+0x544>
movsd  -0x40(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1135 <hessian_edmiston+0x35>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<hessian_edmiston>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %rcx,0x8(%rsp)
lea    -0x1(%rdi),%eax
mov    %eax,-0x68(%rsp)
imul   %edi,%eax
cmp    $0x2,%eax
jl     16f3 <hessian_edmiston+0x5f3>
mov    %rsi,%r10
mov    %eax,%ecx
mov    %eax,%ebp
shr    $0x1f,%eax
add    %ecx,%eax
sar    %eax
movslq %edi,%r12
mov    %eax,-0x1c(%rsp)
mov    %eax,%eax
mov    %rax,(%rsp)
add    $0xffffffff,%ebp
mov    %edi,%ebx
not    %ebx
lea    -0x2(%rdi),%r8d
mov    %edi,%esi
imul   %edi,%esi
mov    %edi,%eax
imul   %esi,%eax
mov    %esi,-0xc(%rsp)
mov    %esi,%ecx
imul   -0x68(%rsp),%ecx
mov    %edi,%esi
mov    %ebp,-0x18(%rsp)
imul   %ebp,%esi
add    $0xfffffffe,%esi
mov    %esi,-0x24(%rsp)
mov    %edi,%esi
imul   %ebx,%esi
add    $0xffffffff,%esi
mov    %rsi,0x38(%rsp)
mov    %eax,-0x7c(%rsp)
mov    %eax,%ebp
mov    %r8d,-0x20(%rsp)
imul   %r8d,%ebp
mov    %ebp,-0x28(%rsp)
lea    0x1(%rdi),%eax
imul   %edi,%eax
mov    %eax,-0x10(%rsp)
mov    $0x1,%eax
mov    %rax,-0x38(%rsp)
xor    %esi,%esi
movsd  0xe51(%rip),%xmm0
movsd  0xe51(%rip),%xmm1
movapd 0xe51(%rip),%xmm2
mov    %edi,%eax
mov    %rax,-0x40(%rsp)
mov    %rdx,-0x50(%rsp)
mov    %r10,-0x48(%rsp)
mov    %rdi,-0x78(%rsp)
mov    %r12,-0x58(%rsp)
mov    %rbx,0x40(%rsp)
jmp    1232 <hessian_edmiston+0x132>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xorpd  %xmm2,%xmm3
mov    0x8(%rsp),%rax
mov    -0x60(%rsp),%rsi
movlpd %xmm3,(%rax,%rsi,8)
add    $0x1,%rsi
mov    -0x38(%rsp),%rax
add    $0x1,%eax
mov    %rax,-0x38(%rsp)
mov    -0x40(%rsp),%rax
add    $0x1,%eax
mov    %rax,-0x40(%rsp)
mov    -0x14(%rsp),%ecx
add    $0x1,%ecx
cmp    (%rsp),%rsi
je     16f3 <hessian_edmiston+0x5f3>
mov    %rsi,-0x60(%rsp)
mov    %esi,%eax
not    %eax
add    -0x1c(%rsp),%eax
mov    %eax,-0x70(%rsp)
mov    $0x2,%r11d
mov    $0xffffffff,%edi
mov    -0x28(%rsp),%ebp
mov    -0x24(%rsp),%esi
mov    -0x68(%rsp),%r12d
mov    %ecx,-0x14(%rsp)
mov    -0x20(%rsp),%r13d
mov    -0x18(%rsp),%r9d
nopl   0x0(%rax)
mov    %r11d,%ebx
mov    %r9d,%edx
mov    %ecx,%r8d
mov    %r12d,%r10d
mov    %esi,%r15d
mov    %ebp,%r14d
lea    -0x1(%rbx),%ecx
imul   %r11d,%ecx
mov    %ecx,%eax
shr    $0x1f,%eax
add    %ecx,%eax
sar    %eax
lea    0x1(%rbx),%r11d
mov    %rdx,0x50(%rsp)
mov    0x40(%rsp),%rcx
lea    (%rdx,%rcx,1),%r9d
mov    %r13,%rdx
add    $0xffffffff,%r13d
mov    %r8,0x48(%rsp)
mov    %r8d,%ecx
sub    -0xc(%rsp),%ecx
mov    %rdi,%r8
add    $0xffffffff,%edi
sub    -0x7c(%rsp),%ebp
lea    -0x1(%r10),%r12d
mov    0x38(%rsp),%rsi
add    %r15d,%esi
cmp    -0x70(%rsp),%eax
jle    1270 <hessian_edmiston+0x170>
mov    %r10,0x10(%rsp)
mov    %r14,0x20(%rsp)
mov    %r15,0x28(%rsp)
mov    %rdx,%rbp
mov    %r12d,-0x64(%rsp)
mov    -0x78(%rsp),%rdi
mov    %edi,%ecx
sub    %ebx,%ecx
mov    %ecx,%esi
imul   %edi,%esi
mov    -0x60(%rsp),%rax
mov    %rsi,0x18(%rsp)
sub    %esi,%eax
lea    0x1(%rcx),%ebx
mov    %rcx,-0x30(%rsp)
add    $0x2,%ecx
imul   %ebx,%ecx
mov    %ecx,%edx
shr    $0x1f,%edx
add    %ecx,%edx
sar    %edx
mov    %rdx,%r14
add    %edx,%eax
xorpd  %xmm3,%xmm3
mov    %rax,-0x70(%rsp)
movslq %eax,%rcx
mov    %rcx,-0x8(%rsp)
cmp    %edi,%ebx
mov    -0x48(%rsp),%rdi
mov    0x50(%rsp),%r10
mov    0x48(%rsp),%r11
jge    1435 <hessian_edmiston+0x335>
mov    %r8,%rdx
mov    -0x70(%rsp),%rcx
mov    -0x78(%rsp),%r9
imul   %r9d,%ecx
mov    %rcx,0x30(%rsp)
movslq %ebx,%r13
mov    %r14,%rcx
sub    %ecx,%r10d
mov    -0x40(%rsp),%rax
lea    (%r14,%rax,1),%ebx
imul   %r9d,%ebx
add    %ebp,%ebx
imul   %r9d,%ebx
add    %ebp,%ebx
add    %ecx,%r11d
imul   %r9d,%r11d
add    %ebp,%r11d
movslq 0x10(%rsp),%rax
imul   -0x58(%rsp),%rax
xorpd  %xmm3,%xmm3
cs nopw 0x0(%rax,%rax,1)
nop
mov    -0x78(%rsp),%r12
lea    (%r12,%rdx,1),%ecx
mov    0x30(%rsp),%rsi
add    %ecx,%esi
imul   %r12d,%esi
add    %ecx,%esi
mov    -0x8(%rsp),%r9
lea    (%rax,%r9,1),%rcx
mov    -0x58(%rsp),%rbp
imul   %rbp,%rcx
add    %r9,%rcx
imul   %rbp,%rcx
add    %r9,%rcx
movslq %r11d,%r11
movslq %ebx,%rbx
movsd  (%rdi,%rbx,8),%xmm4
addsd  %xmm4,%xmm4
movsd  (%rdi,%r11,8),%xmm5
mulsd  %xmm0,%xmm5
addsd  %xmm4,%xmm5
subsd  (%rdi,%rcx,8),%xmm5
imul   %r12d,%esi
movslq %esi,%rcx
mov    -0x50(%rsp),%rsi
add    %r13,%rcx
subsd  (%rdi,%rcx,8),%xmm5
addsd  %xmm5,%xmm5
movslq %r10d,%r10
mulsd  (%rsi,%r10,8),%xmm5
addsd  %xmm5,%xmm3
add    $0x1,%r13
add    $0x1,%r10d
mov    -0x7c(%rsp),%ecx
add    %ecx,%ebx
add    %ecx,%r11d
add    %rbp,%rax
inc    %edx
jne    13a0 <hessian_edmiston+0x2a0>
mov    -0x70(%rsp),%r11
test   %r11d,%r11d
mov    0x28(%rsp),%rdx
mov    0x20(%rsp),%rbp
jle    1510 <hessian_edmiston+0x410>
mov    %r11d,%r13d
mov    -0x60(%rsp),%rax
mov    %r14,%r8
lea    (%r14,%rax,1),%r9d
mov    -0x78(%rsp),%rdi
imul   %edi,%r9d
xor    %r15d,%r15d
xor    %ebx,%ebx
mov    -0x48(%rsp),%r10
mov    -0x58(%rsp),%r12
mov    -0x7c(%rsp),%ecx
nopl   0x0(%rax,%rax,1)
movslq %ebp,%rbp
lea    (%r15,%r11,1),%eax
imul   %edi,%eax
add    -0x30(%rsp),%eax
imul   %edi,%eax
add    %r11d,%eax
mov    %ecx,%r14d
lea    (%r9,%r11,1),%ecx
movslq %ecx,%rcx
movsd  (%r10,%rcx,8),%xmm4
lea    (%rbx,%rbp,1),%rcx
cltq
addsd  %xmm4,%xmm4
movsd  (%r10,%rax,8),%xmm5
mulsd  %xmm0,%xmm5
addsd  %xmm4,%xmm5
movslq %edx,%rdx
subsd  (%r10,%rdx,8),%xmm5
subsd  (%r10,%rcx,8),%xmm5
lea    0x1(%rbx),%rax
add    $0x2,%ebx
mov    %eax,%ecx
imul   %ebx,%ecx
lea    (%r15,%r13,1),%rsi
shr    %ecx
sub    %rcx,%rsi
mov    %r14d,%ecx
addsd  %xmm5,%xmm5
mov    -0x50(%rsp),%rbx
mulsd  (%rbx,%rsi,8),%xmm5
addsd  %xmm5,%xmm3
add    %r14d,%edx
add    %r14d,%r9d
add    %r12,%r15
add    -0x10(%rsp),%ebp
mov    %rax,%rbx
cmp    %r13,%rax
jb     1480 <hessian_edmiston+0x380>
jmp    1522 <hessian_edmiston+0x422>
mov    -0x48(%rsp),%r10
mov    -0x78(%rsp),%rdi
mov    -0x58(%rsp),%r12
mov    %r14,%r8
mov    -0x50(%rsp),%r9
lea    0x1(%r11),%eax
cmp    %edi,%eax
mov    %r11,%r15
jge    1610 <hessian_edmiston+0x510>
mov    %r15d,%r13d
imul   %edi,%r13d
lea    0x2(%r15),%esi
imul   %esi,%eax
mov    %eax,%esi
shr    $0x1f,%esi
add    %eax,%esi
sar    %esi
sub    %esi,%r13d
mov    -0x38(%rsp),%rax
add    %r8d,%eax
mov    0x18(%rsp),%r11
sub    %r11d,%eax
movslq %eax,%rbx
mov    -0x30(%rsp),%rdx
nopl   0x0(%rax,%rax,1)
mov    %ebx,%eax
imul   %edi,%eax
lea    (%rax,%r15,1),%ebp
imul   %edi,%ebp
add    %edx,%ebp
imul   %edi,%ebp
add    %edx,%eax
imul   %edi,%eax
lea    (%rax,%r15,1),%esi
add    %r15d,%ebp
imul   %edi,%esi
add    %edx,%eax
imul   %edi,%eax
add    %r15d,%esi
lea    (%r11,%rbx,1),%ecx
imul   %edi,%ecx
add    %ebx,%ecx
add    %edx,%eax
imul   %edi,%ecx
movslq %ebp,%rbp
movslq %esi,%rsi
movsd  (%r10,%rsi,8),%xmm4
addsd  %xmm4,%xmm4
movsd  (%r10,%rbp,8),%xmm5
mulsd  %xmm0,%xmm5
addsd  %xmm4,%xmm5
cltq
subsd  (%r10,%rax,8),%xmm5
movslq %ecx,%rax
add    %rbx,%rax
subsd  (%r10,%rax,8),%xmm5
lea    (%rbx,%r13,1),%eax
cltq
mulsd  %xmm1,%xmm5
mulsd  (%r9,%rax,8),%xmm5
addsd  %xmm5,%xmm3
add    $0x1,%rbx
cmp    %r12,%rbx
jl     1570 <hessian_edmiston+0x470>
test   %edx,%edx
jg     161d <hessian_edmiston+0x51d>
jmp    11f0 <hessian_edmiston+0xf0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    -0x30(%rsp),%rdx
test   %edx,%edx
jle    11f0 <hessian_edmiston+0xf0>
mov    %r15,%r11
imul   %edi,%r11d
add    -0x60(%rsp),%r8d
imul   %edi,%r8d
mov    -0x64(%rsp),%eax
add    %eax,%r8d
mov    %eax,%r13d
mov    $0x4,%r9d
mov    $0x2,%r14d
xor    %ebx,%ebx
xor    %eax,%eax
mov    -0x8(%rsp),%r15
nopl   0x0(%rax)
mov    %r14d,%esi
shr    %esi
mov    -0x64(%rsp),%ecx
sub    %esi,%ecx
mov    -0x70(%rsp),%rsi
lea    (%rbx,%rsi,1),%ebp
imul   %edi,%ebp
add    %edx,%ebp
imul   %edi,%ebp
add    %edx,%ebp
movslq %ebp,%rsi
movsd  (%r10,%rsi,8),%xmm4
lea    (%rbx,%r15,1),%rsi
imul   %r12,%rsi
add    %r15,%rsi
imul   %r12,%rsi
add    %r15,%rsi
movslq %r8d,%r8
movsd  (%r10,%r8,8),%xmm5
mulsd  %xmm0,%xmm5
addsd  %xmm4,%xmm4
addsd  %xmm4,%xmm5
subsd  (%r10,%rsi,8),%xmm5
lea    (%r11,%rax,1),%esi
imul   %edi,%esi
add    %eax,%esi
imul   %edi,%esi
movslq %esi,%rsi
add    %rax,%rsi
subsd  (%r10,%rsi,8),%xmm5
add    %ebx,%ecx
movslq %ecx,%rcx
mulsd  %xmm1,%xmm5
mov    -0x50(%rsp),%rsi
mulsd  (%rsi,%rcx,8),%xmm5
addsd  %xmm5,%xmm3
add    %r12,%rbx
add    $0x1,%rax
add    %r9d,%r14d
add    $0x2,%r9d
add    -0x7c(%rsp),%r8d
cmp    %rax,%r13
jne    1650 <hessian_edmiston+0x550>
jmp    11f0 <hessian_edmiston+0xf0>
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @hessian_edmiston(i32 %0, double* %1, double* %2, double* %3) {
%5 = alloca i32
%6 = alloca double*
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca double
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca double
%21 = alloca i32
%22 = alloca i32
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca double
%27 = alloca i32
%28 = alloca i32
%29 = alloca i32
%30 = alloca i32
%31 = alloca i32
%32 = alloca double
%33 = alloca i32
%34 = alloca i32
%35 = alloca i32
%36 = alloca i32
%37 = alloca i32
%38 = alloca double
%39 = alloca i32
store i32 %0, i32* %5
store double* %1, double** %6
store double* %2, double** %7
store double* %3, double** %8
%40 = load i32, i32* %5
%41 = load i32, i32* %5
%42 = sub nsw i32 %41, 1
%43 = mul nsw i32 %40, %42
%44 = sdiv i32 %43, 2
store i32 %44, i32* %14
store i32 0, i32* %9
br label %45
45: ; preds = %510, %4
%46 = load i32, i32* %9
%47 = load i32, i32* %14
%48 = icmp slt i32 %46, %47
br i1 %48, label %49, label %513
49: ; preds = %45
%50 = load i32, i32* %14
%51 = sub nsw i32 %50, 1
%52 = load i32, i32* %9
%53 = sub nsw i32 %51, %52
store i32 %53, i32* %13
store i32 2, i32* %12
br label %54
54: ; preds = %62, %49
%55 = load i32, i32* %12
%56 = load i32, i32* %12
%57 = sub nsw i32 %56, 1
%58 = mul nsw i32 %55, %57
%59 = sdiv i32 %58, 2
%60 = load i32, i32* %13
%61 = icmp sle i32 %59, %60
br i1 %61, label %62, label %65
62: ; preds = %54
%63 = load i32, i32* %12
%64 = add nsw i32 %63, 1
store i32 %64, i32* %12
br label %54
65: ; preds = %54
%66 = load i32, i32* %5
%67 = load i32, i32* %12
%68 = sub nsw i32 %66, %67
store i32 %68, i32* %10
%69 = load i32, i32* %9
%70 = load i32, i32* %10
%71 = load i32, i32* %5
%72 = mul nsw i32 %70, %71
%73 = sub nsw i32 %69, %72
%74 = load i32, i32* %10
%75 = add nsw i32 %74, 1
%76 = load i32, i32* %10
%77 = add nsw i32 %76, 2
%78 = mul nsw i32 %75, %77
%79 = sdiv i32 %78, 2
%80 = add nsw i32 %73, %79
store i32 %80, i32* %11
store double 0.000000e+00, double* %15
%81 = load i32, i32* %10
%82 = add nsw i32 %81, 1
store i32 %82, i32* %13
br label %83
83: ; preds = %178, %65
%84 = load i32, i32* %13
%85 = load i32, i32* %5
%86 = icmp slt i32 %84, %85
br i1 %86, label %87, label %181
87: ; preds = %83
%88 = load i32, i32* %10
%89 = load i32, i32* %5
%90 = load i32, i32* %11
%91 = load i32, i32* %5
%92 = load i32, i32* %10
%93 = load i32, i32* %5
%94 = load i32, i32* %13
%95 = mul nsw i32 %93, %94
%96 = add nsw i32 %92, %95
%97 = mul nsw i32 %91, %96
%98 = add nsw i32 %90, %97
%99 = mul nsw i32 %89, %98
%100 = add nsw i32 %88, %99
store i32 %100, i32* %16
%101 = load i32, i32* %10
%102 = load i32, i32* %5
%103 = load i32, i32* %10
%104 = load i32, i32* %5
%105 = load i32, i32* %11
%106 = load i32, i32* %5
%107 = load i32, i32* %13
%108 = mul nsw i32 %106, %107
%109 = add nsw i32 %105, %108
%110 = mul nsw i32 %104, %109
%111 = add nsw i32 %103, %110
%112 = mul nsw i32 %102, %111
%113 = add nsw i32 %101, %112
store i32 %113, i32* %17
%114 = load i32, i32* %11
%115 = load i32, i32* %5
%116 = load i32, i32* %11
%117 = load i32, i32* %5
%118 = load i32, i32* %11
%119 = load i32, i32* %5
%120 = load i32, i32* %13
%121 = mul nsw i32 %119, %120
%122 = add nsw i32 %118, %121
%123 = mul nsw i32 %117, %122
%124 = add nsw i32 %116, %123
%125 = mul nsw i32 %115, %124
%126 = add nsw i32 %114, %125
store i32 %126, i32* %18
%127 = load i32, i32* %13
%128 = load i32, i32* %5
%129 = load i32, i32* %13
%130 = load i32, i32* %5
%131 = load i32, i32* %13
%132 = load i32, i32* %5
%133 = load i32, i32* %11
%134 = mul nsw i32 %132, %133
%135 = add nsw i32 %131, %134
%136 = mul nsw i32 %130, %135
%137 = add nsw i32 %129, %136
%138 = mul nsw i32 %128, %137
%139 = add nsw i32 %127, %138
store i32 %139, i32* %19
%140 = load double*, double** %6
%141 = load i32, i32* %16
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds double, double* %140, i64 %142
%144 = load double, double* %143
%145 = load double*, double** %6
%146 = load i32, i32* %17
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds double, double* %145, i64 %147
%149 = load double, double* %148
%150 = fmul double 2.000000e+00, %149
%151 = call double @llvm.fmuladd.f64(double 4.000000e+00, double %144, double %150)
%152 = load double*, double** %6
%153 = load i32, i32* %18
%154 = sext i32 %153 to i64
%155 = getelementptr inbounds double, double* %152, i64 %154
%156 = load double, double* %155
%157 = fsub double %151, %156
%158 = load double*, double** %6
%159 = load i32, i32* %19
%160 = sext i32 %159 to i64
%161 = getelementptr inbounds double, double* %158, i64 %160
%162 = load double, double* %161
%163 = fsub double %157, %162
%164 = fmul double 2.000000e+00, %163
store double %164, double* %20
%165 = load i32, i32* %9
%166 = load i32, i32* %13
%167 = add nsw i32 %165, %166
%168 = load i32, i32* %11
%169 = sub nsw i32 %167, %168
store i32 %169, i32* %21
%170 = load double, double* %20
%171 = load double*, double** %7
%172 = load i32, i32* %21
%173 = sext i32 %172 to i64
%174 = getelementptr inbounds double, double* %171, i64 %173
%175 = load double, double* %174
%176 = load double, double* %15
%177 = call double @llvm.fmuladd.f64(double %170, double %175, double %176)
store double %177, double* %15
br label %178
178: ; preds = %87
%179 = load i32, i32* %13
%180 = add nsw i32 %179, 1
store i32 %180, i32* %13
br label %83
181: ; preds = %83
store i32 0, i32* %12
br label %182
182: ; preds = %284, %181
%183 = load i32, i32* %12
%184 = load i32, i32* %11
%185 = icmp slt i32 %183, %184
br i1 %185, label %186, label %287
186: ; preds = %182
%187 = load i32, i32* %11
%188 = load i32, i32* %5
%189 = load i32, i32* %10
%190 = load i32, i32* %5
%191 = load i32, i32* %11
%192 = load i32, i32* %5
%193 = load i32, i32* %12
%194 = mul nsw i32 %192, %193
%195 = add nsw i32 %191, %194
%196 = mul nsw i32 %190, %195
%197 = add nsw i32 %189, %196
%198 = mul nsw i32 %188, %197
%199 = add nsw i32 %187, %198
store i32 %199, i32* %22
%200 = load i32, i32* %11
%201 = load i32, i32* %5
%202 = load i32, i32* %11
%203 = load i32, i32* %5
%204 = load i32, i32* %10
%205 = load i32, i32* %5
%206 = load i32, i32* %12
%207 = mul nsw i32 %205, %206
%208 = add nsw i32 %204, %207
%209 = mul nsw i32 %203, %208
%210 = add nsw i32 %202, %209
%211 = mul nsw i32 %201, %210
%212 = add nsw i32 %200, %211
store i32 %212, i32* %23
%213 = load i32, i32* %10
%214 = load i32, i32* %5
%215 = load i32, i32* %10
%216 = load i32, i32* %5
%217 = load i32, i32* %10
%218 = load i32, i32* %5
%219 = load i32, i32* %12
%220 = mul nsw i32 %218, %219
%221 = add nsw i32 %217, %220
%222 = mul nsw i32 %216, %221
%223 = add nsw i32 %215, %222
%224 = mul nsw i32 %214, %223
%225 = add nsw i32 %213, %224
store i32 %225, i32* %24
%226 = load i32, i32* %12
%227 = load i32, i32* %5
%228 = load i32, i32* %12
%229 = load i32, i32* %5
%230 = load i32, i32* %12
%231 = load i32, i32* %5
%232 = load i32, i32* %10
%233 = mul nsw i32 %231, %232
%234 = add nsw i32 %230, %233
%235 = mul nsw i32 %229, %234
%236 = add nsw i32 %228, %235
%237 = mul nsw i32 %227, %236
%238 = add nsw i32 %226, %237
store i32 %238, i32* %25
%239 = load double*, double** %6
%240 = load i32, i32* %22
%241 = sext i32 %240 to i64
%242 = getelementptr inbounds double, double* %239, i64 %241
%243 = load double, double* %242
%244 = load double*, double** %6
%245 = load i32, i32* %23
%246 = sext i32 %245 to i64
%247 = getelementptr inbounds double, double* %244, i64 %246
%248 = load double, double* %247
%249 = fmul double 2.000000e+00, %248
%250 = call double @llvm.fmuladd.f64(double 4.000000e+00, double %243, double %249)
%251 = load double*, double** %6
%252 = load i32, i32* %24
%253 = sext i32 %252 to i64
%254 = getelementptr inbounds double, double* %251, i64 %253
%255 = load double, double* %254
%256 = fsub double %250, %255
%257 = load double*, double** %6
%258 = load i32, i32* %25
%259 = sext i32 %258 to i64
%260 = getelementptr inbounds double, double* %257, i64 %259
%261 = load double, double* %260
%262 = fsub double %256, %261
%263 = fmul double 2.000000e+00, %262
store double %263, double* %26
%264 = load i32, i32* %11
%265 = load i32, i32* %12
%266 = load i32, i32* %5
%267 = mul nsw i32 %265, %266
%268 = add nsw i32 %264, %267
%269 = load i32, i32* %12
%270 = add nsw i32 %269, 1
%271 = load i32, i32* %12
%272 = add nsw i32 %271, 2
%273 = mul nsw i32 %270, %272
%274 = sdiv i32 %273, 2
%275 = sub nsw i32 %268, %274
store i32 %275, i32* %27
%276 = load double, double* %26
%277 = load double*, double** %7
%278 = load i32, i32* %27
%279 = sext i32 %278 to i64
%280 = getelementptr inbounds double, double* %277, i64 %279
%281 = load double, double* %280
%282 = load double, double* %15
%283 = call double @llvm.fmuladd.f64(double %276, double %281, double %282)
store double %283, double* %15
br label %284
284: ; preds = %186
%285 = load i32, i32* %12
%286 = add nsw i32 %285, 1
store i32 %286, i32* %12
br label %182
287: ; preds = %182
%288 = load i32, i32* %11
%289 = add nsw i32 %288, 1
store i32 %289, i32* %13
br label %290
290: ; preds = %393, %287
%291 = load i32, i32* %13
%292 = load i32, i32* %5
%293 = icmp slt i32 %291, %292
br i1 %293, label %294, label %396
294: ; preds = %290
%295 = load i32, i32* %11
%296 = load i32, i32* %5
%297 = load i32, i32* %10
%298 = load i32, i32* %5
%299 = load i32, i32* %11
%300 = load i32, i32* %5
%301 = load i32, i32* %13
%302 = mul nsw i32 %300, %301
%303 = add nsw i32 %299, %302
%304 = mul nsw i32 %298, %303
%305 = add nsw i32 %297, %304
%306 = mul nsw i32 %296, %305
%307 = add nsw i32 %295, %306
store i32 %307, i32* %28
%308 = load i32, i32* %11
%309 = load i32, i32* %5
%310 = load i32, i32* %11
%311 = load i32, i32* %5
%312 = load i32, i32* %10
%313 = load i32, i32* %5
%314 = load i32, i32* %13
%315 = mul nsw i32 %313, %314
%316 = add nsw i32 %312, %315
%317 = mul nsw i32 %311, %316
%318 = add nsw i32 %310, %317
%319 = mul nsw i32 %309, %318
%320 = add nsw i32 %308, %319
store i32 %320, i32* %29
%321 = load i32, i32* %10
%322 = load i32, i32* %5
%323 = load i32, i32* %10
%324 = load i32, i32* %5
%325 = load i32, i32* %10
%326 = load i32, i32* %5
%327 = load i32, i32* %13
%328 = mul nsw i32 %326, %327
%329 = add nsw i32 %325, %328
%330 = mul nsw i32 %324, %329
%331 = add nsw i32 %323, %330
%332 = mul nsw i32 %322, %331
%333 = add nsw i32 %321, %332
store i32 %333, i32* %30
%334 = load i32, i32* %13
%335 = load i32, i32* %5
%336 = load i32, i32* %13
%337 = load i32, i32* %5
%338 = load i32, i32* %13
%339 = load i32, i32* %5
%340 = load i32, i32* %10
%341 = mul nsw i32 %339, %340
%342 = add nsw i32 %338, %341
%343 = mul nsw i32 %337, %342
%344 = add nsw i32 %336, %343
%345 = mul nsw i32 %335, %344
%346 = add nsw i32 %334, %345
store i32 %346, i32* %31
%347 = load double*, double** %6
%348 = load i32, i32* %28
%349 = sext i32 %348 to i64
%350 = getelementptr inbounds double, double* %347, i64 %349
%351 = load double, double* %350
%352 = load double*, double** %6
%353 = load i32, i32* %29
%354 = sext i32 %353 to i64
%355 = getelementptr inbounds double, double* %352, i64 %354
%356 = load double, double* %355
%357 = fmul double 2.000000e+00, %356
%358 = call double @llvm.fmuladd.f64(double 4.000000e+00, double %351, double %357)
%359 = load double*, double** %6
%360 = load i32, i32* %30
%361 = sext i32 %360 to i64
%362 = getelementptr inbounds double, double* %359, i64 %361
%363 = load double, double* %362
%364 = fsub double %358, %363
%365 = load double*, double** %6
%366 = load i32, i32* %31
%367 = sext i32 %366 to i64
%368 = getelementptr inbounds double, double* %365, i64 %367
%369 = load double, double* %368
%370 = fsub double %364, %369
%371 = fmul double 2.000000e+00, %370
store double %371, double* %32
%372 = load i32, i32* %13
%373 = load i32, i32* %11
%374 = load i32, i32* %5
%375 = mul nsw i32 %373, %374
%376 = add nsw i32 %372, %375
%377 = load i32, i32* %11
%378 = add nsw i32 %377, 1
%379 = load i32, i32* %11
%380 = add nsw i32 %379, 2
%381 = mul nsw i32 %378, %380
%382 = sdiv i32 %381, 2
%383 = sub nsw i32 %376, %382
store i32 %383, i32* %33
%384 = load double, double* %32
%385 = load double*, double** %7
%386 = load i32, i32* %33
%387 = sext i32 %386 to i64
%388 = getelementptr inbounds double, double* %385, i64 %387
%389 = load double, double* %388
%390 = load double, double* %15
%391 = fneg double %384
%392 = call double @llvm.fmuladd.f64(double %391, double %389, double %390)
store double %392, double* %15
br label %393
393: ; preds = %294
%394 = load i32, i32* %13
%395 = add nsw i32 %394, 1
store i32 %395, i32* %13
br label %290
396: ; preds = %290
store i32 0, i32* %12
br label %397
397: ; preds = %500, %396
%398 = load i32, i32* %12
%399 = load i32, i32* %10
%400 = icmp slt i32 %398, %399
br i1 %400, label %401, label %503
401: ; preds = %397
%402 = load i32, i32* %10
%403 = load i32, i32* %5
%404 = load i32, i32* %11
%405 = load i32, i32* %5
%406 = load i32, i32* %10
%407 = load i32, i32* %5
%408 = load i32, i32* %12
%409 = mul nsw i32 %407, %408
%410 = add nsw i32 %406, %409
%411 = mul nsw i32 %405, %410
%412 = add nsw i32 %404, %411
%413 = mul nsw i32 %403, %412
%414 = add nsw i32 %402, %413
store i32 %414, i32* %34
%415 = load i32, i32* %10
%416 = load i32, i32* %5
%417 = load i32, i32* %10
%418 = load i32, i32* %5
%419 = load i32, i32* %11
%420 = load i32, i32* %5
%421 = load i32, i32* %12
%422 = mul nsw i32 %420, %421
%423 = add nsw i32 %419, %422
%424 = mul nsw i32 %418, %423
%425 = add nsw i32 %417, %424
%426 = mul nsw i32 %416, %425
%427 = add nsw i32 %415, %426
store i32 %427, i32* %35
%428 = load i32, i32* %11
%429 = load i32, i32* %5
%430 = load i32, i32* %11
%431 = load i32, i32* %5
%432 = load i32, i32* %11
%433 = load i32, i32* %5
%434 = load i32, i32* %12
%435 = mul nsw i32 %433, %434
%436 = add nsw i32 %432, %435
%437 = mul nsw i32 %431, %436
%438 = add nsw i32 %430, %437
%439 = mul nsw i32 %429, %438
%440 = add nsw i32 %428, %439
store i32 %440, i32* %36
%441 = load i32, i32* %12
%442 = load i32, i32* %5
%443 = load i32, i32* %12
%444 = load i32, i32* %5
%445 = load i32, i32* %12
%446 = load i32, i32* %5
%447 = load i32, i32* %11
%448 = mul nsw i32 %446, %447
%449 = add nsw i32 %445, %448
%450 = mul nsw i32 %444, %449
%451 = add nsw i32 %443, %450
%452 = mul nsw i32 %442, %451
%453 = add nsw i32 %441, %452
store i32 %453, i32* %37
%454 = load double*, double** %6
%455 = load i32, i32* %34
%456 = sext i32 %455 to i64
%457 = getelementptr inbounds double, double* %454, i64 %456
%458 = load double, double* %457
%459 = load double*, double** %6
%460 = load i32, i32* %35
%461 = sext i32 %460 to i64
%462 = getelementptr inbounds double, double* %459, i64 %461
%463 = load double, double* %462
%464 = fmul double 2.000000e+00, %463
%465 = call double @llvm.fmuladd.f64(double 4.000000e+00, double %458, double %464)
%466 = load double*, double** %6
%467 = load i32, i32* %36
%468 = sext i32 %467 to i64
%469 = getelementptr inbounds double, double* %466, i64 %468
%470 = load double, double* %469
%471 = fsub double %465, %470
%472 = load double*, double** %6
%473 = load i32, i32* %37
%474 = sext i32 %473 to i64
%475 = getelementptr inbounds double, double* %472, i64 %474
%476 = load double, double* %475
%477 = fsub double %471, %476
%478 = fmul double 2.000000e+00, %477
store double %478, double* %38
%479 = load i32, i32* %10
%480 = load i32, i32* %12
%481 = load i32, i32* %5
%482 = mul nsw i32 %480, %481
%483 = add nsw i32 %479, %482
%484 = load i32, i32* %12
%485 = add nsw i32 %484, 1
%486 = load i32, i32* %12
%487 = add nsw i32 %486, 2
%488 = mul nsw i32 %485, %487
%489 = sdiv i32 %488, 2
%490 = sub nsw i32 %483, %489
store i32 %490, i32* %39
%491 = load double, double* %38
%492 = load double*, double** %7
%493 = load i32, i32* %39
%494 = sext i32 %493 to i64
%495 = getelementptr inbounds double, double* %492, i64 %494
%496 = load double, double* %495
%497 = load double, double* %15
%498 = fneg double %491
%499 = call double @llvm.fmuladd.f64(double %498, double %496, double %497)
store double %499, double* %15
br label %500
500: ; preds = %401
%501 = load i32, i32* %12
%502 = add nsw i32 %501, 1
store i32 %502, i32* %12
br label %397
503: ; preds = %397
%504 = load double, double* %15
%505 = fneg double %504
%506 = load double*, double** %8
%507 = load i32, i32* %9
%508 = sext i32 %507 to i64
%509 = getelementptr inbounds double, double* %506, i64 %508
store double %505, double* %509
br label %510
510: ; preds = %503
%511 = load i32, i32* %9
%512 = add nsw i32 %511, 1
store i32 %512, i32* %9
br label %45
513: ; preds = %45
ret void
}

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