#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 */
typedef size_t m;
typedef int const lmax;

/* Variables and functions */
double M_PI;
double *calloc(int const, int);
double cos(double const);
int /*<<< orphan*/ free(double *);
double sin(double const);
double sqrt(double);
double *realSphericalHarmonics(const int lmax, const double theta,
                               const double phi, double *fac, double *dfac) {
  int l, m;
  double x;
  double *P;
  double *s, *c;
  double norm0, norm;
  double *rsh;
  s = calloc(lmax + 2, sizeof(double));
  c = calloc(lmax + 2, sizeof(double));
  rsh = calloc(((lmax) + 1) * ((lmax) + 1), sizeof(double));
  P = calloc(((lmax) + 1) * ((lmax) + 1), sizeof(double));
  x = cos(theta);
  if (1.0 == x) {
    for (l = 0; l <= lmax; l++)
      P[(((l)) * ((l)) + ((0) - (-(l))))] = 1.0;
  } else if (-1.0 == x) {
    P[(((0)) * ((0)) + ((0) - (-(0))))] = 1.0;
    for (l = 1; l <= lmax; l++)
      P[(((l)) * ((l)) + ((0) - (-(l))))] =
          -P[(((l - 1)) * ((l - 1)) + ((0) - (-(l - 1))))];
  } else {
    s[1] = sqrt(1.0 - x * x);
    for (l = 2; l <= lmax; l++)
      s[l] = s[l - 1] * s[1];
    for (l = 0; l <= lmax; l++) {
      m = l;
      if (0 == m) {
        P[(((l)) * ((l)) + ((0) - (-(l))))] = 1.0;
      } else {
        P[(((l)) * ((l)) + ((m) - (-(l))))] = s[m] * dfac[2 * m - 1];
        m = l - 1;
        P[(((l)) * ((l)) + ((m) - (-(l))))] =
            x * (2 * m + 1) * P[(((l - 1)) * ((l - 1)) + ((m) - (-(l - 1))))];
        if (l > 1) {
          for (m = 0; m <= l - 2; m++) {
            P[(((l)) * ((l)) + ((m) - (-(l))))] =
                (x * (2 * l - 1) *
                     P[(((l - 1)) * ((l - 1)) + ((m) - (-(l - 1))))] -
                 (l + m - 1) *
                     P[(((l - 2)) * ((l - 2)) + ((m) - (-(l - 2))))]) /
                (l - m);
          }
        }
      }
    }
  }
  for (l = 0; l <= lmax; l++) {
    norm0 = sqrt((2.0 * l + 1.0) / (2.0 * M_PI));
    P[(((l)) * ((l)) + ((0) - (-(l))))] =
        norm0 * P[(((l)) * ((l)) + ((0) - (-(l))))];
    for (m = 1; m <= l; m++) {
      norm = sqrt(fac[l - m] / fac[l + m]) * norm0;
      P[(((l)) * ((l)) + ((m) - (-(l))))] =
          norm * P[(((l)) * ((l)) + ((m) - (-(l))))];
    }
  }
  if (lmax > 0) {
    if (0.0 == phi) {
      for (m = 0; m <= lmax; m++) {
        s[m] = 0.0;
        c[m] = 1.0;
      }
    } else {
      s[1] = sin(phi);
      c[1] = cos(phi);
      for (m = 2; m <= lmax; m++) {
        s[m] = s[1] * c[m - 1] + c[1] * s[m - 1];
        c[m] = c[1] * c[m - 1] - s[1] * s[m - 1];
      }
    }
  }
  for (l = 0; l <= lmax; l++) {
    rsh[(((l)) * ((l)) + ((0) - (-(l))))] =
        P[(((l)) * ((l)) + ((0) - (-(l))))] / sqrt(2.0);
    for (m = 1; m <= l; m++) {
      rsh[(((l)) * ((l)) + ((-m) - (-(l))))] =
          P[(((l)) * ((l)) + ((m) - (-(l))))] * s[m];
      rsh[(((l)) * ((l)) + ((+m) - (-(l))))] =
          P[(((l)) * ((l)) + ((m) - (-(l))))] * c[m];
    }
  }
  free(P);
  free(s);
  free(c);
  return rsh;
}


