#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double exp(double);
int /*<<< orphan*/ free(double *);
scalar_t__ malloc(int);
double sqrt(double);
double *getP0(double *x, double *y, double *z, double *r2, double *alphas,
              double *betas, int Asize, int Nsize) {
  double *P0 = (double *)malloc(Nsize * Nsize * sizeof(double));
  if (Asize == 0) {
    for (int n = 0; n < Nsize * Nsize; n++) {
      P0[n] = 0.0;
    }
    return P0;
  }
  double sumsInner = 0;
  double sumsOuter = 0;
  double oneO1Palpha;
  double oneO1PalphaSqrt;
  double *alphaO1Palpha = (double *)malloc(Nsize * sizeof(double));
  double *oneO1PalphaSqrtCubed = (double *)malloc(Nsize * sizeof(double));
  double *kkdij =
      (double *)malloc(Nsize * Nsize * Asize * Asize * sizeof(double));
  for (int n = 0; n < Nsize; n++) {
    oneO1Palpha = 1.0 / (1.00000000 + alphas[n]);
    oneO1PalphaSqrt = sqrt(oneO1Palpha);
    oneO1PalphaSqrtCubed[n] = oneO1PalphaSqrt * oneO1Palpha;
    alphaO1Palpha[n] = alphas[n] * oneO1Palpha;
  }
  for (int k = 0; k < Nsize; k++) {
    for (int kd = 0; kd < Nsize; kd++) {
      for (int i = 0; i < Asize; i++) {
        for (int j = 0; j < Asize; j++) {
          kkdij[k * Nsize * Asize * Asize + kd * Asize * Asize + i * Asize +
                j] = exp(-alphaO1Palpha[k] * r2[i] - alphaO1Palpha[kd] * r2[j]);
        }
      }
    }
  }
  for (int n = 0; n < Nsize; n++) {
    for (int nd = n; nd < Nsize; nd++) {
      sumsOuter = 0;
      for (int k = 0; k < Nsize; k++) {
        for (int kd = 0; kd < Nsize; kd++) {
          sumsInner = 0;
          for (int i = 0; i < Asize; i++) {
            for (int j = 0; j < Asize; j++) {
              sumsInner += kkdij[k * Nsize * Asize * Asize +
                                 kd * Asize * Asize + i * Asize + j];
            }
          }
          sumsOuter += betas[n * Nsize + k] * betas[nd * Nsize + kd] *
                       oneO1PalphaSqrtCubed[k] * oneO1PalphaSqrtCubed[kd] *
                       sumsInner;
        }
      }
      P0[Nsize * n + nd] = 9.86960440108936 * 0.25 * sumsOuter;
    }
  }
  for (int n = 0; n < Nsize; n++) {
    for (int nd = n + 1; nd < Nsize; nd++) {
      P0[Nsize * nd + n] = P0[Nsize * n + nd];
    }
  }
  free(oneO1PalphaSqrtCubed);
  free(alphaO1Palpha);
  free(kkdij);
  return P0;
}


