#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 atan(double);
scalar_t__ calloc(int, int);
int cos(double);
int /*<<< orphan*/ createGauss(int, double, double *);
int fabs(double);
scalar_t__ floor(int);
int /*<<< orphan*/ free(double *);
double g(double *, int, int, int, int, double *, int);
int pow(double, int);
int /*<<< orphan*/ printf(char *, ...);
int sin(double);
double sqrt(double);
void anisoGauss(double img1[], double img2[], int rows, int cols, int x, int y,
                double sig_u, double sig_v, double a, double angle,
                double *out1, double *out2) {
  double den, t_phi, phi, sig_x, sig_phi, den2;
  double *w, *w1, *g1, *g2;
  int i, k, xmfk, xpfk;
  int g1Dim, g2Dim;
  if (0 == angle) {
    t_phi = 0;
    phi = 0;
    sig_phi = sig_v;
    sig_x = sig_u;
  } else {
    den =
        pow(sig_v, 2) * pow(cos(angle), 2) + pow(sig_u, 2) * pow(sin(angle), 2);
    den2 = (pow(sig_u, 2) - pow(sig_v, 2)) * cos(angle) * sin(angle);
    if (fabs(den2) < 1e-6)
      t_phi = 0;
    else
      t_phi = den / den2;
    phi = atan(t_phi);
    if (fabs(phi) < 1e-6)
      sig_phi = sig_u;
    else
      sig_phi = fabs(1 / sin(phi) * sqrt(den));
    sig_x = fabs((sig_u * sig_v) / sqrt(den));
    sig_x = ((sig_x) < (50) ? (sig_x) : (50));
    sig_phi = ((sig_phi) < (50) ? (sig_phi) : (50));
    if (0) {
      printf("i: %d, j: %d\n", x, y);
      printf("sig_u: %2.4f, sig_v: %2.4f, sig_x: %2.4f\n", sig_u, sig_v, sig_x);
      printf("sig_phi: %2.4f, phi: %2.4f, t_phi: %2.4f\n", sig_phi, phi, t_phi);
    }
  }
  g1Dim = 3 * sig_x + 1;
  g1 = (double *)calloc(g1Dim, sizeof(double));
  createGauss(g1Dim, sig_x, g1);
  g2Dim = 3 * sig_phi + 1;
  g2 = (double *)calloc(g2Dim, sizeof(double));
  createGauss(g2Dim, sig_phi, g2);
  if (x - g1Dim < 0 || y - g2Dim < 0) {
    return;
  }
  w = (double *)calloc(g2Dim, sizeof(double));
  w1 = (double *)calloc(g2Dim, sizeof(double));
  for (i = 0; i < g2Dim; i++) {
    w[i] = g2[i] * a;
    w1[i] = g2[i] * (1 - a);
  }
  out1[0] = g2[0] * g(img1, x, y, cols, rows, g1, g1Dim);
  out2[0] = g2[0] * g(img2, x, y, cols, rows, g1, g1Dim);
  for (k = 1; k < g2Dim; k++) {
    if (fabs(t_phi) < 1e-6) {
      xpfk = x;
      xmfk = x;
      if (0)
        printf("Set interpolation factor to 0\n");
    } else {
      xpfk = (int)floor(x + (double)k / t_phi);
      xmfk = (int)floor(x - (double)k / t_phi);
    }
    out1[0] += w[k] * (g(img1, xmfk, y - k, cols, rows, g1, g1Dim) +
                       g(img1, xpfk, y + k, cols, rows, g1, g1Dim)) +
               w1[k] * (g(img1, xmfk, y - k, cols, rows, g1, g1Dim) +
                        g(img1, xpfk, y + k, cols, rows, g1, g1Dim));
    out2[0] += w[k] * (g(img2, xmfk, y - k, cols, rows, g1, g1Dim) +
                       g(img2, xpfk, y + k, cols, rows, g1, g1Dim)) +
               w1[k] * (g(img2, xmfk, y - k, cols, rows, g1, g1Dim) +
                        g(img2, xpfk, y + k, cols, rows, g1, g1Dim));
  }
  if (0)
    printf("out: %2.4f\n", out1[0]);
  if (0)
    printf("out: %2.4f\n", out2[0]);
  free(g1);
  free(g2);
  free(w);
  free(w1);
}


