int foo(int *__restrict__ ia, int *__restrict__ ib, int *__restrict__ oxa,
        int *__restrict__ oxb) {
  int i;
  int oya[52], oyb[52];
  for (i = 0; i < 52; i++) {
    oya[i] = (ia[i] * oxa[i] + ib[i] * oxb[i]) >> 10;
    oyb[i] = (-ia[i] * oxb[i] + ib[i] * oxa[i]) >> 10;
  }
  return oya[22] + oyb[21];
}


x86 assembly (clang, -O0):
<foo>:
push   %rbp
mov    %rsp,%rbp
sub    $0x150,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
cmpl   $0x34,-0x24(%rbp)
jge    11b8 <foo+0xb8>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%ecx
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rdx
imul   (%rax,%rdx,4),%ecx
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rdx
mov    (%rax,%rdx,4),%eax
mov    -0x20(%rbp),%rdx
movslq -0x24(%rbp),%rsi
imul   (%rdx,%rsi,4),%eax
add    %eax,%ecx
sar    $0xa,%ecx
movslq -0x24(%rbp),%rax
mov    %ecx,-0x100(%rbp,%rax,4)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rdx
xor    %ecx,%ecx
sub    (%rax,%rdx,4),%ecx
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rdx
imul   (%rax,%rdx,4),%ecx
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rdx
mov    (%rax,%rdx,4),%eax
mov    -0x18(%rbp),%rdx
movslq -0x24(%rbp),%rsi
imul   (%rdx,%rsi,4),%eax
add    %eax,%ecx
sar    $0xa,%ecx
movslq -0x24(%rbp),%rax
mov    %ecx,-0x1d0(%rbp,%rax,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1122 <foo+0x22>
mov    -0xa8(%rbp),%eax
add    -0x17c(%rbp),%eax
add    $0x150,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<foo>:
sub    $0x128,%rsp
xor    %eax,%eax
nopl   0x0(%rax)
movdqu (%rdi,%rax,1),%xmm9
movdqu (%rdx,%rax,1),%xmm2
pshufd $0xf5,%xmm2,%xmm8
movdqu (%rsi,%rax,1),%xmm1
movdqu (%rcx,%rax,1),%xmm4
movdqa %xmm4,%xmm10
pmuludq %xmm1,%xmm10
pshufd $0xf5,%xmm1,%xmm5
pmuludq %xmm2,%xmm1
pmuludq %xmm9,%xmm2
pshufd $0xe8,%xmm2,%xmm2
pshufd $0xf5,%xmm9,%xmm7
pshufd $0xf5,%xmm4,%xmm3
movdqa %xmm3,%xmm0
pmuludq %xmm5,%xmm0
pmuludq %xmm8,%xmm5
pmuludq %xmm7,%xmm8
pshufd $0xe8,%xmm8,%xmm6
punpckldq %xmm6,%xmm2
pshufd $0xe8,%xmm10,%xmm6
pshufd $0xe8,%xmm0,%xmm0
punpckldq %xmm0,%xmm6
paddd  %xmm2,%xmm6
psrad  $0xa,%xmm6
movdqa %xmm6,0x50(%rsp,%rax,1)
pmuludq %xmm9,%xmm4
pshufd $0xe8,%xmm4,%xmm0
pmuludq %xmm7,%xmm3
pshufd $0xe8,%xmm3,%xmm2
punpckldq %xmm2,%xmm0
pshufd $0xe8,%xmm1,%xmm1
pshufd $0xe8,%xmm5,%xmm2
punpckldq %xmm2,%xmm1
psubd  %xmm0,%xmm1
psrad  $0xa,%xmm1
movdqa %xmm1,-0x80(%rsp,%rax,1)
add    $0x10,%rax
cmp    $0xd0,%rax
jne    1110 <foo+0x10>
mov    -0x2c(%rsp),%eax
add    0xa8(%rsp),%eax
add    $0x128,%rsp
ret

LLVM IR (-O0):
define i32 @foo(i32* noalias %0, i32* noalias %1, i32* noalias %2, i32* noalias %3) {
%5 = alloca i32*
%6 = alloca i32*
%7 = alloca i32*
%8 = alloca i32*
%9 = alloca i32
%10 = alloca [52 x i32]
%11 = alloca [52 x i32]
store i32* %0, i32** %5
store i32* %1, i32** %6
store i32* %2, i32** %7
store i32* %3, i32** %8
store i32 0, i32* %9
br label %12
12: ; preds = %71, %4
%13 = load i32, i32* %9
%14 = icmp slt i32 %13, 52
br i1 %14, label %15, label %74
15: ; preds = %12
%16 = load i32*, i32** %5
%17 = load i32, i32* %9
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds i32, i32* %16, i64 %18
%20 = load i32, i32* %19
%21 = load i32*, i32** %7
%22 = load i32, i32* %9
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds i32, i32* %21, i64 %23
%25 = load i32, i32* %24
%26 = mul nsw i32 %20, %25
%27 = load i32*, i32** %6
%28 = load i32, i32* %9
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds i32, i32* %27, i64 %29
%31 = load i32, i32* %30
%32 = load i32*, i32** %8
%33 = load i32, i32* %9
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds i32, i32* %32, i64 %34
%36 = load i32, i32* %35
%37 = mul nsw i32 %31, %36
%38 = add nsw i32 %26, %37
%39 = ashr i32 %38, 10
%40 = load i32, i32* %9
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds [52 x i32], [52 x i32]* %10, i64 0, i64 %41
store i32 %39, i32* %42
%43 = load i32*, i32** %5
%44 = load i32, i32* %9
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds i32, i32* %43, i64 %45
%47 = load i32, i32* %46
%48 = sub nsw i32 0, %47
%49 = load i32*, i32** %8
%50 = load i32, i32* %9
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds i32, i32* %49, i64 %51
%53 = load i32, i32* %52
%54 = mul nsw i32 %48, %53
%55 = load i32*, i32** %6
%56 = load i32, i32* %9
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds i32, i32* %55, i64 %57
%59 = load i32, i32* %58
%60 = load i32*, i32** %7
%61 = load i32, i32* %9
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds i32, i32* %60, i64 %62
%64 = load i32, i32* %63
%65 = mul nsw i32 %59, %64
%66 = add nsw i32 %54, %65
%67 = ashr i32 %66, 10
%68 = load i32, i32* %9
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds [52 x i32], [52 x i32]* %11, i64 0, i64 %69
store i32 %67, i32* %70
br label %71
71: ; preds = %15
%72 = load i32, i32* %9
%73 = add nsw i32 %72, 1
store i32 %73, i32* %9
br label %12
74: ; preds = %12
%75 = getelementptr inbounds [52 x i32], [52 x i32]* %10, i64 0, i64 22
%76 = load i32, i32* %75
%77 = getelementptr inbounds [52 x i32], [52 x i32]* %11, i64 0, i64 21
%78 = load i32, i32* %77
%79 = add nsw i32 %76, %78
ret i32 %79
}

LLVM IR (-O2):
define i32 @foo(i32* noalias nocapture %0, i32* noalias nocapture %1, i32* noalias nocapture %2, i32* noalias nocapture %3) {
%5 = alloca [52 x i32]
%6 = alloca [52 x i32]
%7 = bitcast [52 x i32]* %5 to i8*
%8 = bitcast [52 x i32]* %6 to i8*
br label %9
9: ; preds = %9, %4
%10 = phi i64 [ 0, %4 ], [ %35, %9 ]
%11 = getelementptr inbounds i32, i32* %0, i64 %10
%12 = bitcast i32* %11 to <4 x i32>*
%13 = load <4 x i32>, <4 x i32>* %12
%14 = getelementptr inbounds i32, i32* %2, i64 %10
%15 = bitcast i32* %14 to <4 x i32>*
%16 = load <4 x i32>, <4 x i32>* %15
%17 = mul nsw <4 x i32> %16, %13
%18 = getelementptr inbounds i32, i32* %1, i64 %10
%19 = bitcast i32* %18 to <4 x i32>*
%20 = load <4 x i32>, <4 x i32>* %19
%21 = getelementptr inbounds i32, i32* %3, i64 %10
%22 = bitcast i32* %21 to <4 x i32>*
%23 = load <4 x i32>, <4 x i32>* %22
%24 = mul nsw <4 x i32> %23, %20
%25 = add nsw <4 x i32> %24, %17
%26 = ashr <4 x i32> %25, <i32 10, i32 10, i32 10, i32 10>
%27 = getelementptr inbounds [52 x i32], [52 x i32]* %5, i64 0, i64 %10
%28 = bitcast i32* %27 to <4 x i32>*
store <4 x i32> %26, <4 x i32>* %28
%29 = mul <4 x i32> %23, %13
%30 = mul nsw <4 x i32> %20, %16
%31 = sub <4 x i32> %30, %29
%32 = ashr <4 x i32> %31, <i32 10, i32 10, i32 10, i32 10>
%33 = getelementptr inbounds [52 x i32], [52 x i32]* %6, i64 0, i64 %10
%34 = bitcast i32* %33 to <4 x i32>*
store <4 x i32> %32, <4 x i32>* %34
%35 = add nuw i64 %10, 4
%36 = icmp eq i64 %35, 52
br i1 %36, label %37, label %9
37: ; preds = %9
%38 = getelementptr inbounds [52 x i32], [52 x i32]* %5, i64 0, i64 22
%39 = load i32, i32* %38
%40 = getelementptr inbounds [52 x i32], [52 x i32]* %6, i64 0, i64 21
%41 = load i32, i32* %40
%42 = add nsw i32 %41, %39
ret i32 %42
}