#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 abs(double);
double pow(double, int);
void calculate_fg(double Re, double GX, double GY, double alpha, double dt,
                  double dx, double dy, int imax, int jmax, double **U,
                  double **V, double **F, double **G) {
  double uij, vij, uuj, udj, uiu, uid, viu, vid, vuj, vdj, udu, vud;
  double Dx = 1 / dx, Dy = 1 / dy;
  int i, j;
  for (i = 1; i < imax + 1; i++) {
    for (j = 1; j < jmax + 1; j++) {
      uij = U[i][j];
      uuj = U[i + 1][j];
      udj = U[i - 1][j];
      uiu = U[i][j + 1];
      uid = U[i][j - 1];
      udu = U[i - 1][j + 1];
      vij = V[i][j];
      viu = V[i][j + 1];
      vid = V[i][j - 1];
      vuj = V[i + 1][j];
      vdj = V[i - 1][j];
      vud = V[i + 1][j - 1];
      F[i][j] = uij + dt * (1 / Re *
                                ((uuj - 2 * uij + udj) * Dx * Dx +
                                 (uiu - 2 * uij + uid) * Dy * Dy) -
                            0.25 * Dx *
                                (pow((uij + uuj), 2) - pow((udj + uij), 2) +
                                 alpha * (abs(uij + uuj) * (uij - uuj) -
                                          abs(udj + uij) * (udj - uij))) -
                            0.25 * Dy *
                                ((vij + vuj) * (uij + uiu) -
                                 (vid + vud) * (uid + uij) +
                                 alpha * (abs(vij + vuj) * (uij - uiu) -
                                          abs(vid + vud) * (uid - uij))) +
                            GX);
      G[i][j] = vij + dt * (1 / Re *
                                ((vuj - 2 * vij + vdj) * Dx * Dx +
                                 (viu - 2 * vij + vid) * Dy * Dy) -
                            0.25 * Dy *
                                (pow((vij + viu), 2) - pow((vid + vij), 2) +
                                 alpha * (abs(vij + viu) * (vij - viu) -
                                          abs(vid + vij) * (vid - vij))) -
                            0.25 * Dx *
                                ((uij + uiu) * (vij + vuj) -
                                 (udj + udu) * (vdj + vij) +
                                 alpha * (abs(uij + uiu) * (vij - vuj) -
                                          abs(udj + udu) * (vdj - vij))) +
                            GY);
    }
    G[i][0] = V[i][0];
    G[i][jmax] = V[i][jmax];
  }
  for (j = 1; j < jmax + 1; j++) {
    F[0][j] = U[0][j];
    F[imax][j] = U[imax][j];
  }
}


