void SetUp1(npatch, loop, coeff, place, size) int npatch, loop[][2];
double coeff[][8192], place[][8192], size[][8192];
{
  int i, j, k, m, n, iface, ipatch, jface, jpatch;
  double d[2][2][2], d2[2][2][2], tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7,
      tmp8;
  for (iface = 0; iface < 3; iface++) {
    jface = iface + 3;
    tmp1 = place[2][loop[jface][0]] * place[2][loop[jface][0]];
    for (ipatch = loop[iface][0]; ipatch <= loop[iface][1]; ipatch++) {
      for (jpatch = loop[jface][0]; jpatch <= loop[jface][1]; jpatch++) {
        for (j = 0; j < 2; j++) {
          d[0][0][j] = place[j][jpatch] - place[j][ipatch];
          d[1][0][j] = d[0][0][j] + size[j][jpatch];
          d[0][1][j] = d[0][0][j] - size[j][ipatch];
          d[1][1][j] = d[1][0][j] - size[j][ipatch];
          d2[0][0][j] = d[0][0][j] * d[0][0][j];
          d2[1][0][j] = d[1][0][j] * d[1][0][j];
          d2[0][1][j] = d[0][1][j] * d[0][1][j];
          d2[1][1][j] = d[1][1][j] * d[1][1][j];
        }
        tmp2 = 0.0;
        for (j = 0; j < 2; j++) {
          for (m = 0; m < 2; m++) {
            for (i = 0; i < 2; i++) {
              tmp3 = d2[m][i][1 - j] + tmp1;
              tmp4 = sqrt(tmp3);
              tmp5 = 1.0 / tmp4;
              tmp8 = 0.0;
              for (k = 0; k < 2; k++) {
                for (n = 0; n < 2; n++) {
                  tmp6 = d2[k][n][j];
                  tmp7 = d[k][n][j];
                  tmp2 = (tmp3 - tmp6) * log(tmp6 + tmp3) - tmp2;
                  tmp8 = -tmp7 * atan(tmp7 * tmp5) - tmp8;
                }
                tmp2 = -tmp2;
                tmp8 = -tmp8;
              }
              tmp2 = -4.0 * tmp4 * tmp8 - tmp2;
            }
            tmp2 = -tmp2;
          }
        }
        coeff[ipatch][jpatch] = tmp2;
        coeff[jpatch][ipatch] = tmp2;
      }
    }
  }
}


