void interpolate_order2_square(double epsilon, double eta,
                               double *out_5124881133433324070) {
  out_5124881133433324070[0] =
      0.25 * epsilon * eta * (epsilon - 1.0) * (eta - 1.0);
  out_5124881133433324070[1] =
      -0.5 * eta * (epsilon - 1.0) * (epsilon + 1.0) * (eta - 1.0);
  out_5124881133433324070[2] =
      0.25 * epsilon * eta * (epsilon + 1.0) * (eta - 1.0);
  out_5124881133433324070[3] =
      -0.5 * epsilon * (epsilon - 1.0) * (eta - 1.0) * (eta + 1.0);
  out_5124881133433324070[4] =
      1.0 * (epsilon - 1.0) * (epsilon + 1.0) * (eta - 1.0) * (eta + 1.0);
  out_5124881133433324070[5] =
      -0.5 * epsilon * (epsilon + 1.0) * (eta - 1.0) * (eta + 1.0);
  out_5124881133433324070[6] =
      0.25 * epsilon * eta * (epsilon - 1.0) * (eta + 1.0);
  out_5124881133433324070[7] =
      -0.5 * eta * (epsilon - 1.0) * (epsilon + 1.0) * (eta + 1.0);
  out_5124881133433324070[8] =
      0.25 * epsilon * eta * (epsilon + 1.0) * (eta + 1.0);
}