x86 assembly (clang, -O0):
<anisoGauss>:
push   %rbp
mov    %rsp,%rbp
sub    $0x150,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9d,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
movsd  %xmm1,-0x30(%rbp)
movsd  %xmm2,-0x38(%rbp)
movsd  %xmm3,-0x40(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x40(%rbp),%xmm0
jne    122a <anisoGauss+0x7a>
jp     122a <anisoGauss+0x7a>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
jmp    14a6 <anisoGauss+0x2f6>
movsd  -0x30(%rbp),%xmm0
mov    $0x2,%edi
call   10b0 <pow@plt>
mov    %eax,-0xc0(%rbp)
movsd  -0x40(%rbp),%xmm0
call   10a0 <cos@plt>
cvtsi2sd %eax,%xmm0
mov    $0x2,%edi
call   10b0 <pow@plt>
mov    %eax,%ecx
mov    -0xc0(%rbp),%eax
imul   %ecx,%eax
mov    %eax,-0xb8(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    $0x2,%edi
call   10b0 <pow@plt>
mov    %eax,-0xbc(%rbp)
movsd  -0x40(%rbp),%xmm0
call   1070 <sin@plt>
cvtsi2sd %eax,%xmm0
mov    $0x2,%edi
call   10b0 <pow@plt>
mov    -0xbc(%rbp),%ecx
mov    %eax,%edx
mov    -0xb8(%rbp),%eax
imul   %edx,%ecx
add    %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    $0x2,%edi
call   10b0 <pow@plt>
mov    %eax,-0xb4(%rbp)
movsd  -0x30(%rbp),%xmm0
mov    $0x2,%edi
call   10b0 <pow@plt>
mov    %eax,%ecx
mov    -0xb4(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0xb0(%rbp)
movsd  -0x40(%rbp),%xmm0
call   10a0 <cos@plt>
mov    %eax,%ecx
mov    -0xb0(%rbp),%eax
imul   %ecx,%eax
mov    %eax,-0xac(%rbp)
movsd  -0x40(%rbp),%xmm0
call   1070 <sin@plt>
mov    %eax,%ecx
mov    -0xac(%rbp),%eax
imul   %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
call   1060 <fabs@plt>
cvtsi2sd %eax,%xmm1
movsd  0xccc(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    134b <anisoGauss+0x19b>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    135a <anisoGauss+0x1aa>
movsd  -0x48(%rbp),%xmm0
divsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
call   1080 <atan@plt>
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
call   1060 <fabs@plt>
cvtsi2sd %eax,%xmm1
movsd  0xc81(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1398 <anisoGauss+0x1e8>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    13df <anisoGauss+0x22f>
movsd  -0x58(%rbp),%xmm0
call   1070 <sin@plt>
mov    %eax,%ecx
mov    $0x1,%eax
cltd
idiv   %ecx
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
call   10c0 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0xc8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
call   1060 <fabs@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x28(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
call   10c0 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0xd0(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
call   1060 <fabs@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  0xbe8(%rip),%xmm0
ucomisd -0x60(%rbp),%xmm0
jbe    143d <anisoGauss+0x28d>
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
jmp    1452 <anisoGauss+0x2a2>
movsd  0xbc3(%rip),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
jmp    1452 <anisoGauss+0x2a2>
movsd  -0xd8(%rbp),%xmm0
ff
movsd  %xmm0,-0x60(%rbp)
movsd  0xba1(%rip),%xmm0
ucomisd -0x68(%rbp),%xmm0
jbe    1484 <anisoGauss+0x2d4>
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
jmp    1499 <anisoGauss+0x2e9>
movsd  0xb7c(%rip),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
jmp    1499 <anisoGauss+0x2e9>
movsd  -0xe0(%rbp),%xmm0
ff
movsd  %xmm0,-0x68(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  0xb5d(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xb59(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0xa4(%rbp)
mov    -0xa4(%rbp),%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x88(%rbp)
mov    -0xa4(%rbp),%edi
movsd  -0x60(%rbp),%xmm0
mov    -0x88(%rbp),%rsi
call   1040 <createGauss@plt>
movsd  -0x68(%rbp),%xmm0
movsd  0xb08(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xb04(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0xa8(%rbp)
mov    -0xa8(%rbp),%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x90(%rbp)
mov    -0xa8(%rbp),%edi
movsd  -0x68(%rbp),%xmm0
mov    -0x90(%rbp),%rsi
call   1040 <createGauss@plt>
mov    -0x1c(%rbp),%eax
sub    -0xa4(%rbp),%eax
cmp    $0x0,%eax
jl     1574 <anisoGauss+0x3c4>
mov    -0x20(%rbp),%eax
sub    -0xa8(%rbp),%eax
cmp    $0x0,%eax
jge    1579 <anisoGauss+0x3c9>
jmp    1a4c <anisoGauss+0x89c>
mov    -0xa8(%rbp),%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x78(%rbp)
mov    -0xa8(%rbp),%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x80(%rbp)
movl   $0x0,-0x94(%rbp)
mov    -0x94(%rbp),%eax
cmp    -0xa8(%rbp),%eax
jge    162d <anisoGauss+0x47d>
mov    -0x90(%rbp),%rax
movslq -0x94(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x38(%rbp),%xmm0
mov    -0x78(%rbp),%rax
movslq -0x94(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x90(%rbp),%rax
movslq -0x94(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xa18(%rip),%xmm1
subsd  -0x38(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x80(%rbp),%rax
movslq -0x94(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    15ab <anisoGauss+0x3fb>
mov    -0x90(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
mov    -0x8(%rbp),%rdi
mov    -0x1c(%rbp),%esi
mov    -0x20(%rbp),%edx
mov    -0x18(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    -0x88(%rbp),%r9
mov    -0xa4(%rbp),%eax
mov    %eax,(%rsp)
call   1090 <g@plt>
movaps %xmm0,%xmm1
movsd  -0xf0(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
mov    0x10(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x90(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
mov    -0x10(%rbp),%rdi
mov    -0x1c(%rbp),%esi
mov    -0x20(%rbp),%edx
mov    -0x18(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    -0x88(%rbp),%r9
mov    -0xa4(%rbp),%eax
mov    %eax,(%rsp)
call   1090 <g@plt>
movaps %xmm0,%xmm1
movsd  -0xe8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
mov    0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x98(%rbp)
mov    -0x98(%rbp),%eax
cmp    -0xa8(%rbp),%eax
jge    1a22 <anisoGauss+0x872>
movsd  -0x50(%rbp),%xmm0
call   1060 <fabs@plt>
cvtsi2sd %eax,%xmm1
movsd  0x901(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1720 <anisoGauss+0x570>
mov    -0x1c(%rbp),%eax
mov    %eax,-0xa0(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x9c(%rbp)
jmp    176a <anisoGauss+0x5ba>
cvtsi2sdl -0x1c(%rbp),%xmm0
cvtsi2sdl -0x98(%rbp),%xmm1
ff
divsd  -0x50(%rbp),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   10d0 <floor@plt>
mov    %eax,-0xa0(%rbp)
cvtsi2sdl -0x1c(%rbp),%xmm0
cvtsi2sdl -0x98(%rbp),%xmm1
ff
divsd  -0x50(%rbp),%xmm1
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   10d0 <floor@plt>
mov    %eax,-0x9c(%rbp)
mov    -0x78(%rbp),%rax
movslq -0x98(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
mov    -0x8(%rbp),%rdi
mov    -0x9c(%rbp),%esi
mov    -0x20(%rbp),%edx
sub    -0x98(%rbp),%edx
mov    -0x18(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    -0x88(%rbp),%r9
mov    -0xa4(%rbp),%eax
mov    %eax,(%rsp)
call   1090 <g@plt>
movsd  %xmm0,-0x140(%rbp)
ff
mov    -0x8(%rbp),%rdi
mov    -0xa0(%rbp),%esi
mov    -0x20(%rbp),%edx
add    -0x98(%rbp),%edx
mov    -0x18(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    -0x88(%rbp),%r9
mov    -0xa4(%rbp),%eax
mov    %eax,(%rsp)
call   1090 <g@plt>
movaps %xmm0,%xmm1
movsd  -0x140(%rbp),%xmm0
ff
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x128(%rbp)
ff
mov    -0x80(%rbp),%rax
movslq -0x98(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
mov    -0x8(%rbp),%rdi
mov    -0x9c(%rbp),%esi
mov    -0x20(%rbp),%edx
sub    -0x98(%rbp),%edx
mov    -0x18(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    -0x88(%rbp),%r9
mov    -0xa4(%rbp),%eax
mov    %eax,(%rsp)
call   1090 <g@plt>
movsd  %xmm0,-0x138(%rbp)
ff
mov    -0x8(%rbp),%rdi
mov    -0xa0(%rbp),%esi
mov    -0x20(%rbp),%edx
add    -0x98(%rbp),%edx
mov    -0x18(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    -0x88(%rbp),%r9
mov    -0xa4(%rbp),%eax
mov    %eax,(%rsp)
call   1090 <g@plt>
movsd  -0x138(%rbp),%xmm3
ff
movsd  -0x130(%rbp),%xmm1
ff
movsd  -0x128(%rbp),%xmm2
ff
movaps %xmm0,%xmm4
movsd  -0x120(%rbp),%xmm0
ff
addsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x10(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x78(%rbp),%rax
movslq -0x98(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
mov    -0x10(%rbp),%rdi
mov    -0x9c(%rbp),%esi
mov    -0x20(%rbp),%edx
sub    -0x98(%rbp),%edx
mov    -0x18(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    -0x88(%rbp),%r9
mov    -0xa4(%rbp),%eax
mov    %eax,(%rsp)
call   1090 <g@plt>
movsd  %xmm0,-0x118(%rbp)
ff
mov    -0x10(%rbp),%rdi
mov    -0xa0(%rbp),%esi
mov    -0x20(%rbp),%edx
add    -0x98(%rbp),%edx
mov    -0x18(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    -0x88(%rbp),%r9
mov    -0xa4(%rbp),%eax
mov    %eax,(%rsp)
call   1090 <g@plt>
movaps %xmm0,%xmm1
movsd  -0x118(%rbp),%xmm0
ff
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
mov    -0x80(%rbp),%rax
movslq -0x98(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
mov    -0x10(%rbp),%rdi
mov    -0x9c(%rbp),%esi
mov    -0x20(%rbp),%edx
sub    -0x98(%rbp),%edx
mov    -0x18(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    -0x88(%rbp),%r9
mov    -0xa4(%rbp),%eax
mov    %eax,(%rsp)
call   1090 <g@plt>
movsd  %xmm0,-0x110(%rbp)
ff
mov    -0x10(%rbp),%rdi
mov    -0xa0(%rbp),%esi
mov    -0x20(%rbp),%edx
add    -0x98(%rbp),%edx
mov    -0x18(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    -0x88(%rbp),%r9
mov    -0xa4(%rbp),%eax
mov    %eax,(%rsp)
call   1090 <g@plt>
movsd  -0x110(%rbp),%xmm3
ff
movsd  -0x108(%rbp),%xmm1
ff
movsd  -0x100(%rbp),%xmm2
ff
movaps %xmm0,%xmm4
movsd  -0xf8(%rbp),%xmm0
ff
addsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x18(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x98(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x98(%rbp)
jmp    16d7 <anisoGauss+0x527>
mov    -0x88(%rbp),%rdi
call   1030 <free@plt>
mov    -0x90(%rbp),%rdi
call   1030 <free@plt>
mov    -0x78(%rbp),%rdi
call   1030 <free@plt>
mov    -0x80(%rbp),%rdi
call   1030 <free@plt>
add    $0x150,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<anisoGauss>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xc8,%rsp
movaps %xmm2,0xa0(%rsp)
movapd %xmm0,%xmm2
mov    %r9d,%r13d
mov    %r8d,0x2c(%rsp)
mov    %ecx,0x4c(%rsp)
mov    %edx,0x48(%rsp)
mov    %rsi,0x90(%rsp)
mov    %rdi,0x88(%rsp)
mov    0x108(%rsp),%rax
mov    %rax,0x80(%rsp)
mov    0x100(%rsp),%rax
mov    %rax,0x78(%rsp)
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm3
movsd  %xmm1,0x20(%rsp)
jne    122b <anisoGauss+0x7b>
jp     122b <anisoGauss+0x7b>
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x58(%rsp)
jmp    13fe <anisoGauss+0x24e>
movapd %xmm1,%xmm0
mov    $0x2,%edi
movsd  %xmm2,0x10(%rsp)
movsd  %xmm3,0x18(%rsp)
call   10b0 <pow@plt>
mov    %eax,%r14d
movsd  0x18(%rsp),%xmm0
call   10a0 <cos@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    $0x2,%edi
call   10b0 <pow@plt>
mov    %eax,%ebp
imul   %r14d,%ebp
movsd  0x10(%rsp),%xmm0
mov    $0x2,%edi
call   10b0 <pow@plt>
mov    %eax,%r14d
movsd  0x18(%rsp),%xmm0
call   1070 <sin@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    $0x2,%edi
call   10b0 <pow@plt>
imul   %r14d,%eax
add    %ebp,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x40(%rsp)
movsd  0x10(%rsp),%xmm0
mov    $0x2,%edi
call   10b0 <pow@plt>
mov    %eax,%ebp
movsd  0x20(%rsp),%xmm0
mov    $0x2,%edi
call   10b0 <pow@plt>
sub    %eax,%ebp
movsd  0x18(%rsp),%xmm0
call   10a0 <cos@plt>
mov    %eax,%r14d
movsd  0x18(%rsp),%xmm0
call   1070 <sin@plt>
imul   %r14d,%eax
imul   %ebp,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x18(%rsp)
call   1060 <fabs@plt>
movsd  0x40(%rsp),%xmm0
test   %eax,%eax
jle    1314 <anisoGauss+0x164>
divsd  0x18(%rsp),%xmm0
jmp    1318 <anisoGauss+0x168>
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x58(%rsp)
call   1080 <atan@plt>
movsd  %xmm0,0x18(%rsp)
call   1060 <fabs@plt>
movsd  0x10(%rsp),%xmm1
movapd %xmm1,%xmm3
test   %eax,%eax
jle    1397 <anisoGauss+0x1e7>
movsd  0x18(%rsp),%xmm0
call   1070 <sin@plt>
lea    0x1(%rax),%ecx
xor    %edx,%edx
cmp    $0x3,%ecx
cmovb  %eax,%edx
cvtsi2sd %edx,%xmm2
xorpd  %xmm0,%xmm0
movsd  0x40(%rsp),%xmm1
ucomisd %xmm0,%xmm1
movapd %xmm1,%xmm0
jb     1370 <anisoGauss+0x1c0>
sqrtsd %xmm0,%xmm0
jmp    1381 <anisoGauss+0x1d1>
movsd  %xmm2,0x18(%rsp)
call   10c0 <sqrt@plt>
movsd  0x18(%rsp),%xmm2
mulsd  %xmm2,%xmm0
call   1060 <fabs@plt>
xorps  %xmm3,%xmm3
cvtsi2sd %eax,%xmm3
movsd  0x10(%rsp),%xmm1
mulsd  0x20(%rsp),%xmm1
xorpd  %xmm0,%xmm0
movsd  0x40(%rsp),%xmm2
ucomisd %xmm0,%xmm2
movapd %xmm2,%xmm0
movsd  %xmm3,0x18(%rsp)
jb     13bd <anisoGauss+0x20d>
sqrtsd %xmm0,%xmm0
jmp    13ce <anisoGauss+0x21e>
movsd  %xmm1,0x10(%rsp)
call   10c0 <sqrt@plt>
movsd  0x10(%rsp),%xmm1
divsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
call   1060 <fabs@plt>
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
movsd  0xc16(%rip),%xmm0
minsd  %xmm0,%xmm2
movsd  0x18(%rsp),%xmm1
minsd  %xmm0,%xmm1
movsd  %xmm1,0x20(%rsp)
movsd  %xmm2,0x10(%rsp)
mulsd  0xbfc(%rip),%xmm2
movsd  0xbfc(%rip),%xmm1
addsd  %xmm1,%xmm2
cvttsd2si %xmm2,%r15d
mov    %r15d,%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %r15d,%edi
movsd  0x10(%rsp),%xmm0
mov    %rax,0x50(%rsp)
mov    %rax,%rsi
call   1040 <createGauss@plt>
movsd  0xbc0(%rip),%xmm0
mulsd  0x20(%rsp),%xmm0
movsd  0xbba(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%r12d
mov    %r12d,%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,%rbx
mov    %r12d,%edi
movsd  0x20(%rsp),%xmm0
mov    %rax,%rsi
call   1040 <createGauss@plt>
cmp    %r13d,%r12d
jg     156b <anisoGauss+0x3bb>
cmp    0x2c(%rsp),%r15d
jg     156b <anisoGauss+0x3bb>
mov    %r12d,%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,%rbp
mov    %r12d,%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %r12d,%edx
test   %r12d,%r12d
movapd 0xa0(%rsp),%xmm4
movsd  0xb48(%rip),%xmm5
jle    16e7 <anisoGauss+0x537>
subsd  %xmm4,%xmm5
cmp    $0x4,%r12d
jae    157d <anisoGauss+0x3cd>
xor    %edi,%edi
mov    %rdi,%rcx
not    %rcx
test   $0x1,%dl
je     150a <anisoGauss+0x35a>
movsd  (%rbx,%rdi,8),%xmm0
mulsd  %xmm4,%xmm0
movsd  %xmm0,0x0(%rbp,%rdi,8)
movsd  (%rbx,%rdi,8),%xmm0
mulsd  %xmm5,%xmm0
movsd  %xmm0,(%rax,%rdi,8)
or     $0x1,%rdi
add    %rdx,%rcx
je     16e7 <anisoGauss+0x537>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rbx,%rdi,8),%xmm0
mulsd  %xmm4,%xmm0
movsd  %xmm0,0x0(%rbp,%rdi,8)
movsd  (%rbx,%rdi,8),%xmm0
mulsd  %xmm5,%xmm0
movsd  %xmm0,(%rax,%rdi,8)
movsd  0x8(%rbx,%rdi,8),%xmm0
mulsd  %xmm4,%xmm0
movsd  %xmm0,0x8(%rbp,%rdi,8)
movsd  0x8(%rbx,%rdi,8),%xmm0
mulsd  %xmm5,%xmm0
movsd  %xmm0,0x8(%rax,%rdi,8)
add    $0x2,%rdi
cmp    %rdi,%rdx
jne    1520 <anisoGauss+0x370>
jmp    16e7 <anisoGauss+0x537>
add    $0xc8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    0x0(,%rdx,8),%rdi
add    %rbp,%rdi
lea    (%rax,%rdx,8),%rsi
lea    (%rbx,%rdx,8),%rcx
cmp    %rsi,%rbp
setb   %r8b
cmp    %rdi,%rax
setb   %r9b
cmp    %rcx,%rbp
setb   %r14b
cmp    %rdi,%rbx
setb   %r10b
cmp    %rcx,%rax
setb   %r11b
cmp    %rsi,%rbx
setb   %sil
xor    %edi,%edi
test   %r9b,%r8b
jne    14de <anisoGauss+0x32e>
and    %r10b,%r14b
jne    14de <anisoGauss+0x32e>
and    %sil,%r11b
jne    14de <anisoGauss+0x32e>
mov    %edx,%edi
and    $0xfffffffc,%edi
movapd %xmm4,%xmm1
unpcklpd %xmm4,%xmm1
movapd %xmm5,%xmm0
unpcklpd %xmm5,%xmm0
lea    -0x4(%rdi),%rcx
mov    %rcx,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rcx,%rcx
je     1a88 <anisoGauss+0x8d8>
mov    %r8,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %ecx,%ecx
nopl   (%rax)
movupd (%rbx,%rcx,8),%xmm2
movupd 0x10(%rbx,%rcx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,0x0(%rbp,%rcx,8)
movupd %xmm3,0x10(%rbp,%rcx,8)
movupd (%rbx,%rcx,8),%xmm2
movupd 0x10(%rbx,%rcx,8),%xmm3
mulpd  %xmm0,%xmm2
mulpd  %xmm0,%xmm3
movupd %xmm2,(%rax,%rcx,8)
movupd %xmm3,0x10(%rax,%rcx,8)
movupd 0x20(%rbx,%rcx,8),%xmm2
movupd 0x30(%rbx,%rcx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,0x20(%rbp,%rcx,8)
movupd %xmm3,0x30(%rbp,%rcx,8)
movupd 0x20(%rbx,%rcx,8),%xmm2
movupd 0x30(%rbx,%rcx,8),%xmm3
mulpd  %xmm0,%xmm2
mulpd  %xmm0,%xmm3
movupd %xmm2,0x20(%rax,%rcx,8)
movupd %xmm3,0x30(%rax,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rsi
jne    1610 <anisoGauss+0x460>
test   $0x1,%r8b
je     16de <anisoGauss+0x52e>
movupd (%rbx,%rcx,8),%xmm2
movupd 0x10(%rbx,%rcx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,0x0(%rbp,%rcx,8)
movupd %xmm3,0x10(%rbp,%rcx,8)
movupd (%rbx,%rcx,8),%xmm1
movupd 0x10(%rbx,%rcx,8),%xmm2
mulpd  %xmm0,%xmm1
mulpd  %xmm0,%xmm2
movupd %xmm1,(%rax,%rcx,8)
movupd %xmm2,0x10(%rax,%rcx,8)
cmp    %rdx,%rdi
jne    14de <anisoGauss+0x32e>
mov    %rdx,0xb8(%rsp)
mov    %rax,0x60(%rsp)
mov    %rbp,0x68(%rsp)
movsd  (%rbx),%xmm0
movsd  %xmm0,0x10(%rsp)
mov    %r15d,(%rsp)
mov    0x88(%rsp),%rdi
mov    0x2c(%rsp),%esi
mov    %r13d,%edx
mov    0x4c(%rsp),%ebp
mov    %ebp,%ecx
mov    0x48(%rsp),%r14d
mov    %r14d,%r8d
mov    %r12d,0x20(%rsp)
mov    0x50(%rsp),%r12
mov    %r12,%r9
call   1090 <g@plt>
mulsd  0x10(%rsp),%xmm0
mov    0x78(%rsp),%rax
movsd  %xmm0,(%rax)
mov    %rbx,0xb0(%rsp)
movsd  (%rbx),%xmm0
movsd  %xmm0,0x10(%rsp)
mov    %r15d,0x74(%rsp)
mov    %r15d,(%rsp)
mov    0x90(%rsp),%rdi
mov    0x2c(%rsp),%esi
mov    %r13,%rbx
mov    %r13d,%edx
mov    %ebp,%ecx
mov    %r14d,%r8d
mov    %r12,%r9
call   1090 <g@plt>
mulsd  0x10(%rsp),%xmm0
mov    0x80(%rsp),%rax
movsd  %xmm0,(%rax)
cmpl   $0x2,0x20(%rsp)
jl     1a4c <anisoGauss+0x89c>
xorps  %xmm0,%xmm0
cvtsi2sdl 0x2c(%rsp),%xmm0
movsd  %xmm0,0x98(%rsp)
mov    %rbx,%r12
lea    0x1(%rbx),%eax
mov    %eax,0x20(%rsp)
add    $0xffffffff,%r12d
mov    $0x1,%ebp
mov    0x48(%rsp),%r15d
jmp    19e9 <anisoGauss+0x839>
nopl   0x0(%rax,%rax,1)
mov    0x68(%rsp),%rax
movsd  (%rax,%rbp,8),%xmm0
movsd  %xmm0,0x38(%rsp)
mov    0x74(%rsp),%r13d
mov    %r13d,(%rsp)
mov    0x88(%rsp),%r14
mov    %r14,%rdi
mov    %rdx,0x18(%rsp)
mov    %edx,%esi
mov    %r12,%rdx
mov    %r12,0x40(%rsp)
mov    %rbp,0x10(%rsp)
mov    %rbx,%rbp
mov    0x4c(%rsp),%ebx
mov    %ebx,%ecx
mov    %r15d,%r8d
mov    0x50(%rsp),%r12
mov    %r12,%r9
call   1090 <g@plt>
movsd  %xmm0,0x30(%rsp)
mov    %r13d,(%rsp)
mov    %r14,%rdi
mov    %ebp,%esi
mov    %rbp,0xa0(%rsp)
mov    0x20(%rsp),%edx
mov    %ebx,%ecx
mov    %r15d,%r8d
mov    %r12,%r9
call   1090 <g@plt>
addsd  0x30(%rsp),%xmm0
mulsd  0x38(%rsp),%xmm0
movsd  %xmm0,0x38(%rsp)
mov    0x10(%rsp),%rax
mov    0x60(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm0
movsd  %xmm0,0x30(%rsp)
mov    %r13d,(%rsp)
mov    %r14,%rdi
mov    0x18(%rsp),%rsi
mov    0x40(%rsp),%rdx
mov    %ebx,%ecx
mov    %r15d,%r8d
mov    %r12,%r9
call   1090 <g@plt>
movsd  %xmm0,0xc0(%rsp)
mov    %r13d,(%rsp)
mov    %r14,%rdi
mov    %ebp,%esi
mov    0x20(%rsp),%edx
mov    %ebx,%ecx
mov    %r15d,%r8d
mov    %r12,%r9
call   1090 <g@plt>
addsd  0xc0(%rsp),%xmm0
mulsd  0x30(%rsp),%xmm0
addsd  0x38(%rsp),%xmm0
mov    0x78(%rsp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    0x10(%rsp),%rax
mov    0x68(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm0
movsd  %xmm0,0x38(%rsp)
mov    %r13d,(%rsp)
mov    0x90(%rsp),%rbp
mov    %rbp,%rdi
mov    0x18(%rsp),%rsi
mov    0x40(%rsp),%r14
mov    %r14d,%edx
mov    %ebx,%ecx
mov    %r15d,%r8d
mov    %r12,%r9
call   1090 <g@plt>
movsd  %xmm0,0x30(%rsp)
mov    %r13d,(%rsp)
mov    %rbp,%rdi
mov    0xa0(%rsp),%rsi
mov    0x20(%rsp),%edx
mov    %ebx,%ecx
mov    %r15d,%r8d
mov    %r12,%r9
call   1090 <g@plt>
addsd  0x30(%rsp),%xmm0
mulsd  0x38(%rsp),%xmm0
movsd  %xmm0,0x38(%rsp)
mov    0x10(%rsp),%rax
mov    0x60(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm0
movsd  %xmm0,0x30(%rsp)
mov    %r13d,(%rsp)
mov    %rbp,%rdi
mov    0x18(%rsp),%rsi
mov    %r14d,%edx
mov    %ebx,%ecx
mov    %r15d,%r8d
mov    %r12,%r9
call   1090 <g@plt>
movsd  %xmm0,0x18(%rsp)
mov    %r13d,(%rsp)
mov    %rbp,%rdi
mov    0xa0(%rsp),%rsi
mov    0x20(%rsp),%r13d
mov    %r13d,%edx
mov    %ebx,%ecx
mov    0x10(%rsp),%rbp
mov    %r15d,%r8d
mov    %r12,%r9
call   1090 <g@plt>
addsd  0x18(%rsp),%xmm0
mulsd  0x30(%rsp),%xmm0
addsd  0x38(%rsp),%xmm0
mov    0x80(%rsp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
add    $0x1,%rbp
add    $0x1,%r13d
mov    %r13d,0x20(%rsp)
mov    %r14,%r12
add    $0xffffffff,%r12d
cmp    %rbp,0xb8(%rsp)
je     1a4c <anisoGauss+0x89c>
movsd  0x58(%rsp),%xmm0
call   1060 <fabs@plt>
mov    0x2c(%rsp),%ecx
mov    %ecx,%edx
mov    %ecx,%ebx
test   %eax,%eax
jle    17d0 <anisoGauss+0x620>
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
divsd  0x58(%rsp),%xmm0
movsd  %xmm0,0x10(%rsp)
addsd  0x98(%rsp),%xmm0
cvttsd2si %xmm0,%edi
call   10d0 <floor@plt>
mov    %rax,%rbx
movsd  0x98(%rsp),%xmm0
subsd  0x10(%rsp),%xmm0
cvttsd2si %xmm0,%edi
call   10d0 <floor@plt>
mov    %rax,%rdx
jmp    17d0 <anisoGauss+0x620>
mov    0x50(%rsp),%rdi
call   1030 <free@plt>
mov    0xb0(%rsp),%rdi
call   1030 <free@plt>
mov    0x68(%rsp),%rdi
call   1030 <free@plt>
mov    0x60(%rsp),%rdi
add    $0xc8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>
xor    %ecx,%ecx
test   $0x1,%r8b
jne    16a1 <anisoGauss+0x4f1>
jmp    16de <anisoGauss+0x52e>

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

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