x86 assembly (clang, -O0):
<calculate_fg>:
push   %rbp
mov    %rsp,%rbp
sub    $0x1e0,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
movsd  %xmm3,-0x20(%rbp)
movsd  %xmm4,-0x28(%rbp)
movsd  %xmm5,-0x30(%rbp)
movsd  %xmm6,-0x38(%rbp)
mov    %edi,-0x3c(%rbp)
mov    %esi,-0x40(%rbp)
mov    %rdx,-0x48(%rbp)
mov    %rcx,-0x50(%rbp)
mov    %r8,-0x58(%rbp)
mov    %r9,-0x60(%rbp)
movsd  0xe94(%rip),%xmm0
divsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  0xe7f(%rip),%xmm0
divsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
movl   $0x1,-0xd4(%rbp)
mov    -0xd4(%rbp),%eax
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1afc <calculate_fg+0x9dc>
movl   $0x1,-0xd8(%rbp)
mov    -0xd8(%rbp),%eax
mov    -0x40(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1a92 <calculate_fg+0x972>
mov    -0x48(%rbp),%rax
movslq -0xd4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xd8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x48(%rbp),%rax
mov    -0xd4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xd8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x48(%rbp),%rax
mov    -0xd4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xd8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x48(%rbp),%rax
movslq -0xd4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0xd8(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x48(%rbp),%rax
movslq -0xd4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0xd8(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x48(%rbp),%rax
mov    -0xd4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0xd8(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0x50(%rbp),%rax
movslq -0xd4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xd8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x50(%rbp),%rax
movslq -0xd4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0xd8(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x50(%rbp),%rax
movslq -0xd4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0xd8(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0x50(%rbp),%rax
mov    -0xd4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xd8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x50(%rbp),%rax
mov    -0xd4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xd8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x50(%rbp),%rax
mov    -0xd4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0xd8(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x168(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x160(%rbp)
ff
movsd  0xc40(%rip),%xmm0
divsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x1a0(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm1
addsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
addsd  -0x80(%rbp),%xmm0
mulsd  -0xc8(%rbp),%xmm0
ff
movsd  -0xc8(%rbp),%xmm2
ff
movsd  -0x88(%rbp),%xmm1
ff
movsd  -0x68(%rbp),%xmm3
addsd  %xmm3,%xmm3
subsd  %xmm3,%xmm1
addsd  -0x90(%rbp),%xmm1
ff
mulsd  -0xd0(%rbp),%xmm1
ff
mulsd  -0xd0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1a8(%rbp)
ff
movsd  0xbcf(%rip),%xmm0
mulsd  -0xc8(%rbp),%xmm0
ff
movsd  %xmm0,-0x1b0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
addsd  -0x78(%rbp),%xmm0
mov    $0x2,%edi
call   1040 <pow@plt>
movsd  %xmm0,-0x1d8(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
addsd  -0x68(%rbp),%xmm0
mov    $0x2,%edi
call   1040 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0x1d8(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x1b8(%rbp)
ff
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x1c0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
addsd  -0x78(%rbp),%xmm0
call   1030 <abs@plt>
movsd  %xmm0,-0x1d0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
subsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x1c8(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
addsd  -0x68(%rbp),%xmm0
call   1030 <abs@plt>
movsd  -0x1d0(%rbp),%xmm5
ff
movsd  -0x1c8(%rbp),%xmm7
ff
movsd  -0x1c0(%rbp),%xmm3
ff
movsd  -0x1b8(%rbp),%xmm4
ff
movsd  -0x1b0(%rbp),%xmm1
ff
movsd  -0x1a8(%rbp),%xmm2
ff
movaps %xmm0,%xmm6
movsd  -0x1a0(%rbp),%xmm0
ff
movsd  -0x80(%rbp),%xmm8
subsd  -0x68(%rbp),%xmm8
mulsd  %xmm8,%xmm6
movq   %xmm6,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm6
mulsd  %xmm7,%xmm5
addsd  %xmm6,%xmm5
mulsd  %xmm5,%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x170(%rbp)
ff
movsd  0xa8d(%rip),%xmm0
mulsd  -0xd0(%rbp),%xmm0
ff
movsd  %xmm0,-0x178(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
addsd  -0xa8(%rbp),%xmm0
ff
movsd  -0x68(%rbp),%xmm2
addsd  -0x88(%rbp),%xmm2
ff
movsd  -0xa0(%rbp),%xmm1
ff
addsd  -0xc0(%rbp),%xmm1
ff
movsd  -0x90(%rbp),%xmm3
ff
addsd  -0x68(%rbp),%xmm3
mulsd  %xmm3,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x180(%rbp)
ff
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x188(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
addsd  -0xa8(%rbp),%xmm0
ff
call   1030 <abs@plt>
movsd  %xmm0,-0x198(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x190(%rbp)
ff
movsd  -0xa0(%rbp),%xmm0
ff
addsd  -0xc0(%rbp),%xmm0
ff
call   1030 <abs@plt>
movsd  -0x198(%rbp),%xmm6
ff
movsd  -0x190(%rbp),%xmm8
ff ff
movsd  -0x188(%rbp),%xmm4
ff
movsd  -0x180(%rbp),%xmm5
ff
movsd  -0x178(%rbp),%xmm2
ff
movsd  -0x170(%rbp),%xmm3
ff
movsd  -0x168(%rbp),%xmm1
ff
movaps %xmm0,%xmm7
movsd  -0x160(%rbp),%xmm0
ff
movsd  -0x90(%rbp),%xmm9
ff ff
subsd  -0x68(%rbp),%xmm9
mulsd  %xmm9,%xmm7
movq   %xmm7,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm7
mulsd  %xmm8,%xmm6
addsd  %xmm7,%xmm6
mulsd  %xmm6,%xmm4
addsd  %xmm5,%xmm4
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
addsd  -0x10(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x58(%rbp),%rax
movslq -0xd4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xd8(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  0x8d9(%rip),%xmm0
divsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
movsd  -0x70(%rbp),%xmm1
addsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
addsd  -0xb0(%rbp),%xmm0
ff
mulsd  -0xc8(%rbp),%xmm0
ff
movsd  -0xc8(%rbp),%xmm2
ff
movsd  -0x98(%rbp),%xmm1
ff
movsd  -0x70(%rbp),%xmm3
addsd  %xmm3,%xmm3
subsd  %xmm3,%xmm1
addsd  -0xa0(%rbp),%xmm1
ff
mulsd  -0xd0(%rbp),%xmm1
ff
mulsd  -0xd0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x128(%rbp)
ff
movsd  0x862(%rip),%xmm0
mulsd  -0xd0(%rbp),%xmm0
ff
movsd  %xmm0,-0x130(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
mov    $0x2,%edi
call   1040 <pow@plt>
movsd  %xmm0,-0x158(%rbp)
ff
movsd  -0xa0(%rbp),%xmm0
ff
addsd  -0x70(%rbp),%xmm0
mov    $0x2,%edi
call   1040 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0x158(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x138(%rbp)
ff
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x140(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
call   1030 <abs@plt>
movsd  %xmm0,-0x150(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
subsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x148(%rbp)
ff
movsd  -0xa0(%rbp),%xmm0
ff
addsd  -0x70(%rbp),%xmm0
call   1030 <abs@plt>
movsd  -0x150(%rbp),%xmm5
ff
movsd  -0x148(%rbp),%xmm7
ff
movsd  -0x140(%rbp),%xmm3
ff
movsd  -0x138(%rbp),%xmm4
ff
movsd  -0x130(%rbp),%xmm1
ff
movsd  -0x128(%rbp),%xmm2
ff
movaps %xmm0,%xmm6
movsd  -0x120(%rbp),%xmm0
ff
movsd  -0xa0(%rbp),%xmm8
ff ff
subsd  -0x70(%rbp),%xmm8
mulsd  %xmm8,%xmm6
movq   %xmm6,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm6
mulsd  %xmm7,%xmm5
addsd  %xmm6,%xmm5
mulsd  %xmm5,%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  0x70e(%rip),%xmm0
mulsd  -0xc8(%rbp),%xmm0
ff
movsd  %xmm0,-0xf8(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
movsd  -0x70(%rbp),%xmm2
addsd  -0xa8(%rbp),%xmm2
ff
movsd  -0x80(%rbp),%xmm1
addsd  -0xb8(%rbp),%xmm1
ff
movsd  -0xb0(%rbp),%xmm3
ff
addsd  -0x70(%rbp),%xmm3
mulsd  %xmm3,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
call   1030 <abs@plt>
movsd  %xmm0,-0x118(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
subsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0x110(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
addsd  -0xb8(%rbp),%xmm0
ff
call   1030 <abs@plt>
movsd  -0x118(%rbp),%xmm6
ff
movsd  -0x110(%rbp),%xmm8
ff ff
movsd  -0x108(%rbp),%xmm4
ff
movsd  -0x100(%rbp),%xmm5
ff
movsd  -0xf8(%rbp),%xmm2
ff
movsd  -0xf0(%rbp),%xmm3
ff
movsd  -0xe8(%rbp),%xmm1
ff
movaps %xmm0,%xmm7
movsd  -0xe0(%rbp),%xmm0
ff
movsd  -0xb0(%rbp),%xmm9
ff ff
subsd  -0x70(%rbp),%xmm9
mulsd  %xmm9,%xmm7
movq   %xmm7,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm7
mulsd  %xmm8,%xmm6
addsd  %xmm7,%xmm6
mulsd  %xmm6,%xmm4
addsd  %xmm5,%xmm4
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
addsd  -0x18(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x60(%rbp),%rax
movslq -0xd4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xd8(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xd8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xd8(%rbp)
jmp    11b6 <calculate_fg+0x96>
mov    -0x50(%rbp),%rax
movslq -0xd4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movsd  (%rax),%xmm0
mov    -0x60(%rbp),%rax
movslq -0xd4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movsd  %xmm0,(%rax)
mov    -0x50(%rbp),%rax
movslq -0xd4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x60(%rbp),%rax
movslq -0xd4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xd4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xd4(%rbp)
jmp    1198 <calculate_fg+0x78>
movl   $0x1,-0xd8(%rbp)
mov    -0xd8(%rbp),%eax
mov    -0x40(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1b84 <calculate_fg+0xa64>
mov    -0x48(%rbp),%rax
mov    (%rax),%rax
movslq -0xd8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x58(%rbp),%rax
mov    (%rax),%rax
movslq -0xd8(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xd8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x58(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xd8(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xd8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xd8(%rbp)
jmp    1b06 <calculate_fg+0x9e6>
add    $0x1e0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<calculate_fg>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x108,%rsp
mov    %rdx,%rbp
mov    %rdi,0x90(%rsp)
test   %edi,%edi
jle    1769 <calculate_fg+0x649>
mov    %r9,%r14
mov    %rcx,%rbx
movsd  0xead(%rip),%xmm8
movapd %xmm8,%xmm12
divsd  %xmm5,%xmm12
movapd %xmm8,%xmm13
divsd  %xmm6,%xmm13
divsd  %xmm0,%xmm8
movsd  0xe93(%rip),%xmm9
movapd %xmm12,%xmm10
mulsd  %xmm9,%xmm10
mulsd  %xmm13,%xmm9
movslq %esi,%rax
mov    %rax,0xd8(%rsp)
lea    0x1(%rax),%eax
mov    0x90(%rsp),%rcx
lea    0x1(%rcx),%ecx
mov    %rcx,0xd0(%rsp)
lea    -0x8(,%rax,8),%r12
ff
mov    $0x1,%r15d
movsd  0xe54(%rip),%xmm11
mov    %esi,0x8c(%rsp)
movsd  %xmm4,0xb8(%rsp)
movsd  %xmm3,0x40(%rsp)
movsd  %xmm2,0xf0(%rsp)
movsd  %xmm12,0x80(%rsp)
movsd  %xmm13,0x78(%rsp)
movsd  %xmm10,0xb0(%rsp)
movsd  %xmm1,0xe8(%rsp)
movsd  %xmm9,0xa8(%rsp)
mov    %r8,0x98(%rsp)
movsd  %xmm8,0xa0(%rsp)
jmp    1257 <calculate_fg+0x137>
nopl   (%rax)
mov    (%r14,%r15,8),%rax
lea    0x1(%r15),%rdi
mov    (%rbx,%r15,8),%rcx
movsd  (%rcx),%xmm0
movsd  %xmm0,(%rax)
mov    0xd8(%rsp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
movsd  %xmm0,(%rax,%rdx,8)
mov    %rdi,%r15
cmp    0xd0(%rsp),%rdi
je     1769 <calculate_fg+0x649>
test   %esi,%esi
jle    1220 <calculate_fg+0x100>
lea    0x1(%r15),%rax
mov    %rax,0xe0(%rsp)
mov    %eax,%eax
mov    %rax,0xf8(%rsp)
xor    %r13d,%r13d
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    -0x8(%rbp,%r15,8),%rax
mov    0x0(%rbp,%r15,8),%rcx
movsd  (%rcx,%r13,1),%xmm2
movsd  %xmm2,0x38(%rsp)
movsd  0x8(%rcx,%r13,1),%xmm3
mov    0xf8(%rsp),%rsi
mov    0x0(%rbp,%rsi,8),%rdx
movsd  0x8(%rdx,%r13,1),%xmm7
movsd  0x8(%rax,%r13,1),%xmm6
movsd  0x10(%rcx,%r13,1),%xmm4
movsd  %xmm4,0x30(%rsp)
mov    -0x8(%rbx,%r15,8),%rcx
mov    (%rbx,%r15,8),%rdx
movsd  (%rdx,%r13,1),%xmm1
movsd  %xmm1,0x58(%rsp)
movsd  0x8(%rdx,%r13,1),%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  0x10(%rdx,%r13,1),%xmm0
movsd  %xmm0,0x70(%rsp)
mov    (%rbx,%rsi,8),%rdx
movsd  0x8(%rdx,%r13,1),%xmm0
movsd  %xmm0,0x28(%rsp)
movsd  0x8(%rcx,%r13,1),%xmm0
movsd  %xmm0,0x60(%rsp)
movapd %xmm3,%xmm0
movapd %xmm3,%xmm5
movsd  %xmm3,0x8(%rsp)
mulsd  %xmm11,%xmm0
movapd %xmm0,%xmm3
addsd  %xmm7,%xmm3
movsd  %xmm7,0x50(%rsp)
addsd  %xmm6,%xmm3
movsd  %xmm6,0x68(%rsp)
mulsd  %xmm12,%xmm3
addsd  %xmm4,%xmm0
addsd  %xmm2,%xmm0
mulsd  %xmm13,%xmm0
mulsd  %xmm13,%xmm0
mulsd  %xmm12,%xmm3
addsd  %xmm0,%xmm3
movsd  %xmm3,0x18(%rsp)
addsd  %xmm7,%xmm5
movsd  %xmm5,0x48(%rsp)
movsd  (%rdx,%r13,1),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0xc0(%rsp)
movsd  0x10(%rax,%r13,1),%xmm0
addsd  %xmm6,%xmm0
movsd  %xmm0,0xc8(%rsp)
movapd %xmm5,%xmm0
mov    $0x2,%edi
call   1040 <pow@plt>
movsd  %xmm0,0x20(%rsp)
movsd  0x8(%rsp),%xmm0
addsd  0x68(%rsp),%xmm0
movsd  %xmm0,0x100(%rsp)
mov    $0x2,%edi
call   1040 <pow@plt>
movsd  0x20(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x20(%rsp)
movsd  0x48(%rsp),%xmm0
call   1030 <abs@plt>
movsd  %xmm0,0x48(%rsp)
movsd  0x8(%rsp),%xmm0
subsd  0x50(%rsp),%xmm0
movsd  %xmm0,0x50(%rsp)
movsd  0x100(%rsp),%xmm0
call   1030 <abs@plt>
movsd  0x8(%rsp),%xmm2
movsd  0x68(%rsp),%xmm1
subsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
movsd  0x50(%rsp),%xmm0
mulsd  0x48(%rsp),%xmm0
subsd  %xmm1,%xmm0
mulsd  0x40(%rsp),%xmm0
addsd  0x20(%rsp),%xmm0
mulsd  0xb0(%rsp),%xmm0
movapd %xmm0,%xmm1
movsd  0x18(%rsp),%xmm0
mulsd  0xa0(%rsp),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x18(%rsp)
movsd  0x10(%rsp),%xmm0
addsd  0x28(%rsp),%xmm0
movapd %xmm2,%xmm3
addsd  0x30(%rsp),%xmm3
movsd  %xmm3,0x68(%rsp)
movapd %xmm2,%xmm1
addsd  0x38(%rsp),%xmm1
mulsd  0xc0(%rsp),%xmm1
movapd %xmm0,%xmm2
mulsd  %xmm3,%xmm2
movsd  %xmm2,0x20(%rsp)
subsd  %xmm1,%xmm2
movsd  %xmm2,0x50(%rsp)
call   1030 <abs@plt>
movsd  %xmm0,0x48(%rsp)
movsd  0x8(%rsp),%xmm0
subsd  0x30(%rsp),%xmm0
movsd  %xmm0,0x30(%rsp)
movsd  0xc0(%rsp),%xmm0
call   1030 <abs@plt>
movsd  0x38(%rsp),%xmm1
movsd  0x8(%rsp),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
movsd  0x30(%rsp),%xmm0
mulsd  0x48(%rsp),%xmm0
subsd  %xmm1,%xmm0
mulsd  0x40(%rsp),%xmm0
addsd  0x50(%rsp),%xmm0
mulsd  0xa8(%rsp),%xmm0
addsd  0x18(%rsp),%xmm0
addsd  0xe8(%rsp),%xmm0
mulsd  0xb8(%rsp),%xmm0
addsd  %xmm2,%xmm0
mov    0x98(%rsp),%rax
mov    (%rax,%r15,8),%rax
movsd  %xmm0,0x8(%rax,%r13,1)
movsd  0x10(%rsp),%xmm1
movapd %xmm1,%xmm0
mulsd  0xac6(%rip),%xmm0
movapd %xmm0,%xmm3
addsd  0x28(%rsp),%xmm3
addsd  0x60(%rsp),%xmm3
mulsd  0x80(%rsp),%xmm3
movsd  0x70(%rsp),%xmm2
addsd  %xmm2,%xmm0
addsd  0x58(%rsp),%xmm0
mulsd  0x78(%rsp),%xmm0
mulsd  0x78(%rsp),%xmm0
mulsd  0x80(%rsp),%xmm3
addsd  %xmm0,%xmm3
movsd  %xmm3,0x8(%rsp)
movapd %xmm1,%xmm0
addsd  %xmm2,%xmm0
movsd  %xmm0,0x38(%rsp)
mov    $0x2,%edi
call   1040 <pow@plt>
movsd  %xmm0,0x18(%rsp)
movsd  0x10(%rsp),%xmm0
addsd  0x58(%rsp),%xmm0
movsd  %xmm0,0x30(%rsp)
mov    $0x2,%edi
call   1040 <pow@plt>
movsd  0x18(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x18(%rsp)
movsd  0x38(%rsp),%xmm0
call   1030 <abs@plt>
movsd  %xmm0,0x38(%rsp)
movsd  0x10(%rsp),%xmm0
subsd  0x70(%rsp),%xmm0
movsd  %xmm0,0x70(%rsp)
movsd  0x30(%rsp),%xmm0
call   1030 <abs@plt>
movsd  0x10(%rsp),%xmm1
movsd  0x58(%rsp),%xmm2
subsd  %xmm1,%xmm2
mulsd  %xmm0,%xmm2
movsd  0x70(%rsp),%xmm0
mulsd  0x38(%rsp),%xmm0
subsd  %xmm2,%xmm0
mulsd  0x40(%rsp),%xmm0
addsd  0x18(%rsp),%xmm0
mulsd  0xa8(%rsp),%xmm0
movapd %xmm0,%xmm2
movsd  0x8(%rsp),%xmm0
mulsd  0xa0(%rsp),%xmm0
addsd  %xmm2,%xmm0
movsd  %xmm0,0x8(%rsp)
movapd %xmm1,%xmm0
addsd  0x60(%rsp),%xmm0
mulsd  0xc8(%rsp),%xmm0
movsd  0x20(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x20(%rsp)
movsd  0x68(%rsp),%xmm0
call   1030 <abs@plt>
movsd  %xmm0,0x58(%rsp)
movsd  0x10(%rsp),%xmm0
subsd  0x28(%rsp),%xmm0
movsd  %xmm0,0x28(%rsp)
movsd  0xc8(%rsp),%xmm0
call   1030 <abs@plt>
movsd  0x953(%rip),%xmm11
movsd  0xb0(%rsp),%xmm10
movsd  0x78(%rsp),%xmm13
movsd  0x80(%rsp),%xmm12
movsd  0xf0(%rsp),%xmm2
movsd  0x40(%rsp),%xmm3
movsd  0xb8(%rsp),%xmm4
movsd  0x60(%rsp),%xmm1
movsd  0x10(%rsp),%xmm5
subsd  %xmm5,%xmm1
mulsd  %xmm0,%xmm1
movsd  0x28(%rsp),%xmm0
mulsd  0x58(%rsp),%xmm0
subsd  %xmm1,%xmm0
mulsd  %xmm3,%xmm0
addsd  0x20(%rsp),%xmm0
mulsd  %xmm10,%xmm0
addsd  0x8(%rsp),%xmm0
addsd  %xmm2,%xmm0
mulsd  %xmm4,%xmm0
addsd  %xmm5,%xmm0
mov    (%r14,%r15,8),%rax
movsd  %xmm0,0x8(%rax,%r13,1)
add    $0x8,%r13
cmp    %r13,%r12
jne    1280 <calculate_fg+0x160>
mov    0x8c(%rsp),%esi
mov    0x98(%rsp),%r8
mov    0xe0(%rsp),%rdi
jmp    1228 <calculate_fg+0x108>
test   %esi,%esi
jle    1973 <calculate_fg+0x853>
mov    0x0(%rbp),%rax
mov    (%r8),%rcx
movslq 0x90(%rsp),%rdi
mov    0x0(%rbp,%rdi,8),%rdx
mov    (%r8,%rdi,8),%r8
mov    %esi,%edi
lea    0x1(%rdi),%r10
mov    $0x1,%ebp
cmp    $0x4,%esi
jae    180c <calculate_fg+0x6ec>
mov    %rbp,%rdi
not    %rdi
test   $0x1,%r10b
jne    17be <calculate_fg+0x69e>
movsd  (%rax,%rbp,8),%xmm0
movsd  %xmm0,(%rcx,%rbp,8)
movsd  (%rdx,%rbp,8),%xmm0
movsd  %xmm0,(%r8,%rbp,8)
add    $0x1,%rbp
add    %r10,%rdi
je     1973 <calculate_fg+0x853>
nopw   0x0(%rax,%rax,1)
movsd  (%rax,%rbp,8),%xmm0
movsd  %xmm0,(%rcx,%rbp,8)
movsd  (%rdx,%rbp,8),%xmm0
movsd  %xmm0,(%r8,%rbp,8)
movsd  0x8(%rax,%rbp,8),%xmm0
movsd  %xmm0,0x8(%rcx,%rbp,8)
movsd  0x8(%rdx,%rbp,8),%xmm0
movsd  %xmm0,0x8(%r8,%rbp,8)
add    $0x2,%rbp
cmp    %rbp,%r10
jne    17d0 <calculate_fg+0x6b0>
jmp    1973 <calculate_fg+0x853>
mov    %rdi,0x28(%rsp)
lea    0x8(%rcx),%rdi
lea    (%rcx,%r10,8),%rbx
lea    0x8(%r8),%r11
lea    (%r8,%r10,8),%r13
lea    0x8(%rax),%r12
lea    (%rax,%r10,8),%r9
lea    0x8(%rdx),%r15
lea    (%rdx,%r10,8),%r14
cmp    %r13,%rdi
setb   0x10(%rsp)
cmp    %rbx,%r11
setb   0x8(%rsp)
cmp    %r9,%rdi
setb   %sil
cmp    %rbx,%r12
setb   0x40(%rsp)
cmp    %r14,%rdi
setb   %dil
cmp    %rbx,%r15
setb   0x60(%rsp)
cmp    %r9,%r11
setb   %r9b
cmp    %r13,%r12
setb   %r12b
cmp    %r14,%r11
setb   %r11b
cmp    %r13,%r15
setb   %r13b
mov    0x8(%rsp),%bl
test   %bl,0x10(%rsp)
jne    1799 <calculate_fg+0x679>
and    0x40(%rsp),%sil
jne    1799 <calculate_fg+0x679>
and    0x60(%rsp),%dil
jne    1799 <calculate_fg+0x679>
and    %r12b,%r9b
jne    1799 <calculate_fg+0x679>
and    %r13b,%r11b
jne    1799 <calculate_fg+0x679>
mov    0x28(%rsp),%rsi
mov    %esi,%ebp
and    $0xfffffffc,%ebp
lea    -0x4(%rbp),%rdi
mov    %rdi,%r9
shr    $0x2,%r9
add    $0x1,%r9
test   %rdi,%rdi
je     198e <calculate_fg+0x86e>
mov    %r9,%rdi
and    $0xfffffffffffffffe,%rdi
mov    $0x1,%ebx
nop
movups (%rax,%rbx,8),%xmm0
movups 0x10(%rax,%rbx,8),%xmm1
movups %xmm0,(%rcx,%rbx,8)
movups %xmm1,0x10(%rcx,%rbx,8)
movups (%rdx,%rbx,8),%xmm0
movups 0x10(%rdx,%rbx,8),%xmm1
movups %xmm0,(%r8,%rbx,8)
movups %xmm1,0x10(%r8,%rbx,8)
movups 0x20(%rax,%rbx,8),%xmm0
movups 0x30(%rax,%rbx,8),%xmm1
movups %xmm0,0x20(%rcx,%rbx,8)
movups %xmm1,0x30(%rcx,%rbx,8)
movupd 0x20(%rdx,%rbx,8),%xmm0
movupd 0x30(%rdx,%rbx,8),%xmm1
movupd %xmm0,0x20(%r8,%rbx,8)
movupd %xmm1,0x30(%r8,%rbx,8)
add    $0x8,%rbx
add    $0xfffffffffffffffe,%rdi
jne    18e0 <calculate_fg+0x7c0>
test   $0x1,%r9b
je     196e <calculate_fg+0x84e>
movups (%rax,%rbx,8),%xmm0
movups 0x10(%rax,%rbx,8),%xmm1
movups %xmm0,(%rcx,%rbx,8)
movups %xmm1,0x10(%rcx,%rbx,8)
movupd (%rdx,%rbx,8),%xmm0
movupd 0x10(%rdx,%rbx,8),%xmm1
movupd %xmm0,(%r8,%rbx,8)
movupd %xmm1,0x10(%r8,%rbx,8)
cmp    %rbp,%rsi
jne    1985 <calculate_fg+0x865>
add    $0x108,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
or     $0x1,%rbp
jmp    1799 <calculate_fg+0x679>
mov    $0x1,%ebx
test   $0x1,%r9b
jne    1944 <calculate_fg+0x824>
jmp    196e <calculate_fg+0x84e>

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

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