x86 assembly (clang, -O0):
<interpolate_order2_square>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
mov    %rdi,-0x18(%rbp)
movsd  0xeee(%rip),%xmm0
mulsd  -0x8(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  -0x8(%rbp),%xmm1
movsd  0xecf(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x10(%rbp),%xmm1
movsd  0xeba(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  0xeb2(%rip),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  -0x8(%rbp),%xmm1
movsd  0xe90(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xe80(%rip),%xmm1
addsd  -0x8(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x10(%rbp),%xmm1
movsd  0xe6a(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  0xe59(%rip),%xmm0
mulsd  -0x8(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  0xe3f(%rip),%xmm1
addsd  -0x8(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x10(%rbp),%xmm1
movsd  0xe29(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x10(%rax)
movsd  0xe20(%rip),%xmm0
mulsd  -0x8(%rbp),%xmm0
movsd  -0x8(%rbp),%xmm1
movsd  0xdfe(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x10(%rbp),%xmm1
movsd  0xde9(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xdd9(%rip),%xmm1
addsd  -0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x18(%rax)
movsd  -0x8(%rbp),%xmm1
movsd  0xdba(%rip),%xmm0
subsd  %xmm0,%xmm1
movsd  0xdae(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  0xda2(%rip),%xmm1
addsd  -0x8(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x10(%rbp),%xmm1
movsd  0xd8c(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xd7c(%rip),%xmm1
addsd  -0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x20(%rax)
movsd  0xd72(%rip),%xmm0
mulsd  -0x8(%rbp),%xmm0
movsd  0xd55(%rip),%xmm1
addsd  -0x8(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x10(%rbp),%xmm1
movsd  0xd3f(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xd2f(%rip),%xmm1
addsd  -0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x28(%rax)
movsd  0xd1d(%rip),%xmm0
mulsd  -0x8(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  -0x8(%rbp),%xmm1
movsd  0xcfe(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xcee(%rip),%xmm1
addsd  -0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x30(%rax)
movsd  0xce4(%rip),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  -0x8(%rbp),%xmm1
movsd  0xcc2(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xcb2(%rip),%xmm1
addsd  -0x8(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xca1(%rip),%xmm1
addsd  -0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x38(%rax)
movsd  0xc8f(%rip),%xmm0
mulsd  -0x8(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  0xc75(%rip),%xmm1
addsd  -0x8(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xc64(%rip),%xmm1
addsd  -0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x40(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<interpolate_order2_square>:
movapd %xmm0,%xmm5
unpcklpd %xmm1,%xmm5
mulpd  0xef0(%rip),%xmm5
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
movapd 0xef0(%rip),%xmm4
addpd  %xmm4,%xmm3
movapd %xmm1,%xmm2
unpcklpd %xmm3,%xmm2
mulpd  %xmm5,%xmm2
mulsd  0xee8(%rip),%xmm0
movapd %xmm3,%xmm5
shufpd $0x1,%xmm0,%xmm5
mulsd  %xmm3,%xmm0
movapd %xmm3,%xmm6
mulpd  %xmm3,%xmm5
mulpd  %xmm2,%xmm3
unpcklpd %xmm1,%xmm1
addpd  %xmm4,%xmm1
movapd %xmm1,%xmm4
unpcklpd %xmm1,%xmm4
mulpd  %xmm4,%xmm5
mulpd  %xmm3,%xmm4
movupd %xmm4,(%rdi)
shufpd $0x1,%xmm0,%xmm6
unpcklpd %xmm1,%xmm2
mulpd  %xmm6,%xmm2
movapd %xmm1,%xmm0
mulpd  %xmm2,%xmm0
movupd %xmm0,0x10(%rdi)
unpckhpd %xmm1,%xmm1
mulpd  %xmm1,%xmm5
movupd %xmm5,0x20(%rdi)
mulpd  %xmm1,%xmm3
movupd %xmm3,0x30(%rdi)
mulsd  %xmm2,%xmm1
movsd  %xmm1,0x40(%rdi)
ret

LLVM IR (-O0):
define void @interpolate_order2_square(double %0, double %1, double* %2) {
%4 = alloca double
%5 = alloca double
%6 = alloca double*
store double %0, double* %4
store double %1, double* %5
store double* %2, double** %6
%7 = load double, double* %4
%8 = fmul double 2.500000e-01, %7
%9 = load double, double* %5
%10 = fmul double %8, %9
%11 = load double, double* %4
%12 = fsub double %11, 1.000000e+00
%13 = fmul double %10, %12
%14 = load double, double* %5
%15 = fsub double %14, 1.000000e+00
%16 = fmul double %13, %15
%17 = load double*, double** %6
%18 = getelementptr inbounds double, double* %17, i64 0
store double %16, double* %18
%19 = load double, double* %5
%20 = fmul double -5.000000e-01, %19
%21 = load double, double* %4
%22 = fsub double %21, 1.000000e+00
%23 = fmul double %20, %22
%24 = load double, double* %4
%25 = fadd double %24, 1.000000e+00
%26 = fmul double %23, %25
%27 = load double, double* %5
%28 = fsub double %27, 1.000000e+00
%29 = fmul double %26, %28
%30 = load double*, double** %6
%31 = getelementptr inbounds double, double* %30, i64 1
store double %29, double* %31
%32 = load double, double* %4
%33 = fmul double 2.500000e-01, %32
%34 = load double, double* %5
%35 = fmul double %33, %34
%36 = load double, double* %4
%37 = fadd double %36, 1.000000e+00
%38 = fmul double %35, %37
%39 = load double, double* %5
%40 = fsub double %39, 1.000000e+00
%41 = fmul double %38, %40
%42 = load double*, double** %6
%43 = getelementptr inbounds double, double* %42, i64 2
store double %41, double* %43
%44 = load double, double* %4
%45 = fmul double -5.000000e-01, %44
%46 = load double, double* %4
%47 = fsub double %46, 1.000000e+00
%48 = fmul double %45, %47
%49 = load double, double* %5
%50 = fsub double %49, 1.000000e+00
%51 = fmul double %48, %50
%52 = load double, double* %5
%53 = fadd double %52, 1.000000e+00
%54 = fmul double %51, %53
%55 = load double*, double** %6
%56 = getelementptr inbounds double, double* %55, i64 3
store double %54, double* %56
%57 = load double, double* %4
%58 = fsub double %57, 1.000000e+00
%59 = fmul double 1.000000e+00, %58
%60 = load double, double* %4
%61 = fadd double %60, 1.000000e+00
%62 = fmul double %59, %61
%63 = load double, double* %5
%64 = fsub double %63, 1.000000e+00
%65 = fmul double %62, %64
%66 = load double, double* %5
%67 = fadd double %66, 1.000000e+00
%68 = fmul double %65, %67
%69 = load double*, double** %6
%70 = getelementptr inbounds double, double* %69, i64 4
store double %68, double* %70
%71 = load double, double* %4
%72 = fmul double -5.000000e-01, %71
%73 = load double, double* %4
%74 = fadd double %73, 1.000000e+00
%75 = fmul double %72, %74
%76 = load double, double* %5
%77 = fsub double %76, 1.000000e+00
%78 = fmul double %75, %77
%79 = load double, double* %5
%80 = fadd double %79, 1.000000e+00
%81 = fmul double %78, %80
%82 = load double*, double** %6
%83 = getelementptr inbounds double, double* %82, i64 5
store double %81, double* %83
%84 = load double, double* %4
%85 = fmul double 2.500000e-01, %84
%86 = load double, double* %5
%87 = fmul double %85, %86
%88 = load double, double* %4
%89 = fsub double %88, 1.000000e+00
%90 = fmul double %87, %89
%91 = load double, double* %5
%92 = fadd double %91, 1.000000e+00
%93 = fmul double %90, %92
%94 = load double*, double** %6
%95 = getelementptr inbounds double, double* %94, i64 6
store double %93, double* %95
%96 = load double, double* %5
%97 = fmul double -5.000000e-01, %96
%98 = load double, double* %4
%99 = fsub double %98, 1.000000e+00
%100 = fmul double %97, %99
%101 = load double, double* %4
%102 = fadd double %101, 1.000000e+00
%103 = fmul double %100, %102
%104 = load double, double* %5
%105 = fadd double %104, 1.000000e+00
%106 = fmul double %103, %105
%107 = load double*, double** %6
%108 = getelementptr inbounds double, double* %107, i64 7
store double %106, double* %108
%109 = load double, double* %4
%110 = fmul double 2.500000e-01, %109
%111 = load double, double* %5
%112 = fmul double %110, %111
%113 = load double, double* %4
%114 = fadd double %113, 1.000000e+00
%115 = fmul double %112, %114
%116 = load double, double* %5
%117 = fadd double %116, 1.000000e+00
%118 = fmul double %115, %117
%119 = load double*, double** %6
%120 = getelementptr inbounds double, double* %119, i64 8
store double %118, double* %120
ret void
}

LLVM IR (-O2):
define void @interpolate_order2_square(double %0, double %1, double* nocapture writeonly %2) {
%4 = insertelement <2 x double> poison, double %0, i64 0
%5 = insertelement <2 x double> %4, double %1, i64 1
%6 = fmul <2 x double> %5, <double 2.500000e-01, double -5.000000e-01>
%7 = shufflevector <2 x double> %4, <2 x double> poison, <2 x i32> zeroinitializer
%8 = fadd <2 x double> %7, <double -1.000000e+00, double 1.000000e+00>
%9 = insertelement <2 x double> poison, double %1, i64 0
%10 = extractelement <2 x double> %8, i64 0
%11 = shufflevector <2 x double> %9, <2 x double> %8, <2 x i32> <i32 0, i32 2>
%12 = fmul <2 x double> %6, %11
%13 = fmul <2 x double> %8, %12
%14 = shufflevector <2 x double> %9, <2 x double> poison, <2 x i32> zeroinitializer
%15 = fadd <2 x double> %14, <double -1.000000e+00, double 1.000000e+00>
%16 = shufflevector <2 x double> %15, <2 x double> undef, <2 x i32> zeroinitializer
%17 = fmul <2 x double> %16, %13
%18 = bitcast double* %2 to <2 x double>*
store <2 x double> %17, <2 x double>* %18
%19 = getelementptr inbounds double, double* %2, i64 2
%20 = fmul double %0, -5.000000e-01
%21 = fmul double %20, %10
%22 = shufflevector <2 x double> %8, <2 x double> undef, <2 x i32> <i32 1, i32 undef>
%23 = insertelement <2 x double> %22, double %21, i64 1
%24 = shufflevector <2 x double> %12, <2 x double> %15, <2 x i32> <i32 0, i32 2>
%25 = fmul <2 x double> %23, %24
%26 = fmul <2 x double> %15, %25
%27 = bitcast double* %19 to <2 x double>*
store <2 x double> %26, <2 x double>* %27
%28 = extractelement <2 x double> %15, i64 1
%29 = getelementptr inbounds double, double* %2, i64 4
%30 = insertelement <2 x double> %22, double %20, i64 1
%31 = fmul <2 x double> %30, %8
%32 = fmul <2 x double> %31, %16
%33 = shufflevector <2 x double> %15, <2 x double> undef, <2 x i32> <i32 1, i32 1>
%34 = fmul <2 x double> %33, %32
%35 = bitcast double* %29 to <2 x double>*
store <2 x double> %34, <2 x double>* %35
%36 = getelementptr inbounds double, double* %2, i64 6
%37 = fmul <2 x double> %33, %13
%38 = bitcast double* %36 to <2 x double>*
store <2 x double> %37, <2 x double>* %38
%39 = extractelement <2 x double> %25, i64 0
%40 = fmul double %28, %39
%41 = getelementptr inbounds double, double* %2, i64 8
store double %40, double* %41
ret void
}