void harmonicDiffuse(float *C, const float *Y, const float *D) {
  C[0] = (float)(0.429043 * Y[0 * 9 + 8] * (D[0] * D[0] - D[1] * D[1]) +
                 0.743125 * Y[0 * 9 + 7] * D[2] * D[2] +
                 0.886227 * Y[0 * 9 + 0] - 0.247708 * Y[0 * 9 + 7] +
                 0.429043 * 2 *
                     (Y[0 * 9 + 4] * D[0] * D[1] + Y[0 * 9 + 5] * D[1] * D[2] +
                      Y[0 * 9 + 6] * D[2] * D[0]) +
                 0.511664 * 2 *
                     (Y[0 * 9 + 1] * D[0] + Y[0 * 9 + 2] * D[1] +
                      Y[0 * 9 + 3] * D[2]));
  C[1] = (float)(0.429043 * Y[1 * 9 + 8] * (D[0] * D[0] - D[1] * D[1]) +
                 0.743125 * Y[1 * 9 + 7] * D[2] * D[2] +
                 0.886227 * Y[1 * 9 + 0] - 0.247708 * Y[1 * 9 + 7] +
                 0.429043 * 2 *
                     (Y[1 * 9 + 4] * D[0] * D[1] + Y[1 * 9 + 5] * D[1] * D[2] +
                      Y[1 * 9 + 6] * D[2] * D[0]) +
                 0.511664 * 2 *
                     (Y[1 * 9 + 1] * D[0] + Y[1 * 9 + 2] * D[1] +
                      Y[1 * 9 + 3] * D[2]));
  C[2] = (float)(0.429043 * Y[2 * 9 + 8] * (D[0] * D[0] - D[1] * D[1]) +
                 0.743125 * Y[2 * 9 + 7] * D[2] * D[2] +
                 0.886227 * Y[2 * 9 + 0] - 0.247708 * Y[2 * 9 + 7] +
                 0.429043 * 2 *
                     (Y[2 * 9 + 4] * D[0] * D[1] + Y[2 * 9 + 5] * D[1] * D[2] +
                      Y[2 * 9 + 6] * D[2] * D[0]) +
                 0.511664 * 2 *
                     (Y[2 * 9 + 1] * D[0] + Y[2 * 9 + 2] * D[1] +
                      Y[2 * 9 + 3] * D[2]));
}