x86 assembly (clang, -O0):
<SetUp1>:
push   %rbp
mov    %rsp,%rbp
sub    $0x120,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movl   $0x0,-0x40(%rbp)
cmpl   $0x3,-0x40(%rbp)
jge    16e1 <SetUp1+0x5b1>
mov    -0x40(%rbp),%eax
add    $0x3,%eax
mov    %eax,-0x48(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x48(%rbp),%rdx
movslq (%rcx,%rdx,8),%rcx
movsd  0x20000(%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x48(%rbp),%rdx
movslq (%rcx,%rdx,8),%rcx
mulsd  0x20000(%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
mov    (%rax,%rcx,8),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    -0x10(%rbp),%rcx
movslq -0x40(%rbp),%rdx
cmp    0x4(%rcx,%rdx,8),%eax
jg     16ce <SetUp1+0x59e>
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,8),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    -0x10(%rbp),%rcx
movslq -0x48(%rbp),%rdx
cmp    0x4(%rcx,%rdx,8),%eax
jg     16bb <SetUp1+0x58b>
movl   $0x0,-0x30(%rbp)
cmpl   $0x2,-0x30(%rbp)
jge    1358 <SetUp1+0x228>
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
shl    $0x10,%rcx
add    %rcx,%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
shl    $0x10,%rcx
add    %rcx,%rax
movslq -0x44(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movslq -0x30(%rbp),%rax
movsd  %xmm0,-0x90(%rbp,%rax,8)
ff ff
movslq -0x30(%rbp),%rax
movsd  -0x90(%rbp,%rax,8),%xmm0
ff ff
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
shl    $0x10,%rcx
add    %rcx,%rax
movslq -0x4c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movslq -0x30(%rbp),%rax
movsd  %xmm0,-0x70(%rbp,%rax,8)
movslq -0x30(%rbp),%rax
movsd  -0x90(%rbp,%rax,8),%xmm0
ff ff
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
shl    $0x10,%rcx
add    %rcx,%rax
movslq -0x44(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movslq -0x30(%rbp),%rax
movsd  %xmm0,-0x80(%rbp,%rax,8)
movslq -0x30(%rbp),%rax
movsd  -0x70(%rbp,%rax,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
shl    $0x10,%rcx
add    %rcx,%rax
movslq -0x44(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movslq -0x30(%rbp),%rax
movsd  %xmm0,-0x60(%rbp,%rax,8)
movslq -0x30(%rbp),%rax
movsd  -0x90(%rbp,%rax,8),%xmm0
ff ff
movslq -0x30(%rbp),%rax
mulsd  -0x90(%rbp,%rax,8),%xmm0
ff ff
movslq -0x30(%rbp),%rax
movsd  %xmm0,-0xd0(%rbp,%rax,8)
ff ff
movslq -0x30(%rbp),%rax
movsd  -0x70(%rbp,%rax,8),%xmm0
movslq -0x30(%rbp),%rax
mulsd  -0x70(%rbp,%rax,8),%xmm0
movslq -0x30(%rbp),%rax
movsd  %xmm0,-0xb0(%rbp,%rax,8)
ff ff
movslq -0x30(%rbp),%rax
movsd  -0x80(%rbp,%rax,8),%xmm0
movslq -0x30(%rbp),%rax
mulsd  -0x80(%rbp,%rax,8),%xmm0
movslq -0x30(%rbp),%rax
movsd  %xmm0,-0xc0(%rbp,%rax,8)
ff ff
movslq -0x30(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm0
movslq -0x30(%rbp),%rax
mulsd  -0x60(%rbp,%rax,8),%xmm0
movslq -0x30(%rbp),%rax
movsd  %xmm0,-0xa0(%rbp,%rax,8)
ff ff
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11ef <SetUp1+0xbf>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movl   $0x0,-0x30(%rbp)
cmpl   $0x2,-0x30(%rbp)
jge    166d <SetUp1+0x53d>
movl   $0x0,-0x38(%rbp)
cmpl   $0x2,-0x38(%rbp)
jge    165a <SetUp1+0x52a>
movl   $0x0,-0x2c(%rbp)
cmpl   $0x2,-0x2c(%rbp)
jge    1625 <SetUp1+0x4f5>
movslq -0x38(%rbp),%rcx
lea    -0xd0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
mov    $0x1,%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0xd8(%rbp),%xmm0
ff
movsd  %xmm0,-0xe8(%rbp)
ff
movsd  -0xe8(%rbp),%xmm0
ff
call   1050 <sqrt@plt>
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  0xc10(%rip),%xmm0
divsd  -0xf0(%rbp),%xmm0
ff
movsd  %xmm0,-0xf8(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x110(%rbp)
ff
movl   $0x0,-0x34(%rbp)
cmpl   $0x2,-0x34(%rbp)
jge    15d0 <SetUp1+0x4a0>
movl   $0x0,-0x3c(%rbp)
cmpl   $0x2,-0x3c(%rbp)
jge    1574 <SetUp1+0x444>
movslq -0x34(%rbp),%rcx
lea    -0xd0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x3c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
movslq -0x34(%rbp),%rcx
lea    -0x90(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x3c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
movsd  -0xe8(%rbp),%xmm0
ff
subsd  -0x100(%rbp),%xmm0
ff
movsd  %xmm0,-0x120(%rbp)
ff
movsd  -0x100(%rbp),%xmm0
ff
addsd  -0xe8(%rbp),%xmm0
ff
call   1040 <log@plt>
movaps %xmm0,%xmm2
movsd  -0x120(%rbp),%xmm0
ff
movsd  -0xe0(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  -0x108(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x118(%rbp)
ff
movsd  -0x108(%rbp),%xmm0
ff
mulsd  -0xf8(%rbp),%xmm0
ff
call   1030 <atan@plt>
movaps %xmm0,%xmm2
movsd  -0x118(%rbp),%xmm0
ff
movsd  -0x110(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x110(%rbp)
ff
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1423 <SetUp1+0x2f3>
movsd  -0xe0(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  -0x110(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x110(%rbp)
ff
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1412 <SetUp1+0x2e2>
movsd  0xa30(%rip),%xmm0
mulsd  -0xf0(%rbp),%xmm0
ff
movsd  -0x110(%rbp),%xmm2
ff
movsd  -0xe0(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    138c <SetUp1+0x25c>
movsd  -0xe0(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    137b <SetUp1+0x24b>
jmp    165f <SetUp1+0x52f>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    136a <SetUp1+0x23a>
movsd  -0xe0(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
shl    $0x10,%rcx
add    %rcx,%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xe0(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
shl    $0x10,%rcx
add    %rcx,%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11d3 <SetUp1+0xa3>
jmp    16c0 <SetUp1+0x590>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    11b0 <SetUp1+0x80>
jmp    16d3 <SetUp1+0x5a3>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1155 <SetUp1+0x25>
add    $0x120,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<SetUp1>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x1a8,%rsp
mov    %r8,0xa8(%rsp)
mov    %rcx,0x50(%rsp)
mov    %rdx,0xa0(%rsp)
mov    %rsi,0x28(%rsp)
xor    %edi,%edi
jmp    116e <SetUp1+0x3e>
nop
add    $0x1,%rdi
cmp    $0x3,%rdi
je     182c <SetUp1+0x6fc>
mov    0x28(%rsp),%rax
movslq (%rax,%rdi,8),%rbp
mov    0x4(%rax,%rdi,8),%eax
cmp    %eax,%ebp
jg     1160 <SetUp1+0x30>
mov    0x28(%rsp),%rsi
movslq 0x18(%rsi,%rdi,8),%rdx
mov    0x50(%rsp),%rcx
movsd  0x20000(%rcx,%rdx,8),%xmm0
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x58(%rsp)
mov    0x1c(%rsi,%rdi,8),%ecx
mov    %rdi,0xb0(%rsp)
cmp    %ecx,%edx
jg     1809 <SetUp1+0x6d9>
cs nopw 0x0(%rax,%rax,1)
nop
movslq %edx,%r12
mov    %rbp,0xb8(%rsp)
jmp    1222 <SetUp1+0xf2>
nopl   (%rax)
mov    0xb8(%rsp),%rbp
mov    %rbp,%rax
shl    $0x10,%rax
mov    0xa0(%rsp),%rcx
add    %rcx,%rax
movsd  %xmm1,(%rax,%r12,8)
mov    %r12,%rax
shl    $0x10,%rax
add    %rcx,%rax
movsd  %xmm1,(%rax,%rbp,8)
lea    0x1(%r12),%rax
mov    0x28(%rsp),%rcx
mov    0xb0(%rsp),%rdi
movslq 0x1c(%rcx,%rdi,8),%rcx
cmp    %rcx,%r12
mov    %rax,%r12
jge    1800 <SetUp1+0x6d0>
mov    0x50(%rsp),%rax
movsd  (%rax,%r12,8),%xmm1
movhpd 0x10000(%rax,%r12,8),%xmm1
movsd  (%rax,%rbp,8),%xmm0
movhpd 0x10000(%rax,%rbp,8),%xmm0
mov    0xa8(%rsp),%rax
movsd  (%rax,%r12,8),%xmm2
movhpd 0x10000(%rax,%r12,8),%xmm2
subpd  %xmm0,%xmm1
addpd  %xmm1,%xmm2
movlpd %xmm2,0x180(%rsp)
movsd  (%rax,%rbp,8),%xmm0
movhpd 0x10000(%rax,%rbp,8),%xmm0
movapd %xmm1,%xmm3
subpd  %xmm0,%xmm3
movapd %xmm3,0x170(%rsp)
movapd %xmm2,%xmm4
subpd  %xmm0,%xmm4
movhpd %xmm1,0x168(%rsp)
movapd %xmm2,%xmm0
shufpd $0x1,%xmm4,%xmm0
movupd %xmm0,0x188(%rsp)
movhpd %xmm4,0x198(%rsp)
movapd %xmm1,0x40(%rsp)
mulpd  %xmm1,%xmm1
movapd %xmm1,0x90(%rsp)
movapd %xmm1,0x120(%rsp)
movapd %xmm2,%xmm0
mulpd  %xmm2,%xmm0
movapd %xmm0,0x80(%rsp)
movapd %xmm0,0x140(%rsp)
movapd %xmm3,%xmm0
mulpd  %xmm3,%xmm0
movapd %xmm0,0x70(%rsp)
movapd %xmm0,0x130(%rsp)
movapd %xmm4,%xmm0
mulpd  %xmm4,%xmm0
movapd %xmm0,0x60(%rsp)
movapd %xmm0,0x150(%rsp)
xorpd  %xmm1,%xmm1
mov    $0x1,%r13b
xor    %r14d,%r14d
nopw   0x0(%rax,%rax,1)
mov    %r14,%r15
xor    $0x1,%r15
movapd %xmm3,0x100(%rsp)
movapd 0xcb8(%rip),%xmm0
xorpd  %xmm0,%xmm3
movapd %xmm3,0xe0(%rsp)
movapd %xmm2,0x110(%rsp)
xorpd  %xmm0,%xmm2
movapd %xmm2,0xd0(%rsp)
movapd %xmm4,0xf0(%rsp)
xorpd  %xmm0,%xmm4
movapd %xmm4,0xc0(%rsp)
mov    $0x1,%bl
xor    %eax,%eax
jmp    1550 <SetUp1+0x420>
nopw   0x0(%rax,%rax,1)
call   1050 <sqrt@plt>
movapd %xmm0,%xmm1
movapd %xmm1,0x30(%rsp)
movsd  0xc69(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x20(%rsp)
movsd  0x8(%rsp),%xmm0
movapd %xmm0,%xmm2
movapd 0x90(%rsp),%xmm1
subsd  %xmm1,%xmm2
movsd  %xmm2,(%rsp)
addsd  %xmm1,%xmm0
call   1040 <log@plt>
mulsd  (%rsp),%xmm0
subsd  0x18(%rsp),%xmm0
movsd  %xmm0,0x18(%rsp)
movsd  0x20(%rsp),%xmm0
mulsd  0x40(%rsp),%xmm0
call   1030 <atan@plt>
movsd  %xmm0,(%rsp)
movsd  0x8(%rsp),%xmm0
movapd %xmm0,%xmm2
movapd 0x70(%rsp),%xmm1
subsd  %xmm1,%xmm2
movsd  %xmm2,0x10(%rsp)
addsd  %xmm1,%xmm0
call   1040 <log@plt>
mulsd  0x10(%rsp),%xmm0
subsd  0x18(%rsp),%xmm0
movsd  %xmm0,0x18(%rsp)
movsd  0x20(%rsp),%xmm0
mulsd  0x100(%rsp),%xmm0
call   1030 <atan@plt>
movsd  (%rsp),%xmm1
mulsd  0x40(%rsp),%xmm1
mulsd  0xe0(%rsp),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rsp)
movsd  0x8(%rsp),%xmm0
movapd %xmm0,%xmm2
movapd 0x80(%rsp),%xmm1
subsd  %xmm1,%xmm2
movsd  %xmm2,0x10(%rsp)
addsd  %xmm1,%xmm0
call   1040 <log@plt>
mulsd  0x10(%rsp),%xmm0
addsd  0x18(%rsp),%xmm0
movsd  %xmm0,0x18(%rsp)
movsd  0x20(%rsp),%xmm0
mulsd  0x110(%rsp),%xmm0
call   1030 <atan@plt>
mulsd  0xd0(%rsp),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x8(%rsp),%xmm1
movapd %xmm1,%xmm2
movapd 0x60(%rsp),%xmm0
subsd  %xmm0,%xmm2
movsd  %xmm2,0x10(%rsp)
addsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
call   1040 <log@plt>
mulsd  0x10(%rsp),%xmm0
subsd  0x18(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x20(%rsp),%xmm0
mulsd  0xf0(%rsp),%xmm0
call   1030 <atan@plt>
movapd 0x30(%rsp),%xmm1
mulsd  0xc0(%rsp),%xmm0
subsd  (%rsp),%xmm0
mulsd  0xaed(%rip),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x8(%rsp),%xmm1
xorpd  0xac3(%rip),%xmm1
mov    $0x1,%eax
test   $0x1,%bl
mov    $0x0,%ebx
je     1770 <SetUp1+0x640>
shl    $0x5,%rax
lea    (%rsp,%rax,1),%rbp
add    $0x120,%rbp
movsd  0x0(%rbp,%r15,8),%xmm0
addsd  0x58(%rsp),%xmm0
ucomisd 0xaac(%rip),%xmm0
movapd %xmm1,0x30(%rsp)
movsd  %xmm0,0x8(%rsp)
jb     1590 <SetUp1+0x460>
xorps  %xmm1,%xmm1
sqrtsd %xmm0,%xmm1
jmp    1599 <SetUp1+0x469>
nopl   0x0(%rax,%rax,1)
call   1050 <sqrt@plt>
movapd %xmm0,%xmm1
movsd  %xmm1,0x18(%rsp)
movsd  0xa69(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x20(%rsp)
movsd  0x8(%rsp),%xmm0
movapd %xmm0,%xmm2
movapd 0x90(%rsp),%xmm1
subsd  %xmm1,%xmm2
movsd  %xmm2,(%rsp)
addsd  %xmm1,%xmm0
call   1040 <log@plt>
mulsd  (%rsp),%xmm0
subsd  0x30(%rsp),%xmm0
movsd  %xmm0,0x30(%rsp)
movsd  0x20(%rsp),%xmm0
mulsd  0x40(%rsp),%xmm0
call   1030 <atan@plt>
movsd  %xmm0,(%rsp)
movsd  0x8(%rsp),%xmm0
movapd %xmm0,%xmm2
movapd 0x70(%rsp),%xmm1
subsd  %xmm1,%xmm2
movsd  %xmm2,0x10(%rsp)
addsd  %xmm1,%xmm0
call   1040 <log@plt>
mulsd  0x10(%rsp),%xmm0
subsd  0x30(%rsp),%xmm0
movsd  %xmm0,0x30(%rsp)
movsd  0x20(%rsp),%xmm0
mulsd  0x100(%rsp),%xmm0
call   1030 <atan@plt>
movsd  (%rsp),%xmm1
mulsd  0x40(%rsp),%xmm1
mulsd  0xe0(%rsp),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rsp)
movsd  0x8(%rsp),%xmm0
movapd %xmm0,%xmm2
movapd 0x80(%rsp),%xmm1
subsd  %xmm1,%xmm2
movsd  %xmm2,0x10(%rsp)
addsd  %xmm1,%xmm0
call   1040 <log@plt>
mulsd  0x10(%rsp),%xmm0
addsd  0x30(%rsp),%xmm0
movsd  %xmm0,0x30(%rsp)
movsd  0x20(%rsp),%xmm0
mulsd  0x110(%rsp),%xmm0
call   1030 <atan@plt>
mulsd  0xd0(%rsp),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x8(%rsp),%xmm1
movapd %xmm1,%xmm2
movapd 0x60(%rsp),%xmm0
subsd  %xmm0,%xmm2
movsd  %xmm2,0x10(%rsp)
addsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
call   1040 <log@plt>
mulsd  0x10(%rsp),%xmm0
subsd  0x30(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x20(%rsp),%xmm0
mulsd  0xf0(%rsp),%xmm0
call   1030 <atan@plt>
movsd  0x18(%rsp),%xmm1
mulsd  0xc0(%rsp),%xmm0
subsd  (%rsp),%xmm0
mulsd  0x8ed(%rip),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x8(%rsp),%xmm1
movsd  0x10(%rbp,%r15,8),%xmm0
addsd  0x58(%rsp),%xmm0
ucomisd 0x8d6(%rip),%xmm0
movsd  %xmm1,0x18(%rsp)
movsd  %xmm0,0x8(%rsp)
jb     1390 <SetUp1+0x260>
xorps  %xmm1,%xmm1
sqrtsd %xmm0,%xmm1
jmp    1399 <SetUp1+0x269>
nopl   0x0(%rax,%rax,1)
test   $0x1,%r13b
je     11d0 <SetUp1+0xa0>
movsd  0x128(%rsp,%r14,8),%xmm0
movaps %xmm0,0x90(%rsp)
movsd  0x138(%rsp,%r14,8),%xmm0
movaps %xmm0,0x70(%rsp)
movsd  0x168(%rsp,%r14,8),%xmm0
movaps %xmm0,0x40(%rsp)
movsd  0x178(%rsp,%r14,8),%xmm3
movsd  0x148(%rsp,%r14,8),%xmm0
movaps %xmm0,0x80(%rsp)
movsd  0x188(%rsp,%r14,8),%xmm2
movsd  0x158(%rsp,%r14,8),%xmm0
movapd %xmm0,0x60(%rsp)
movsd  0x198(%rsp,%r14,8),%xmm4
mov    $0x1,%r14d
xor    %r13d,%r13d
jmp    1330 <SetUp1+0x200>
nopl   0x0(%rax,%rax,1)
mov    0x28(%rsp),%rax
mov    0x4(%rax,%rdi,8),%eax
movslq %eax,%rdx
cmp    %rdx,%rbp
jge    1160 <SetUp1+0x30>
add    $0x1,%rbp
mov    0x28(%rsp),%rdx
mov    0x18(%rdx,%rdi,8),%edx
cmp    %ecx,%edx
jle    11c0 <SetUp1+0x90>
jmp    1809 <SetUp1+0x6d9>
add    $0x1a8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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