#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 */
void NPdsymm_sum(int n, double *a, double *out, int hermi) {
  size_t i, j, j0, j1;
  double tmp;
  if (hermi == 1 || hermi == 3) {
    for (j0 = 0; j0 < n; j0 += 120)
      for (i = 0, j1 = ((j0 + 120) < (n) ? (j0 + 120) : (n)); i < j1; i++)
        for (j = ((i) > (j0) ? (i) : (j0)); j < j1; j++) {
          tmp = a[i * n + j] + a[j * n + i];
          out[i * n + j] = tmp;
          out[j * n + i] = tmp;
        }
  } else {
    for (j0 = 0; j0 < n; j0 += 120)
      for (i = 0, j1 = ((j0 + 120) < (n) ? (j0 + 120) : (n)); i < j1; i++)
        for (j = ((i) > (j0) ? (i) : (j0)); j < j1; j++) {
          tmp = a[i * n + j] - a[j * n + i];
          out[i * n + j] = tmp;
          out[j * n + i] = -tmp;
        }
  }
}


x86 assembly (clang, -O0):
<NPdsymm_sum>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
cmpl   $0x1,-0x1c(%rbp)
je     1126 <NPdsymm_sum+0x26>
cmpl   $0x3,-0x1c(%rbp)
jne    1279 <NPdsymm_sum+0x179>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x4(%rbp),%rcx
cmp    %rcx,%rax
jae    1274 <NPdsymm_sum+0x174>
movq   $0x0,-0x28(%rbp)
mov    -0x38(%rbp),%rax
add    $0x78,%rax
movslq -0x4(%rbp),%rcx
cmp    %rcx,%rax
jae    116d <NPdsymm_sum+0x6d>
mov    -0x38(%rbp),%rax
add    $0x78,%rax
mov    %rax,-0x50(%rbp)
jmp    1175 <NPdsymm_sum+0x75>
movslq -0x4(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x40(%rbp),%rax
jae    125e <NPdsymm_sum+0x15e>
mov    -0x28(%rbp),%rax
cmp    -0x38(%rbp),%rax
jbe    11a6 <NPdsymm_sum+0xa6>
mov    -0x28(%rbp),%rax
mov    %rax,-0x58(%rbp)
jmp    11ae <NPdsymm_sum+0xae>
mov    -0x38(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x58(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x40(%rbp),%rax
jae    1248 <NPdsymm_sum+0x148>
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
movslq -0x4(%rbp),%rdx
imul   %rdx,%rcx
add    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movslq -0x4(%rbp),%rdx
imul   %rdx,%rcx
add    -0x28(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%rcx
movslq -0x4(%rbp),%rdx
imul   %rdx,%rcx
add    -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
movslq -0x4(%rbp),%rdx
imul   %rdx,%rcx
add    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    11b6 <NPdsymm_sum+0xb6>
jmp    124d <NPdsymm_sum+0x14d>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    117d <NPdsymm_sum+0x7d>
jmp    1263 <NPdsymm_sum+0x163>
mov    -0x38(%rbp),%rax
add    $0x78,%rax
mov    %rax,-0x38(%rbp)
jmp    112e <NPdsymm_sum+0x2e>
jmp    13e3 <NPdsymm_sum+0x2e3>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x4(%rbp),%rcx
cmp    %rcx,%rax
jae    13de <NPdsymm_sum+0x2de>
movq   $0x0,-0x28(%rbp)
mov    -0x38(%rbp),%rax
add    $0x78,%rax
movslq -0x4(%rbp),%rcx
cmp    %rcx,%rax
jae    12c0 <NPdsymm_sum+0x1c0>
mov    -0x38(%rbp),%rax
add    $0x78,%rax
mov    %rax,-0x60(%rbp)
jmp    12c8 <NPdsymm_sum+0x1c8>
movslq -0x4(%rbp),%rax
mov    %rax,-0x60(%rbp)
mov    -0x60(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x40(%rbp),%rax
jae    13c8 <NPdsymm_sum+0x2c8>
mov    -0x28(%rbp),%rax
cmp    -0x38(%rbp),%rax
jbe    12f9 <NPdsymm_sum+0x1f9>
mov    -0x28(%rbp),%rax
mov    %rax,-0x68(%rbp)
jmp    1301 <NPdsymm_sum+0x201>
mov    -0x38(%rbp),%rax
mov    %rax,-0x68(%rbp)
mov    -0x68(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x40(%rbp),%rax
jae    13b2 <NPdsymm_sum+0x2b2>
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
movslq -0x4(%rbp),%rdx
imul   %rdx,%rcx
add    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movslq -0x4(%rbp),%rdx
imul   %rdx,%rcx
add    -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%rcx
movslq -0x4(%rbp),%rdx
imul   %rdx,%rcx
add    -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
movslq -0x4(%rbp),%rdx
imul   %rdx,%rcx
add    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1309 <NPdsymm_sum+0x209>
jmp    13b7 <NPdsymm_sum+0x2b7>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    12d0 <NPdsymm_sum+0x1d0>
jmp    13cd <NPdsymm_sum+0x2cd>
mov    -0x38(%rbp),%rax
add    $0x78,%rax
mov    %rax,-0x38(%rbp)
jmp    1281 <NPdsymm_sum+0x181>
jmp    13e3 <NPdsymm_sum+0x2e3>
pop    %rbp
ret

x86 assembly (clang, -O2):
<NPdsymm_sum>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rdx,-0x68(%rsp)
mov    %rsi,-0x70(%rsp)
and    $0xfffffffd,%ecx
movslq %edi,%r14
cmp    $0x1,%ecx
mov    %edi,-0x7c(%rsp)
jne    13d9 <NPdsymm_sum+0x2d9>
test   %edi,%edi
je     172c <NPdsymm_sum+0x62c>
mov    -0x68(%rsp),%rax
add    $0x10,%rax
mov    %rax,-0x18(%rsp)
lea    0x0(,%r14,8),%r11
mov    -0x70(%rsp),%rax
add    $0x10,%rax
mov    %rax,-0x20(%rsp)
lea    0x0(,%r14,4),%r13
mov    $0x78,%eax
xor    %ecx,%ecx
mov    %r11,-0x40(%rsp)
mov    %r13,(%rsp)
jmp    1197 <NPdsymm_sum+0x97>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    -0x8(%rsp),%rax
add    $0x78,%rax
mov    -0x10(%rsp),%rcx
cmp    %r14,%rcx
jae    172c <NPdsymm_sum+0x62c>
mov    %rcx,%rsi
cmp    %r14,%rax
mov    %r14,%rcx
mov    %rax,-0x8(%rsp)
cmovb  %rax,%rcx
mov    %rcx,-0x78(%rsp)
lea    0x78(%rsi),%rcx
cmp    %r14,%rcx
mov    %r14,%rax
mov    %rcx,-0x10(%rsp)
cmovb  %rcx,%rax
test   %rax,%rax
je     1180 <NPdsymm_sum+0x80>
mov    -0x70(%rsp),%rdx
mov    %rdx,%rcx
mov    -0x68(%rsp),%r8
mov    %r8,%rbp
mov    -0x20(%rsp),%rbx
mov    %rbx,%r9
mov    %rsi,0x10(%rsp)
mov    -0x18(%rsp),%rsi
mov    %rsi,%r10
mov    %rsi,-0x50(%rsp)
mov    %rbx,%rsi
xor    %ebx,%ebx
jmp    1232 <NPdsymm_sum+0x132>
nopw   0x0(%rax,%rax,1)
mov    -0x48(%rsp),%rbx
add    $0x1,%rbx
add    %r11,-0x50(%rsp)
add    %r11,%rsi
add    $0x8,%r10
add    $0x8,%r9
add    $0x8,%r8
add    $0x8,%rdx
add    %r11,%rbp
add    %r11,%rcx
cmp    -0x78(%rsp),%rbx
je     1180 <NPdsymm_sum+0x80>
mov    0x10(%rsp),%r12
cmp    %rbx,%r12
mov    %rbx,%r15
cmova  %r12,%r15
mov    %rbx,-0x48(%rsp)
cmovb  %rbx,%r12
cmp    %rax,%r12
jae    1200 <NPdsymm_sum+0x100>
mov    -0x78(%rsp),%rbx
sub    %r15,%rbx
cmp    $0x4,%rbx
jb     13a0 <NPdsymm_sum+0x2a0>
cmp    $0x1,%edi
jne    13a0 <NPdsymm_sum+0x2a0>
mov    %rbx,-0x58(%rsp)
mov    %rsi,-0x60(%rsp)
mov    -0x48(%rsp),%rdi
lea    (%r15,%rdi,1),%rbx
mov    -0x68(%rsp),%r11
lea    (%r11,%rbx,8),%rsi
mov    %rsi,-0x28(%rsp)
mov    -0x78(%rsp),%rsi
lea    (%rsi,%rdi,1),%r13
lea    (%r11,%r13,8),%rsi
mov    %rsi,0x8(%rsp)
mov    -0x70(%rsp),%rsi
lea    (%rsi,%rbx,8),%rbx
mov    %rbx,-0x30(%rsp)
lea    (%rsi,%r13,8),%rbx
imul   %r14,%rdi
lea    (%r15,%rdi,1),%r13
lea    (%r11,%r13,8),%rsi
mov    %rsi,-0x38(%rsp)
add    -0x78(%rsp),%rdi
lea    (%r11,%rdi,8),%r11
cmp    %rbx,-0x28(%rsp)
setb   %r13b
mov    -0x30(%rsp),%rsi
cmp    0x8(%rsp),%rsi
setb   %dil
cmp    %rbx,-0x38(%rsp)
setb   %bl
cmp    %r11,%rsi
setb   %r11b
test   %dil,%r13b
jne    1392 <NPdsymm_sum+0x292>
and    %r11b,%bl
jne    1392 <NPdsymm_sum+0x292>
mov    -0x50(%rsp),%rsi
lea    (%rsi,%r15,8),%rbx
mov    -0x60(%rsp),%rsi
lea    (%rsi,%r15,8),%r13
imul   %r14,%r15
mov    -0x58(%rsp),%r11
and    $0xfffffffffffffffc,%r11
add    %r11,%r12
xor    %edi,%edi
mov    (%rsp),%rsi
nopl   0x0(%rax)
movupd -0x10(%r13,%rdi,8),%xmm0
movupd 0x0(%r13,%rdi,8),%xmm1
movupd -0x10(%r9,%r15,8),%xmm2
addpd  %xmm0,%xmm2
movupd (%r9,%r15,8),%xmm0
addpd  %xmm1,%xmm0
movupd %xmm2,-0x10(%rbx,%rdi,8)
movupd %xmm0,(%rbx,%rdi,8)
movupd %xmm2,-0x10(%r10,%r15,8)
movupd %xmm0,(%r10,%r15,8)
add    $0x4,%rdi
add    %rsi,%r15
cmp    %rdi,%r11
jne    1330 <NPdsymm_sum+0x230>
cmp    %r11,-0x58(%rsp)
mov    -0x7c(%rsp),%edi
mov    -0x40(%rsp),%r11
mov    -0x60(%rsp),%rsi
je     1200 <NPdsymm_sum+0x100>
jmp    13a0 <NPdsymm_sum+0x2a0>
mov    -0x7c(%rsp),%edi
mov    -0x40(%rsp),%r11
mov    -0x60(%rsp),%rsi
mov    %r14,%rbx
imul   %r12,%rbx
nopw   0x0(%rax,%rax,1)
movsd  (%rcx,%r12,8),%xmm0
addsd  (%rdx,%rbx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%r12,8)
movsd  %xmm0,(%r8,%rbx,8)
add    $0x1,%r12
add    %r14,%rbx
cmp    %rax,%r12
jb     13b0 <NPdsymm_sum+0x2b0>
jmp    1200 <NPdsymm_sum+0x100>
test   %edi,%edi
je     172c <NPdsymm_sum+0x62c>
mov    %r14,%rax
shl    $0x5,%rax
mov    %rax,-0x20(%rsp)
lea    0x0(,%r14,8),%r11
mov    $0x78,%eax
xor    %ecx,%ecx
movapd 0xbfc(%rip),%xmm0
mov    %r11,-0x58(%rsp)
jmp    142a <NPdsymm_sum+0x32a>
nopl   0x0(%rax,%rax,1)
mov    -0x8(%rsp),%rax
add    $0x78,%rax
mov    -0x10(%rsp),%rcx
cmp    %r14,%rcx
mov    %r10d,%edi
jae    172c <NPdsymm_sum+0x62c>
mov    %edi,%r10d
cmp    %r14,%rax
mov    %r14,%r15
mov    %rax,-0x8(%rsp)
cmovb  %rax,%r15
lea    0x78(%rcx),%rax
cmp    %r14,%rax
mov    %r14,%r13
mov    %rax,-0x10(%rsp)
cmovb  %rax,%r13
test   %r13,%r13
je     1410 <NPdsymm_sum+0x310>
mov    %rcx,%rsi
lea    -0x2(%r15),%rax
mov    %rax,-0x18(%rsp)
mov    -0x70(%rsp),%rbp
mov    %rbp,%rcx
mov    -0x68(%rsp),%rdx
mov    %rdx,%rbx
xor    %eax,%eax
mov    %rax,-0x60(%rsp)
xor    %edi,%edi
mov    %rsi,-0x30(%rsp)
mov    %r15,-0x38(%rsp)
jmp    14b1 <NPdsymm_sum+0x3b1>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rdi
addq   $0x8,-0x60(%rsp)
add    %r11,%rdx
add    %r11,%rbp
add    $0x8,%rbx
add    $0x8,%rcx
cmp    %r15,%rdi
je     1410 <NPdsymm_sum+0x310>
cmp    %rdi,%rsi
mov    %rdi,%r9
cmova  %rsi,%r9
mov    %rsi,%r12
cmovb  %rdi,%r12
cmp    %r13,%r12
jae    1490 <NPdsymm_sum+0x390>
mov    %r15,%rax
sub    %r9,%rax
cmp    $0x2,%rax
jb     16f0 <NPdsymm_sum+0x5f0>
mov    %r10d,%r8d
cmp    $0x1,%r10d
jne    16f0 <NPdsymm_sum+0x5f0>
mov    %rax,-0x28(%rsp)
lea    (%r9,%rdi,1),%r11
mov    -0x68(%rsp),%rax
lea    (%rax,%r11,8),%rsi
mov    %rsi,-0x48(%rsp)
lea    (%r15,%rdi,1),%r8
lea    (%rax,%r8,8),%rsi
mov    %rsi,-0x78(%rsp)
mov    -0x70(%rsp),%rsi
lea    (%rsi,%r11,8),%r11
lea    (%rsi,%r8,8),%r10
mov    %rdi,%r8
imul   %r14,%r8
lea    (%r9,%r8,1),%rsi
lea    (%rax,%rsi,8),%rsi
mov    %rsi,-0x50(%rsp)
mov    %r8,0x8(%rsp)
lea    (%r15,%r8,1),%rsi
lea    (%rax,%rsi,8),%rsi
cmp    %r10,-0x48(%rsp)
setb   %r15b
cmp    -0x78(%rsp),%r11
setb   %r8b
cmp    %r10,-0x50(%rsp)
setb   %al
cmp    %rsi,%r11
setb   %sil
test   %r8b,%r15b
jne    158a <NPdsymm_sum+0x48a>
and    %sil,%al
jne    158a <NPdsymm_sum+0x48a>
mov    -0x28(%rsp),%rax
add    $0xfffffffffffffffe,%rax
mov    %rax,%rsi
shr    %rsi
add    $0x1,%rsi
cmp    $0x2,%rax
mov    %rsi,-0x40(%rsp)
mov    %rdi,(%rsp)
jae    15a3 <NPdsymm_sum+0x4a3>
xor    %eax,%eax
jmp    1671 <NPdsymm_sum+0x571>
mov    -0x7c(%rsp),%r10d
mov    -0x58(%rsp),%r11
mov    -0x30(%rsp),%rsi
mov    -0x38(%rsp),%r15
jmp    16f0 <NPdsymm_sum+0x5f0>
lea    0x2(%r9),%rax
imul   %r14,%rax
mov    -0x68(%rsp),%r11
lea    (%r11,%rax,8),%rsi
mov    %rsi,-0x48(%rsp)
mov    -0x70(%rsp),%rsi
lea    (%rsi,%rax,8),%rax
mov    %rax,-0x78(%rsp)
mov    -0x18(%rsp),%r15
sub    %r9,%r15
shr    %r15
add    $0x1,%r15
and    $0xfffffffffffffffe,%r15
mov    %r14,%rax
imul   %r9,%rax
lea    (%rsi,%rax,8),%rsi
mov    %rsi,-0x50(%rsp)
lea    (%rdx,%r9,8),%r10
lea    0x0(,%r9,8),%rsi
add    %rbp,%rsi
lea    (%r11,%rax,8),%rax
mov    %rax,0x10(%rsp)
mov    -0x60(%rsp),%r11
xor    %eax,%eax
mov    -0x20(%rsp),%r8
xchg   %ax,%ax
movupd (%rsi,%rax,8),%xmm1
mov    -0x50(%rsp),%rdi
movupd (%rdi,%r11,1),%xmm2
subpd  %xmm2,%xmm1
movupd %xmm1,(%r10,%rax,8)
xorpd  %xmm0,%xmm1
mov    0x10(%rsp),%rdi
movupd %xmm1,(%rdi,%r11,1)
movupd 0x10(%rsi,%rax,8),%xmm1
mov    -0x78(%rsp),%rdi
movupd (%rdi,%r11,1),%xmm2
subpd  %xmm2,%xmm1
movupd %xmm1,0x10(%r10,%rax,8)
xorpd  %xmm0,%xmm1
mov    -0x48(%rsp),%r9
movupd %xmm1,(%r9,%r11,1)
add    $0x4,%rax
add    %r8,%r11
add    $0xfffffffffffffffe,%r15
jne    1610 <NPdsymm_sum+0x510>
mov    -0x28(%rsp),%r8
mov    %r8,%r9
and    $0xfffffffffffffffe,%r9
testb  $0x1,-0x40(%rsp)
mov    (%rsp),%rdi
je     16c4 <NPdsymm_sum+0x5c4>
add    %r12,%rax
mov    0x8(%rsp),%rsi
add    %rax,%rsi
mov    -0x70(%rsp),%r10
movupd (%r10,%rsi,8),%xmm1
imul   %r14,%rax
add    %rdi,%rax
movupd (%r10,%rax,8),%xmm2
subpd  %xmm2,%xmm1
mov    -0x68(%rsp),%r10
movupd %xmm1,(%r10,%rsi,8)
xorpd  %xmm0,%xmm1
movupd %xmm1,(%r10,%rax,8)
cmp    %r9,%r8
mov    -0x7c(%rsp),%r10d
mov    -0x58(%rsp),%r11
mov    -0x30(%rsp),%rsi
mov    -0x38(%rsp),%r15
je     1490 <NPdsymm_sum+0x390>
add    %r9,%r12
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %r14,%rax
imul   %r12,%rax
nopw   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%r12,8),%xmm1
subsd  (%rcx,%rax,8),%xmm1
movsd  %xmm1,(%rdx,%r12,8)
xorpd  %xmm0,%xmm1
movlpd %xmm1,(%rbx,%rax,8)
add    $0x1,%r12
add    %r14,%rax
cmp    %r13,%r12
jb     1700 <NPdsymm_sum+0x600>
jmp    1490 <NPdsymm_sum+0x390>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @NPdsymm_sum(i32 %0, double* %1, double* %2, i32 %3) {
%5 = alloca i32
%6 = alloca double*
%7 = alloca double*
%8 = alloca i32
%9 = alloca i64
%10 = alloca i64
%11 = alloca i64
%12 = alloca i64
%13 = alloca double
store i32 %0, i32* %5
store double* %1, double** %6
store double* %2, double** %7
store i32 %3, i32* %8
%14 = load i32, i32* %8
%15 = icmp eq i32 %14, 1
br i1 %15, label %19, label %16
16: ; preds = %4
%17 = load i32, i32* %8
%18 = icmp eq i32 %17, 3
br i1 %18, label %19, label %107
19: ; preds = %16, %4
store i64 0, i64* %11
br label %20
20: ; preds = %103, %19
%21 = load i64, i64* %11
%22 = load i32, i32* %5
%23 = sext i32 %22 to i64
%24 = icmp ult i64 %21, %23
br i1 %24, label %25, label %106
25: ; preds = %20
store i64 0, i64* %9
%26 = load i64, i64* %11
%27 = add i64 %26, 120
%28 = load i32, i32* %5
%29 = sext i32 %28 to i64
%30 = icmp ult i64 %27, %29
br i1 %30, label %31, label %34
31: ; preds = %25
%32 = load i64, i64* %11
%33 = add i64 %32, 120
br label %37
34: ; preds = %25
%35 = load i32, i32* %5
%36 = sext i32 %35 to i64
br label %37
37: ; preds = %34, %31
%38 = phi i64 [ %33, %31 ], [ %36, %34 ]
store i64 %38, i64* %12
br label %39
39: ; preds = %99, %37
%40 = load i64, i64* %9
%41 = load i64, i64* %12
%42 = icmp ult i64 %40, %41
br i1 %42, label %43, label %102
43: ; preds = %39
%44 = load i64, i64* %9
%45 = load i64, i64* %11
%46 = icmp ugt i64 %44, %45
br i1 %46, label %47, label %49
47: ; preds = %43
%48 = load i64, i64* %9
br label %51
49: ; preds = %43
%50 = load i64, i64* %11
br label %51
51: ; preds = %49, %47
%52 = phi i64 [ %48, %47 ], [ %50, %49 ]
store i64 %52, i64* %10
br label %53
53: ; preds = %95, %51
%54 = load i64, i64* %10
%55 = load i64, i64* %12
%56 = icmp ult i64 %54, %55
br i1 %56, label %57, label %98
57: ; preds = %53
%58 = load double*, double** %6
%59 = load i64, i64* %9
%60 = load i32, i32* %5
%61 = sext i32 %60 to i64
%62 = mul i64 %59, %61
%63 = load i64, i64* %10
%64 = add i64 %62, %63
%65 = getelementptr inbounds double, double* %58, i64 %64
%66 = load double, double* %65
%67 = load double*, double** %6
%68 = load i64, i64* %10
%69 = load i32, i32* %5
%70 = sext i32 %69 to i64
%71 = mul i64 %68, %70
%72 = load i64, i64* %9
%73 = add i64 %71, %72
%74 = getelementptr inbounds double, double* %67, i64 %73
%75 = load double, double* %74
%76 = fadd double %66, %75
store double %76, double* %13
%77 = load double, double* %13
%78 = load double*, double** %7
%79 = load i64, i64* %9
%80 = load i32, i32* %5
%81 = sext i32 %80 to i64
%82 = mul i64 %79, %81
%83 = load i64, i64* %10
%84 = add i64 %82, %83
%85 = getelementptr inbounds double, double* %78, i64 %84
store double %77, double* %85
%86 = load double, double* %13
%87 = load double*, double** %7
%88 = load i64, i64* %10
%89 = load i32, i32* %5
%90 = sext i32 %89 to i64
%91 = mul i64 %88, %90
%92 = load i64, i64* %9
%93 = add i64 %91, %92
%94 = getelementptr inbounds double, double* %87, i64 %93
store double %86, double* %94
br label %95
95: ; preds = %57
%96 = load i64, i64* %10
%97 = add i64 %96, 1
store i64 %97, i64* %10
br label %53
98: ; preds = %53
br label %99
99: ; preds = %98
%100 = load i64, i64* %9
%101 = add i64 %100, 1
store i64 %101, i64* %9
br label %39
102: ; preds = %39
br label %103
103: ; preds = %102
%104 = load i64, i64* %11
%105 = add i64 %104, 120
store i64 %105, i64* %11
br label %20
106: ; preds = %20
br label %196
107: ; preds = %16
store i64 0, i64* %11
br label %108
108: ; preds = %192, %107
%109 = load i64, i64* %11
%110 = load i32, i32* %5
%111 = sext i32 %110 to i64
%112 = icmp ult i64 %109, %111
br i1 %112, label %113, label %195
113: ; preds = %108
store i64 0, i64* %9
%114 = load i64, i64* %11
%115 = add i64 %114, 120
%116 = load i32, i32* %5
%117 = sext i32 %116 to i64
%118 = icmp ult i64 %115, %117
br i1 %118, label %119, label %122
119: ; preds = %113
%120 = load i64, i64* %11
%121 = add i64 %120, 120
br label %125
122: ; preds = %113
%123 = load i32, i32* %5
%124 = sext i32 %123 to i64
br label %125
125: ; preds = %122, %119
%126 = phi i64 [ %121, %119 ], [ %124, %122 ]
store i64 %126, i64* %12
br label %127
127: ; preds = %188, %125
%128 = load i64, i64* %9
%129 = load i64, i64* %12
%130 = icmp ult i64 %128, %129
br i1 %130, label %131, label %191
131: ; preds = %127
%132 = load i64, i64* %9
%133 = load i64, i64* %11
%134 = icmp ugt i64 %132, %133
br i1 %134, label %135, label %137
135: ; preds = %131
%136 = load i64, i64* %9
br label %139
137: ; preds = %131
%138 = load i64, i64* %11
br label %139
139: ; preds = %137, %135
%140 = phi i64 [ %136, %135 ], [ %138, %137 ]
store i64 %140, i64* %10
br label %141
141: ; preds = %184, %139
%142 = load i64, i64* %10
%143 = load i64, i64* %12
%144 = icmp ult i64 %142, %143
br i1 %144, label %145, label %187
145: ; preds = %141
%146 = load double*, double** %6
%147 = load i64, i64* %9
%148 = load i32, i32* %5
%149 = sext i32 %148 to i64
%150 = mul i64 %147, %149
%151 = load i64, i64* %10
%152 = add i64 %150, %151
%153 = getelementptr inbounds double, double* %146, i64 %152
%154 = load double, double* %153
%155 = load double*, double** %6
%156 = load i64, i64* %10
%157 = load i32, i32* %5
%158 = sext i32 %157 to i64
%159 = mul i64 %156, %158
%160 = load i64, i64* %9
%161 = add i64 %159, %160
%162 = getelementptr inbounds double, double* %155, i64 %161
%163 = load double, double* %162
%164 = fsub double %154, %163
store double %164, double* %13
%165 = load double, double* %13
%166 = load double*, double** %7
%167 = load i64, i64* %9
%168 = load i32, i32* %5
%169 = sext i32 %168 to i64
%170 = mul i64 %167, %169
%171 = load i64, i64* %10
%172 = add i64 %170, %171
%173 = getelementptr inbounds double, double* %166, i64 %172
store double %165, double* %173
%174 = load double, double* %13
%175 = fneg double %174
%176 = load double*, double** %7
%177 = load i64, i64* %10
%178 = load i32, i32* %5
%179 = sext i32 %178 to i64
%180 = mul i64 %177, %179
%181 = load i64, i64* %9
%182 = add i64 %180, %181
%183 = getelementptr inbounds double, double* %176, i64 %182
store double %175, double* %183
br label %184
184: ; preds = %145
%185 = load i64, i64* %10
%186 = add i64 %185, 1
store i64 %186, i64* %10
br label %141
187: ; preds = %141
br label %188
188: ; preds = %187
%189 = load i64, i64* %9
%190 = add i64 %189, 1
store i64 %190, i64* %9
br label %127
191: ; preds = %127
br label %192
192: ; preds = %191
%193 = load i64, i64* %11
%194 = add i64 %193, 120
store i64 %194, i64* %11
br label %108
195: ; preds = %108
br label %196
196: ; preds = %195, %106
ret void
}

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