x86 assembly (clang, -O0):
<realSphericalHarmonics>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
mov    -0x4(%rbp),%edi
add    $0x2,%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
mov    %rax,-0x48(%rbp)
mov    -0x4(%rbp),%edi
add    $0x2,%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
mov    %rax,-0x50(%rbp)
mov    -0x4(%rbp),%edi
add    $0x1,%edi
mov    -0x4(%rbp),%eax
add    $0x1,%eax
imul   %eax,%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
mov    %rax,-0x68(%rbp)
mov    -0x4(%rbp),%edi
add    $0x1,%edi
mov    -0x4(%rbp),%eax
add    $0x1,%eax
imul   %eax,%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
mov    %rax,-0x40(%rbp)
movsd  -0x10(%rbp),%xmm0
call   1060 <cos@plt>
movsd  %xmm0,-0x38(%rbp)
movsd  0xe1a(%rip),%xmm0
ucomisd -0x38(%rbp),%xmm0
jne    1243 <realSphericalHarmonics+0xf3>
jp     1243 <realSphericalHarmonics+0xf3>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     123e <realSphericalHarmonics+0xee>
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
xor    %edx,%edx
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  0xdd5(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11fe <realSphericalHarmonics+0xae>
jmp    1592 <realSphericalHarmonics+0x442>
movsd  0xdbd(%rip),%xmm0
ucomisd -0x38(%rbp),%xmm0
jne    12f2 <realSphericalHarmonics+0x1a2>
jp     12f2 <realSphericalHarmonics+0x1a2>
mov    -0x40(%rbp),%rax
movsd  0xd98(%rip),%xmm0
movsd  %xmm0,(%rax)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     12ed <realSphericalHarmonics+0x19d>
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x2c(%rbp),%edx
sub    $0x1,%edx
imul   %edx,%ecx
mov    -0x2c(%rbp),%edx
sub    $0x1,%edx
xor    %esi,%esi
sub    %edx,%esi
xor    %edx,%edx
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
xor    %edx,%edx
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1273 <realSphericalHarmonics+0x123>
jmp    158d <realSphericalHarmonics+0x43d>
movsd  -0x38(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xce1(%rip),%xmm1
addsd  %xmm1,%xmm0
call   1070 <sqrt@plt>
mov    -0x48(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movl   $0x2,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     137a <realSphericalHarmonics+0x22a>
mov    -0x48(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
mulsd  0x8(%rax),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1338 <realSphericalHarmonics+0x1e8>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1588 <realSphericalHarmonics+0x438>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
xor    %eax,%eax
cmp    -0x30(%rbp),%eax
jne    13c9 <realSphericalHarmonics+0x279>
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
xor    %edx,%edx
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  0xc41(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1575 <realSphericalHarmonics+0x425>
mov    -0x48(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%ecx
shl    $0x1,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x30(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x2c(%rbp),%edx
sub    $0x1,%edx
imul   %edx,%ecx
mov    -0x30(%rbp),%edx
mov    -0x2c(%rbp),%edi
sub    $0x1,%edi
xor    %esi,%esi
sub    %edi,%esi
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
cmpl   $0x1,-0x2c(%rbp)
jle    1570 <realSphericalHarmonics+0x420>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jg     156b <realSphericalHarmonics+0x41b>
movsd  -0x38(%rbp),%xmm0
mov    -0x2c(%rbp),%eax
shl    $0x1,%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x2c(%rbp),%edx
sub    $0x1,%edx
imul   %edx,%ecx
mov    -0x30(%rbp),%edx
mov    -0x2c(%rbp),%edi
sub    $0x1,%edi
xor    %esi,%esi
sub    %edi,%esi
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x2c(%rbp),%eax
add    -0x30(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
mov    -0x2c(%rbp),%edx
sub    $0x2,%edx
imul   %edx,%ecx
mov    -0x30(%rbp),%edx
mov    -0x2c(%rbp),%edi
sub    $0x2,%edi
xor    %esi,%esi
sub    %edi,%esi
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x2c(%rbp),%eax
sub    -0x30(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1485 <realSphericalHarmonics+0x335>
jmp    1570 <realSphericalHarmonics+0x420>
jmp    1575 <realSphericalHarmonics+0x425>
jmp    157a <realSphericalHarmonics+0x42a>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1381 <realSphericalHarmonics+0x231>
jmp    158d <realSphericalHarmonics+0x43d>
jmp    1592 <realSphericalHarmonics+0x442>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     16c6 <realSphericalHarmonics+0x576>
cvtsi2sdl -0x2c(%rbp),%xmm0
addsd  %xmm0,%xmm0
movsd  0xa4a(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  0xa4e(%rip),%xmm1
mov    0x2a1f(%rip),%rax
mulsd  (%rax),%xmm1
divsd  %xmm1,%xmm0
call   1070 <sqrt@plt>
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
xor    %edx,%edx
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
xor    %edx,%edx
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     16b3 <realSphericalHarmonics+0x563>
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
call   1070 <sqrt@plt>
mulsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1623 <realSphericalHarmonics+0x4d3>
jmp    16b8 <realSphericalHarmonics+0x568>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1599 <realSphericalHarmonics+0x449>
cmpl   $0x0,-0x4(%rbp)
jle    182d <realSphericalHarmonics+0x6dd>
xorps  %xmm0,%xmm0
ucomisd -0x18(%rbp),%xmm0
jne    172f <realSphericalHarmonics+0x5df>
jp     172f <realSphericalHarmonics+0x5df>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     172a <realSphericalHarmonics+0x5da>
mov    -0x48(%rbp),%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  0x8e9(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    16eb <realSphericalHarmonics+0x59b>
jmp    1828 <realSphericalHarmonics+0x6d8>
movsd  -0x18(%rbp),%xmm0
call   1050 <sin@plt>
mov    -0x48(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0x18(%rbp),%xmm0
call   1060 <cos@plt>
mov    -0x50(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movl   $0x2,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1823 <realSphericalHarmonics+0x6d3>
mov    -0x48(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x50(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x50(%rbp),%rax
movsd  0x8(%rax),%xmm1
mov    -0x48(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x48(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x50(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x48(%rbp),%rax
movsd  0x8(%rax),%xmm1
mov    -0x48(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x50(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    175c <realSphericalHarmonics+0x60c>
jmp    1828 <realSphericalHarmonics+0x6d8>
jmp    182d <realSphericalHarmonics+0x6dd>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1966 <realSphericalHarmonics+0x816>
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
xor    %edx,%edx
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  0x7a5(%rip),%xmm0
call   1070 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x70(%rbp),%xmm0
divsd  %xmm1,%xmm0
mov    -0x68(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
xor    %edx,%edx
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     1953 <realSphericalHarmonics+0x803>
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x30(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x68(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
xor    %edx,%edx
sub    -0x30(%rbp),%edx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x30(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x68(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
xor    %esi,%esi
sub    -0x2c(%rbp),%esi
sub    %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    18a1 <realSphericalHarmonics+0x751>
jmp    1958 <realSphericalHarmonics+0x808>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1834 <realSphericalHarmonics+0x6e4>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
mov    -0x50(%rbp),%rdi
call   1030 <free@plt>
mov    -0x68(%rbp),%rax
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<realSphericalHarmonics>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %rdx,%r12
mov    %rsi,0x58(%rsp)
movsd  %xmm1,0x38(%rsp)
movsd  %xmm0,0x18(%rsp)
mov    %edi,%ebp
lea    0x2(%rbp),%ebx
mov    %ebx,%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,%r13
mov    %ebx,%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,%r14
lea    0x1(%rbp),%ebx
mov    %ebx,0xc(%rsp)
imul   %ebx,%ebx
mov    %ebx,%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,%r15
mov    %ebx,%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,(%rsp)
movsd  0x18(%rsp),%xmm0
call   1070 <cos@plt>
ucomisd 0xe1f(%rip),%xmm0
mov    %r14,0x48(%rsp)
mov    %r15,0x40(%rsp)
mov    %rbp,0x28(%rsp)
jne    122d <realSphericalHarmonics+0xcd>
jp     122d <realSphericalHarmonics+0xcd>
test   %ebp,%ebp
mov    (%rsp),%rdi
js     1d73 <realSphericalHarmonics+0xc13>
mov    0xc(%rsp),%r9d
lea    -0x1(%r9),%rax
mov    %r9d,%r10d
and    $0x3,%r10d
cmp    $0x3,%rax
jae    1319 <realSphericalHarmonics+0x1b9>
mov    %rbp,%rbx
xor    %ecx,%ecx
test   %r10,%r10
jne    13e3 <realSphericalHarmonics+0x283>
jmp    141d <realSphericalHarmonics+0x2bd>
movapd %xmm0,%xmm2
ucomisd 0xdcf(%rip),%xmm0
jne    15c4 <realSphericalHarmonics+0x464>
jp     15c4 <realSphericalHarmonics+0x464>
movabs $0x3ff0000000000000,%rax
mov    (%rsp),%rdi
mov    %rax,(%rdi)
mov    0x28(%rsp),%rbx
test   %ebx,%ebx
jle    141d <realSphericalHarmonics+0x2bd>
mov    0xc(%rsp),%edx
mov    %edx,%ecx
add    $0xffffffffffffffff,%rcx
mov    $0x1,%eax
cmp    $0x2,%edx
je     12e9 <realSphericalHarmonics+0x189>
mov    %rcx,%rdx
and    $0xfffffffffffffffe,%rdx
xor    %eax,%eax
movapd 0xdb8(%rip),%xmm0
nopl   0x0(%rax,%rax,1)
lea    0x1(%rax),%esi
mov    %eax,%edi
imul   %esi,%edi
movslq %edi,%rdi
mov    (%rsp),%rbp
movsd  0x0(%rbp,%rdi,8),%xmm1
xorpd  %xmm0,%xmm1
mov    %esi,%edi
imul   %edi,%edi
add    %eax,%edi
add    $0x1,%edi
mov    (%rsp),%rbp
movlpd %xmm1,0x0(%rbp,%rdi,8)
mov    (%rsp),%rdi
add    $0x2,%rax
imul   %eax,%esi
movslq %esi,%rsi
movsd  (%rdi,%rsi,8),%xmm1
xorpd  %xmm0,%xmm1
mov    %eax,%esi
imul   %esi,%esi
add    %eax,%esi
movlpd %xmm1,(%rdi,%rsi,8)
cmp    %rax,%rdx
jne    1290 <realSphericalHarmonics+0x130>
add    $0x1,%rax
test   $0x1,%cl
je     141d <realSphericalHarmonics+0x2bd>
lea    -0x1(%rax),%ecx
imul   %eax,%ecx
movslq %ecx,%rcx
movsd  (%rdi,%rcx,8),%xmm0
xorpd  0xd38(%rip),%xmm0
mov    %eax,%ecx
imul   %ecx,%ecx
add    %eax,%ecx
movlpd %xmm0,(%rdi,%rcx,8)
jmp    141d <realSphericalHarmonics+0x2bd>
mov    %r13,0x20(%rsp)
lea    0x10(%rdi),%r8
and    $0xfffffffc,%r9d
neg    %r9
mov    $0x1,%ecx
mov    $0x20,%r11d
mov    $0x4,%r15d
mov    $0x18,%r14d
mov    $0x1,%r12d
mov    $0x28,%eax
mov    %rdi,%rbx
mov    $0x9,%edi
mov    $0x10,%esi
xor    %edx,%edx
movabs $0x3ff0000000000000,%r13
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %edx,%ebp
and    $0xfffffff0,%ebp
lea    (%rbx,%rbp,8),%rbp
mov    %r13,-0x8(%rbp,%rcx,8)
mov    %r12d,%ebp
add    %rcx,%rbp
mov    %r13,(%rbx,%rbp,8)
mov    %r15d,%ebp
and    $0xfffffffc,%ebp
lea    (%r8,%rbp,8),%rbp
mov    %r13,-0x8(%rbp,%rcx,8)
mov    %edi,%ebp
add    %rcx,%rbp
mov    %r13,0x10(%rbx,%rbp,8)
add    %r11,%r15
add    $0x20,%r11
add    %r14d,%r12d
add    $0x20,%r14d
add    %eax,%edi
add    $0x20,%eax
add    %rsi,%rdx
add    $0x20,%rsi
lea    (%r9,%rcx,1),%rbp
add    $0x4,%rbp
add    $0x4,%rcx
cmp    $0x1,%rbp
jne    1370 <realSphericalHarmonics+0x210>
add    $0xffffffffffffffff,%rcx
mov    0x20(%rsp),%r13
mov    %rbx,%rdi
mov    0x28(%rsp),%rbx
test   %r10,%r10
je     141d <realSphericalHarmonics+0x2bd>
mov    %ecx,%eax
imul   %ecx,%eax
lea    (%rcx,%rcx,1),%edx
add    $0x1,%edx
movabs $0x3ff0000000000000,%rsi
mov    (%rsp),%rbp
nopl   0x0(%rax)
mov    %eax,%edi
add    %rcx,%rdi
mov    %rsi,0x0(%rbp,%rdi,8)
mov    (%rsp),%rdi
add    $0x1,%rcx
add    %edx,%eax
add    $0x2,%edx
add    $0xffffffffffffffff,%r10
jne    1400 <realSphericalHarmonics+0x2a0>
test   %ebx,%ebx
mov    0x48(%rsp),%r14
mov    0x40(%rsp),%r15
js     1d73 <realSphericalHarmonics+0xc13>
mov    %r13,0x20(%rsp)
mov    0xc(%rsp),%eax
mov    %rax,0x10(%rsp)
mov    0x58(%rsp),%rbx
lea    0x8(%rbx),%r13
add    $0xfffffffffffffff8,%rbx
mov    $0x2,%eax
mov    $0x1,%r14d
xor    %r15d,%r15d
xorpd  %xmm2,%xmm2
jmp    1486 <realSphericalHarmonics+0x326>
xchg   %ax,%ax
add    $0x1,%r15
add    $0x8,%r13
mov    0x30(%rsp),%rax
add    %rax,%r14
add    $0x2,%rax
add    $0x8,%rbx
mov    0x10(%rsp),%rbp
cmp    %rbp,%r15
je     155f <realSphericalHarmonics+0x3ff>
mov    %rax,0x30(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %r15d,%xmm0
addsd  %xmm0,%xmm0
addsd  0xb61(%rip),%xmm0
mov    0x2b42(%rip),%rax
movsd  (%rax),%xmm1
addsd  %xmm1,%xmm1
divsd  %xmm1,%xmm0
ucomisd %xmm2,%xmm0
jb     14d0 <realSphericalHarmonics+0x370>
xorps  %xmm1,%xmm1
sqrtsd %xmm0,%xmm1
jmp    14e1 <realSphericalHarmonics+0x381>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
call   1080 <sqrt@plt>
xorpd  %xmm2,%xmm2
mov    (%rsp),%rdi
movapd %xmm0,%xmm1
mov    %r15d,%eax
imul   %eax,%eax
add    %r15d,%eax
movsd  (%rdi,%rax,8),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,(%rdi,%rax,8)
test   %r15,%r15
je     1460 <realSphericalHarmonics+0x300>
mov    %rbx,%r12
xor    %ebp,%ebp
movsd  %xmm1,0x18(%rsp)
jmp    1546 <realSphericalHarmonics+0x3e6>
xchg   %ax,%ax
call   1080 <sqrt@plt>
movsd  0x18(%rsp),%xmm1
xorpd  %xmm2,%xmm2
mov    (%rsp),%rdi
mulsd  %xmm1,%xmm0
lea    (%r14,%rbp,1),%eax
mulsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,(%rdi,%rax,8)
add    $0x1,%rbp
add    $0xfffffffffffffff8,%r12
cmp    %rbp,%r15
je     1460 <realSphericalHarmonics+0x300>
movsd  (%r12),%xmm0
divsd  0x0(%r13,%rbp,8),%xmm0
ucomisd %xmm2,%xmm0
jb     1510 <realSphericalHarmonics+0x3b0>
sqrtsd %xmm0,%xmm0
jmp    1523 <realSphericalHarmonics+0x3c3>
mov    0x28(%rsp),%rbx
test   %ebx,%ebx
mov    0x20(%rsp),%r13
mov    0x48(%rsp),%r14
mov    0x40(%rsp),%r15
jle    1c2d <realSphericalHarmonics+0xacd>
xorpd  %xmm0,%xmm0
movsd  0x38(%rsp),%xmm1
ucomisd %xmm0,%xmm1
jne    15f4 <realSphericalHarmonics+0x494>
jp     15f4 <realSphericalHarmonics+0x494>
test   %ebx,%ebx
js     1d73 <realSphericalHarmonics+0xc13>
mov    %ebx,%eax
lea    0x8(,%rax,8),%rdx
xor    %ebp,%ebp
mov    %r13,%rdi
xor    %esi,%esi
call   1040 <memset@plt>
cmpl   $0x4,0xc(%rsp)
jae    1af2 <realSphericalHarmonics+0x992>
mov    (%rsp),%rdi
mov    0x10(%rsp),%rsi
jmp    1c09 <realSphericalHarmonics+0xaa9>
movapd %xmm2,%xmm1
mulsd  %xmm2,%xmm1
movsd  0xa2c(%rip),%xmm0
subsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
mov    0x28(%rsp),%rbx
jb     16bf <realSphericalHarmonics+0x55f>
sqrtsd %xmm0,%xmm0
jmp    16d0 <realSphericalHarmonics+0x570>
movapd %xmm1,%xmm0
call   1060 <sin@plt>
movsd  %xmm0,0x18(%rsp)
movsd  %xmm0,0x8(%r13)
movsd  0x38(%rsp),%xmm0
call   1070 <cos@plt>
mov    (%rsp),%rdi
movsd  %xmm0,0x8(%r14)
cmp    $0x2,%ebx
jb     1c2d <realSphericalHarmonics+0xacd>
movapd %xmm0,%xmm2
movsd  0x18(%rsp),%xmm3
mulsd  %xmm3,%xmm2
addsd  %xmm2,%xmm2
movsd  %xmm2,0x10(%r13)
movsd  0x8(%r13),%xmm1
mulsd  %xmm1,%xmm3
mulsd  %xmm0,%xmm0
subsd  %xmm3,%xmm0
movsd  %xmm0,0x10(%r14)
cmpl   $0x3,0xc(%rsp)
je     1c2d <realSphericalHarmonics+0xacd>
mov    $0x3,%eax
nopw   0x0(%rax,%rax,1)
movsd  0x8(%r14),%xmm3
movapd %xmm3,%xmm4
mulsd  %xmm2,%xmm4
mulsd  %xmm0,%xmm1
movapd %xmm1,%xmm5
addsd  %xmm4,%xmm5
movsd  %xmm5,0x0(%r13,%rax,8)
movsd  0x8(%r13),%xmm1
mulsd  %xmm1,%xmm2
mulsd  %xmm0,%xmm3
subsd  %xmm2,%xmm3
movsd  %xmm3,(%r14,%rax,8)
add    $0x1,%rax
movapd %xmm3,%xmm0
movapd %xmm5,%xmm2
cmp    %rax,%rbp
jne    1670 <realSphericalHarmonics+0x510>
jmp    1c2d <realSphericalHarmonics+0xacd>
movsd  %xmm2,0x18(%rsp)
call   1080 <sqrt@plt>
movsd  0x18(%rsp),%xmm2
mov    (%rsp),%rdi
movsd  %xmm0,0x8(%r13)
mov    0xc(%rsp),%eax
mov    %rax,0x50(%rsp)
cmp    $0x2,%ebx
jl     17a4 <realSphericalHarmonics+0x644>
mulsd  %xmm0,%xmm0
movsd  %xmm0,0x10(%r13)
cmpl   $0x3,0xc(%rsp)
je     17a4 <realSphericalHarmonics+0x644>
mov    0x50(%rsp),%rax
lea    -0x3(%rax),%rcx
lea    -0x4(%rax),%rsi
mov    %ecx,%eax
and    $0x3,%eax
mov    $0x3,%edx
cmp    $0x3,%rsi
jb     1778 <realSphericalHarmonics+0x618>
and    $0xfffffffffffffffc,%rcx
neg    %rcx
mov    $0x3,%edx
nopw   0x0(%rax,%rax,1)
mulsd  0x8(%r13),%xmm0
movsd  %xmm0,0x0(%r13,%rdx,8)
movsd  0x8(%r13),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%r13,%rdx,8)
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x10(%r13,%rdx,8)
mulsd  0x8(%r13),%xmm0
movsd  %xmm0,0x18(%r13,%rdx,8)
lea    (%rcx,%rdx,1),%rsi
add    $0x4,%rsi
add    $0x4,%rdx
cmp    $0x3,%rsi
jne    1730 <realSphericalHarmonics+0x5d0>
test   %rax,%rax
je     17a4 <realSphericalHarmonics+0x644>
lea    0x0(,%rdx,8),%rcx
add    %r13,%rcx
xor    %edx,%edx
nopw   0x0(%rax,%rax,1)
mulsd  0x8(%r13),%xmm0
movsd  %xmm0,(%rcx,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1790 <realSphericalHarmonics+0x630>
test   %ebx,%ebx
js     1d73 <realSphericalHarmonics+0xc13>
mov    %r12,0x60(%rsp)
mov    $0x1,%esi
mov    $0xffffffffffffffff,%r15
mov    $0xfffffffffffffffe,%rax
mov    $0x2,%r12d
xor    %ebp,%ebp
movdqa 0x84b(%rip),%xmm8
pcmpeqd %xmm9,%xmm9
movdqa 0x84d(%rip),%xmm10
xor    %ecx,%ecx
xor    %r11d,%r11d
mov    $0x2,%r14d
xor    %r9d,%r9d
mov    %r13,0x20(%rsp)
jmp    1843 <realSphericalHarmonics+0x6e3>
nopl   0x0(%rax,%rax,1)
movabs $0x3ff0000000000000,%rax
mov    %rax,(%rdi)
add    $0x1,%r9
add    $0x1,%r15
add    %r14,%r11
mov    0x30(%rsp),%rcx
add    %rcx,%rbp
mov    0x18(%rsp),%rax
add    %rax,%r12
add    $0x1,%rsi
add    $0x2,%r14
add    $0x2,%rcx
add    $0x2,%rax
cmp    0x50(%rsp),%r9
je     141d <realSphericalHarmonics+0x2bd>
test   %r9,%r9
mov    %rax,0x18(%rsp)
mov    %rcx,0x30(%rsp)
je     1800 <realSphericalHarmonics+0x6a0>
movsd  0x0(%r13,%r9,8),%xmm3
lea    (%r9,%r9,1),%rax
movslq %eax,%rdx
mov    0x60(%rsp),%rcx
mulsd  -0x8(%rcx,%rdx,8),%xmm3
mov    %r9,%rdx
imul   %r9,%rdx
mov    %rsi,0x10(%rsp)
lea    -0x1(%r9),%rsi
lea    (%rsi,%rsi,1),%edi
add    $0x1,%edi
xorps  %xmm4,%xmm4
cvtsi2sd %edi,%xmm4
mov    %edx,%edi
lea    (%rdi,%r9,2),%rdi
mov    (%rsp),%rcx
movsd  %xmm3,(%rcx,%rdi,8)
mulsd  %xmm2,%xmm4
lea    (%r9,%rsi,1),%edi
add    $0xffffffff,%edi
imul   %rsi,%rsi
movslq %edi,%rdi
add    %rsi,%rdi
mov    (%rsp),%rcx
mulsd  (%rcx,%rdi,8),%xmm4
mov    (%rsp),%rdi
lea    (%r9,%r9,1),%rsi
add    $0xffffffffffffffff,%rsi
add    %rdx,%rsi
movsd  %xmm4,(%rdi,%rsi,8)
mov    0x10(%rsp),%rsi
cmp    $0x1,%r9
je     180d <realSphericalHarmonics+0x6ad>
mov    %r15d,%r10d
add    $0xffffffff,%eax
xorps  %xmm3,%xmm3
cvtsi2sd %eax,%xmm3
mulsd  %xmm2,%xmm3
cmp    $0x2,%r10d
jae    1950 <realSphericalHarmonics+0x7f0>
xor    %r8d,%r8d
mov    %r9d,%eax
sub    %r8d,%eax
nopl   0x0(%rax,%rax,1)
lea    (%r8,%rbp,1),%ecx
lea    (%r15,%r8,1),%edx
xorps  %xmm0,%xmm0
cvtsi2sd %edx,%xmm0
lea    (%r12,%r8,1),%edx
movslq %edx,%rdx
mulsd  (%rdi,%rdx,8),%xmm0
movslq %ecx,%rcx
movsd  (%rdi,%rcx,8),%xmm1
mulsd  %xmm3,%xmm1
subsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
divsd  %xmm0,%xmm1
lea    (%r11,%r8,1),%ecx
movsd  %xmm1,(%rdi,%rcx,8)
add    $0x1,%r8
add    $0xffffffff,%eax
cmp    %r8,%r10
jne    1900 <realSphericalHarmonics+0x7a0>
jmp    180d <realSphericalHarmonics+0x6ad>
lea    -0x1(%r9),%eax
add    $0xffffffffffffffff,%rax
mov    %r11d,%edx
add    %eax,%edx
jb     18ef <realSphericalHarmonics+0x78f>
mov    %rax,%rdx
shr    $0x20,%rdx
jne    18ef <realSphericalHarmonics+0x78f>
lea    (%rax,%rbp,1),%esi
cmp    %ebp,%esi
jl     1ae5 <realSphericalHarmonics+0x985>
test   %rdx,%rdx
mov    0x10(%rsp),%rsi
jne    18ef <realSphericalHarmonics+0x78f>
add    %r12d,%eax
cmp    %r12d,%eax
jl     18ef <realSphericalHarmonics+0x78f>
mov    $0x0,%r8d
test   %rdx,%rdx
jne    18f2 <realSphericalHarmonics+0x792>
mov    %r11d,%eax
lea    (%rdi,%rax,8),%rdx
add    %r10,%rax
lea    (%rdi,%rax,8),%rsi
movslq %ebp,%rax
lea    (%rdi,%rax,8),%r8
add    %r10,%rax
lea    (%rdi,%rax,8),%rax
mov    %rdi,%rbx
movslq %r12d,%rdi
lea    (%rbx,%rdi,8),%rcx
add    %r10,%rdi
lea    (%rbx,%rdi,8),%rdi
cmp    %rax,%rdx
setb   %r13b
cmp    %rsi,%r8
setb   %r8b
cmp    %rdi,%rdx
setb   %al
cmp    %rsi,%rcx
setb   %dl
test   %r8b,%r13b
jne    1aca <realSphericalHarmonics+0x96a>
mov    $0x0,%r8d
and    %dl,%al
mov    0x20(%rsp),%r13
mov    (%rsp),%rdi
mov    0x28(%rsp),%rbx
mov    0x10(%rsp),%rsi
jne    18f2 <realSphericalHarmonics+0x792>
mov    %r15d,%edx
and    $0xfffffffe,%edx
mov    %r10d,%r8d
and    $0xfffffffe,%r8d
movapd %xmm3,%xmm4
unpcklpd %xmm3,%xmm4
movq   %r9,%xmm5
pshufd $0x44,%xmm5,%xmm5
mov    %rdi,%rbx
mov    %r12d,%edi
mov    %ebp,%eax
mov    %rsi,%r13
mov    %r11,%rsi
movdqa %xmm8,%xmm6
nop
cltq
movupd (%rbx,%rax,8),%xmm7
mulpd  %xmm4,%xmm7
movd   %r13d,%xmm0
lea    -0x1(%r13),%ecx
movd   %ecx,%xmm1
punpckldq %xmm0,%xmm1
paddd  %xmm9,%xmm1
cvtdq2pd %xmm1,%xmm0
movslq %edi,%rdi
movupd (%rbx,%rdi,8),%xmm1
mulpd  %xmm0,%xmm1
subpd  %xmm1,%xmm7
movdqa %xmm5,%xmm0
psubq  %xmm6,%xmm0
pshufd $0xe8,%xmm0,%xmm0
cvtdq2pd %xmm0,%xmm0
divpd  %xmm0,%xmm7
mov    %esi,%ecx
movupd %xmm7,(%rbx,%rcx,8)
paddq  %xmm10,%xmm6
add    $0x2,%r13
add    $0x2,%rsi
add    $0x2,%eax
add    $0x2,%edi
add    $0xfffffffffffffffe,%rdx
jne    1a40 <realSphericalHarmonics+0x8e0>
cmp    %r10,%r8
mov    0x20(%rsp),%r13
mov    %rbx,%rdi
mov    0x28(%rsp),%rbx
mov    0x10(%rsp),%rsi
je     180d <realSphericalHarmonics+0x6ad>
jmp    18f2 <realSphericalHarmonics+0x792>
xor    %r8d,%r8d
mov    0x20(%rsp),%r13
mov    (%rsp),%rdi
mov    0x28(%rsp),%rbx
mov    0x10(%rsp),%rsi
jmp    18f2 <realSphericalHarmonics+0x792>
xor    %r8d,%r8d
mov    0x10(%rsp),%rsi
jmp    18f2 <realSphericalHarmonics+0x792>
mov    0x10(%rsp),%rsi
mov    %esi,%ebp
and    $0xfffffffc,%ebp
lea    -0x4(%rbp),%rdx
mov    %rdx,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %ecx,%eax
and    $0x7,%eax
cmp    $0x1c,%rdx
jae    1b21 <realSphericalHarmonics+0x9c1>
xor    %edx,%edx
mov    (%rsp),%rdi
jmp    1bd5 <realSphericalHarmonics+0xa75>
and    $0xfffffffffffffff8,%rcx
xor    %edx,%edx
movapd 0x521(%rip),%xmm0
mov    (%rsp),%rdi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd %xmm0,(%r14,%rdx,8)
movupd %xmm0,0x10(%r14,%rdx,8)
movupd %xmm0,0x20(%r14,%rdx,8)
movupd %xmm0,0x30(%r14,%rdx,8)
movupd %xmm0,0x40(%r14,%rdx,8)
movupd %xmm0,0x50(%r14,%rdx,8)
movupd %xmm0,0x60(%r14,%rdx,8)
movupd %xmm0,0x70(%r14,%rdx,8)
movupd %xmm0,0x80(%r14,%rdx,8)
movupd %xmm0,0x90(%r14,%rdx,8)
movupd %xmm0,0xa0(%r14,%rdx,8)
movupd %xmm0,0xb0(%r14,%rdx,8)
movupd %xmm0,0xc0(%r14,%rdx,8)
movupd %xmm0,0xd0(%r14,%rdx,8)
movupd %xmm0,0xe0(%r14,%rdx,8)
movupd %xmm0,0xf0(%r14,%rdx,8)
add    $0x20,%rdx
add    $0xfffffffffffffff8,%rcx
jne    1b40 <realSphericalHarmonics+0x9e0>
test   %rax,%rax
je     1c04 <realSphericalHarmonics+0xaa4>
lea    (%r14,%rdx,8),%rcx
add    $0x10,%rcx
shl    $0x5,%rax
xor    %edx,%edx
movapd 0x460(%rip),%xmm0
movupd %xmm0,-0x10(%rcx,%rdx,1)
movupd %xmm0,(%rcx,%rdx,1)
add    $0x20,%rdx
cmp    %rdx,%rax
jne    1bf0 <realSphericalHarmonics+0xa90>
cmp    %rsi,%rbp
je     1c2d <realSphericalHarmonics+0xacd>
movabs $0x3ff0000000000000,%rax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %rax,(%r14,%rbp,8)
add    $0x1,%rbp
cmp    %rbp,%rsi
jne    1c20 <realSphericalHarmonics+0xac0>
test   %ebx,%ebx
js     1d73 <realSphericalHarmonics+0xc13>
mov    $0xffffffff,%r8d
mov    $0x1,%r12d
mov    $0x2,%edx
xor    %r11d,%r11d
movsd  0x3bf(%rip),%xmm0
mov    $0x2,%r9d
jmp    1c7c <realSphericalHarmonics+0xb1c>
nopl   0x0(%rax)
add    $0x1,%r11
add    %r9d,%r8d
add    %r9,%r12
add    %r9,%rdx
add    $0x2,%r9
cmp    0x10(%rsp),%r11
je     1d73 <realSphericalHarmonics+0xc13>
mov    %r11d,%r10d
imul   %r10d,%r10d
lea    (%r10,%r11,1),%rax
movsd  (%rdi,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r15,%rax,8)
test   %r11,%r11
je     1c60 <realSphericalHarmonics+0xb00>
mov    $0x1,%eax
cmp    $0x1,%r11
mov    %rdi,%rcx
je     1d2e <realSphericalHarmonics+0xbce>
mov    %r11,%rsi
and    $0xfffffffffffffffe,%rsi
neg    %rsi
mov    $0x1,%eax
mov    %r8d,%edi
nop
lea    (%r12,%rax,1),%ebp
add    $0xffffffff,%ebp
movsd  (%rcx,%rbp,8),%xmm1
movsd  0x0(%r13,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
movslq %edi,%rdi
mulsd  (%r14,%rax,8),%xmm1
movsd  %xmm2,(%r15,%rdi,8)
movsd  %xmm1,(%r15,%rbp,8)
lea    -0x1(%rdx,%rax,1),%ebp
movsd  (%rcx,%rbp,8),%xmm1
movsd  0x8(%r13,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
lea    -0x1(%rdi),%ebx
movslq %ebx,%rbx
movsd  %xmm2,(%r15,%rbx,8)
mulsd  0x8(%r14,%rax,8),%xmm1
movsd  %xmm1,(%r15,%rbp,8)
add    $0xfffffffe,%edi
lea    (%rsi,%rax,1),%rbx
add    $0x2,%rbx
add    $0x2,%rax
cmp    $0x1,%rbx
jne    1cc0 <realSphericalHarmonics+0xb60>
test   $0x1,%r11b
mov    %rcx,%rdi
je     1c60 <realSphericalHarmonics+0xb00>
add    %r11d,%r10d
lea    (%r10,%rax,1),%esi
movsd  (%rdi,%rsi,8),%xmm1
movsd  0x0(%r13,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
sub    %eax,%r10d
movslq %r10d,%rdi
movsd  %xmm2,(%r15,%rdi,8)
mov    (%rsp),%rdi
mulsd  (%r14,%rax,8),%xmm1
movsd  %xmm1,(%r15,%rsi,8)
jmp    1c60 <realSphericalHarmonics+0xb00>
call   1030 <free@plt>
mov    %r13,%rdi
call   1030 <free@plt>
mov    %r14,%rdi
call   1030 <free@plt>
mov    %r15,%rax
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@M_PI = global double 0.000000e+00
define double* @realSphericalHarmonics(i32 %0, double %1, double %2, double* %3, double* %4) {
%6 = alloca i32
%7 = alloca double
%8 = alloca double
%9 = alloca double*
%10 = alloca double*
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
%14 = alloca double*
%15 = alloca double*
%16 = alloca double*
%17 = alloca double
%18 = alloca double
%19 = alloca double*
store i32 %0, i32* %6
store double %1, double* %7
store double %2, double* %8
store double* %3, double** %9
store double* %4, double** %10
%20 = load i32, i32* %6
%21 = add nsw i32 %20, 2
%22 = call double* @calloc(i32 %21, i32 8)
store double* %22, double** %15
%23 = load i32, i32* %6
%24 = add nsw i32 %23, 2
%25 = call double* @calloc(i32 %24, i32 8)
store double* %25, double** %16
%26 = load i32, i32* %6
%27 = add nsw i32 %26, 1
%28 = load i32, i32* %6
%29 = add nsw i32 %28, 1
%30 = mul nsw i32 %27, %29
%31 = call double* @calloc(i32 %30, i32 8)
store double* %31, double** %19
%32 = load i32, i32* %6
%33 = add nsw i32 %32, 1
%34 = load i32, i32* %6
%35 = add nsw i32 %34, 1
%36 = mul nsw i32 %33, %35
%37 = call double* @calloc(i32 %36, i32 8)
store double* %37, double** %14
%38 = load double, double* %7
%39 = call double @cos(double %38)
store double %39, double* %13
%40 = load double, double* %13
%41 = fcmp oeq double 1.000000e+00, %40
br i1 %41, label %42, label %62
42: ; preds = %5
store i32 0, i32* %11
br label %43
43: ; preds = %58, %42
%44 = load i32, i32* %11
%45 = load i32, i32* %6
%46 = icmp sle i32 %44, %45
br i1 %46, label %47, label %61
47: ; preds = %43
%48 = load double*, double** %14
%49 = load i32, i32* %11
%50 = load i32, i32* %11
%51 = mul nsw i32 %49, %50
%52 = load i32, i32* %11
%53 = sub nsw i32 0, %52
%54 = sub nsw i32 0, %53
%55 = add nsw i32 %51, %54
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds double, double* %48, i64 %56
store double 1.000000e+00, double* %57
br label %58
58: ; preds = %47
%59 = load i32, i32* %11
%60 = add nsw i32 %59, 1
store i32 %60, i32* %11
br label %43
61: ; preds = %43
br label %292
62: ; preds = %5
%63 = load double, double* %13
%64 = fcmp oeq double -1.000000e+00, %63
br i1 %64, label %65, label %102
65: ; preds = %62
%66 = load double*, double** %14
%67 = getelementptr inbounds double, double* %66, i64 0
store double 1.000000e+00, double* %67
store i32 1, i32* %11
br label %68
68: ; preds = %98, %65
%69 = load i32, i32* %11
%70 = load i32, i32* %6
%71 = icmp sle i32 %69, %70
br i1 %71, label %72, label %101
72: ; preds = %68
%73 = load double*, double** %14
%74 = load i32, i32* %11
%75 = sub nsw i32 %74, 1
%76 = load i32, i32* %11
%77 = sub nsw i32 %76, 1
%78 = mul nsw i32 %75, %77
%79 = load i32, i32* %11
%80 = sub nsw i32 %79, 1
%81 = sub nsw i32 0, %80
%82 = sub nsw i32 0, %81
%83 = add nsw i32 %78, %82
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %73, i64 %84
%86 = load double, double* %85
%87 = fneg double %86
%88 = load double*, double** %14
%89 = load i32, i32* %11
%90 = load i32, i32* %11
%91 = mul nsw i32 %89, %90
%92 = load i32, i32* %11
%93 = sub nsw i32 0, %92
%94 = sub nsw i32 0, %93
%95 = add nsw i32 %91, %94
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds double, double* %88, i64 %96
store double %87, double* %97
br label %98
98: ; preds = %72
%99 = load i32, i32* %11
%100 = add nsw i32 %99, 1
store i32 %100, i32* %11
br label %68
101: ; preds = %68
br label %291
102: ; preds = %62
%103 = load double, double* %13
%104 = load double, double* %13
%105 = fneg double %103
%106 = call double @llvm.fmuladd.f64(double %105, double %104, double 1.000000e+00)
%107 = call double @sqrt(double %106)
%108 = load double*, double** %15
%109 = getelementptr inbounds double, double* %108, i64 1
store double %107, double* %109
store i32 2, i32* %11
br label %110
110: ; preds = %129, %102
%111 = load i32, i32* %11
%112 = load i32, i32* %6
%113 = icmp sle i32 %111, %112
br i1 %113, label %114, label %132
114: ; preds = %110
%115 = load double*, double** %15
%116 = load i32, i32* %11
%117 = sub nsw i32 %116, 1
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds double, double* %115, i64 %118
%120 = load double, double* %119
%121 = load double*, double** %15
%122 = getelementptr inbounds double, double* %121, i64 1
%123 = load double, double* %122
%124 = fmul double %120, %123
%125 = load double*, double** %15
%126 = load i32, i32* %11
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds double, double* %125, i64 %127
store double %124, double* %128
br label %129
129: ; preds = %114
%130 = load i32, i32* %11
%131 = add nsw i32 %130, 1
store i32 %131, i32* %11
br label %110
132: ; preds = %110
store i32 0, i32* %11
br label %133
133: ; preds = %287, %132
%134 = load i32, i32* %11
%135 = load i32, i32* %6
%136 = icmp sle i32 %134, %135
br i1 %136, label %137, label %290
137: ; preds = %133
%138 = load i32, i32* %11
store i32 %138, i32* %12
%139 = load i32, i32* %12
%140 = icmp eq i32 0, %139
br i1 %140, label %141, label %152
141: ; preds = %137
%142 = load double*, double** %14
%143 = load i32, i32* %11
%144 = load i32, i32* %11
%145 = mul nsw i32 %143, %144
%146 = load i32, i32* %11
%147 = sub nsw i32 0, %146
%148 = sub nsw i32 0, %147
%149 = add nsw i32 %145, %148
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds double, double* %142, i64 %150
store double 1.000000e+00, double* %151
br label %286
152: ; preds = %137
%153 = load double*, double** %15
%154 = load i32, i32* %12
%155 = sext i32 %154 to i64
%156 = getelementptr inbounds double, double* %153, i64 %155
%157 = load double, double* %156
%158 = load double*, double** %10
%159 = load i32, i32* %12
%160 = mul nsw i32 2, %159
%161 = sub nsw i32 %160, 1
%162 = sext i32 %161 to i64
%163 = getelementptr inbounds double, double* %158, i64 %162
%164 = load double, double* %163
%165 = fmul double %157, %164
%166 = load double*, double** %14
%167 = load i32, i32* %11
%168 = load i32, i32* %11
%169 = mul nsw i32 %167, %168
%170 = load i32, i32* %12
%171 = load i32, i32* %11
%172 = sub nsw i32 0, %171
%173 = sub nsw i32 %170, %172
%174 = add nsw i32 %169, %173
%175 = sext i32 %174 to i64
%176 = getelementptr inbounds double, double* %166, i64 %175
store double %165, double* %176
%177 = load i32, i32* %11
%178 = sub nsw i32 %177, 1
store i32 %178, i32* %12
%179 = load double, double* %13
%180 = load i32, i32* %12
%181 = mul nsw i32 2, %180
%182 = add nsw i32 %181, 1
%183 = sitofp i32 %182 to double
%184 = fmul double %179, %183
%185 = load double*, double** %14
%186 = load i32, i32* %11
%187 = sub nsw i32 %186, 1
%188 = load i32, i32* %11
%189 = sub nsw i32 %188, 1
%190 = mul nsw i32 %187, %189
%191 = load i32, i32* %12
%192 = load i32, i32* %11
%193 = sub nsw i32 %192, 1
%194 = sub nsw i32 0, %193
%195 = sub nsw i32 %191, %194
%196 = add nsw i32 %190, %195
%197 = sext i32 %196 to i64
%198 = getelementptr inbounds double, double* %185, i64 %197
%199 = load double, double* %198
%200 = fmul double %184, %199
%201 = load double*, double** %14
%202 = load i32, i32* %11
%203 = load i32, i32* %11
%204 = mul nsw i32 %202, %203
%205 = load i32, i32* %12
%206 = load i32, i32* %11
%207 = sub nsw i32 0, %206
%208 = sub nsw i32 %205, %207
%209 = add nsw i32 %204, %208
%210 = sext i32 %209 to i64
%211 = getelementptr inbounds double, double* %201, i64 %210
store double %200, double* %211
%212 = load i32, i32* %11
%213 = icmp sgt i32 %212, 1
br i1 %213, label %214, label %285
214: ; preds = %152
store i32 0, i32* %12
br label %215
215: ; preds = %281, %214
%216 = load i32, i32* %12
%217 = load i32, i32* %11
%218 = sub nsw i32 %217, 2
%219 = icmp sle i32 %216, %218
br i1 %219, label %220, label %284
220: ; preds = %215
%221 = load double, double* %13
%222 = load i32, i32* %11
%223 = mul nsw i32 2, %222
%224 = sub nsw i32 %223, 1
%225 = sitofp i32 %224 to double
%226 = fmul double %221, %225
%227 = load double*, double** %14
%228 = load i32, i32* %11
%229 = sub nsw i32 %228, 1
%230 = load i32, i32* %11
%231 = sub nsw i32 %230, 1
%232 = mul nsw i32 %229, %231
%233 = load i32, i32* %12
%234 = load i32, i32* %11
%235 = sub nsw i32 %234, 1
%236 = sub nsw i32 0, %235
%237 = sub nsw i32 %233, %236
%238 = add nsw i32 %232, %237
%239 = sext i32 %238 to i64
%240 = getelementptr inbounds double, double* %227, i64 %239
%241 = load double, double* %240
%242 = load i32, i32* %11
%243 = load i32, i32* %12
%244 = add nsw i32 %242, %243
%245 = sub nsw i32 %244, 1
%246 = sitofp i32 %245 to double
%247 = load double*, double** %14
%248 = load i32, i32* %11
%249 = sub nsw i32 %248, 2
%250 = load i32, i32* %11
%251 = sub nsw i32 %250, 2
%252 = mul nsw i32 %249, %251
%253 = load i32, i32* %12
%254 = load i32, i32* %11
%255 = sub nsw i32 %254, 2
%256 = sub nsw i32 0, %255
%257 = sub nsw i32 %253, %256
%258 = add nsw i32 %252, %257
%259 = sext i32 %258 to i64
%260 = getelementptr inbounds double, double* %247, i64 %259
%261 = load double, double* %260
%262 = fmul double %246, %261
%263 = fneg double %262
%264 = call double @llvm.fmuladd.f64(double %226, double %241, double %263)
%265 = load i32, i32* %11
%266 = load i32, i32* %12
%267 = sub nsw i32 %265, %266
%268 = sitofp i32 %267 to double
%269 = fdiv double %264, %268
%270 = load double*, double** %14
%271 = load i32, i32* %11
%272 = load i32, i32* %11
%273 = mul nsw i32 %271, %272
%274 = load i32, i32* %12
%275 = load i32, i32* %11
%276 = sub nsw i32 0, %275
%277 = sub nsw i32 %274, %276
%278 = add nsw i32 %273, %277
%279 = sext i32 %278 to i64
%280 = getelementptr inbounds double, double* %270, i64 %279
store double %269, double* %280
br label %281
281: ; preds = %220
%282 = load i32, i32* %12
%283 = add nsw i32 %282, 1
store i32 %283, i32* %12
br label %215
284: ; preds = %215
br label %285
285: ; preds = %284, %152
br label %286
286: ; preds = %285, %141
br label %287
287: ; preds = %286
%288 = load i32, i32* %11
%289 = add nsw i32 %288, 1
store i32 %289, i32* %11
br label %133
290: ; preds = %133
br label %291
291: ; preds = %290, %101
br label %292
292: ; preds = %291, %61
store i32 0, i32* %11
br label %293
293: ; preds = %380, %292
%294 = load i32, i32* %11
%295 = load i32, i32* %6
%296 = icmp sle i32 %294, %295
br i1 %296, label %297, label %383
297: ; preds = %293
%298 = load i32, i32* %11
%299 = sitofp i32 %298 to double
%300 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %299, double 1.000000e+00)
%301 = load double, double* @M_PI
%302 = fmul double 2.000000e+00, %301
%303 = fdiv double %300, %302
%304 = call double @sqrt(double %303)
store double %304, double* %17
%305 = load double, double* %17
%306 = load double*, double** %14
%307 = load i32, i32* %11
%308 = load i32, i32* %11
%309 = mul nsw i32 %307, %308
%310 = load i32, i32* %11
%311 = sub nsw i32 0, %310
%312 = sub nsw i32 0, %311
%313 = add nsw i32 %309, %312
%314 = sext i32 %313 to i64
%315 = getelementptr inbounds double, double* %306, i64 %314
%316 = load double, double* %315
%317 = fmul double %305, %316
%318 = load double*, double** %14
%319 = load i32, i32* %11
%320 = load i32, i32* %11
%321 = mul nsw i32 %319, %320
%322 = load i32, i32* %11
%323 = sub nsw i32 0, %322
%324 = sub nsw i32 0, %323
%325 = add nsw i32 %321, %324
%326 = sext i32 %325 to i64
%327 = getelementptr inbounds double, double* %318, i64 %326
store double %317, double* %327
store i32 1, i32* %12
br label %328
328: ; preds = %376, %297
%329 = load i32, i32* %12
%330 = load i32, i32* %11
%331 = icmp sle i32 %329, %330
br i1 %331, label %332, label %379
332: ; preds = %328
%333 = load double*, double** %9
%334 = load i32, i32* %11
%335 = load i32, i32* %12
%336 = sub nsw i32 %334, %335
%337 = sext i32 %336 to i64
%338 = getelementptr inbounds double, double* %333, i64 %337
%339 = load double, double* %338
%340 = load double*, double** %9
%341 = load i32, i32* %11
%342 = load i32, i32* %12
%343 = add nsw i32 %341, %342
%344 = sext i32 %343 to i64
%345 = getelementptr inbounds double, double* %340, i64 %344
%346 = load double, double* %345
%347 = fdiv double %339, %346
%348 = call double @sqrt(double %347)
%349 = load double, double* %17
%350 = fmul double %348, %349
store double %350, double* %18
%351 = load double, double* %18
%352 = load double*, double** %14
%353 = load i32, i32* %11
%354 = load i32, i32* %11
%355 = mul nsw i32 %353, %354
%356 = load i32, i32* %12
%357 = load i32, i32* %11
%358 = sub nsw i32 0, %357
%359 = sub nsw i32 %356, %358
%360 = add nsw i32 %355, %359
%361 = sext i32 %360 to i64
%362 = getelementptr inbounds double, double* %352, i64 %361
%363 = load double, double* %362
%364 = fmul double %351, %363
%365 = load double*, double** %14
%366 = load i32, i32* %11
%367 = load i32, i32* %11
%368 = mul nsw i32 %366, %367
%369 = load i32, i32* %12
%370 = load i32, i32* %11
%371 = sub nsw i32 0, %370
%372 = sub nsw i32 %369, %371
%373 = add nsw i32 %368, %372
%374 = sext i32 %373 to i64
%375 = getelementptr inbounds double, double* %365, i64 %374
store double %364, double* %375
br label %376
376: ; preds = %332
%377 = load i32, i32* %12
%378 = add nsw i32 %377, 1
store i32 %378, i32* %12
br label %328
379: ; preds = %328
br label %380
380: ; preds = %379
%381 = load i32, i32* %11
%382 = add nsw i32 %381, 1
store i32 %382, i32* %11
br label %293
383: ; preds = %293
%384 = load i32, i32* %6
%385 = icmp sgt i32 %384, 0
br i1 %385, label %386, label %475
386: ; preds = %383
%387 = load double, double* %8
%388 = fcmp oeq double 0.000000e+00, %387
br i1 %388, label %389, label %407
389: ; preds = %386
store i32 0, i32* %12
br label %390
390: ; preds = %403, %389
%391 = load i32, i32* %12
%392 = load i32, i32* %6
%393 = icmp sle i32 %391, %392
br i1 %393, label %394, label %406
394: ; preds = %390
%395 = load double*, double** %15
%396 = load i32, i32* %12
%397 = sext i32 %396 to i64
%398 = getelementptr inbounds double, double* %395, i64 %397
store double 0.000000e+00, double* %398
%399 = load double*, double** %16
%400 = load i32, i32* %12
%401 = sext i32 %400 to i64
%402 = getelementptr inbounds double, double* %399, i64 %401
store double 1.000000e+00, double* %402
br label %403
403: ; preds = %394
%404 = load i32, i32* %12
%405 = add nsw i32 %404, 1
store i32 %405, i32* %12
br label %390
406: ; preds = %390
br label %474
407: ; preds = %386
%408 = load double, double* %8
%409 = call double @sin(double %408)
%410 = load double*, double** %15
%411 = getelementptr inbounds double, double* %410, i64 1
store double %409, double* %411
%412 = load double, double* %8
%413 = call double @cos(double %412)
%414 = load double*, double** %16
%415 = getelementptr inbounds double, double* %414, i64 1
store double %413, double* %415
store i32 2, i32* %12
br label %416
416: ; preds = %470, %407
%417 = load i32, i32* %12
%418 = load i32, i32* %6
%419 = icmp sle i32 %417, %418
br i1 %419, label %420, label %473
420: ; preds = %416
%421 = load double*, double** %15
%422 = getelementptr inbounds double, double* %421, i64 1
%423 = load double, double* %422
%424 = load double*, double** %16
%425 = load i32, i32* %12
%426 = sub nsw i32 %425, 1
%427 = sext i32 %426 to i64
%428 = getelementptr inbounds double, double* %424, i64 %427
%429 = load double, double* %428
%430 = load double*, double** %16
%431 = getelementptr inbounds double, double* %430, i64 1
%432 = load double, double* %431
%433 = load double*, double** %15
%434 = load i32, i32* %12
%435 = sub nsw i32 %434, 1
%436 = sext i32 %435 to i64
%437 = getelementptr inbounds double, double* %433, i64 %436
%438 = load double, double* %437
%439 = fmul double %432, %438
%440 = call double @llvm.fmuladd.f64(double %423, double %429, double %439)
%441 = load double*, double** %15
%442 = load i32, i32* %12
%443 = sext i32 %442 to i64
%444 = getelementptr inbounds double, double* %441, i64 %443
store double %440, double* %444
%445 = load double*, double** %16
%446 = getelementptr inbounds double, double* %445, i64 1
%447 = load double, double* %446
%448 = load double*, double** %16
%449 = load i32, i32* %12
%450 = sub nsw i32 %449, 1
%451 = sext i32 %450 to i64
%452 = getelementptr inbounds double, double* %448, i64 %451
%453 = load double, double* %452
%454 = load double*, double** %15
%455 = getelementptr inbounds double, double* %454, i64 1
%456 = load double, double* %455
%457 = load double*, double** %15
%458 = load i32, i32* %12
%459 = sub nsw i32 %458, 1
%460 = sext i32 %459 to i64
%461 = getelementptr inbounds double, double* %457, i64 %460
%462 = load double, double* %461
%463 = fmul double %456, %462
%464 = fneg double %463
%465 = call double @llvm.fmuladd.f64(double %447, double %453, double %464)
%466 = load double*, double** %16
%467 = load i32, i32* %12
%468 = sext i32 %467 to i64
%469 = getelementptr inbounds double, double* %466, i64 %468
store double %465, double* %469
br label %470
470: ; preds = %420
%471 = load i32, i32* %12
%472 = add nsw i32 %471, 1
store i32 %472, i32* %12
br label %416
473: ; preds = %416
br label %474
474: ; preds = %473, %406
br label %475
475: ; preds = %474, %383
store i32 0, i32* %11
br label %476
476: ; preds = %572, %475
%477 = load i32, i32* %11
%478 = load i32, i32* %6
%479 = icmp sle i32 %477, %478
br i1 %479, label %480, label %575
480: ; preds = %476
%481 = load double*, double** %14
%482 = load i32, i32* %11
%483 = load i32, i32* %11
%484 = mul nsw i32 %482, %483
%485 = load i32, i32* %11
%486 = sub nsw i32 0, %485
%487 = sub nsw i32 0, %486
%488 = add nsw i32 %484, %487
%489 = sext i32 %488 to i64
%490 = getelementptr inbounds double, double* %481, i64 %489
%491 = load double, double* %490
%492 = call double @sqrt(double 2.000000e+00)
%493 = fdiv double %491, %492
%494 = load double*, double** %19
%495 = load i32, i32* %11
%496 = load i32, i32* %11
%497 = mul nsw i32 %495, %496
%498 = load i32, i32* %11
%499 = sub nsw i32 0, %498
%500 = sub nsw i32 0, %499
%501 = add nsw i32 %497, %500
%502 = sext i32 %501 to i64
%503 = getelementptr inbounds double, double* %494, i64 %502
store double %493, double* %503
store i32 1, i32* %12
br label %504
504: ; preds = %568, %480
%505 = load i32, i32* %12
%506 = load i32, i32* %11
%507 = icmp sle i32 %505, %506
br i1 %507, label %508, label %571
508: ; preds = %504
%509 = load double*, double** %14
%510 = load i32, i32* %11
%511 = load i32, i32* %11
%512 = mul nsw i32 %510, %511
%513 = load i32, i32* %12
%514 = load i32, i32* %11
%515 = sub nsw i32 0, %514
%516 = sub nsw i32 %513, %515
%517 = add nsw i32 %512, %516
%518 = sext i32 %517 to i64
%519 = getelementptr inbounds double, double* %509, i64 %518
%520 = load double, double* %519
%521 = load double*, double** %15
%522 = load i32, i32* %12
%523 = sext i32 %522 to i64
%524 = getelementptr inbounds double, double* %521, i64 %523
%525 = load double, double* %524
%526 = fmul double %520, %525
%527 = load double*, double** %19
%528 = load i32, i32* %11
%529 = load i32, i32* %11
%530 = mul nsw i32 %528, %529
%531 = load i32, i32* %12
%532 = sub nsw i32 0, %531
%533 = load i32, i32* %11
%534 = sub nsw i32 0, %533
%535 = sub nsw i32 %532, %534
%536 = add nsw i32 %530, %535
%537 = sext i32 %536 to i64
%538 = getelementptr inbounds double, double* %527, i64 %537
store double %526, double* %538
%539 = load double*, double** %14
%540 = load i32, i32* %11
%541 = load i32, i32* %11
%542 = mul nsw i32 %540, %541
%543 = load i32, i32* %12
%544 = load i32, i32* %11
%545 = sub nsw i32 0, %544
%546 = sub nsw i32 %543, %545
%547 = add nsw i32 %542, %546
%548 = sext i32 %547 to i64
%549 = getelementptr inbounds double, double* %539, i64 %548
%550 = load double, double* %549
%551 = load double*, double** %16
%552 = load i32, i32* %12
%553 = sext i32 %552 to i64
%554 = getelementptr inbounds double, double* %551, i64 %553
%555 = load double, double* %554
%556 = fmul double %550, %555
%557 = load double*, double** %19
%558 = load i32, i32* %11
%559 = load i32, i32* %11
%560 = mul nsw i32 %558, %559
%561 = load i32, i32* %12
%562 = load i32, i32* %11
%563 = sub nsw i32 0, %562
%564 = sub nsw i32 %561, %563
%565 = add nsw i32 %560, %564
%566 = sext i32 %565 to i64
%567 = getelementptr inbounds double, double* %557, i64 %566
store double %556, double* %567
br label %568
568: ; preds = %508
%569 = load i32, i32* %12
%570 = add nsw i32 %569, 1
store i32 %570, i32* %12
br label %504
571: ; preds = %504
br label %572
572: ; preds = %571
%573 = load i32, i32* %11
%574 = add nsw i32 %573, 1
store i32 %574, i32* %11
br label %476
575: ; preds = %476
%576 = load double*, double** %14
%577 = call i32 @free(double* %576)
%578 = load double*, double** %15
%579 = call i32 @free(double* %578)
%580 = load double*, double** %16
%581 = call i32 @free(double* %580)
%582 = load double*, double** %19
ret double* %582
}

LLVM IR (-O2):
@M_PI = global double 0.000000e+00
define noalias double* @realSphericalHarmonics(i32 %0, double %1, double %2, double* nocapture %3, double* nocapture %4) {
%6 = add nsw i32 %0, 2
%7 = tail call double* @calloc(i32 %6, i32 8)
%8 = bitcast double* %7 to i8*
%9 = tail call double* @calloc(i32 %6, i32 8)
%10 = add i32 %0, 1
%11 = mul nsw i32 %10, %10
%12 = tail call double* @calloc(i32 %11, i32 8)
%13 = tail call double* @calloc(i32 %11, i32 8)
%14 = tail call double @cos(double %1)
%15 = fcmp oeq double %14, 1.000000e+00
br i1 %15, label %16, label %50
16: ; preds = %5
%17 = icmp slt i32 %0, 0
br i1 %17, label %602, label %18
18: ; preds = %16
%19 = zext i32 %10 to i64
%20 = add nsw i64 %19, -1
%21 = and i64 %19, 3
%22 = icmp ult i64 %20, 3
br i1 %22, label %330, label %23
23: ; preds = %18
%24 = and i64 %19, 4294967292
br label %25
25: ; preds = %25, %23
%26 = phi i64 [ 0, %23 ], [ %47, %25 ]
%27 = phi i64 [ 0, %23 ], [ %48, %25 ]
%28 = mul i64 %26, %26
%29 = and i64 %28, 4294967280
%30 = add nuw nsw i64 %29, %26
%31 = getelementptr inbounds double, double* %13, i64 %30
store double 1.000000e+00, double* %31
%32 = or i64 %26, 1
%33 = mul i64 %32, %32
%34 = and i64 %33, 4294967293
%35 = add nuw nsw i64 %34, %32
%36 = getelementptr inbounds double, double* %13, i64 %35
store double 1.000000e+00, double* %36
%37 = or i64 %26, 2
%38 = mul i64 %37, %37
%39 = and i64 %38, 4294967292
%40 = add nuw nsw i64 %39, %37
%41 = getelementptr inbounds double, double* %13, i64 %40
store double 1.000000e+00, double* %41
%42 = or i64 %26, 3
%43 = mul i64 %42, %42
%44 = and i64 %43, 4294967293
%45 = add nuw nsw i64 %44, %42
%46 = getelementptr inbounds double, double* %13, i64 %45
store double 1.000000e+00, double* %46
%47 = add nuw nsw i64 %26, 4
%48 = add i64 %27, 4
%49 = icmp eq i64 %48, %24
br i1 %49, label %330, label %25
50: ; preds = %5
%51 = fcmp oeq double %14, -1.000000e+00
br i1 %51, label %52, label %89
52: ; preds = %50
store double 1.000000e+00, double* %13
%53 = icmp slt i32 %0, 1
br i1 %53, label %358, label %54
54: ; preds = %52
%55 = zext i32 %10 to i64
%56 = add nsw i64 %55, -1
%57 = and i64 %56, 1
%58 = icmp eq i32 %10, 2
br i1 %58, label %343, label %59
59: ; preds = %54
%60 = and i64 %56, -2
br label %61
61: ; preds = %61, %59
%62 = phi i64 [ 1, %59 ], [ %86, %61 ]
%63 = phi i64 [ 0, %59 ], [ %87, %61 ]
%64 = add nsw i64 %62, -1
%65 = mul i64 %64, %62
%66 = shl i64 %65, 32
%67 = ashr exact i64 %66, 32
%68 = getelementptr inbounds double, double* %13, i64 %67
%69 = load double, double* %68
%70 = fneg double %69
%71 = mul i64 %62, %62
%72 = add i64 %71, %62
%73 = and i64 %72, 4294967295
%74 = getelementptr inbounds double, double* %13, i64 %73
store double %70, double* %74
%75 = add nuw nsw i64 %62, 1
%76 = mul i64 %62, %75
%77 = shl i64 %76, 32
%78 = ashr exact i64 %77, 32
%79 = getelementptr inbounds double, double* %13, i64 %78
%80 = load double, double* %79
%81 = fneg double %80
%82 = mul i64 %75, %75
%83 = add i64 %82, %75
%84 = and i64 %83, 4294967295
%85 = getelementptr inbounds double, double* %13, i64 %84
store double %81, double* %85
%86 = add nuw nsw i64 %62, 2
%87 = add i64 %63, 2
%88 = icmp eq i64 %87, %60
br i1 %88, label %343, label %61
89: ; preds = %50
%90 = fneg double %14
%91 = tail call double @llvm.fmuladd.f64(double %90, double %14, double 1.000000e+00)
%92 = tail call double @sqrt(double %91)
%93 = getelementptr inbounds double, double* %7, i64 1
store double %92, double* %93
%94 = icmp slt i32 %0, 2
br i1 %94, label %121, label %95
95: ; preds = %89
%96 = zext i32 %10 to i64
%97 = fmul double %92, %92
%98 = getelementptr inbounds double, double* %7, i64 2
store double %97, double* %98
%99 = icmp eq i32 %10, 3
br i1 %99, label %121, label %100
100: ; preds = %95
%101 = add nsw i64 %96, -3
%102 = add nsw i64 %96, -4
%103 = and i64 %101, 3
%104 = icmp ult i64 %102, 3
br i1 %104, label %107, label %105
105: ; preds = %100
%106 = and i64 %101, -4
br label %125
107: ; preds = %125, %100
%108 = phi i64 [ 3, %100 ], [ %143, %125 ]
%109 = phi double [ %97, %100 ], [ %141, %125 ]
%110 = icmp eq i64 %103, 0
br i1 %110, label %121, label %111
111: ; preds = %107, %111
%112 = phi i64 [ %118, %111 ], [ %108, %107 ]
%113 = phi double [ %116, %111 ], [ %109, %107 ]
%114 = phi i64 [ %119, %111 ], [ 0, %107 ]
%115 = load double, double* %93
%116 = fmul double %113, %115
%117 = getelementptr inbounds double, double* %7, i64 %112
store double %116, double* %117
%118 = add nuw nsw i64 %112, 1
%119 = add i64 %114, 1
%120 = icmp eq i64 %119, %103
br i1 %120, label %121, label %111
121: ; preds = %107, %111, %95, %89
%122 = icmp slt i32 %0, 0
br i1 %122, label %602, label %123
123: ; preds = %121
%124 = zext i32 %10 to i64
br label %146
125: ; preds = %125, %105
%126 = phi i64 [ 3, %105 ], [ %143, %125 ]
%127 = phi double [ %97, %105 ], [ %141, %125 ]
%128 = phi i64 [ 0, %105 ], [ %144, %125 ]
%129 = load double, double* %93
%130 = fmul double %127, %129
%131 = getelementptr inbounds double, double* %7, i64 %126
store double %130, double* %131
%132 = add nuw nsw i64 %126, 1
%133 = load double, double* %93
%134 = fmul double %130, %133
%135 = getelementptr inbounds double, double* %7, i64 %132
store double %134, double* %135
%136 = add nuw nsw i64 %126, 2
%137 = fmul double %134, %133
%138 = getelementptr inbounds double, double* %7, i64 %136
store double %137, double* %138
%139 = add nuw nsw i64 %126, 3
%140 = load double, double* %93
%141 = fmul double %137, %140
%142 = getelementptr inbounds double, double* %7, i64 %139
store double %141, double* %142
%143 = add nuw nsw i64 %126, 4
%144 = add i64 %128, 4
%145 = icmp eq i64 %144, %106
br i1 %145, label %107, label %125
146: ; preds = %123, %323
%147 = phi i32 [ 2, %123 ], [ %329, %323 ]
%148 = phi i32 [ 0, %123 ], [ %328, %323 ]
%149 = phi i32 [ 0, %123 ], [ %327, %323 ]
%150 = phi i64 [ 0, %123 ], [ %324, %323 ]
%151 = phi i32 [ -1, %123 ], [ %325, %323 ]
%152 = zext i32 %151 to i64
%153 = zext i32 %149 to i64
%154 = getelementptr double, double* %13, i64 %153
%155 = zext i32 %151 to i64
%156 = add nuw nsw i64 %153, %155
%157 = getelementptr double, double* %13, i64 %156
%158 = sext i32 %148 to i64
%159 = getelementptr double, double* %13, i64 %158
%160 = add nsw i64 %158, %155
%161 = getelementptr double, double* %13, i64 %160
%162 = sext i32 %147 to i64
%163 = getelementptr double, double* %13, i64 %162
%164 = add nsw i64 %162, %155
%165 = getelementptr double, double* %13, i64 %164
%166 = trunc i64 %150 to i32
%167 = shl i32 %166, 1
%168 = add i32 %167, -2
%169 = trunc i64 %150 to i32
%170 = shl i32 %169, 1
%171 = trunc i64 %150 to i32
%172 = shl i32 %171, 1
%173 = add i32 %172, 2
%174 = add nuw i64 %150, 4294967295
%175 = and i64 %174, 4294967295
%176 = add nsw i64 %175, -1
%177 = trunc i64 %150 to i32
%178 = icmp eq i64 %150, 0
br i1 %178, label %179, label %180
179: ; preds = %146
store double 1.000000e+00, double* %13
br label %323
180: ; preds = %146
%181 = getelementptr inbounds double, double* %7, i64 %150
%182 = load double, double* %181
%183 = shl nuw nsw i64 %150, 1
%184 = trunc i64 %183 to i32
%185 = add nsw i32 %184, -1
%186 = sext i32 %185 to i64
%187 = getelementptr inbounds double, double* %4, i64 %186
%188 = load double, double* %187
%189 = fmul double %182, %188
%190 = mul nsw i64 %150, %150
%191 = and i64 %190, 4294967295
%192 = add nuw nsw i64 %191, %183
%193 = getelementptr inbounds double, double* %13, i64 %192
store double %189, double* %193
%194 = add nsw i64 %150, -1
%195 = trunc i64 %194 to i32
%196 = shl nsw i32 %195, 1
%197 = or i32 %196, 1
%198 = sitofp i32 %197 to double
%199 = fmul double %14, %198
%200 = mul nsw i64 %194, %194
%201 = add nsw i32 %177, -1
%202 = add i32 %201, %195
%203 = sext i32 %202 to i64
%204 = add nsw i64 %200, %203
%205 = getelementptr inbounds double, double* %13, i64 %204
%206 = load double, double* %205
%207 = fmul double %199, %206
%208 = add nsw i64 %194, %150
%209 = add nsw i64 %208, %190
%210 = getelementptr inbounds double, double* %13, i64 %209
store double %207, double* %210
%211 = icmp eq i64 %150, 1
br i1 %211, label %323, label %212
212: ; preds = %180
%213 = add nsw i32 %177, -2
%214 = sitofp i32 %185 to double
%215 = fmul double %14, %214
%216 = trunc i64 %200 to i32
%217 = add i32 %201, %216
%218 = mul nsw i32 %213, %213
%219 = trunc i64 %150 to i32
%220 = add i32 %219, -2
%221 = add i32 %218, %220
%222 = zext i32 %151 to i64
%223 = icmp ult i32 %151, 2
br i1 %223, label %294, label %224
224: ; preds = %212
%225 = trunc i64 %176 to i32
%226 = xor i32 %149, -1
%227 = icmp ult i32 %226, %225
%228 = icmp ugt i64 %176, 4294967295
%229 = or i1 %227, %228
%230 = trunc i64 %176 to i32
%231 = add i32 %148, %230
%232 = icmp slt i32 %231, %148
%233 = icmp ugt i64 %176, 4294967295
%234 = or i1 %232, %233
%235 = trunc i64 %176 to i32
%236 = add i32 %147, %235
%237 = icmp slt i32 %236, %147
%238 = icmp ugt i64 %176, 4294967295
%239 = or i1 %237, %238
%240 = or i1 %229, %234
%241 = or i1 %240, %239
br i1 %241, label %294, label %242
242: ; preds = %224
%243 = icmp ult double* %154, %161
%244 = icmp ult double* %159, %157
%245 = and i1 %243, %244
%246 = icmp ult double* %154, %165
%247 = icmp ult double* %163, %157
%248 = and i1 %246, %247
%249 = or i1 %245, %248
br i1 %249, label %294, label %250
250: ; preds = %242
%251 = and i64 %152, 4294967294
%252 = insertelement <2 x double> poison, double %215, i64 0
%253 = shufflevector <2 x double> %252, <2 x double> poison, <2 x i32> zeroinitializer
%254 = insertelement <2 x i64> poison, i64 %150, i64 0
%255 = shufflevector <2 x i64> %254, <2 x i64> poison, <2 x i32> zeroinitializer
br label %256
256: ; preds = %256, %250
%257 = phi i64 [ 0, %250 ], [ %289, %256 ]
%258 = phi <2 x i64> [ <i64 0, i64 1>, %250 ], [ %290, %256 ]
%259 = or i64 %257, 1
%260 = trunc i64 %257 to i32
%261 = add i32 %217, %260
%262 = sext i32 %261 to i64
%263 = getelementptr inbounds double, double* %13, i64 %262
%264 = bitcast double* %263 to <2 x double>*
%265 = load <2 x double>, <2 x double>* %264
%266 = add nuw nsw i64 %257, %150
%267 = add nuw nsw i64 %259, %150
%268 = insertelement <2 x i64> poison, i64 %266, i64 0
%269 = insertelement <2 x i64> %268, i64 %267, i64 1
%270 = trunc <2 x i64> %269 to <2 x i32>
%271 = add <2 x i32> %270, <i32 -1, i32 -1>
%272 = sitofp <2 x i32> %271 to <2 x double>
%273 = add i32 %221, %260
%274 = sext i32 %273 to i64
%275 = getelementptr inbounds double, double* %13, i64 %274
%276 = bitcast double* %275 to <2 x double>*
%277 = load <2 x double>, <2 x double>* %276
%278 = fneg <2 x double> %277
%279 = fmul <2 x double> %272, %278
%280 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %253, <2 x double> %265, <2 x double> %279)
%281 = sub nsw <2 x i64> %255, %258
%282 = trunc <2 x i64> %281 to <2 x i32>
%283 = sitofp <2 x i32> %282 to <2 x double>
%284 = fdiv <2 x double> %280, %283
%285 = add nuw nsw i64 %266, %190
%286 = and i64 %285, 4294967295
%287 = getelementptr inbounds double, double* %13, i64 %286
%288 = bitcast double* %287 to <2 x double>*
store <2 x double> %284, <2 x double>* %288
%289 = add nuw i64 %257, 2
%290 = add <2 x i64> %258, <i64 2, i64 2>
%291 = icmp eq i64 %289, %251
br i1 %291, label %292, label %256
292: ; preds = %256
%293 = icmp eq i64 %251, %152
br i1 %293, label %323, label %294
294: ; preds = %242, %224, %212, %292
%295 = phi i64 [ 0, %242 ], [ 0, %224 ], [ 0, %212 ], [ %251, %292 ]
br label %296
296: ; preds = %294, %296
%297 = phi i64 [ %321, %296 ], [ %295, %294 ]
%298 = trunc i64 %297 to i32
%299 = add i32 %217, %298
%300 = sext i32 %299 to i64
%301 = getelementptr inbounds double, double* %13, i64 %300
%302 = load double, double* %301
%303 = add nuw nsw i64 %297, %150
%304 = trunc i64 %303 to i32
%305 = add i32 %304, -1
%306 = sitofp i32 %305 to double
%307 = add i32 %221, %298
%308 = sext i32 %307 to i64
%309 = getelementptr inbounds double, double* %13, i64 %308
%310 = load double, double* %309
%311 = fneg double %310
%312 = fmul double %306, %311
%313 = tail call double @llvm.fmuladd.f64(double %215, double %302, double %312)
%314 = sub nsw i64 %150, %297
%315 = trunc i64 %314 to i32
%316 = sitofp i32 %315 to double
%317 = fdiv double %313, %316
%318 = add nuw nsw i64 %303, %190
%319 = and i64 %318, 4294967295
%320 = getelementptr inbounds double, double* %13, i64 %319
store double %317, double* %320
%321 = add nuw nsw i64 %297, 1
%322 = icmp eq i64 %321, %222
br i1 %322, label %323, label %296
323: ; preds = %296, %292, %179, %180
%324 = add nuw nsw i64 %150, 1
%325 = add nsw i32 %151, 1
%326 = icmp eq i64 %324, %124
%327 = add i32 %149, %173
%328 = add i32 %148, %170
%329 = add i32 %147, %168
br i1 %326, label %358, label %146
330: ; preds = %25, %18
%331 = phi i64 [ 0, %18 ], [ %47, %25 ]
%332 = icmp eq i64 %21, 0
br i1 %332, label %358, label %333
333: ; preds = %330, %333
%334 = phi i64 [ %340, %333 ], [ %331, %330 ]
%335 = phi i64 [ %341, %333 ], [ 0, %330 ]
%336 = mul i64 %334, %334
%337 = and i64 %336, 4294967295
%338 = add nuw nsw i64 %337, %334
%339 = getelementptr inbounds double, double* %13, i64 %338
store double 1.000000e+00, double* %339
%340 = add nuw nsw i64 %334, 1
%341 = add i64 %335, 1
%342 = icmp eq i64 %341, %21
br i1 %342, label %358, label %333
343: ; preds = %61, %54
%344 = phi i64 [ 1, %54 ], [ %86, %61 ]
%345 = icmp eq i64 %57, 0
br i1 %345, label %358, label %346
346: ; preds = %343
%347 = add nsw i64 %344, -1
%348 = mul i64 %347, %344
%349 = shl i64 %348, 32
%350 = ashr exact i64 %349, 32
%351 = getelementptr inbounds double, double* %13, i64 %350
%352 = load double, double* %351
%353 = fneg double %352
%354 = mul i64 %344, %344
%355 = add i64 %354, %344
%356 = and i64 %355, 4294967295
%357 = getelementptr inbounds double, double* %13, i64 %356
store double %353, double* %357
br label %358
358: ; preds = %323, %346, %343, %330, %333, %52
%359 = icmp slt i32 %0, 0
br i1 %359, label %602, label %360
360: ; preds = %358
%361 = zext i32 %10 to i64
br label %362
362: ; preds = %360, %399
%363 = phi i64 [ 0, %360 ], [ %400, %399 ]
%364 = phi i64 [ 1, %360 ], [ %401, %399 ]
%365 = trunc i64 %363 to i32
%366 = sitofp i32 %365 to double
%367 = tail call double @llvm.fmuladd.f64(double %366, double 2.000000e+00, double 1.000000e+00)
%368 = load double, double* @M_PI
%369 = fmul double %368, 2.000000e+00
%370 = fdiv double %367, %369
%371 = tail call double @sqrt(double %370)
%372 = mul nsw i64 %363, %363
%373 = mul i64 %363, %363
%374 = add i64 %373, %363
%375 = and i64 %374, 4294967295
%376 = getelementptr inbounds double, double* %13, i64 %375
%377 = load double, double* %376
%378 = fmul double %371, %377
store double %378, double* %376
%379 = icmp eq i64 %363, 0
br i1 %379, label %399, label %380
380: ; preds = %362, %380
%381 = phi i64 [ %397, %380 ], [ 1, %362 ]
%382 = sub nsw i64 %363, %381
%383 = getelementptr inbounds double, double* %3, i64 %382
%384 = load double, double* %383
%385 = add nuw nsw i64 %381, %363
%386 = and i64 %385, 4294967295
%387 = getelementptr inbounds double, double* %3, i64 %386
%388 = load double, double* %387
%389 = fdiv double %384, %388
%390 = tail call double @sqrt(double %389)
%391 = fmul double %371, %390
%392 = add nuw nsw i64 %385, %372
%393 = and i64 %392, 4294967295
%394 = getelementptr inbounds double, double* %13, i64 %393
%395 = load double, double* %394
%396 = fmul double %391, %395
store double %396, double* %394
%397 = add nuw nsw i64 %381, 1
%398 = icmp eq i64 %397, %364
br i1 %398, label %399, label %380
399: ; preds = %380, %362
%400 = add nuw nsw i64 %363, 1
%401 = add nuw nsw i64 %364, 1
%402 = icmp eq i64 %400, %361
br i1 %402, label %403, label %362
403: ; preds = %399
%404 = icmp sgt i32 %0, 0
br i1 %404, label %405, label %523
405: ; preds = %403
%406 = fcmp oeq double %2, 0.000000e+00
br i1 %406, label %407, label %490
407: ; preds = %405
br i1 %359, label %602, label %408
408: ; preds = %407
%409 = zext i32 %0 to i64
%410 = shl nuw nsw i64 %409, 3
%411 = add nuw nsw i64 %410, 8
%412 = zext i32 %10 to i64
%413 = icmp ult i32 %10, 4
br i1 %413, label %483, label %414
414: ; preds = %408
%415 = and i64 %361, 4294967292
%416 = add nsw i64 %415, -4
%417 = lshr exact i64 %416, 2
%418 = add nuw nsw i64 %417, 1
%419 = and i64 %418, 7
%420 = icmp ult i64 %416, 28
br i1 %420, label %468, label %421
421: ; preds = %414
%422 = and i64 %418, 9223372036854775800
br label %423
423: ; preds = %423, %421
%424 = phi i64 [ 0, %421 ], [ %465, %423 ]
%425 = phi i64 [ 0, %421 ], [ %466, %423 ]
%426 = getelementptr inbounds double, double* %9, i64 %424
%427 = bitcast double* %426 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %427
%428 = getelementptr inbounds double, double* %426, i64 2
%429 = bitcast double* %428 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %429
%430 = or i64 %424, 4
%431 = getelementptr inbounds double, double* %9, i64 %430
%432 = bitcast double* %431 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %432
%433 = getelementptr inbounds double, double* %431, i64 2
%434 = bitcast double* %433 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %434
%435 = or i64 %424, 8
%436 = getelementptr inbounds double, double* %9, i64 %435
%437 = bitcast double* %436 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %437
%438 = getelementptr inbounds double, double* %436, i64 2
%439 = bitcast double* %438 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %439
%440 = or i64 %424, 12
%441 = getelementptr inbounds double, double* %9, i64 %440
%442 = bitcast double* %441 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %442
%443 = getelementptr inbounds double, double* %441, i64 2
%444 = bitcast double* %443 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %444
%445 = or i64 %424, 16
%446 = getelementptr inbounds double, double* %9, i64 %445
%447 = bitcast double* %446 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %447
%448 = getelementptr inbounds double, double* %446, i64 2
%449 = bitcast double* %448 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %449
%450 = or i64 %424, 20
%451 = getelementptr inbounds double, double* %9, i64 %450
%452 = bitcast double* %451 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %452
%453 = getelementptr inbounds double, double* %451, i64 2
%454 = bitcast double* %453 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %454
%455 = or i64 %424, 24
%456 = getelementptr inbounds double, double* %9, i64 %455
%457 = bitcast double* %456 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %457
%458 = getelementptr inbounds double, double* %456, i64 2
%459 = bitcast double* %458 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %459
%460 = or i64 %424, 28
%461 = getelementptr inbounds double, double* %9, i64 %460
%462 = bitcast double* %461 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %462
%463 = getelementptr inbounds double, double* %461, i64 2
%464 = bitcast double* %463 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %464
%465 = add nuw i64 %424, 32
%466 = add i64 %425, 8
%467 = icmp eq i64 %466, %422
br i1 %467, label %468, label %423
468: ; preds = %423, %414
%469 = phi i64 [ 0, %414 ], [ %465, %423 ]
%470 = icmp eq i64 %419, 0
br i1 %470, label %481, label %471
471: ; preds = %468, %471
%472 = phi i64 [ %478, %471 ], [ %469, %468 ]
%473 = phi i64 [ %479, %471 ], [ 0, %468 ]
%474 = getelementptr inbounds double, double* %9, i64 %472
%475 = bitcast double* %474 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %475
%476 = getelementptr inbounds double, double* %474, i64 2
%477 = bitcast double* %476 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %477
%478 = add nuw i64 %472, 4
%479 = add i64 %473, 1
%480 = icmp eq i64 %479, %419
br i1 %480, label %481, label %471
481: ; preds = %471, %468
%482 = icmp eq i64 %415, %361
br i1 %482, label %523, label %483
483: ; preds = %408, %481
%484 = phi i64 [ 0, %408 ], [ %415, %481 ]
br label %485
485: ; preds = %483, %485
%486 = phi i64 [ %488, %485 ], [ %484, %483 ]
%487 = getelementptr inbounds double, double* %9, i64 %486
store double 1.000000e+00, double* %487
%488 = add nuw nsw i64 %486, 1
%489 = icmp eq i64 %488, %412
br i1 %489, label %523, label %485
490: ; preds = %405
%491 = tail call double @sin(double %2)
%492 = getelementptr inbounds double, double* %7, i64 1
store double %491, double* %492
%493 = tail call double @cos(double %2)
%494 = getelementptr inbounds double, double* %9, i64 1
store double %493, double* %494
%495 = icmp ult i32 %0, 2
br i1 %495, label %523, label %496
496: ; preds = %490
%497 = zext i32 %10 to i64
%498 = fmul double %493, %491
%499 = tail call double @llvm.fmuladd.f64(double %491, double %493, double %498)
%500 = getelementptr inbounds double, double* %7, i64 2
store double %499, double* %500
%501 = load double, double* %492
%502 = fneg double %501
%503 = fmul double %491, %502
%504 = tail call double @llvm.fmuladd.f64(double %493, double %493, double %503)
%505 = getelementptr inbounds double, double* %9, i64 2
store double %504, double* %505
%506 = icmp eq i32 %10, 3
br i1 %506, label %523, label %507
507: ; preds = %496, %507
%508 = phi i64 [ %521, %507 ], [ 3, %496 ]
%509 = phi double [ %519, %507 ], [ %504, %496 ]
%510 = phi double [ %516, %507 ], [ %501, %496 ]
%511 = phi double [ %514, %507 ], [ %499, %496 ]
%512 = load double, double* %494
%513 = fmul double %512, %511
%514 = tail call double @llvm.fmuladd.f64(double %510, double %509, double %513)
%515 = getelementptr inbounds double, double* %7, i64 %508
store double %514, double* %515
%516 = load double, double* %492
%517 = fneg double %516
%518 = fmul double %511, %517
%519 = tail call double @llvm.fmuladd.f64(double %512, double %509, double %518)
%520 = getelementptr inbounds double, double* %9, i64 %508
store double %519, double* %520
%521 = add nuw nsw i64 %508, 1
%522 = icmp eq i64 %521, %497
br i1 %522, label %523, label %507
523: ; preds = %507, %485, %496, %481, %490, %403
br i1 %359, label %602, label %524
524: ; preds = %523
%525 = zext i32 %10 to i64
br label %526
526: ; preds = %524, %599
%527 = phi i64 [ 0, %524 ], [ %600, %599 ]
%528 = trunc i64 %527 to i32
%529 = mul nsw i32 %528, %528
%530 = zext i32 %529 to i64
%531 = add nuw nsw i64 %527, %530
%532 = getelementptr inbounds double, double* %13, i64 %531
%533 = load double, double* %532
%534 = fdiv double %533, 0x3FF6A09E667F3BCD
%535 = getelementptr inbounds double, double* %12, i64 %531
store double %534, double* %535
%536 = icmp eq i64 %527, 0
br i1 %536, label %599, label %537
537: ; preds = %526
%538 = add nuw i32 %529, %528
%539 = and i64 %527, 1
%540 = icmp eq i64 %527, 1
br i1 %540, label %580, label %541
541: ; preds = %537
%542 = and i64 %527, 9223372036854775806
br label %543
543: ; preds = %543, %541
%544 = phi i64 [ 1, %541 ], [ %577, %543 ]
%545 = phi i64 [ 0, %541 ], [ %578, %543 ]
%546 = trunc i64 %544 to i32
%547 = add i32 %538, %546
%548 = zext i32 %547 to i64
%549 = getelementptr inbounds double, double* %13, i64 %548
%550 = load double, double* %549
%551 = getelementptr inbounds double, double* %7, i64 %544
%552 = load double, double* %551
%553 = fmul double %550, %552
%554 = sub i32 %538, %546
%555 = sext i32 %554 to i64
%556 = getelementptr inbounds double, double* %12, i64 %555
store double %553, double* %556
%557 = getelementptr inbounds double, double* %9, i64 %544
%558 = load double, double* %557
%559 = fmul double %550, %558
%560 = getelementptr inbounds double, double* %12, i64 %548
store double %559, double* %560
%561 = add nuw nsw i64 %544, 1
%562 = trunc i64 %561 to i32
%563 = add i32 %538, %562
%564 = zext i32 %563 to i64
%565 = getelementptr inbounds double, double* %13, i64 %564
%566 = load double, double* %565
%567 = getelementptr inbounds double, double* %7, i64 %561
%568 = load double, double* %567
%569 = fmul double %566, %568
%570 = sub i32 %538, %562
%571 = sext i32 %570 to i64
%572 = getelementptr inbounds double, double* %12, i64 %571
store double %569, double* %572
%573 = getelementptr inbounds double, double* %9, i64 %561
%574 = load double, double* %573
%575 = fmul double %566, %574
%576 = getelementptr inbounds double, double* %12, i64 %564
store double %575, double* %576
%577 = add nuw nsw i64 %544, 2
%578 = add i64 %545, 2
%579 = icmp eq i64 %578, %542
br i1 %579, label %580, label %543
580: ; preds = %543, %537
%581 = phi i64 [ 1, %537 ], [ %577, %543 ]
%582 = icmp eq i64 %539, 0
br i1 %582, label %599, label %583
583: ; preds = %580
%584 = trunc i64 %581 to i32
%585 = add i32 %538, %584
%586 = zext i32 %585 to i64
%587 = getelementptr inbounds double, double* %13, i64 %586
%588 = load double, double* %587
%589 = getelementptr inbounds double, double* %7, i64 %581
%590 = load double, double* %589
%591 = fmul double %588, %590
%592 = sub i32 %538, %584
%593 = sext i32 %592 to i64
%594 = getelementptr inbounds double, double* %12, i64 %593
store double %591, double* %594
%595 = getelementptr inbounds double, double* %9, i64 %581
%596 = load double, double* %595
%597 = fmul double %588, %596
%598 = getelementptr inbounds double, double* %12, i64 %586
store double %597, double* %598
br label %599
599: ; preds = %583, %580, %526
%600 = add nuw nsw i64 %527, 1
%601 = icmp eq i64 %600, %525
br i1 %601, label %602, label %526
602: ; preds = %599, %358, %407, %121, %16, %523
%603 = tail call i32 @free(double* %13)
%604 = tail call i32 @free(double* %7)
%605 = tail call i32 @free(double* %9)
ret double* %12
}