void syndrome(unsigned char *s, const unsigned char *pk, unsigned char *e) {
  unsigned char b, row[8192 / 8];
  const unsigned char *pk_ptr = pk;
  int i, j;
  for (i = 0; i < (((128 * 13) + 7) / 8); i++)
    s[i] = 0;
  for (i = 0; i < (128 * 13); i++) {
    for (j = 0; j < 8192 / 8; j++)
      row[j] = 0;
    for (j = 0; j < (((8192 - (128 * 13)) + 7) / 8); j++)
      row[8192 / 8 - (((8192 - (128 * 13)) + 7) / 8) + j] = pk_ptr[j];
    row[i / 8] |= 1 << (i % 8);
    b = 0;
    for (j = 0; j < 8192 / 8; j++)
      b ^= row[j] & e[j];
    b ^= b >> 4;
    b ^= b >> 2;
    b ^= b >> 1;
    b &= 1;
    s[i / 8] |= (b << (i % 8));
    pk_ptr += (((8192 - (128 * 13)) + 7) / 8);
  }
}


x86 assembly (clang, -O0):
<syndrome>:
push   %rbp
mov    %rsp,%rbp
sub    $0x3c0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x428(%rbp)
movl   $0x0,-0x42c(%rbp)
cmpl   $0xd0,-0x42c(%rbp)
jge    115f <syndrome+0x5f>
mov    -0x8(%rbp),%rax
movslq -0x42c(%rbp),%rcx
movb   $0x0,(%rax,%rcx,1)
mov    -0x42c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x42c(%rbp)
jmp    112c <syndrome+0x2c>
movl   $0x0,-0x42c(%rbp)
cmpl   $0x680,-0x42c(%rbp)
jge    134e <syndrome+0x24e>
movl   $0x0,-0x430(%rbp)
cmpl   $0x400,-0x430(%rbp)
jge    11b6 <syndrome+0xb6>
movslq -0x430(%rbp),%rax
movb   $0x0,-0x420(%rbp,%rax,1)
mov    -0x430(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x430(%rbp)
jmp    1183 <syndrome+0x83>
movl   $0x0,-0x430(%rbp)
cmpl   $0x330,-0x430(%rbp)
jge    1209 <syndrome+0x109>
mov    -0x428(%rbp),%rax
movslq -0x430(%rbp),%rcx
mov    (%rax,%rcx,1),%cl
mov    -0x430(%rbp),%eax
add    $0xd0,%eax
cltq
mov    %cl,-0x420(%rbp,%rax,1)
mov    -0x430(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x430(%rbp)
jmp    11c0 <syndrome+0xc0>
mov    -0x42c(%rbp),%eax
mov    $0x8,%ecx
cltd
idiv   %ecx
mov    %edx,%ecx
mov    $0x1,%eax
shl    %cl,%eax
mov    %eax,-0x434(%rbp)
mov    -0x42c(%rbp),%eax
mov    $0x8,%ecx
cltd
idiv   %ecx
mov    -0x434(%rbp),%edx
cltq
movzbl -0x420(%rbp,%rax,1),%ecx
ff
or     %edx,%ecx
mov    %cl,-0x420(%rbp,%rax,1)
movb   $0x0,-0x19(%rbp)
movl   $0x0,-0x430(%rbp)
cmpl   $0x400,-0x430(%rbp)
jge    12a8 <syndrome+0x1a8>
movslq -0x430(%rbp),%rax
movzbl -0x420(%rbp,%rax,1),%ecx
ff
mov    -0x18(%rbp),%rax
movslq -0x430(%rbp),%rdx
movzbl (%rax,%rdx,1),%eax
and    %eax,%ecx
movzbl -0x19(%rbp),%eax
xor    %ecx,%eax
mov    %al,-0x19(%rbp)
mov    -0x430(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x430(%rbp)
jmp    125b <syndrome+0x15b>
movzbl -0x19(%rbp),%ecx
sar    $0x4,%ecx
movzbl -0x19(%rbp),%eax
xor    %ecx,%eax
mov    %al,-0x19(%rbp)
movzbl -0x19(%rbp),%ecx
sar    $0x2,%ecx
movzbl -0x19(%rbp),%eax
xor    %ecx,%eax
mov    %al,-0x19(%rbp)
movzbl -0x19(%rbp),%ecx
sar    $0x1,%ecx
movzbl -0x19(%rbp),%eax
xor    %ecx,%eax
mov    %al,-0x19(%rbp)
movzbl -0x19(%rbp),%eax
and    $0x1,%eax
mov    %al,-0x19(%rbp)
movzbl -0x19(%rbp),%esi
mov    -0x42c(%rbp),%eax
mov    $0x8,%ecx
cltd
idiv   %ecx
mov    %edx,%ecx
shl    %cl,%esi
mov    -0x8(%rbp),%rax
mov    %rax,-0x440(%rbp)
mov    -0x42c(%rbp),%eax
mov    $0x8,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x440(%rbp),%rax
movslq %ecx,%rcx
movzbl (%rax,%rcx,1),%edx
or     %esi,%edx
mov    %dl,(%rax,%rcx,1)
mov    -0x428(%rbp),%rax
add    $0x330,%rax
mov    %rax,-0x428(%rbp)
mov    -0x42c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x42c(%rbp)
jmp    1169 <syndrome+0x69>
add    $0x3c0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<syndrome>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x400,%rsp
mov    %rdx,%rbx
mov    %rsi,%r14
mov    %rdi,%r15
pxor   %xmm0,%xmm0
movdqu %xmm0,0xc0(%rdi)
movdqu %xmm0,0xb0(%rdi)
movdqu %xmm0,0xa0(%rdi)
movdqu %xmm0,0x90(%rdi)
movdqu %xmm0,0x80(%rdi)
movdqu %xmm0,0x70(%rdi)
movdqu %xmm0,0x60(%rdi)
movdqu %xmm0,0x50(%rdi)
movdqu %xmm0,0x40(%rdi)
movdqu %xmm0,0x30(%rdi)
movdqu %xmm0,0x20(%rdi)
movdqu %xmm0,0x10(%rdi)
movdqu %xmm0,(%rdi)
xor    %r13d,%r13d
lea    0xd0(%rsp),%r12
nopw   0x0(%rax,%rax,1)
pxor   %xmm0,%xmm0
movdqa %xmm0,0xc0(%rsp)
movdqa %xmm0,0xb0(%rsp)
movdqa %xmm0,0xa0(%rsp)
movdqa %xmm0,0x90(%rsp)
movdqa %xmm0,0x80(%rsp)
movdqa %xmm0,0x70(%rsp)
movdqa %xmm0,0x60(%rsp)
movdqa %xmm0,0x50(%rsp)
movdqa %xmm0,0x40(%rsp)
movdqa %xmm0,0x30(%rsp)
movdqa %xmm0,0x20(%rsp)
movdqa %xmm0,0x10(%rsp)
movdqa %xmm0,(%rsp)
mov    $0x330,%edx
mov    %r12,%rdi
mov    %r14,%rsi
call   1030 <memcpy@plt>
mov    %r13d,%ecx
and    $0x7,%ecx
mov    $0x1,%edx
shl    %cl,%edx
mov    %r13d,%esi
shr    $0x3,%rsi
or     %dl,(%rsp,%rsi,1)
xor    %edx,%edx
pxor   %xmm1,%xmm1
pxor   %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movdqu (%rbx,%rdx,1),%xmm2
movdqu 0x10(%rbx,%rdx,1),%xmm3
pand   (%rsp,%rdx,1),%xmm2
pxor   %xmm1,%xmm2
movdqu 0x20(%rbx,%rdx,1),%xmm1
pand   0x10(%rsp,%rdx,1),%xmm3
pxor   %xmm0,%xmm3
movdqu 0x30(%rbx,%rdx,1),%xmm0
pand   0x20(%rsp,%rdx,1),%xmm1
pand   0x30(%rsp,%rdx,1),%xmm0
pxor   %xmm2,%xmm1
pxor   %xmm3,%xmm0
add    $0x40,%rdx
cmp    $0x400,%rdx
jne    1230 <syndrome+0x120>
pxor   %xmm1,%xmm0
pshufd $0xee,%xmm0,%xmm1
pxor   %xmm0,%xmm1
pshufd $0x55,%xmm1,%xmm0
pxor   %xmm1,%xmm0
movdqa %xmm0,%xmm1
psrld  $0x10,%xmm1
pxor   %xmm0,%xmm1
movdqa %xmm1,%xmm0
psrlw  $0x8,%xmm0
pxor   %xmm1,%xmm0
movd   %xmm0,%edx
mov    %edx,%eax
shr    $0x4,%al
xor    %dl,%al
mov    %eax,%edx
shr    $0x2,%dl
xor    %al,%dl
mov    %edx,%eax
shr    %al
xor    %dl,%al
and    $0x1,%al
movzbl %al,%eax
shl    %cl,%eax
or     %al,(%r15,%rsi,1)
add    $0x330,%r14
add    $0x1,%r13d
cmp    $0x680,%r13d
jne    1190 <syndrome+0x80>
add    $0x400,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @syndrome(i8* %0, i8* %1, i8* %2) {
%4 = alloca i8*
%5 = alloca i8*
%6 = alloca i8*
%7 = alloca i8
%8 = alloca [1024 x i8]
%9 = alloca i8*
%10 = alloca i32
%11 = alloca i32
store i8* %0, i8** %4
store i8* %1, i8** %5
store i8* %2, i8** %6
%12 = load i8*, i8** %5
store i8* %12, i8** %9
store i32 0, i32* %10
br label %13
13: ; preds = %21, %3
%14 = load i32, i32* %10
%15 = icmp slt i32 %14, 208
br i1 %15, label %16, label %24
16: ; preds = %13
%17 = load i8*, i8** %4
%18 = load i32, i32* %10
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds i8, i8* %17, i64 %19
store i8 0, i8* %20
br label %21
21: ; preds = %16
%22 = load i32, i32* %10
%23 = add nsw i32 %22, 1
store i32 %23, i32* %10
br label %13
24: ; preds = %13
store i32 0, i32* %10
br label %25
25: ; preds = %133, %24
%26 = load i32, i32* %10
%27 = icmp slt i32 %26, 1664
br i1 %27, label %28, label %136
28: ; preds = %25
store i32 0, i32* %11
br label %29
29: ; preds = %36, %28
%30 = load i32, i32* %11
%31 = icmp slt i32 %30, 1024
br i1 %31, label %32, label %39
32: ; preds = %29
%33 = load i32, i32* %11
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds [1024 x i8], [1024 x i8]* %8, i64 0, i64 %34
store i8 0, i8* %35
br label %36
36: ; preds = %32
%37 = load i32, i32* %11
%38 = add nsw i32 %37, 1
store i32 %38, i32* %11
br label %29
39: ; preds = %29
store i32 0, i32* %11
br label %40
40: ; preds = %53, %39
%41 = load i32, i32* %11
%42 = icmp slt i32 %41, 816
br i1 %42, label %43, label %56
43: ; preds = %40
%44 = load i8*, i8** %9
%45 = load i32, i32* %11
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds i8, i8* %44, i64 %46
%48 = load i8, i8* %47
%49 = load i32, i32* %11
%50 = add nsw i32 208, %49
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds [1024 x i8], [1024 x i8]* %8, i64 0, i64 %51
store i8 %48, i8* %52
br label %53
53: ; preds = %43
%54 = load i32, i32* %11
%55 = add nsw i32 %54, 1
store i32 %55, i32* %11
br label %40
56: ; preds = %40
%57 = load i32, i32* %10
%58 = srem i32 %57, 8
%59 = shl i32 1, %58
%60 = load i32, i32* %10
%61 = sdiv i32 %60, 8
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds [1024 x i8], [1024 x i8]* %8, i64 0, i64 %62
%64 = load i8, i8* %63
%65 = zext i8 %64 to i32
%66 = or i32 %65, %59
%67 = trunc i32 %66 to i8
store i8 %67, i8* %63
store i8 0, i8* %7
store i32 0, i32* %11
br label %68
68: ; preds = %88, %56
%69 = load i32, i32* %11
%70 = icmp slt i32 %69, 1024
br i1 %70, label %71, label %91
71: ; preds = %68
%72 = load i32, i32* %11
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds [1024 x i8], [1024 x i8]* %8, i64 0, i64 %73
%75 = load i8, i8* %74
%76 = zext i8 %75 to i32
%77 = load i8*, i8** %6
%78 = load i32, i32* %11
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds i8, i8* %77, i64 %79
%81 = load i8, i8* %80
%82 = zext i8 %81 to i32
%83 = and i32 %76, %82
%84 = load i8, i8* %7
%85 = zext i8 %84 to i32
%86 = xor i32 %85, %83
%87 = trunc i32 %86 to i8
store i8 %87, i8* %7
br label %88
88: ; preds = %71
%89 = load i32, i32* %11
%90 = add nsw i32 %89, 1
store i32 %90, i32* %11
br label %68
91: ; preds = %68
%92 = load i8, i8* %7
%93 = zext i8 %92 to i32
%94 = ashr i32 %93, 4
%95 = load i8, i8* %7
%96 = zext i8 %95 to i32
%97 = xor i32 %96, %94
%98 = trunc i32 %97 to i8
store i8 %98, i8* %7
%99 = load i8, i8* %7
%100 = zext i8 %99 to i32
%101 = ashr i32 %100, 2
%102 = load i8, i8* %7
%103 = zext i8 %102 to i32
%104 = xor i32 %103, %101
%105 = trunc i32 %104 to i8
store i8 %105, i8* %7
%106 = load i8, i8* %7
%107 = zext i8 %106 to i32
%108 = ashr i32 %107, 1
%109 = load i8, i8* %7
%110 = zext i8 %109 to i32
%111 = xor i32 %110, %108
%112 = trunc i32 %111 to i8
store i8 %112, i8* %7
%113 = load i8, i8* %7
%114 = zext i8 %113 to i32
%115 = and i32 %114, 1
%116 = trunc i32 %115 to i8
store i8 %116, i8* %7
%117 = load i8, i8* %7
%118 = zext i8 %117 to i32
%119 = load i32, i32* %10
%120 = srem i32 %119, 8
%121 = shl i32 %118, %120
%122 = load i8*, i8** %4
%123 = load i32, i32* %10
%124 = sdiv i32 %123, 8
%125 = sext i32 %124 to i64
%126 = getelementptr inbounds i8, i8* %122, i64 %125
%127 = load i8, i8* %126
%128 = zext i8 %127 to i32
%129 = or i32 %128, %121
%130 = trunc i32 %129 to i8
store i8 %130, i8* %126
%131 = load i8*, i8** %9
%132 = getelementptr inbounds i8, i8* %131, i64 816
store i8* %132, i8** %9
br label %133
133: ; preds = %91
%134 = load i32, i32* %10
%135 = add nsw i32 %134, 1
store i32 %135, i32* %10
br label %25
136: ; preds = %25
ret void
}

LLVM IR (-O2):
define void @syndrome(i8* nocapture %0, i8* nocapture %1, i8* nocapture %2) {
%4 = alloca [1024 x i8]
%5 = getelementptr inbounds [1024 x i8], [1024 x i8]* %4, i64 0, i64 0
%6 = getelementptr inbounds [1024 x i8], [1024 x i8]* %4, i64 0, i64 0
%7 = getelementptr inbounds [1024 x i8], [1024 x i8]* %4, i64 0, i64 208
br label %8
8: ; preds = %3, %58
%9 = phi i32 [ %75, %58 ], [ 0, %3 ]
%10 = phi i8* [ %74, %58 ], [ %1, %3 ]
%11 = and i32 %9, 7
%12 = shl nuw nsw i32 1, %11
%13 = lshr i32 %9, 3
%14 = zext i32 %13 to i64
%15 = getelementptr inbounds [1024 x i8], [1024 x i8]* %4, i64 0, i64 %14
%16 = load i8, i8* %15
%17 = trunc i32 %12 to i8
%18 = or i8 %16, %17
store i8 %18, i8* %15
br label %19
19: ; preds = %19, %8
%20 = phi i64 [ 0, %8 ], [ %56, %19 ]
%21 = phi <16 x i8> [ zeroinitializer, %8 ], [ %54, %19 ]
%22 = phi <16 x i8> [ zeroinitializer, %8 ], [ %55, %19 ]
%23 = getelementptr inbounds [1024 x i8], [1024 x i8]* %4, i64 0, i64 %20
%24 = bitcast i8* %23 to <16 x i8>*
%25 = load <16 x i8>, <16 x i8>* %24
%26 = getelementptr inbounds i8, i8* %23, i64 16
%27 = bitcast i8* %26 to <16 x i8>*
%28 = load <16 x i8>, <16 x i8>* %27
%29 = getelementptr inbounds i8, i8* %2, i64 %20
%30 = bitcast i8* %29 to <16 x i8>*
%31 = load <16 x i8>, <16 x i8>* %30
%32 = getelementptr inbounds i8, i8* %29, i64 16
%33 = bitcast i8* %32 to <16 x i8>*
%34 = load <16 x i8>, <16 x i8>* %33
%35 = and <16 x i8> %31, %25
%36 = and <16 x i8> %34, %28
%37 = xor <16 x i8> %35, %21
%38 = xor <16 x i8> %36, %22
%39 = or i64 %20, 32
%40 = getelementptr inbounds [1024 x i8], [1024 x i8]* %4, i64 0, i64 %39
%41 = bitcast i8* %40 to <16 x i8>*
%42 = load <16 x i8>, <16 x i8>* %41
%43 = getelementptr inbounds i8, i8* %40, i64 16
%44 = bitcast i8* %43 to <16 x i8>*
%45 = load <16 x i8>, <16 x i8>* %44
%46 = getelementptr inbounds i8, i8* %2, i64 %39
%47 = bitcast i8* %46 to <16 x i8>*
%48 = load <16 x i8>, <16 x i8>* %47
%49 = getelementptr inbounds i8, i8* %46, i64 16
%50 = bitcast i8* %49 to <16 x i8>*
%51 = load <16 x i8>, <16 x i8>* %50
%52 = and <16 x i8> %48, %42
%53 = and <16 x i8> %51, %45
%54 = xor <16 x i8> %52, %37
%55 = xor <16 x i8> %53, %38
%56 = add nuw nsw i64 %20, 64
%57 = icmp eq i64 %56, 1024
br i1 %57, label %58, label %19
58: ; preds = %19
%59 = xor <16 x i8> %55, %54
%60 = call i8 @llvm.vector.reduce.xor.v16i8(<16 x i8> %59)
%61 = lshr i8 %60, 4
%62 = xor i8 %61, %60
%63 = lshr i8 %62, 2
%64 = xor i8 %63, %62
%65 = lshr i8 %64, 1
%66 = xor i8 %65, %64
%67 = and i8 %66, 1
%68 = zext i8 %67 to i32
%69 = shl nuw i32 %68, %11
%70 = getelementptr inbounds i8, i8* %0, i64 %14
%71 = load i8, i8* %70
%72 = trunc i32 %69 to i8
%73 = or i8 %71, %72
store i8 %73, i8* %70
%74 = getelementptr inbounds i8, i8* %10, i64 816
%75 = add nuw nsw i32 %9, 1
%76 = icmp eq i32 %75, 1664
br i1 %76, label %77, label %8
77: ; preds = %58
ret void
}