void tpenc4(unsigned char *in, unsigned n, unsigned char *out) {
  unsigned char *op, *ip, *e;
  unsigned stride = n / 4;
  e = in + (n & ~(4 - 1));
  for (ip = in, op = out; ip < e; op++, ip += 4) {
    unsigned char *p = op;
    p[0] = ip[0];
    *(p += stride) = ip[1];
    *(p += stride) = ip[2];
    *(p += stride) = ip[3];
  }
  for (op = out + stride * 4; ip < in + n;)
    *op++ = *ip++;
}


x86 assembly (clang, -O0):
<tpenc4>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0xc(%rbp),%eax
shr    $0x2,%eax
mov    %eax,-0x34(%rbp)
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
and    $0xfffffffc,%ecx
mov    %ecx,%ecx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x20(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x30(%rbp),%rax
jae    11d1 <tpenc4+0xd1>
mov    -0x20(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%rax
mov    (%rax),%cl
mov    -0x40(%rbp),%rax
mov    %cl,(%rax)
mov    -0x28(%rbp),%rax
mov    0x1(%rax),%dl
mov    -0x34(%rbp),%ecx
mov    -0x40(%rbp),%rax
mov    %ecx,%ecx
mov    %rax,%rsi
add    %rcx,%rsi
mov    %rsi,-0x40(%rbp)
mov    %dl,(%rax,%rcx,1)
mov    -0x28(%rbp),%rax
mov    0x2(%rax),%dl
mov    -0x34(%rbp),%ecx
mov    -0x40(%rbp),%rax
mov    %ecx,%ecx
mov    %rax,%rsi
add    %rcx,%rsi
mov    %rsi,-0x40(%rbp)
mov    %dl,(%rax,%rcx,1)
mov    -0x28(%rbp),%rax
mov    0x3(%rax),%dl
mov    -0x34(%rbp),%ecx
mov    -0x40(%rbp),%rax
mov    %ecx,%ecx
mov    %rax,%rsi
add    %rcx,%rsi
mov    %rsi,-0x40(%rbp)
mov    %dl,(%rax,%rcx,1)
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
mov    -0x28(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x28(%rbp)
jmp    113b <tpenc4+0x3b>
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
mov    %ecx,%ecx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x8(%rbp),%rcx
mov    -0xc(%rbp),%edx
add    %rdx,%rcx
cmp    %rcx,%rax
jae    1222 <tpenc4+0x122>
mov    -0x28(%rbp),%rax
mov    %rax,%rcx
add    $0x1,%rcx
mov    %rcx,-0x28(%rbp)
mov    (%rax),%cl
mov    -0x20(%rbp),%rax
mov    %rax,%rdx
add    $0x1,%rdx
mov    %rdx,-0x20(%rbp)
mov    %cl,(%rax)
jmp    11e4 <tpenc4+0xe4>
pop    %rbp
ret

x86 assembly (clang, -O2):
<tpenc4>:
push   %r14
push   %rbx
mov    %esi,%r10d
and    $0xfffffffc,%esi
lea    (%rdi,%rsi,1),%r11
mov    %rdi,%rax
test   %esi,%esi
je     115a <tpenc4+0x5a>
mov    %r10,%r9
shr    $0x2,%r9
lea    (%r9,%r9,2),%r8
mov    %rdx,%rcx
mov    %rdi,%rax
cs nopw 0x0(%rax,%rax,1)
nop
movzbl (%rax),%ebx
mov    %bl,(%rcx)
movzbl 0x1(%rax),%ebx
mov    %bl,(%rcx,%r9,1)
movzbl 0x2(%rax),%ebx
mov    %bl,(%rcx,%r9,2)
movzbl 0x3(%rax),%ebx
mov    %bl,(%rcx,%r8,1)
add    $0x1,%rcx
add    $0x4,%rax
cmp    %r11,%rax
jb     1130 <tpenc4+0x30>
add    %r10,%rdi
cmp    %rdi,%rax
jae    11a2 <tpenc4+0xa2>
lea    (%rdx,%rsi,1),%r9
mov    %rdi,%r8
sub    %rax,%r8
cmp    $0x8,%r8
jb     1185 <tpenc4+0x85>
cmp    %rdi,%r9
jae    11a6 <tpenc4+0xa6>
add    %r10,%r11
sub    %rax,%r11
add    %rdx,%r11
cmp    %r11,%rax
jae    11a6 <tpenc4+0xa6>
mov    %r9,%rdx
mov    %rax,%rsi
nopl   0x0(%rax,%rax,1)
movzbl (%rsi),%eax
add    $0x1,%rsi
mov    %al,(%rdx)
add    $0x1,%rdx
cmp    %rdi,%rsi
jb     1190 <tpenc4+0x90>
pop    %rbx
pop    %r14
ret
cmp    $0x20,%r8
jae    11b3 <tpenc4+0xb3>
xor    %ecx,%ecx
jmp    129c <tpenc4+0x19c>
mov    %r8,%rcx
and    $0xffffffffffffffe0,%rcx
lea    -0x20(%rcx),%rbx
mov    %rbx,%r11
shr    $0x5,%r11
add    $0x1,%r11
mov    %r11d,%r10d
and    $0x3,%r10d
cmp    $0x60,%rbx
jae    11da <tpenc4+0xda>
xor    %ebx,%ebx
jmp    1250 <tpenc4+0x150>
and    $0xfffffffffffffffc,%r11
lea    (%rsi,%rdx,1),%r14
add    $0x70,%r14
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
movups (%rax,%rbx,1),%xmm0
movups 0x10(%rax,%rbx,1),%xmm1
movups %xmm0,-0x70(%r14,%rbx,1)
movups %xmm1,-0x60(%r14,%rbx,1)
movups 0x20(%rax,%rbx,1),%xmm0
movups 0x30(%rax,%rbx,1),%xmm1
movups %xmm0,-0x50(%r14,%rbx,1)
movups %xmm1,-0x40(%r14,%rbx,1)
movups 0x40(%rax,%rbx,1),%xmm0
movups 0x50(%rax,%rbx,1),%xmm1
movups %xmm0,-0x30(%r14,%rbx,1)
movups %xmm1,-0x20(%r14,%rbx,1)
movups 0x60(%rax,%rbx,1),%xmm0
movups 0x70(%rax,%rbx,1),%xmm1
movups %xmm0,-0x10(%r14,%rbx,1)
movups %xmm1,(%r14,%rbx,1)
sub    $0xffffffffffffff80,%rbx
add    $0xfffffffffffffffc,%r11
jne    11f0 <tpenc4+0xf0>
test   %r10,%r10
je     128d <tpenc4+0x18d>
lea    (%rbx,%rax,1),%r11
add    $0x10,%r11
add    %rbx,%rsi
add    %rsi,%rdx
add    $0x10,%rdx
shl    $0x5,%r10
xor    %esi,%esi
nopl   (%rax)
movups -0x10(%r11,%rsi,1),%xmm0
movups (%r11,%rsi,1),%xmm1
movups %xmm0,-0x10(%rdx,%rsi,1)
movups %xmm1,(%rdx,%rsi,1)
add    $0x20,%rsi
cmp    %rsi,%r10
jne    1270 <tpenc4+0x170>
cmp    %rcx,%r8
je     11a2 <tpenc4+0xa2>
test   $0x18,%r8b
je     12cf <tpenc4+0x1cf>
mov    %r8,%r10
and    $0xfffffffffffffff8,%r10
lea    (%r9,%r10,1),%rdx
lea    (%rax,%r10,1),%rsi
nopl   0x0(%rax,%rax,1)
mov    (%rax,%rcx,1),%rbx
mov    %rbx,(%r9,%rcx,1)
add    $0x8,%rcx
cmp    %rcx,%r10
jne    12b0 <tpenc4+0x1b0>
cmp    %r10,%r8
jne    1190 <tpenc4+0x90>
jmp    11a2 <tpenc4+0xa2>
add    %rcx,%rax
add    %rcx,%r9
jmp    1185 <tpenc4+0x85>

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

LLVM IR (-O2):
define void @tpenc4(i8* %0, i32 %1, i8* nocapture writeonly %2) {
%4 = ptrtoint i8* %0 to i64
%5 = and i32 %1, -4
%6 = zext i32 %5 to i64
%7 = getelementptr inbounds i8, i8* %0, i64 %6
%8 = icmp eq i32 %5, 0
br i1 %8, label %28, label %9
9: ; preds = %3
%10 = lshr i32 %1, 2
%11 = zext i32 %10 to i64
br label %12
12: ; preds = %9, %12
%13 = phi i8* [ %2, %9 ], [ %25, %12 ]
%14 = phi i8* [ %0, %9 ], [ %26, %12 ]
%15 = load i8, i8* %14
store i8 %15, i8* %13
%16 = getelementptr inbounds i8, i8* %14, i64 1
%17 = load i8, i8* %16
%18 = getelementptr inbounds i8, i8* %13, i64 %11
store i8 %17, i8* %18
%19 = getelementptr inbounds i8, i8* %14, i64 2
%20 = load i8, i8* %19
%21 = getelementptr inbounds i8, i8* %18, i64 %11
store i8 %20, i8* %21
%22 = getelementptr inbounds i8, i8* %14, i64 3
%23 = load i8, i8* %22
%24 = getelementptr inbounds i8, i8* %21, i64 %11
store i8 %23, i8* %24
%25 = getelementptr inbounds i8, i8* %13, i64 1
%26 = getelementptr inbounds i8, i8* %14, i64 4
%27 = icmp ult i8* %26, %7
br i1 %27, label %12, label %28
28: ; preds = %12, %3
%29 = phi i8* [ %0, %3 ], [ %26, %12 ]
%30 = ptrtoint i8* %29 to i64
%31 = zext i32 %1 to i64
%32 = getelementptr inbounds i8, i8* %0, i64 %31
%33 = icmp ult i8* %29, %32
br i1 %33, label %34, label %164
34: ; preds = %28
%35 = getelementptr i8, i8* %2, i64 %6
%36 = add i64 %4, %31
%37 = sub i64 %36, %30
%38 = icmp ult i64 %37, 8
br i1 %38, label %154, label %39
39: ; preds = %34
%40 = add i64 %4, %6
%41 = add i64 %40, %31
%42 = sub i64 %41, %30
%43 = getelementptr i8, i8* %2, i64 %42
%44 = add i64 %4, %31
%45 = sub i64 %44, %30
%46 = getelementptr i8, i8* %29, i64 %45
%47 = icmp ult i8* %35, %46
%48 = icmp ult i8* %29, %43
%49 = and i1 %47, %48
br i1 %49, label %154, label %50
50: ; preds = %39
%51 = icmp ult i64 %37, 32
br i1 %51, label %136, label %52
52: ; preds = %50
%53 = and i64 %37, -32
%54 = add i64 %53, -32
%55 = lshr exact i64 %54, 5
%56 = add nuw nsw i64 %55, 1
%57 = and i64 %56, 3
%58 = icmp ult i64 %54, 96
br i1 %58, label %110, label %59
59: ; preds = %52
%60 = and i64 %56, 1152921504606846972
br label %61
61: ; preds = %61, %59
%62 = phi i64 [ 0, %59 ], [ %107, %61 ]
%63 = phi i64 [ 0, %59 ], [ %108, %61 ]
%64 = getelementptr i8, i8* %35, i64 %62
%65 = getelementptr i8, i8* %29, i64 %62
%66 = bitcast i8* %65 to <16 x i8>*
%67 = load <16 x i8>, <16 x i8>* %66
%68 = getelementptr i8, i8* %65, i64 16
%69 = bitcast i8* %68 to <16 x i8>*
%70 = load <16 x i8>, <16 x i8>* %69
%71 = bitcast i8* %64 to <16 x i8>*
store <16 x i8> %67, <16 x i8>* %71
%72 = getelementptr i8, i8* %64, i64 16
%73 = bitcast i8* %72 to <16 x i8>*
store <16 x i8> %70, <16 x i8>* %73
%74 = or i64 %62, 32
%75 = getelementptr i8, i8* %35, i64 %74
%76 = getelementptr i8, i8* %29, i64 %74
%77 = bitcast i8* %76 to <16 x i8>*
%78 = load <16 x i8>, <16 x i8>* %77
%79 = getelementptr i8, i8* %76, i64 16
%80 = bitcast i8* %79 to <16 x i8>*
%81 = load <16 x i8>, <16 x i8>* %80
%82 = bitcast i8* %75 to <16 x i8>*
store <16 x i8> %78, <16 x i8>* %82
%83 = getelementptr i8, i8* %75, i64 16
%84 = bitcast i8* %83 to <16 x i8>*
store <16 x i8> %81, <16 x i8>* %84
%85 = or i64 %62, 64
%86 = getelementptr i8, i8* %35, i64 %85
%87 = getelementptr i8, i8* %29, i64 %85
%88 = bitcast i8* %87 to <16 x i8>*
%89 = load <16 x i8>, <16 x i8>* %88
%90 = getelementptr i8, i8* %87, i64 16
%91 = bitcast i8* %90 to <16 x i8>*
%92 = load <16 x i8>, <16 x i8>* %91
%93 = bitcast i8* %86 to <16 x i8>*
store <16 x i8> %89, <16 x i8>* %93
%94 = getelementptr i8, i8* %86, i64 16
%95 = bitcast i8* %94 to <16 x i8>*
store <16 x i8> %92, <16 x i8>* %95
%96 = or i64 %62, 96
%97 = getelementptr i8, i8* %35, i64 %96
%98 = getelementptr i8, i8* %29, i64 %96
%99 = bitcast i8* %98 to <16 x i8>*
%100 = load <16 x i8>, <16 x i8>* %99
%101 = getelementptr i8, i8* %98, i64 16
%102 = bitcast i8* %101 to <16 x i8>*
%103 = load <16 x i8>, <16 x i8>* %102
%104 = bitcast i8* %97 to <16 x i8>*
store <16 x i8> %100, <16 x i8>* %104
%105 = getelementptr i8, i8* %97, i64 16
%106 = bitcast i8* %105 to <16 x i8>*
store <16 x i8> %103, <16 x i8>* %106
%107 = add nuw i64 %62, 128
%108 = add i64 %63, 4
%109 = icmp eq i64 %108, %60
br i1 %109, label %110, label %61
110: ; preds = %61, %52
%111 = phi i64 [ 0, %52 ], [ %107, %61 ]
%112 = icmp eq i64 %57, 0
br i1 %112, label %129, label %113
113: ; preds = %110, %113
%114 = phi i64 [ %126, %113 ], [ %111, %110 ]
%115 = phi i64 [ %127, %113 ], [ 0, %110 ]
%116 = getelementptr i8, i8* %35, i64 %114
%117 = getelementptr i8, i8* %29, i64 %114
%118 = bitcast i8* %117 to <16 x i8>*
%119 = load <16 x i8>, <16 x i8>* %118
%120 = getelementptr i8, i8* %117, i64 16
%121 = bitcast i8* %120 to <16 x i8>*
%122 = load <16 x i8>, <16 x i8>* %121
%123 = bitcast i8* %116 to <16 x i8>*
store <16 x i8> %119, <16 x i8>* %123
%124 = getelementptr i8, i8* %116, i64 16
%125 = bitcast i8* %124 to <16 x i8>*
store <16 x i8> %122, <16 x i8>* %125
%126 = add nuw i64 %114, 32
%127 = add i64 %115, 1
%128 = icmp eq i64 %127, %57
br i1 %128, label %129, label %113
129: ; preds = %113, %110
%130 = icmp eq i64 %37, %53
br i1 %130, label %164, label %131
131: ; preds = %129
%132 = getelementptr i8, i8* %29, i64 %53
%133 = getelementptr i8, i8* %35, i64 %53
%134 = and i64 %37, 24
%135 = icmp eq i64 %134, 0
br i1 %135, label %154, label %136
136: ; preds = %50, %131
%137 = phi i64 [ %53, %131 ], [ 0, %50 ]
%138 = add i64 %4, %31
%139 = sub i64 %138, %30
%140 = and i64 %139, -8
%141 = getelementptr i8, i8* %35, i64 %140
%142 = getelementptr i8, i8* %29, i64 %140
br label %143
143: ; preds = %143, %136
%144 = phi i64 [ %137, %136 ], [ %150, %143 ]
%145 = getelementptr i8, i8* %35, i64 %144
%146 = getelementptr i8, i8* %29, i64 %144
%147 = bitcast i8* %146 to <8 x i8>*
%148 = load <8 x i8>, <8 x i8>* %147
%149 = bitcast i8* %145 to <8 x i8>*
store <8 x i8> %148, <8 x i8>* %149
%150 = add nuw i64 %144, 8
%151 = icmp eq i64 %150, %140
br i1 %151, label %152, label %143
152: ; preds = %143
%153 = icmp eq i64 %139, %140
br i1 %153, label %164, label %154
154: ; preds = %39, %34, %131, %152
%155 = phi i8* [ %35, %34 ], [ %35, %39 ], [ %133, %131 ], [ %141, %152 ]
%156 = phi i8* [ %29, %34 ], [ %29, %39 ], [ %132, %131 ], [ %142, %152 ]
br label %157
157: ; preds = %154, %157
%158 = phi i8* [ %162, %157 ], [ %155, %154 ]
%159 = phi i8* [ %160, %157 ], [ %156, %154 ]
%160 = getelementptr inbounds i8, i8* %159, i64 1
%161 = load i8, i8* %159
%162 = getelementptr inbounds i8, i8* %158, i64 1
store i8 %161, i8* %158
%163 = icmp ult i8* %160, %32
br i1 %163, label %157, label %164
164: ; preds = %157, %129, %152, %28
ret void
}