double inspherefast(double *pa, double *pb, double *pc, double *pd,
                    double *pe) {
  double aex, bex, cex, dex;
  double aey, bey, cey, dey;
  double aez, bez, cez, dez;
  double alift, blift, clift, dlift;
  double ab, bc, cd, da, ac, bd;
  double abc, bcd, cda, dab;
  aex = pa[0] - pe[0];
  bex = pb[0] - pe[0];
  cex = pc[0] - pe[0];
  dex = pd[0] - pe[0];
  aey = pa[1] - pe[1];
  bey = pb[1] - pe[1];
  cey = pc[1] - pe[1];
  dey = pd[1] - pe[1];
  aez = pa[2] - pe[2];
  bez = pb[2] - pe[2];
  cez = pc[2] - pe[2];
  dez = pd[2] - pe[2];
  ab = aex * bey - bex * aey;
  bc = bex * cey - cex * bey;
  cd = cex * dey - dex * cey;
  da = dex * aey - aex * dey;
  ac = aex * cey - cex * aey;
  bd = bex * dey - dex * bey;
  abc = aez * bc - bez * ac + cez * ab;
  bcd = bez * cd - cez * bd + dez * bc;
  cda = cez * da + dez * ac + aez * cd;
  dab = dez * ab + aez * bd + bez * da;
  alift = aex * aex + aey * aey + aez * aez;
  blift = bex * bex + bey * bey + bez * bez;
  clift = cex * cex + cey * cey + cez * cez;
  dlift = dex * dex + dey * dey + dez * dez;
  return (dlift * abc - clift * dab) + (blift * cda - alift * bcd);
}