x86 assembly (clang, -O0):
<harmonicDiffuse>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
movss  0x20(%rax),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xf03(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm1
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm3
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm2
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm2
movd   %xmm2,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm2
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
cvtss2sd %xmm1,%xmm2
mov    -0x10(%rbp),%rax
movss  0x1c(%rax),%xmm1
cvtss2sd %xmm1,%xmm3
movsd  0xea7(%rip),%xmm1
mulsd  %xmm3,%xmm1
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm3,%xmm1
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xe5d(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movss  0x1c(%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe38(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  0x10(%rax),%xmm2
mov    -0x18(%rbp),%rax
mulss  (%rax),%xmm2
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm3
mov    -0x10(%rbp),%rax
movss  0x14(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulss  0x8(%rax),%xmm1
mulss  %xmm3,%xmm2
addss  %xmm1,%xmm2
mov    -0x10(%rbp),%rax
movss  0x18(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulss  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm3
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xdbd(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm2
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm3
mov    -0x10(%rbp),%rax
movss  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm0
mulss  %xmm3,%xmm2
addss  %xmm0,%xmm2
mov    -0x10(%rbp),%rax
movss  0xc(%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm3
mulss  %xmm3,%xmm0
addss  %xmm2,%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xd5c(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
movss  0x44(%rax),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd5b(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm1
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm3
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm2
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm2
movd   %xmm2,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm2
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
cvtss2sd %xmm1,%xmm2
mov    -0x10(%rbp),%rax
movss  0x40(%rax),%xmm1
cvtss2sd %xmm1,%xmm3
movsd  0xcff(%rip),%xmm1
mulsd  %xmm3,%xmm1
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm3,%xmm1
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  0x24(%rax),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xcb4(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movss  0x40(%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xc8f(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  0x34(%rax),%xmm2
mov    -0x18(%rbp),%rax
mulss  (%rax),%xmm2
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm3
mov    -0x10(%rbp),%rax
movss  0x38(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulss  0x8(%rax),%xmm1
mulss  %xmm3,%xmm2
addss  %xmm1,%xmm2
mov    -0x10(%rbp),%rax
movss  0x3c(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulss  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm3
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xc14(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movss  0x28(%rax),%xmm2
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm3
mov    -0x10(%rbp),%rax
movss  0x2c(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm0
mulss  %xmm3,%xmm2
addss  %xmm0,%xmm2
mov    -0x10(%rbp),%rax
movss  0x30(%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm3
mulss  %xmm3,%xmm0
addss  %xmm2,%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xbb3(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x4(%rax)
mov    -0x10(%rbp),%rax
movss  0x68(%rax),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xbb1(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm1
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm3
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm2
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm2
movd   %xmm2,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm2
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
cvtss2sd %xmm1,%xmm2
mov    -0x10(%rbp),%rax
movss  0x64(%rax),%xmm1
cvtss2sd %xmm1,%xmm3
movsd  0xb55(%rip),%xmm1
mulsd  %xmm3,%xmm1
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm3,%xmm1
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  0x48(%rax),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xb0a(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movss  0x64(%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xae5(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  0x58(%rax),%xmm2
mov    -0x18(%rbp),%rax
mulss  (%rax),%xmm2
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm3
mov    -0x10(%rbp),%rax
movss  0x5c(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulss  0x8(%rax),%xmm1
mulss  %xmm3,%xmm2
addss  %xmm1,%xmm2
mov    -0x10(%rbp),%rax
movss  0x60(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulss  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm3
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xa6a(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movss  0x4c(%rax),%xmm2
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm3
mov    -0x10(%rbp),%rax
movss  0x50(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm0
mulss  %xmm3,%xmm2
addss  %xmm0,%xmm2
mov    -0x10(%rbp),%rax
movss  0x54(%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm3
mulss  %xmm3,%xmm0
addss  %xmm2,%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xa09(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x8(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<harmonicDiffuse>:
movss  0x20(%rsi),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xeee(%rip),%xmm8
mulsd  %xmm8,%xmm0
movss  (%rdx),%xmm5
movss  0x4(%rdx),%xmm6
movaps %xmm6,%xmm1
mulss  %xmm6,%xmm1
movaps %xmm5,%xmm2
mulss  %xmm5,%xmm2
subss  %xmm1,%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm0,%xmm2
movss  (%rsi),%xmm9
movss  0x1c(%rsi),%xmm0
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
movsd  0xeb5(%rip),%xmm0
movaps %xmm1,%xmm4
mulsd  %xmm0,%xmm4
movss  0x8(%rdx),%xmm7
cvtss2sd %xmm7,%xmm3
mulsd  %xmm3,%xmm4
mulsd  %xmm3,%xmm4
addsd  %xmm2,%xmm4
xorps  %xmm2,%xmm2
cvtss2sd %xmm9,%xmm2
movsd  0xe90(%rip),%xmm9
mulsd  %xmm9,%xmm2
addsd  %xmm4,%xmm2
movsd  0xe86(%rip),%xmm10
mulsd  %xmm10,%xmm1
addsd  %xmm2,%xmm1
movss  0x10(%rsi),%xmm2
mulss  %xmm5,%xmm2
movss  0x14(%rsi),%xmm4
mulss  %xmm6,%xmm4
mulss  %xmm7,%xmm4
mulss  %xmm6,%xmm2
addss  %xmm4,%xmm2
movss  0x18(%rsi),%xmm4
mulss  %xmm7,%xmm4
mulss  %xmm5,%xmm4
addss  %xmm2,%xmm4
xorps  %xmm2,%xmm2
cvtss2sd %xmm4,%xmm2
movsd  0xe46(%rip),%xmm11
mulsd  %xmm11,%xmm2
mulss  0x8(%rsi),%xmm6
mulss  0x4(%rsi),%xmm5
addsd  %xmm1,%xmm2
addss  %xmm6,%xmm5
mulss  0xc(%rsi),%xmm7
addss  %xmm5,%xmm7
xorps  %xmm1,%xmm1
cvtss2sd %xmm7,%xmm1
movsd  0xe1e(%rip),%xmm12
mulsd  %xmm12,%xmm1
addsd  %xmm2,%xmm1
cvtsd2ss %xmm1,%xmm1
movss  %xmm1,(%rdi)
movss  0x44(%rsi),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm8,%xmm1
movss  (%rdx),%xmm3
movss  0x4(%rdx),%xmm7
movaps %xmm7,%xmm2
mulss  %xmm7,%xmm2
movaps %xmm3,%xmm6
mulss  %xmm3,%xmm6
subss  %xmm2,%xmm6
cvtss2sd %xmm6,%xmm6
mulsd  %xmm1,%xmm6
movss  0x24(%rsi),%xmm13
movss  0x40(%rsi),%xmm1
xorps  %xmm5,%xmm5
cvtss2sd %xmm1,%xmm5
movaps %xmm5,%xmm2
mulsd  %xmm0,%xmm2
movss  0x8(%rdx),%xmm1
xorps  %xmm4,%xmm4
cvtss2sd %xmm1,%xmm4
mulsd  %xmm4,%xmm2
mulsd  %xmm4,%xmm2
addsd  %xmm6,%xmm2
xorps  %xmm4,%xmm4
cvtss2sd %xmm13,%xmm4
mulsd  %xmm9,%xmm4
addsd  %xmm2,%xmm4
mulsd  %xmm10,%xmm5
addsd  %xmm4,%xmm5
movss  0x34(%rsi),%xmm2
mulss  %xmm3,%xmm2
movss  0x38(%rsi),%xmm4
mulss  %xmm7,%xmm4
mulss  %xmm1,%xmm4
mulss  %xmm7,%xmm2
addss  %xmm4,%xmm2
movss  0x3c(%rsi),%xmm4
mulss  %xmm1,%xmm4
mulss  %xmm3,%xmm4
addss  %xmm2,%xmm4
xorps  %xmm2,%xmm2
cvtss2sd %xmm4,%xmm2
mulsd  %xmm11,%xmm2
addsd  %xmm5,%xmm2
mulss  0x2c(%rsi),%xmm7
mulss  0x28(%rsi),%xmm3
mulss  0x30(%rsi),%xmm1
addss  %xmm7,%xmm3
addss  %xmm3,%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm12,%xmm1
addsd  %xmm2,%xmm1
cvtsd2ss %xmm1,%xmm1
movss  %xmm1,0x4(%rdi)
movss  0x68(%rsi),%xmm1
xorps  %xmm2,%xmm2
cvtss2sd %xmm1,%xmm2
mulsd  %xmm8,%xmm2
movss  (%rdx),%xmm1
movss  0x4(%rdx),%xmm6
movaps %xmm6,%xmm3
mulss  %xmm6,%xmm3
movaps %xmm1,%xmm4
mulss  %xmm1,%xmm4
subss  %xmm3,%xmm4
xorps  %xmm3,%xmm3
cvtss2sd %xmm4,%xmm3
mulsd  %xmm2,%xmm3
movss  0x48(%rsi),%xmm4
movss  0x64(%rsi),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm2,%xmm0
movss  0x8(%rdx),%xmm7
xorps  %xmm5,%xmm5
cvtss2sd %xmm7,%xmm5
mulsd  %xmm5,%xmm0
mulsd  %xmm5,%xmm0
addsd  %xmm3,%xmm0
xorps  %xmm3,%xmm3
cvtss2sd %xmm4,%xmm3
mulsd  %xmm9,%xmm3
addsd  %xmm0,%xmm3
mulsd  %xmm10,%xmm2
addsd  %xmm3,%xmm2
movss  0x58(%rsi),%xmm0
mulss  %xmm1,%xmm0
movss  0x5c(%rsi),%xmm3
mulss  %xmm6,%xmm3
mulss  %xmm7,%xmm3
mulss  %xmm6,%xmm0
addss  %xmm3,%xmm0
movss  0x60(%rsi),%xmm3
mulss  %xmm7,%xmm3
mulss  %xmm1,%xmm3
addss  %xmm0,%xmm3
xorps  %xmm0,%xmm0
cvtss2sd %xmm3,%xmm0
mulsd  %xmm11,%xmm0
addsd  %xmm2,%xmm0
mulss  0x50(%rsi),%xmm6
mulss  0x4c(%rsi),%xmm1
mulss  0x54(%rsi),%xmm7
addss  %xmm6,%xmm1
addss  %xmm1,%xmm7
xorps  %xmm1,%xmm1
cvtss2sd %xmm7,%xmm1
mulsd  %xmm12,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm0,0x8(%rdi)
ret

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

LLVM IR (-O2):
define void @harmonicDiffuse(float* nocapture writeonly %0, float* %1, float* nocapture %2) {
%4 = getelementptr inbounds float, float* %1, i64 8
%5 = load float, float* %4
%6 = fpext float %5 to double
%7 = fmul double %6, 4.290430e-01
%8 = load float, float* %2
%9 = getelementptr inbounds float, float* %2, i64 1
%10 = load float, float* %9
%11 = fneg float %10
%12 = fmul float %10, %11
%13 = tail call float @llvm.fmuladd.f32(float %8, float %8, float %12)
%14 = fpext float %13 to double
%15 = getelementptr inbounds float, float* %1, i64 7
%16 = load float, float* %15
%17 = fpext float %16 to double
%18 = fmul double %17, 7.431250e-01
%19 = getelementptr inbounds float, float* %2, i64 2
%20 = load float, float* %19
%21 = fpext float %20 to double
%22 = fmul double %18, %21
%23 = fmul double %22, %21
%24 = tail call double @llvm.fmuladd.f64(double %7, double %14, double %23)
%25 = load float, float* %1
%26 = fpext float %25 to double
%27 = tail call double @llvm.fmuladd.f64(double %26, double 8.862270e-01, double %24)
%28 = tail call double @llvm.fmuladd.f64(double %17, double -2.477080e-01, double %27)
%29 = getelementptr inbounds float, float* %1, i64 4
%30 = load float, float* %29
%31 = fmul float %8, %30
%32 = getelementptr inbounds float, float* %1, i64 5
%33 = load float, float* %32
%34 = fmul float %10, %33
%35 = fmul float %20, %34
%36 = tail call float @llvm.fmuladd.f32(float %31, float %10, float %35)
%37 = getelementptr inbounds float, float* %1, i64 6
%38 = load float, float* %37
%39 = fmul float %20, %38
%40 = tail call float @llvm.fmuladd.f32(float %39, float %8, float %36)
%41 = fpext float %40 to double
%42 = tail call double @llvm.fmuladd.f64(double %41, double 8.580860e-01, double %28)
%43 = getelementptr inbounds float, float* %1, i64 1
%44 = load float, float* %43
%45 = getelementptr inbounds float, float* %1, i64 2
%46 = load float, float* %45
%47 = fmul float %10, %46
%48 = tail call float @llvm.fmuladd.f32(float %44, float %8, float %47)
%49 = getelementptr inbounds float, float* %1, i64 3
%50 = load float, float* %49
%51 = tail call float @llvm.fmuladd.f32(float %50, float %20, float %48)
%52 = fpext float %51 to double
%53 = tail call double @llvm.fmuladd.f64(double %52, double 0x3FF05F8D2E514C23, double %42)
%54 = fptrunc double %53 to float
store float %54, float* %0
%55 = getelementptr inbounds float, float* %1, i64 17
%56 = load float, float* %55
%57 = fpext float %56 to double
%58 = fmul double %57, 4.290430e-01
%59 = load float, float* %2
%60 = load float, float* %9
%61 = fneg float %60
%62 = fmul float %60, %61
%63 = tail call float @llvm.fmuladd.f32(float %59, float %59, float %62)
%64 = fpext float %63 to double
%65 = getelementptr inbounds float, float* %1, i64 16
%66 = load float, float* %65
%67 = fpext float %66 to double
%68 = fmul double %67, 7.431250e-01
%69 = load float, float* %19
%70 = fpext float %69 to double
%71 = fmul double %68, %70
%72 = fmul double %71, %70
%73 = tail call double @llvm.fmuladd.f64(double %58, double %64, double %72)
%74 = getelementptr inbounds float, float* %1, i64 9
%75 = load float, float* %74
%76 = fpext float %75 to double
%77 = tail call double @llvm.fmuladd.f64(double %76, double 8.862270e-01, double %73)
%78 = tail call double @llvm.fmuladd.f64(double %67, double -2.477080e-01, double %77)
%79 = getelementptr inbounds float, float* %1, i64 13
%80 = load float, float* %79
%81 = fmul float %59, %80
%82 = getelementptr inbounds float, float* %1, i64 14
%83 = load float, float* %82
%84 = fmul float %60, %83
%85 = fmul float %69, %84
%86 = tail call float @llvm.fmuladd.f32(float %81, float %60, float %85)
%87 = getelementptr inbounds float, float* %1, i64 15
%88 = load float, float* %87
%89 = fmul float %69, %88
%90 = tail call float @llvm.fmuladd.f32(float %89, float %59, float %86)
%91 = fpext float %90 to double
%92 = tail call double @llvm.fmuladd.f64(double %91, double 8.580860e-01, double %78)
%93 = getelementptr inbounds float, float* %1, i64 10
%94 = load float, float* %93
%95 = getelementptr inbounds float, float* %1, i64 11
%96 = load float, float* %95
%97 = fmul float %60, %96
%98 = tail call float @llvm.fmuladd.f32(float %94, float %59, float %97)
%99 = getelementptr inbounds float, float* %1, i64 12
%100 = load float, float* %99
%101 = tail call float @llvm.fmuladd.f32(float %100, float %69, float %98)
%102 = fpext float %101 to double
%103 = tail call double @llvm.fmuladd.f64(double %102, double 0x3FF05F8D2E514C23, double %92)
%104 = fptrunc double %103 to float
%105 = getelementptr inbounds float, float* %0, i64 1
store float %104, float* %105
%106 = getelementptr inbounds float, float* %1, i64 26
%107 = load float, float* %106
%108 = fpext float %107 to double
%109 = fmul double %108, 4.290430e-01
%110 = load float, float* %2
%111 = load float, float* %9
%112 = fneg float %111
%113 = fmul float %111, %112
%114 = tail call float @llvm.fmuladd.f32(float %110, float %110, float %113)
%115 = fpext float %114 to double
%116 = getelementptr inbounds float, float* %1, i64 25
%117 = load float, float* %116
%118 = fpext float %117 to double
%119 = fmul double %118, 7.431250e-01
%120 = load float, float* %19
%121 = fpext float %120 to double
%122 = fmul double %119, %121
%123 = fmul double %122, %121
%124 = tail call double @llvm.fmuladd.f64(double %109, double %115, double %123)
%125 = getelementptr inbounds float, float* %1, i64 18
%126 = load float, float* %125
%127 = fpext float %126 to double
%128 = tail call double @llvm.fmuladd.f64(double %127, double 8.862270e-01, double %124)
%129 = tail call double @llvm.fmuladd.f64(double %118, double -2.477080e-01, double %128)
%130 = getelementptr inbounds float, float* %1, i64 22
%131 = load float, float* %130
%132 = fmul float %110, %131
%133 = getelementptr inbounds float, float* %1, i64 23
%134 = load float, float* %133
%135 = fmul float %111, %134
%136 = fmul float %120, %135
%137 = tail call float @llvm.fmuladd.f32(float %132, float %111, float %136)
%138 = getelementptr inbounds float, float* %1, i64 24
%139 = load float, float* %138
%140 = fmul float %120, %139
%141 = tail call float @llvm.fmuladd.f32(float %140, float %110, float %137)
%142 = fpext float %141 to double
%143 = tail call double @llvm.fmuladd.f64(double %142, double 8.580860e-01, double %129)
%144 = getelementptr inbounds float, float* %1, i64 19
%145 = load float, float* %144
%146 = getelementptr inbounds float, float* %1, i64 20
%147 = load float, float* %146
%148 = fmul float %111, %147
%149 = tail call float @llvm.fmuladd.f32(float %145, float %110, float %148)
%150 = getelementptr inbounds float, float* %1, i64 21
%151 = load float, float* %150
%152 = tail call float @llvm.fmuladd.f32(float %151, float %120, float %149)
%153 = fpext float %152 to double
%154 = tail call double @llvm.fmuladd.f64(double %153, double 0x3FF05F8D2E514C23, double %143)
%155 = fptrunc double %154 to float
%156 = getelementptr inbounds float, float* %0, i64 2
store float %155, float* %156
ret void
}