void mix3(double *ct, int sig) {
  double *p1 = ct + sig, *p2 = p1 + sig;
  double a, b, c;
  for (int i = 0; i < sig; i++) {
    a = ct[i];
    b = p1[i];
    c = p2[i];
    ct[i] = (a + b + c) * 0.57735026918962576451;
    p1[i] = (a - c) * 0.70710678118654752438;
    p2[i] = (2 * b - a - c) * 0.40824829046386301636;
  }
  sig = ((sig)*UINT64_C(1431655766) >> 32);
  if (sig) {
    mix3(ct, sig);
    mix3(p1, sig);
    mix3(p2, sig);
  }
}


x86 assembly (clang, -O0):
<mix3>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x18(%rbp),%rax
movslq -0xc(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    123e <mix3+0x11e>
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x28(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  0xe5b(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x28(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
movsd  0xe30(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
addsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
subsd  -0x38(%rbp),%xmm0
movsd  0xde1(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    115c <mix3+0x3c>
mov    -0xc(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    $0x55555556,%edi
mov    $0x0,%al
call   1040 <UINT64_C@plt>
mov    %eax,%ecx
mov    -0x40(%rbp),%eax
imul   %ecx,%eax
mov    $0x20,%ecx
sar    %cl,%eax
mov    %eax,-0xc(%rbp)
cmpl   $0x0,-0xc(%rbp)
je     1290 <mix3+0x170>
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
call   1030 <mix3@plt>
mov    -0x18(%rbp),%rdi
mov    -0xc(%rbp),%esi
call   1030 <mix3@plt>
mov    -0x20(%rbp),%rdi
mov    -0xc(%rbp),%esi
call   1030 <mix3@plt>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<mix3>:
push   %rbp
push   %r14
push   %rbx
test   %esi,%esi
jle    1281 <mix3+0x171>
movslq %esi,%r8
lea    (%rdi,%r8,8),%rax
lea    (%rax,%r8,8),%r11
mov    %esi,%r9d
cmp    $0x1,%esi
jne    1136 <mix3+0x26>
xor    %esi,%esi
jmp    1217 <mix3+0x107>
lea    (%rdi,%r9,8),%rdx
lea    (%r8,%r9,1),%rcx
lea    (%rdi,%rcx,8),%rsi
lea    (%r9,%r8,2),%rcx
lea    (%rdi,%rcx,8),%rbx
cmp    %rdi,%rsi
seta   %bpl
cmp    %rdx,%rax
setb   %r14b
cmp    %rdi,%rbx
seta   %cl
cmp    %rdx,%r11
setb   %r10b
cmp    %rbx,%rax
setb   %dl
cmp    %rsi,%r11
setb   %bl
xor    %esi,%esi
test   %r14b,%bpl
jne    1217 <mix3+0x107>
and    %r10b,%cl
jne    1217 <mix3+0x107>
and    %bl,%dl
jne    1217 <mix3+0x107>
mov    %r9d,%esi
and    $0xfffffffe,%esi
mov    %r8,%r10
shl    $0x4,%r10
movapd 0xe5e(%rip),%xmm0
movapd 0xe66(%rip),%xmm1
movapd 0xe6e(%rip),%xmm2
mov    %rsi,%rdx
mov    %rdi,%rcx
nopl   0x0(%rax,%rax,1)
movupd (%rcx),%xmm3
movupd (%rcx,%r8,8),%xmm4
movupd (%rcx,%r10,1),%xmm5
movapd %xmm3,%xmm6
addpd  %xmm4,%xmm6
addpd  %xmm5,%xmm6
mulpd  %xmm0,%xmm6
movupd %xmm6,(%rcx)
addpd  %xmm4,%xmm4
subpd  %xmm3,%xmm4
subpd  %xmm5,%xmm3
mulpd  %xmm1,%xmm3
movupd %xmm3,(%rcx,%r8,8)
subpd  %xmm5,%xmm4
mulpd  %xmm2,%xmm4
movupd %xmm4,(%rcx,%r10,1)
add    $0x10,%rcx
add    $0xfffffffffffffffe,%rdx
jne    11c0 <mix3+0xb0>
cmp    %r9,%rsi
je     1281 <mix3+0x171>
movsd  0xe11(%rip),%xmm0
movsd  0xe11(%rip),%xmm1
movsd  0xe11(%rip),%xmm2
nop
movsd  (%rdi,%rsi,8),%xmm3
movsd  (%rax,%rsi,8),%xmm4
movsd  (%r11,%rsi,8),%xmm5
movapd %xmm3,%xmm6
addsd  %xmm4,%xmm6
addsd  %xmm5,%xmm6
mulsd  %xmm0,%xmm6
movsd  %xmm6,(%rdi,%rsi,8)
addsd  %xmm4,%xmm4
subsd  %xmm3,%xmm4
subsd  %xmm5,%xmm3
mulsd  %xmm1,%xmm3
movsd  %xmm3,(%rax,%rsi,8)
subsd  %xmm5,%xmm4
mulsd  %xmm2,%xmm4
movsd  %xmm4,(%r11,%rsi,8)
add    $0x1,%rsi
cmp    %rsi,%r9
jne    1230 <mix3+0x120>
mov    $0x55555556,%edi
xor    %eax,%eax
pop    %rbx
pop    %r14
pop    %rbp
jmp    1030 <UINT64_C@plt>

LLVM IR (-O0):
define void @mix3(double* %0, i32 %1) {
%3 = alloca double*
%4 = alloca i32
%5 = alloca double*
%6 = alloca double*
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca i32
store double* %0, double** %3
store i32 %1, i32* %4
%11 = load double*, double** %3
%12 = load i32, i32* %4
%13 = sext i32 %12 to i64
%14 = getelementptr inbounds double, double* %11, i64 %13
store double* %14, double** %5
%15 = load double*, double** %5
%16 = load i32, i32* %4
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds double, double* %15, i64 %17
store double* %18, double** %6
store i32 0, i32* %10
br label %19
19: ; preds = %68, %2
%20 = load i32, i32* %10
%21 = load i32, i32* %4
%22 = icmp slt i32 %20, %21
br i1 %22, label %23, label %71
23: ; preds = %19
%24 = load double*, double** %3
%25 = load i32, i32* %10
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds double, double* %24, i64 %26
%28 = load double, double* %27
store double %28, double* %7
%29 = load double*, double** %5
%30 = load i32, i32* %10
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds double, double* %29, i64 %31
%33 = load double, double* %32
store double %33, double* %8
%34 = load double*, double** %6
%35 = load i32, i32* %10
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %34, i64 %36
%38 = load double, double* %37
store double %38, double* %9
%39 = load double, double* %7
%40 = load double, double* %8
%41 = fadd double %39, %40
%42 = load double, double* %9
%43 = fadd double %41, %42
%44 = fmul double %43, 0x3FE279A74590331C
%45 = load double*, double** %3
%46 = load i32, i32* %10
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %45, i64 %47
store double %44, double* %48
%49 = load double, double* %7
%50 = load double, double* %9
%51 = fsub double %49, %50
%52 = fmul double %51, 0x3FE6A09E667F3BCD
%53 = load double*, double** %5
%54 = load i32, i32* %10
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %53, i64 %55
store double %52, double* %56
%57 = load double, double* %8
%58 = load double, double* %7
%59 = fneg double %58
%60 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %57, double %59)
%61 = load double, double* %9
%62 = fsub double %60, %61
%63 = fmul double %62, 0x3FDA20BD700C2C3E
%64 = load double*, double** %6
%65 = load i32, i32* %10
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double, double* %64, i64 %66
store double %63, double* %67
br label %68
68: ; preds = %23
%69 = load i32, i32* %10
%70 = add nsw i32 %69, 1
store i32 %70, i32* %10
br label %19
71: ; preds = %19
%72 = load i32, i32* %4
%73 = call i32 (i32, ...) bitcast (i32 (...)* @UINT64_C to i32 (i32, ...)*)(i32 1431655766)
%74 = mul nsw i32 %72, %73
%75 = ashr i32 %74, 32
store i32 %75, i32* %4
%76 = load i32, i32* %4
%77 = icmp ne i32 %76, 0
br i1 %77, label %78, label %85
78: ; preds = %71
%79 = load double*, double** %3
%80 = load i32, i32* %4
%81 = load double*, double** %5
%82 = load i32, i32* %4
%83 = load double*, double** %6
%84 = load i32, i32* %4
br label %85
85: ; preds = %78, %71
ret void
}

LLVM IR (-O2):
define void @mix3(double* nocapture %0, i32 %1) {
%3 = sext i32 %1 to i64
%4 = getelementptr double, double* %0, i64 %3
%5 = getelementptr double, double* %4, i64 %3
%6 = icmp sgt i32 %1, 0
br i1 %6, label %7, label %59
7: ; preds = %2
%8 = zext i32 %1 to i64
%9 = icmp eq i32 %1, 1
br i1 %9, label %57, label %10
10: ; preds = %7
%11 = getelementptr double, double* %0, i64 %8
%12 = add nsw i64 %3, %8
%13 = getelementptr double, double* %0, i64 %12
%14 = shl nsw i64 %3, 1
%15 = add nsw i64 %14, %8
%16 = getelementptr double, double* %0, i64 %15
%17 = icmp ugt double* %13, %0
%18 = icmp ult double* %4, %11
%19 = and i1 %17, %18
%20 = icmp ugt double* %16, %0
%21 = icmp ult double* %5, %11
%22 = and i1 %20, %21
%23 = or i1 %19, %22
%24 = icmp ult double* %4, %16
%25 = icmp ult double* %5, %13
%26 = and i1 %24, %25
%27 = or i1 %23, %26
br i1 %27, label %57, label %28
28: ; preds = %10
%29 = and i64 %8, 4294967294
br label %30
30: ; preds = %30, %28
%31 = phi i64 [ 0, %28 ], [ %53, %30 ]
%32 = getelementptr inbounds double, double* %0, i64 %31
%33 = bitcast double* %32 to <2 x double>*
%34 = load <2 x double>, <2 x double>* %33
%35 = getelementptr inbounds double, double* %4, i64 %31
%36 = bitcast double* %35 to <2 x double>*
%37 = load <2 x double>, <2 x double>* %36
%38 = getelementptr inbounds double, double* %5, i64 %31
%39 = bitcast double* %38 to <2 x double>*
%40 = load <2 x double>, <2 x double>* %39
%41 = fadd <2 x double> %34, %37
%42 = fadd <2 x double> %41, %40
%43 = fmul <2 x double> %42, <double 0x3FE279A74590331C, double 0x3FE279A74590331C>
%44 = bitcast double* %32 to <2 x double>*
store <2 x double> %43, <2 x double>* %44
%45 = fsub <2 x double> %34, %40
%46 = fmul <2 x double> %45, <double 0x3FE6A09E667F3BCD, double 0x3FE6A09E667F3BCD>
%47 = bitcast double* %35 to <2 x double>*
store <2 x double> %46, <2 x double>* %47
%48 = fneg <2 x double> %34
%49 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %37, <2 x double> <double 2.000000e+00, double 2.000000e+00>, <2 x double> %48)
%50 = fsub <2 x double> %49, %40
%51 = fmul <2 x double> %50, <double 0x3FDA20BD700C2C3E, double 0x3FDA20BD700C2C3E>
%52 = bitcast double* %38 to <2 x double>*
store <2 x double> %51, <2 x double>* %52
%53 = add nuw i64 %31, 2
%54 = icmp eq i64 %53, %29
br i1 %54, label %55, label %30
55: ; preds = %30
%56 = icmp eq i64 %29, %8
br i1 %56, label %59, label %57
57: ; preds = %10, %7, %55
%58 = phi i64 [ 0, %10 ], [ 0, %7 ], [ %29, %55 ]
br label %61
59: ; preds = %61, %55, %2
%60 = tail call i32 (i32, ...) bitcast (i32 (...)* @UINT64_C to i32 (i32, ...)*)(i32 1431655766)
ret void
61: ; preds = %57, %61
%62 = phi i64 [ %78, %61 ], [ %58, %57 ]
%63 = getelementptr inbounds double, double* %0, i64 %62
%64 = load double, double* %63
%65 = getelementptr inbounds double, double* %4, i64 %62
%66 = load double, double* %65
%67 = getelementptr inbounds double, double* %5, i64 %62
%68 = load double, double* %67
%69 = fadd double %64, %66
%70 = fadd double %69, %68
%71 = fmul double %70, 0x3FE279A74590331C
store double %71, double* %63
%72 = fsub double %64, %68
%73 = fmul double %72, 0x3FE6A09E667F3BCD
store double %73, double* %65
%74 = fneg double %64
%75 = tail call double @llvm.fmuladd.f64(double %66, double 2.000000e+00, double %74)
%76 = fsub double %75, %68
%77 = fmul double %76, 0x3FDA20BD700C2C3E
store double %77, double* %67
%78 = add nuw nsw i64 %62, 1
%79 = icmp eq i64 %78, %8
br i1 %79, label %59, label %61
}