x86 assembly (clang, -O0):
<inspherefast>:
push   %rbp
mov    %rsp,%rbp
sub    $0x78,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x28(%rbp),%rax
subsd  (%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x28(%rbp),%rax
subsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x28(%rbp),%rax
subsd  (%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x28(%rbp),%rax
subsd  (%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x28(%rbp),%rax
subsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x28(%rbp),%rax
subsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x28(%rbp),%rax
subsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x28(%rbp),%rax
subsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x28(%rbp),%rax
subsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x28(%rbp),%rax
subsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x28(%rbp),%rax
subsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x20(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x28(%rbp),%rax
subsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  -0x58(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x40(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  -0x58(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movsd  -0x70(%rbp),%xmm1
movsd  -0xb8(%rbp),%xmm2
ff
movsd  -0x78(%rbp),%xmm0
mulsd  -0xd0(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x80(%rbp),%xmm0
movsd  -0xb0(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  -0x78(%rbp),%xmm1
movsd  -0xc0(%rbp),%xmm2
ff
movsd  -0x80(%rbp),%xmm0
mulsd  -0xd8(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0xb8(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
movsd  -0x80(%rbp),%xmm1
movsd  -0xc8(%rbp),%xmm2
ff
movsd  -0x88(%rbp),%xmm0
ff
mulsd  -0xd0(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x70(%rbp),%xmm0
movsd  -0xc0(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  -0x88(%rbp),%xmm1
ff
movsd  -0xb0(%rbp),%xmm2
ff
movsd  -0x70(%rbp),%xmm0
mulsd  -0xd8(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x78(%rbp),%xmm0
movsd  -0xc8(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
movsd  -0x30(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x70(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x38(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x78(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x40(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm2
movsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x80(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x48(%rbp),%xmm1
movsd  -0x48(%rbp),%xmm2
movsd  -0x68(%rbp),%xmm0
mulsd  -0x68(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0x88(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
movsd  -0xe0(%rbp),%xmm2
ff
movsd  -0xa0(%rbp),%xmm1
ff
mulsd  -0xf8(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x98(%rbp),%xmm1
ff
movsd  -0xf0(%rbp),%xmm3
ff
movsd  -0x90(%rbp),%xmm2
ff
mulsd  -0xe8(%rbp),%xmm2
ff
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
add    $0x78,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<inspherefast>:
movsd  (%rdi),%xmm12
movsd  0x8(%rdi),%xmm10
movhpd (%rdx),%xmm12
movsd  (%r8),%xmm0
movsd  0x8(%r8),%xmm2
movsd  0x10(%r8),%xmm4
unpcklpd %xmm0,%xmm0
subpd  %xmm0,%xmm12
movsd  (%rcx),%xmm11
movsd  0x8(%rcx),%xmm8
movsd  0x10(%rcx),%xmm1
movhpd (%rsi),%xmm11
movhpd 0x8(%rdx),%xmm10
subpd  %xmm0,%xmm11
unpcklpd %xmm2,%xmm2
movhpd 0x8(%rsi),%xmm8
subpd  %xmm2,%xmm10
movsd  0x10(%rdx),%xmm5
movhpd 0x10(%rdi),%xmm5
subpd  %xmm2,%xmm8
unpcklpd %xmm4,%xmm4
movhpd 0x10(%rsi),%xmm1
subpd  %xmm4,%xmm5
subpd  %xmm4,%xmm1
movapd %xmm11,%xmm13
unpckhpd %xmm11,%xmm13
movapd %xmm11,%xmm4
shufpd $0x1,%xmm11,%xmm4
movapd 0xe69(%rip),%xmm9
movapd %xmm10,%xmm0
mulpd  %xmm4,%xmm0
xorpd  %xmm9,%xmm4
movapd %xmm8,%xmm3
shufpd $0x1,%xmm8,%xmm3
movapd %xmm12,%xmm2
mulpd  %xmm3,%xmm2
subpd  %xmm0,%xmm2
movapd %xmm12,%xmm7
shufpd $0x1,%xmm12,%xmm7
mulpd  %xmm7,%xmm3
xorpd  %xmm9,%xmm7
movapd %xmm10,%xmm0
unpckhpd %xmm10,%xmm0
movapd %xmm11,%xmm6
mulpd  %xmm10,%xmm6
shufpd $0x1,%xmm6,%xmm6
subpd  %xmm3,%xmm6
mulsd  %xmm10,%xmm7
mulsd  %xmm12,%xmm0
addsd  %xmm7,%xmm0
mulpd  %xmm8,%xmm4
unpckhpd %xmm4,%xmm4
mulsd  %xmm8,%xmm13
addsd  %xmm4,%xmm13
movapd %xmm1,%xmm3
unpckhpd %xmm1,%xmm3
xorpd  %xmm9,%xmm3
xorpd  %xmm5,%xmm9
unpcklpd %xmm0,%xmm0
unpcklpd %xmm1,%xmm3
mulpd  %xmm0,%xmm3
movapd %xmm5,%xmm0
unpckhpd %xmm5,%xmm0
unpcklpd %xmm5,%xmm0
mulpd  %xmm6,%xmm0
addpd  %xmm3,%xmm0
movapd %xmm5,%xmm4
mulpd  %xmm2,%xmm4
addpd  %xmm0,%xmm4
unpcklpd %xmm13,%xmm13
movapd %xmm5,%xmm0
shufpd $0x1,%xmm9,%xmm0
mulpd  %xmm13,%xmm0
mulpd  %xmm1,%xmm2
addpd  %xmm0,%xmm2
mulpd  %xmm1,%xmm6
shufpd $0x1,%xmm6,%xmm6
addpd  %xmm2,%xmm6
mulpd  %xmm10,%xmm10
mulpd  %xmm12,%xmm12
addpd  %xmm10,%xmm12
shufpd $0x1,%xmm12,%xmm12
mulpd  %xmm5,%xmm5
addpd  %xmm12,%xmm5
mulpd  %xmm6,%xmm5
mulpd  %xmm8,%xmm8
mulpd  %xmm11,%xmm11
addpd  %xmm8,%xmm11
mulpd  %xmm1,%xmm1
addpd  %xmm11,%xmm1
mulpd  %xmm4,%xmm1
subpd  %xmm5,%xmm1
movapd %xmm1,%xmm0
unpckhpd %xmm1,%xmm0
addsd  %xmm1,%xmm0
ret

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

LLVM IR (-O2):
define double @inspherefast(double* nocapture %0, double* nocapture %1, double* nocapture %2, double* nocapture %3, double* nocapture %4) {
%6 = load double, double* %0
%7 = load double, double* %4
%8 = load double, double* %1
%9 = load double, double* %2
%10 = insertelement <2 x double> poison, double %6, i64 0
%11 = insertelement <2 x double> %10, double %9, i64 1
%12 = insertelement <2 x double> poison, double %7, i64 0
%13 = shufflevector <2 x double> %12, <2 x double> poison, <2 x i32> zeroinitializer
%14 = fsub <2 x double> %11, %13
%15 = load double, double* %3
%16 = insertelement <2 x double> poison, double %15, i64 0
%17 = insertelement <2 x double> %16, double %8, i64 1
%18 = fsub <2 x double> %17, %13
%19 = getelementptr inbounds double, double* %0, i64 1
%20 = load double, double* %19
%21 = getelementptr inbounds double, double* %4, i64 1
%22 = load double, double* %21
%23 = getelementptr inbounds double, double* %1, i64 1
%24 = load double, double* %23
%25 = getelementptr inbounds double, double* %2, i64 1
%26 = load double, double* %25
%27 = insertelement <2 x double> poison, double %20, i64 0
%28 = insertelement <2 x double> %27, double %26, i64 1
%29 = insertelement <2 x double> poison, double %22, i64 0
%30 = shufflevector <2 x double> %29, <2 x double> poison, <2 x i32> zeroinitializer
%31 = fsub <2 x double> %28, %30
%32 = getelementptr inbounds double, double* %3, i64 1
%33 = load double, double* %32
%34 = insertelement <2 x double> poison, double %33, i64 0
%35 = insertelement <2 x double> %34, double %24, i64 1
%36 = fsub <2 x double> %35, %30
%37 = getelementptr inbounds double, double* %0, i64 2
%38 = load double, double* %37
%39 = getelementptr inbounds double, double* %4, i64 2
%40 = load double, double* %39
%41 = getelementptr inbounds double, double* %1, i64 2
%42 = load double, double* %41
%43 = getelementptr inbounds double, double* %2, i64 2
%44 = load double, double* %43
%45 = insertelement <2 x double> poison, double %44, i64 0
%46 = insertelement <2 x double> %45, double %38, i64 1
%47 = insertelement <2 x double> poison, double %40, i64 0
%48 = shufflevector <2 x double> %47, <2 x double> poison, <2 x i32> zeroinitializer
%49 = fsub <2 x double> %46, %48
%50 = getelementptr inbounds double, double* %3, i64 2
%51 = load double, double* %50
%52 = insertelement <2 x double> poison, double %51, i64 0
%53 = insertelement <2 x double> %52, double %42, i64 1
%54 = fsub <2 x double> %53, %48
%55 = extractelement <2 x double> %18, i64 1
%56 = shufflevector <2 x double> %18, <2 x double> undef, <2 x i32> <i32 1, i32 0>
%57 = fneg <2 x double> %56
%58 = fmul <2 x double> %31, %57
%59 = extractelement <2 x double> %36, i64 0
%60 = shufflevector <2 x double> %36, <2 x double> undef, <2 x i32> <i32 1, i32 0>
%61 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %14, <2 x double> %60, <2 x double> %58)
%62 = extractelement <2 x double> %14, i64 0
%63 = shufflevector <2 x double> %14, <2 x double> undef, <2 x i32> <i32 1, i32 0>
%64 = fneg <2 x double> %63
%65 = fmul <2 x double> %60, %64
%66 = extractelement <2 x double> %31, i64 1
%67 = shufflevector <2 x double> %31, <2 x double> undef, <2 x i32> <i32 1, i32 0>
%68 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %56, <2 x double> %67, <2 x double> %65)
%69 = fmul <2 x double> %31, %64
%70 = extractelement <2 x double> %69, i64 0
%71 = tail call double @llvm.fmuladd.f64(double %62, double %66, double %70)
%72 = fmul <2 x double> %36, %57
%73 = extractelement <2 x double> %72, i64 1
%74 = tail call double @llvm.fmuladd.f64(double %55, double %59, double %73)
%75 = extractelement <2 x double> %54, i64 1
%76 = fneg double %75
%77 = extractelement <2 x double> %49, i64 0
%78 = fneg double %77
%79 = insertelement <2 x double> poison, double %71, i64 0
%80 = shufflevector <2 x double> %79, <2 x double> poison, <2 x i32> zeroinitializer
%81 = insertelement <2 x double> poison, double %76, i64 0
%82 = shufflevector <2 x double> %81, <2 x double> %54, <2 x i32> <i32 0, i32 2>
%83 = fmul <2 x double> %80, %82
%84 = shufflevector <2 x double> %49, <2 x double> undef, <2 x i32> <i32 1, i32 undef>
%85 = shufflevector <2 x double> %84, <2 x double> %49, <2 x i32> <i32 0, i32 2>
%86 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %85, <2 x double> %68, <2 x double> %83)
%87 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %49, <2 x double> %61, <2 x double> %86)
%88 = insertelement <2 x double> poison, double %74, i64 0
%89 = shufflevector <2 x double> %88, <2 x double> poison, <2 x i32> zeroinitializer
%90 = insertelement <2 x double> %84, double %78, i64 1
%91 = fmul <2 x double> %89, %90
%92 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %54, <2 x double> %61, <2 x double> %91)
%93 = shufflevector <2 x double> %54, <2 x double> undef, <2 x i32> <i32 1, i32 0>
%94 = shufflevector <2 x double> %68, <2 x double> undef, <2 x i32> <i32 1, i32 0>
%95 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %93, <2 x double> %94, <2 x double> %92)
%96 = fmul <2 x double> %31, %31
%97 = shufflevector <2 x double> %96, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%98 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %63, <2 x double> %63, <2 x double> %97)
%99 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %49, <2 x double> %49, <2 x double> %98)
%100 = fmul <2 x double> %36, %36
%101 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %18, <2 x double> %18, <2 x double> %100)
%102 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %54, <2 x double> %54, <2 x double> %101)
%103 = fneg <2 x double> %99
%104 = fmul <2 x double> %95, %103
%105 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %102, <2 x double> %87, <2 x double> %104)
%106 = shufflevector <2 x double> %105, <2 x double> poison, <2 x i32> <i32 1, i32 undef>
%107 = fadd <2 x double> %105, %106
%108 = extractelement <2 x double> %107, i64 0
ret double %108
}