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


x86 assembly (clang, -O0):
<foo>:
push   %rbp
mov    %rsp,%rbp
sub    $0x300,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
cmpl   $0x34,-0x34(%rbp)
jge    11fc <foo+0xfc>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x18(%rbp),%rcx
movslq -0x34(%rbp),%rdx
imul   (%rcx,%rdx,4),%eax
mov    -0x10(%rbp),%rcx
movslq -0x34(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
mov    -0x20(%rbp),%rdx
movslq -0x34(%rbp),%rsi
imul   (%rdx,%rsi,4),%ecx
add    %ecx,%eax
movslq %eax,%rcx
movslq -0x34(%rbp),%rax
mov    %rcx,-0x1e0(%rbp,%rax,8)
ff
mov    -0x8(%rbp),%rcx
movslq -0x34(%rbp),%rdx
xor    %eax,%eax
sub    (%rcx,%rdx,4),%eax
mov    -0x20(%rbp),%rcx
movslq -0x34(%rbp),%rdx
imul   (%rcx,%rdx,4),%eax
mov    -0x10(%rbp),%rcx
movslq -0x34(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
mov    -0x18(%rbp),%rdx
movslq -0x34(%rbp),%rsi
imul   (%rdx,%rsi,4),%ecx
add    %ecx,%eax
movslq %eax,%rcx
movslq -0x34(%rbp),%rax
mov    %rcx,-0x380(%rbp,%rax,8)
ff
movslq -0x34(%rbp),%rax
mov    -0x1e0(%rbp,%rax,8),%rax
ff
sar    $0xa,%rax
mov    %eax,%edx
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
movslq -0x34(%rbp),%rax
mov    -0x380(%rbp,%rax,8),%rax
ff
sar    $0xa,%rax
mov    %eax,%edx
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    112a <foo+0x2a>
add    $0x300,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<foo>:
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
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,%xmm6
pmuludq %xmm2,%xmm1
pmuludq %xmm9,%xmm2
pshufd $0xe8,%xmm2,%xmm2
pshufd $0xf5,%xmm9,%xmm7
pshufd $0xf5,%xmm4,%xmm3
movdqa %xmm3,%xmm0
pmuludq %xmm6,%xmm0
pmuludq %xmm8,%xmm6
pmuludq %xmm7,%xmm8
pshufd $0xe8,%xmm8,%xmm5
punpckldq %xmm5,%xmm2
pshufd $0xe8,%xmm10,%xmm5
pshufd $0xe8,%xmm0,%xmm0
punpckldq %xmm0,%xmm5
paddd  %xmm2,%xmm5
pmuludq %xmm9,%xmm4
pshufd $0xe8,%xmm4,%xmm0
pmuludq %xmm7,%xmm3
pshufd $0xe8,%xmm3,%xmm2
punpckldq %xmm2,%xmm0
pshufd $0xe8,%xmm1,%xmm1
pshufd $0xe8,%xmm6,%xmm2
punpckldq %xmm2,%xmm1
psubd  %xmm0,%xmm1
psrad  $0xa,%xmm5
movdqu %xmm5,(%r8,%rax,1)
psrad  $0xa,%xmm1
movdqu %xmm1,(%r9,%rax,1)
add    $0x10,%rax
cmp    $0xd0,%rax
jne    1110 <foo+0x10>
ret

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

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