x86 assembly (clang, -O0):
<getP0>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
mov    0x18(%rbp),%eax
imul   0x18(%rbp),%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x40(%rbp)
cmpl   $0x0,0x10(%rbp)
jne    11cf <getP0+0x8f>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    0x18(%rbp),%ecx
imul   0x18(%rbp),%ecx
cmp    %ecx,%eax
jge    11c2 <getP0+0x82>
mov    -0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1192 <getP0+0x52>
mov    -0x40(%rbp),%rax
mov    %rax,-0x8(%rbp)
jmp    1717 <getP0+0x5d7>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movslq 0x18(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x70(%rbp)
movslq 0x18(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x78(%rbp)
mov    0x18(%rbp),%eax
imul   0x18(%rbp),%eax
imul   0x10(%rbp),%eax
imul   0x10(%rbp),%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x80(%rbp)
movl   $0x0,-0x84(%rbp)
mov    -0x84(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    12c9 <getP0+0x189>
mov    -0x30(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  0xdb7(%rip),%xmm1
addsd  (%rax,%rcx,8),%xmm1
movsd  0xdaa(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
call   1060 <sqrt@plt>
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
mov    -0x78(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x60(%rbp),%xmm0
mov    -0x70(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    122f <getP0+0xef>
movl   $0x0,-0x88(%rbp)
mov    -0x88(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    1447 <getP0+0x307>
movl   $0x0,-0x8c(%rbp)
mov    -0x8c(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    142e <getP0+0x2ee>
movl   $0x0,-0x90(%rbp)
mov    -0x90(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    1415 <getP0+0x2d5>
movl   $0x0,-0x94(%rbp)
mov    -0x94(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    13fc <getP0+0x2bc>
mov    -0x70(%rbp),%rax
movslq -0x88(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x90(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x70(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x28(%rbp),%rax
movslq -0x94(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1050 <exp@plt>
mov    -0x80(%rbp),%rax
mov    -0x88(%rbp),%ecx
imul   0x18(%rbp),%ecx
imul   0x10(%rbp),%ecx
imul   0x10(%rbp),%ecx
mov    -0x8c(%rbp),%edx
imul   0x10(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
mov    -0x90(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x94(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    131e <getP0+0x1de>
jmp    1401 <getP0+0x2c1>
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    1305 <getP0+0x1c5>
jmp    141a <getP0+0x2da>
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
jmp    12ec <getP0+0x1ac>
jmp    1433 <getP0+0x2f3>
mov    -0x88(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x88(%rbp)
jmp    12d3 <getP0+0x193>
movl   $0x0,-0x98(%rbp)
mov    -0x98(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    1658 <getP0+0x518>
mov    -0x98(%rbp),%eax
mov    %eax,-0x9c(%rbp)
mov    -0x9c(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    163f <getP0+0x4ff>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0xa0(%rbp)
mov    -0xa0(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    1602 <getP0+0x4c2>
movl   $0x0,-0xa4(%rbp)
mov    -0xa4(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    15e9 <getP0+0x4a9>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0xa8(%rbp)
mov    -0xa8(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    1566 <getP0+0x426>
movl   $0x0,-0xac(%rbp)
mov    -0xac(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    154d <getP0+0x40d>
mov    -0x80(%rbp),%rax
mov    -0xa0(%rbp),%ecx
imul   0x18(%rbp),%ecx
imul   0x10(%rbp),%ecx
imul   0x10(%rbp),%ecx
mov    -0xa4(%rbp),%edx
imul   0x10(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
mov    -0xa8(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0xac(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0xac(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xac(%rbp)
jmp    14e0 <getP0+0x3a0>
jmp    1552 <getP0+0x412>
mov    -0xa8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa8(%rbp)
jmp    14c7 <getP0+0x387>
mov    -0x38(%rbp),%rax
mov    -0x98(%rbp),%ecx
imul   0x18(%rbp),%ecx
add    -0xa0(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x9c(%rbp),%ecx
imul   0x18(%rbp),%ecx
add    -0xa4(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x78(%rbp),%rax
movslq -0xa0(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x78(%rbp),%rax
movslq -0xa4(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0xa4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa4(%rbp)
jmp    14a6 <getP0+0x366>
jmp    15ee <getP0+0x4ae>
mov    -0xa0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa0(%rbp)
jmp    148d <getP0+0x34d>
movsd  0x9f6(%rip),%xmm0
mulsd  -0x58(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    0x18(%rbp),%ecx
imul   -0x98(%rbp),%ecx
add    -0x9c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x9c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x9c(%rbp)
jmp    146c <getP0+0x32c>
jmp    1644 <getP0+0x504>
mov    -0x98(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x98(%rbp)
jmp    1451 <getP0+0x311>
movl   $0x0,-0xb0(%rbp)
mov    -0xb0(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    16f4 <getP0+0x5b4>
mov    -0xb0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xb4(%rbp)
mov    -0xb4(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    16db <getP0+0x59b>
mov    -0x40(%rbp),%rax
mov    0x18(%rbp),%ecx
imul   -0xb0(%rbp),%ecx
add    -0xb4(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
mov    0x18(%rbp),%ecx
imul   -0xb4(%rbp),%ecx
add    -0xb0(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xb4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xb4(%rbp)
jmp    1680 <getP0+0x540>
jmp    16e0 <getP0+0x5a0>
mov    -0xb0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xb0(%rbp)
jmp    1662 <getP0+0x522>
mov    -0x78(%rbp),%rdi
call   1030 <free@plt>
mov    -0x70(%rbp),%rdi
call   1030 <free@plt>
mov    -0x80(%rbp),%rdi
call   1030 <free@plt>
mov    -0x40(%rbp),%rax
mov    %rax,-0x8(%rbp)
mov    -0x8(%rbp),%rax
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<getP0>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %r9,0x58(%rsp)
mov    %r8,%r13
mov    %rcx,%rbp
mov    0xa8(%rsp),%ebx
mov    0xa0(%rsp),%r15d
imul   %ebx,%ebx
lea    0x0(,%rbx,8),%edi
call   1050 <malloc@plt>
mov    %rax,0x20(%rsp)
test   %r15d,%r15d
je     13b8 <getP0+0x268>
mov    %rbp,0x28(%rsp)
mov    0xa8(%rsp),%eax
lea    0x0(,%rax,8),%r14d
mov    %r14d,%edi
call   1050 <malloc@plt>
mov    %rax,%r12
mov    %r14d,%edi
mov    0xa8(%rsp),%ebp
call   1050 <malloc@plt>
mov    %rax,0x30(%rsp)
mov    %r15d,%eax
imul   %r15d,%eax
mov    %eax,0x4(%rsp)
imul   %eax,%ebx
shl    $0x3,%ebx
mov    %ebx,%edi
call   1050 <malloc@plt>
mov    %rax,%rbx
test   %ebp,%ebp
jle    1628 <getP0+0x4d8>
mov    %ebp,%ecx
xor    %r14d,%r14d
movsd  0xe08(%rip),%xmm1
xorpd  %xmm2,%xmm2
mov    %rcx,0x50(%rsp)
mov    %r13,%rdx
jmp    125d <getP0+0x10d>
cs nopw 0x0(%rax,%rax,1)
movapd %xmm3,%xmm0
movsd  %xmm3,0x10(%rsp)
call   1070 <sqrt@plt>
movsd  0x10(%rsp),%xmm3
xorpd  %xmm2,%xmm2
movsd  0xdcf(%rip),%xmm1
mov    %r13,%rdx
mov    0x50(%rsp),%rcx
mulsd  %xmm3,%xmm0
mov    0x30(%rsp),%rax
movsd  %xmm0,(%rax,%r14,8)
mulsd  (%rdx,%r14,8),%xmm3
movsd  %xmm3,(%r12,%r14,8)
add    $0x1,%r14
cmp    %r14,%rcx
je     127e <getP0+0x12e>
movsd  (%rdx,%r14,8),%xmm0
addsd  %xmm1,%xmm0
movapd %xmm1,%xmm3
divsd  %xmm0,%xmm3
ucomisd %xmm2,%xmm3
jb     1210 <getP0+0xc0>
xorps  %xmm0,%xmm0
sqrtsd %xmm3,%xmm0
jmp    1239 <getP0+0xe9>
mov    0xa8(%rsp),%eax
test   %eax,%eax
jle    1628 <getP0+0x4d8>
mov    %r12,0x48(%rsp)
mov    %r15d,%r13d
imul   0x4(%rsp),%eax
mov    %eax,0x1c(%rsp)
xor    %eax,%eax
movapd 0xd68(%rip),%xmm2
xor    %ecx,%ecx
mov    %rbx,0x60(%rsp)
mov    0x28(%rsp),%rbp
jmp    12d6 <getP0+0x186>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%rcx
mov    0x38(%rsp),%rax
add    0x1c(%rsp),%eax
cmp    %r11,%rcx
je     13d7 <getP0+0x287>
mov    %rax,0x38(%rsp)
mov    %eax,0x8(%rsp)
xor    %edx,%edx
jmp    130a <getP0+0x1ba>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rdx
mov    0x8(%rsp),%eax
add    0x4(%rsp),%eax
mov    %eax,0x8(%rsp)
mov    0x50(%rsp),%r11
cmp    %r11,%rdx
je     12c0 <getP0+0x170>
test   %r15d,%r15d
jle    12f0 <getP0+0x1a0>
mov    0x8(%rsp),%eax
xor    %r14d,%r14d
cs nopw 0x0(%rax,%rax,1)
mov    %eax,0x40(%rsp)
cltq
lea    (%rbx,%rax,8),%rax
mov    %rax,0x10(%rsp)
xor    %r12d,%r12d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x48(%rsp),%rax
movsd  (%rax,%rcx,8),%xmm0
movsd  0x0(%rbp,%r12,8),%xmm1
mulsd  (%rax,%rdx,8),%xmm1
xorpd  %xmm2,%xmm0
mulsd  0x0(%rbp,%r14,8),%xmm0
subsd  %xmm1,%xmm0
mov    %rcx,%r15
mov    %rdx,%rbx
call   1060 <exp@plt>
mov    %rbx,%rdx
mov    %r15,%rcx
movapd 0xc92(%rip),%xmm2
mov    0x10(%rsp),%rax
movsd  %xmm0,(%rax,%r12,8)
add    $0x1,%r12
cmp    %r12,%r13
jne    1340 <getP0+0x1f0>
add    $0x1,%r14
mov    0xa0(%rsp),%r15d
mov    0x40(%rsp),%eax
add    %r15d,%eax
cmp    %r13,%r14
mov    0x60(%rsp),%rbx
jne    1320 <getP0+0x1d0>
jmp    12f0 <getP0+0x1a0>
test   %ebx,%ebx
je     1642 <getP0+0x4f2>
mov    %ebx,%edx
shl    $0x3,%rdx
mov    0x20(%rsp),%rdi
xor    %esi,%esi
call   1040 <memset@plt>
jmp    1642 <getP0+0x4f2>
mov    0xa8(%rsp),%eax
test   %eax,%eax
mov    0x48(%rsp),%r12
jle    1628 <getP0+0x4d8>
lea    -0x1(%r13),%rax
mov    %r13d,%ecx
and    $0x3,%ecx
mov    %r13d,%edx
and    $0xfffffffc,%edx
lea    0x18(%rbx),%rsi
xor    %edi,%edi
movsd  0xbff(%rip),%xmm0
jmp    1422 <getP0+0x2d2>
nopl   0x0(%rax,%rax,1)
mov    0x38(%rsp),%rdi
add    $0x1,%rdi
cmp    %r11,%rdi
je     15b0 <getP0+0x460>
mov    %rdi,%rbp
imul   %r11,%rbp
mov    %rbp,0x28(%rsp)
mov    %rdi,0x38(%rsp)
mov    %rdi,%r12
jmp    1465 <getP0+0x315>
nopl   0x0(%rax,%rax,1)
mulsd  %xmm0,%xmm1
mov    0x28(%rsp),%rdi
mov    0x8(%rsp),%r12
add    %r12,%rdi
mov    0x20(%rsp),%rbp
movsd  %xmm1,0x0(%rbp,%rdi,8)
add    $0x1,%r12
cmp    %r11,%r12
je     1410 <getP0+0x2c0>
mov    %r12,0x8(%rsp)
imul   %r11,%r12
xorpd  %xmm1,%xmm1
xor    %r8d,%r8d
xor    %r9d,%r9d
jmp    1498 <getP0+0x348>
nopw   0x0(%rax,%rax,1)
mov    0x40(%rsp),%r9
add    $0x1,%r9
mov    0x10(%rsp),%r8
add    0x1c(%rsp),%r8d
cmp    %r11,%r9
je     1440 <getP0+0x2f0>
mov    0x28(%rsp),%rdi
add    %r9,%rdi
mov    0x58(%rsp),%rbp
movsd  0x0(%rbp,%rdi,8),%xmm2
mov    0x30(%rsp),%rdi
mov    %r9,0x40(%rsp)
movsd  (%rdi,%r9,8),%xmm3
mov    %r8,0x10(%rsp)
mov    %r8d,%r10d
xor    %r9d,%r9d
jmp    1511 <getP0+0x3c1>
nopl   0x0(%rax,%rax,1)
lea    (%r9,%r12,1),%rdi
mov    0x58(%rsp),%rbp
movsd  0x0(%rbp,%rdi,8),%xmm5
mulsd  %xmm2,%xmm5
mulsd  %xmm3,%xmm5
mov    0x30(%rsp),%rdi
mulsd  (%rdi,%r9,8),%xmm5
mulsd  %xmm4,%xmm5
addsd  %xmm5,%xmm1
add    $0x1,%r9
add    0x4(%rsp),%r10d
mov    0x50(%rsp),%r11
cmp    %r11,%r9
je     1480 <getP0+0x330>
xorpd  %xmm4,%xmm4
test   %r15d,%r15d
jle    14d0 <getP0+0x380>
mov    %r10d,%r8d
xor    %r11d,%r11d
jmp    153c <getP0+0x3ec>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r11
add    %r15d,%r8d
cmp    %r13,%r11
je     14d0 <getP0+0x380>
movslq %r8d,%rdi
cmp    $0x3,%rax
jae    1550 <getP0+0x400>
xor    %ebp,%ebp
jmp    1584 <getP0+0x434>
nopl   0x0(%rax)
lea    (%rsi,%rdi,8),%r14
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
addsd  -0x18(%r14,%rbp,8),%xmm4
addsd  -0x10(%r14,%rbp,8),%xmm4
addsd  -0x8(%r14,%rbp,8),%xmm4
addsd  (%r14,%rbp,8),%xmm4
add    $0x4,%rbp
cmp    %rbp,%rdx
jne    1560 <getP0+0x410>
test   %rcx,%rcx
je     1530 <getP0+0x3e0>
add    %rdi,%rbp
lea    (%rbx,%rbp,8),%rdi
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
addsd  (%rdi,%rbp,8),%xmm4
add    $0x1,%rbp
cmp    %rbp,%rcx
jne    15a0 <getP0+0x450>
jmp    1530 <getP0+0x3e0>
mov    0xa8(%rsp),%eax
test   %eax,%eax
mov    0x48(%rsp),%r12
jle    1628 <getP0+0x4d8>
mov    0x20(%rsp),%rbp
lea    0x0(,%r11,8),%rsi
add    %rbp,%rsi
lea    0x0(,%r11,8),%rcx
lea    0x8(,%r11,8),%r8
xor    %r9d,%r9d
mov    %r11,%rdi
jmp    15ff <getP0+0x4af>
nopl   0x0(%rax,%rax,1)
add    %r8,%rsi
add    %r8,%rbp
add    $0xffffffffffffffff,%rdi
cmp    %r11,%r9
je     1628 <getP0+0x4d8>
add    $0x1,%r9
cmp    %r11,%r9
jae    15f0 <getP0+0x4a0>
mov    $0x1,%edx
mov    %rsi,%rax
movsd  0x0(%rbp,%rdx,8),%xmm0
movsd  %xmm0,(%rax)
add    %rcx,%rax
add    $0x1,%rdx
cmp    %rdx,%rdi
jne    1610 <getP0+0x4c0>
jmp    15f0 <getP0+0x4a0>
mov    0x30(%rsp),%rdi
call   1030 <free@plt>
mov    %r12,%rdi
call   1030 <free@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
mov    0x20(%rsp),%rax
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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