yencode 1.0.8 → 1.1.2

package/src/decoder_neon.cc

@@ -0,0 +1,474 @@
+#include "common.h"
+#ifdef __ARM_NEON
+
+#ifndef __aarch64__
+#define YENC_DEC_USE_THINTABLE 1
+#endif
+#include "decoder_common.h"
+
+
+#if defined(_MSC_VER) && !defined(__clang__)
+# define vld1_u8_align(p, a) vld1_u8_ex(p, a*8)
+# define vld1q_u8_align(p, a) vld1q_u8_ex(p, a*8)
+#elif defined(__GNUC__)
+# define vld1_u8_align(p, n) vld1_u8((uint8_t*)__builtin_assume_aligned(p, n))
+# define vld1q_u8_align(p, n) vld1q_u8((uint8_t*)__builtin_assume_aligned(p, n))
+#else
+# define vld1_u8_align(p, n) vld1_u8(p)
+# define vld1q_u8_align(p, n) vld1q_u8(p)
+#endif
+
+
+// for compilers that lack these functions
+#if defined(__clang__) || (defined(__GNUC__) && (defined(__aarch64__) && __GNUC__ >= 8))
+# define vld1q_u8_x2_align(p, n) vld1q_u8_x2((uint8_t*)__builtin_assume_aligned(p, n))
+#else
+# define vld1q_u8_x2_align(p, n) vcreate2_u8(vld1q_u8_align(p, (n)/2), vld1q_u8_align((p)+16, (n)/2))
+#endif
+// Clang wrongly assumes alignment on vld1q_u8_x2, and ARMv7 GCC doesn't support the function, so effectively, it can only be used in ARMv8 compilers
+#if defined(__aarch64__) && (defined(__clang__) || (defined(__GNUC__) && __GNUC__ >= 9))
+# define vst1q_u8_x2_unaligned vst1q_u8_x2
+#else
+static HEDLEY_ALWAYS_INLINE void vst1q_u8_x2_unaligned(uint8_t* p, uint8x16x2_t data) {
+	vst1q_u8(p, data.val[0]);
+	vst1q_u8(p+16, data.val[1]);
+}
+#endif
+
+#ifdef YENC_DEC_USE_THINTABLE
+static uint64_t ALIGN_TO(8, compactLUT[256]);
+#else
+# pragma pack(16)
+static struct { char bytes[16]; } ALIGN_TO(16, compactLUT[32768]);
+# pragma pack()
+#endif
+
+static uint8_t eqFixLUT[256];
+
+
+
+static bool neon_vect_is_nonzero(uint8x16_t v) {
+# ifdef __aarch64__
+	return !!(vget_lane_u64(vreinterpret_u64_u32(vqmovn_u64(vreinterpretq_u64_u8(v))), 0));
+# else
+	uint32x4_t tmp1 = vreinterpretq_u32_u8(v);
+	uint32x2_t tmp2 = vorr_u32(vget_low_u32(tmp1), vget_high_u32(tmp1));
+	return !!(vget_lane_u32(vpmax_u32(tmp2, tmp2), 0));
+# endif
+}
+
+
+template<bool isRaw, bool searchEnd>
+HEDLEY_ALWAYS_INLINE void do_decode_neon(const uint8_t* HEDLEY_RESTRICT src, long& len, unsigned char* HEDLEY_RESTRICT & p, unsigned char& escFirst, uint16_t& nextMask) {
+	HEDLEY_ASSUME(escFirst == 0 || escFirst == 1);
+	HEDLEY_ASSUME(nextMask == 0 || nextMask == 1 || nextMask == 2);
+	uint8x16_t yencOffset = escFirst ? vmakeq_u8(42+64,42,42,42,42,42,42,42,42,42,42,42,42,42,42,42) : vdupq_n_u8(42);
+#ifdef __aarch64__
+	uint8x16_t nextMaskMix = vdupq_n_u8(0);
+	if(nextMask == 1)
+		nextMaskMix = vsetq_lane_u8(1, nextMaskMix, 0);
+	if(nextMask == 2)
+		nextMaskMix = vsetq_lane_u8(2, nextMaskMix, 1);
+#else
+	uint8x16_t lfCompare = vdupq_n_u8('\n');
+	if(isRaw) {
+		if(nextMask == 1)
+			lfCompare = vsetq_lane_u8('.', lfCompare, 0);
+		if(nextMask == 2)
+			lfCompare = vsetq_lane_u8('.', lfCompare, 1);
+	}
+#endif
+	long i;
+	for(i = -len; i; i += sizeof(uint8x16_t)*2) {
+		uint8x16x2_t data = vld1q_u8_x2_align(src+i, 32);
+		uint8x16_t dataA = data.val[0];
+		uint8x16_t dataB = data.val[1];
+		
+		// search for special chars
+		uint8x16_t cmpEqA = vceqq_u8(dataA, vdupq_n_u8('=')),
+		cmpEqB = vceqq_u8(dataB, vdupq_n_u8('=')),
+#ifdef __aarch64__
+		cmpA = vqtbx1q_u8(
+			cmpEqA,
+			//                             \n      \r
+			vmakeq_u8(0,0,0,0,0,0,0,0,0,0,255,0,0,255,0,0),
+			dataA
+		),
+		cmpB = vqtbx1q_u8(
+			cmpEqB,
+			vmakeq_u8(0,0,0,0,0,0,0,0,0,0,255,0,0,255,0,0),
+			dataB
+		);
+		if(isRaw) cmpA = vorrq_u8(cmpA, nextMaskMix);
+#else
+		cmpCrA = vceqq_u8(dataA, vdupq_n_u8('\r')),
+		cmpCrB = vceqq_u8(dataB, vdupq_n_u8('\r')),
+		cmpA = vorrq_u8(
+			vorrq_u8(
+				cmpCrA,
+				vceqq_u8(dataA, lfCompare)
+			),
+			cmpEqA
+		),
+		cmpB = vorrq_u8(
+			vorrq_u8(
+				cmpCrB,
+				vceqq_u8(dataB, vdupq_n_u8('\n'))
+			),
+			cmpEqB
+		);
+#endif
+		
+		
+#ifdef __aarch64__
+		if (LIKELIHOOD(0.42 /*guess*/, neon_vect_is_nonzero(vorrq_u8(cmpA, cmpB)))) {
+			cmpA = vandq_u8(cmpA, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128));
+			cmpB = vandq_u8(cmpB, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128));
+			uint8x16_t cmpMerge = vpaddq_u8(cmpA, cmpB);
+			uint8x16_t cmpEqMerge = vpaddq_u8(
+				vandq_u8(cmpEqA, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128)),
+				vandq_u8(cmpEqB, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128))
+			);
+			
+			uint8x16_t cmpCombined = vpaddq_u8(cmpMerge, cmpEqMerge);
+			cmpCombined = vpaddq_u8(cmpCombined, cmpCombined);
+			uint8x8_t cmpPacked = vget_low_u8(cmpCombined);
+			uint32_t mask = vgetq_lane_u32(vreinterpretq_u32_u8(cmpCombined), 0);
+			uint32_t maskEq = vgetq_lane_u32(vreinterpretq_u32_u8(cmpCombined), 1);
+#else
+		cmpA = vandq_u8(cmpA, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128));
+		cmpB = vandq_u8(cmpB, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128));
+		// no vpaddq_u8 in ARMv7, so need extra 64-bit VPADD
+		uint8x8_t cmpPacked = vpadd_u8(
+			vpadd_u8(
+				vget_low_u8(cmpA), vget_high_u8(cmpA)
+			),
+			vpadd_u8(
+				vget_low_u8(cmpB), vget_high_u8(cmpB)
+			)
+		);
+		cmpPacked = vpadd_u8(cmpPacked, cmpPacked);
+		uint32_t mask = vget_lane_u32(vreinterpret_u32_u8(cmpPacked), 0);
+		if(LIKELIHOOD(0.42, mask != 0)) {
+			uint8x16_t cmpEqMaskedA = vandq_u8(cmpEqA, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128));
+			uint8x16_t cmpEqMaskedB = vandq_u8(cmpEqB, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128));
+			uint8x8_t cmpEqPacked = vpadd_u8(
+				vpadd_u8(
+					vget_low_u8(cmpEqMaskedA), vget_high_u8(cmpEqMaskedA)
+				),
+				vpadd_u8(
+					vget_low_u8(cmpEqMaskedB), vget_high_u8(cmpEqMaskedB)
+				)
+			);
+			cmpEqPacked = vpadd_u8(cmpEqPacked, cmpEqPacked);
+			uint32_t maskEq = vget_lane_u32(vreinterpret_u32_u8(cmpEqPacked), 0);
+#endif
+			
+			// handle \r\n. sequences
+			// RFC3977 requires the first dot on a line to be stripped, due to dot-stuffing
+			if((isRaw || searchEnd) && LIKELIHOOD(0.15, mask != maskEq)) {
+				// vext seems to be a cheap operation on ARM, relative to loads, so only avoid it if there's only one load (isRaw only)
+				uint8x16_t tmpData2, nextData;
+				if(isRaw && !searchEnd) {
+					tmpData2 = vld1q_u8(src+i + 2 + sizeof(uint8x16_t));
+				} else {
+					nextData = vld1q_u8_align(src+i + sizeof(uint8x16_t)*2, 16); // only 32-bits needed, but there doesn't appear a nice way to do this via intrinsics: https://stackoverflow.com/questions/46910799/arm-neon-intrinsics-convert-d-64-bit-register-to-low-half-of-q-128-bit-regis
+					tmpData2 = vextq_u8(dataB, nextData, 2);
+				}
+#ifdef __aarch64__
+				uint8x16_t cmpCrA = vceqq_u8(dataA, vdupq_n_u8('\r'));
+				uint8x16_t cmpCrB = vceqq_u8(dataB, vdupq_n_u8('\r'));
+# define NEXT_DATA(n) vextq_u8(dataB, nextData, n)
+#else
+// on ARMv7, prefer loading over VEXT to avoid holding onto nextData reference; this reduces register spills. Shouldn't be an issue on ARMv8 due to 32x 128-bit registers
+# define NEXT_DATA(n) vld1q_u8(src+i + n+sizeof(uint8x16_t))
+#endif
+				uint8x16_t match2EqA, match2Cr_DotA;
+				uint8x16_t match2EqB, match2Cr_DotB;
+				if(searchEnd) {
+					match2EqB = vceqq_u8(tmpData2, vdupq_n_u8('='));
+				}
+				if(isRaw) {
+					match2Cr_DotA = vandq_u8(cmpCrA, vceqq_u8(vextq_u8(dataA, dataB, 2), vdupq_n_u8('.')));
+					match2Cr_DotB = vandq_u8(cmpCrB, vceqq_u8(tmpData2, vdupq_n_u8('.')));
+				}
+				
+				// find patterns of \r_.
+				if(isRaw && LIKELIHOOD(0.001, neon_vect_is_nonzero(
+					vorrq_u8(match2Cr_DotA, match2Cr_DotB)
+				))) {
+					uint8x16_t match1LfA = vceqq_u8(vextq_u8(dataA, dataB, 1), vdupq_n_u8('\n'));
+					uint8x16_t match1LfB;
+					if(searchEnd)
+						match1LfB = vceqq_u8(NEXT_DATA(1), vdupq_n_u8('\n'));
+					else
+						match1LfB = vceqq_u8(vld1q_u8(src+i + 1+sizeof(uint8x16_t)), vdupq_n_u8('\n'));
+					// merge matches of \r_. with those for \n
+					uint8x16_t match2NlDotA = vandq_u8(match2Cr_DotA, match1LfA);
+					uint8x16_t match2NlDotB = vandq_u8(match2Cr_DotB, match1LfB);
+					if(searchEnd) {
+						uint8x16_t match1NlA = vandq_u8(match1LfA, cmpCrA);
+						uint8x16_t match1NlB = vandq_u8(match1LfB, cmpCrB);
+						
+						uint8x16_t tmpData3 = NEXT_DATA(3);
+						uint8x16_t tmpData4 = NEXT_DATA(4);
+						// match instances of \r\n.\r\n and \r\n.=y
+						uint8x16_t match3CrB = vceqq_u8(tmpData3, vdupq_n_u8('\r'));
+						uint8x16_t match4LfB = vceqq_u8(tmpData4, vdupq_n_u8('\n'));
+						uint8x16_t match4Nl = vbslq_u8(vdupq_n_u8('\r'), // exact VBSL vector doesn't matter, so reuse the '\r' vector
+							vextq_u8(match1NlA, match1NlB, 3),
+							vandq_u8(match3CrB, match4LfB)
+						);
+						uint8x16_t match4EqY = vbslq_u8(vdupq_n_u8('\r'),
+							// match =y
+							vreinterpretq_u8_u16(vceqq_u16(vreinterpretq_u16_u8(vextq_u8(dataA, dataB, 4)), vdupq_n_u16(0x793d))),
+							vreinterpretq_u8_u16(vceqq_u16(vreinterpretq_u16_u8(tmpData4), vdupq_n_u16(0x793d)))
+						);
+						
+						match2EqA = vextq_u8(cmpEqA, cmpEqB, 2);
+						uint8x16_t match3EqY = vbslq_u8(vdupq_n_u8('\r'),
+							vandq_u8(
+								vceqq_u8(vextq_u8(dataA, dataB, 3), vdupq_n_u8('y')),
+								match2EqA
+							), vandq_u8(
+								vceqq_u8(tmpData3, vdupq_n_u8('y')),
+								match2EqB
+							)
+						);
+						// merge \r\n and =y matches for tmpData4
+						uint8x16_t match4End = vorrq_u8(
+							match4Nl,
+							vreinterpretq_u8_u16(vsriq_n_u16(vreinterpretq_u16_u8(match4EqY), vreinterpretq_u16_u8(match3EqY), 8))
+						);
+						// merge with \r\n.
+						uint8x16_t match2NlDot = vbslq_u8(vdupq_n_u8('\r'), match2NlDotA, match2NlDotB);
+						match4End = vandq_u8(match4End, match2NlDot);
+						// match \r\n=y
+						uint8x16_t match1Nl = vbslq_u8(vdupq_n_u8('\r'), match1NlA, match1NlB);
+						uint8x16_t match3End = vandq_u8(match3EqY, match1Nl);
+						// combine match sequences
+						if(LIKELIHOOD(0.001, neon_vect_is_nonzero(vorrq_u8(match4End, match3End)))) {
+							// terminator found
+							// there's probably faster ways to do this, but reverting to scalar code should be good enough
+							len += i;
+							break;
+						}
+					}
+#ifdef __aarch64__
+					uint8x16_t match2NlDotBMasked = vandq_u8(match2NlDotB, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128));
+					uint8x16_t mergeKillDots = vpaddq_u8(
+						vandq_u8(match2NlDotA, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128)),
+						match2NlDotBMasked
+					);
+					uint8x8_t mergeKillDots2 = vget_low_u8(vpaddq_u8(mergeKillDots, mergeKillDots));
+#else
+					uint8x16_t match2NlDotMaskedA = vandq_u8(match2NlDotA, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128));
+					uint8x16_t match2NlDotMaskedB = vandq_u8(match2NlDotB, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128));
+					uint8x8_t mergeKillDots2 = vpadd_u8(
+						vpadd_u8(
+							vget_low_u8(match2NlDotMaskedA), vget_high_u8(match2NlDotMaskedA)
+						),
+						vpadd_u8(
+							vget_low_u8(match2NlDotMaskedB), vget_high_u8(match2NlDotMaskedB)
+						)
+					);
+#endif
+					mergeKillDots2 = vpadd_u8(mergeKillDots2, mergeKillDots2);
+					uint32x2_t mergeKillDotsShifted = vshl_n_u32(vreinterpret_u32_u8(mergeKillDots2), 2);
+					mask |= vget_lane_u32(mergeKillDotsShifted, 0);
+					cmpPacked = vorr_u8(cmpPacked, vreinterpret_u8_u32(mergeKillDotsShifted));
+#ifdef __aarch64__
+					nextMaskMix = vextq_u8(match2NlDotB, vdupq_n_u8(0), 14);
+#else
+					lfCompare = vcombine_u8(vbsl_u8(
+						vext_u8(vget_high_u8(match2NlDotB), vdup_n_u8('\n'), 6),
+						vdup_n_u8('.'),
+						vget_high_u8(lfCompare)
+					), vget_high_u8(lfCompare));
+#endif
+				} else if(searchEnd) {
+					match2EqA = vextq_u8(cmpEqA, cmpEqB, 2);
+					uint8x16_t match3EqYA = vandq_u8(match2EqA, vceqq_u8(vextq_u8(dataA, dataB, 3), vdupq_n_u8('y')));
+					uint8x16_t match3EqYB = vandq_u8(match2EqB, vceqq_u8(NEXT_DATA(3), vdupq_n_u8('y')));
+					if(LIKELIHOOD(0.001, neon_vect_is_nonzero(vorrq_u8(
+						match3EqYA, match3EqYB
+					)))) {
+						uint8x16_t match1LfA = vceqq_u8(vextq_u8(dataA, dataB, 1), vdupq_n_u8('\n'));
+						uint8x16_t match1LfB = vceqq_u8(NEXT_DATA(1), vdupq_n_u8('\n'));
+						uint8x16_t matchEnd = vorrq_u8(
+							vandq_u8(match3EqYA, vandq_u8(match1LfA, cmpCrA)),
+							vandq_u8(match3EqYB, vandq_u8(match1LfB, cmpCrB))
+						);
+						if(LIKELIHOOD(0.001, neon_vect_is_nonzero(matchEnd))) {
+							len += i;
+							break;
+						}
+					}
+#undef NEXT_DATA
+					if(isRaw)
+#ifdef __aarch64__
+						nextMaskMix = vdupq_n_u8(0);
+#else
+						lfCompare = vcombine_u8(vget_high_u8(lfCompare), vget_high_u8(lfCompare));
+#endif
+				} else if(isRaw) // no \r_. found
+#ifdef __aarch64__
+					nextMaskMix = vdupq_n_u8(0);
+#else
+					lfCompare = vcombine_u8(vget_high_u8(lfCompare), vget_high_u8(lfCompare));
+#endif
+			}
+			
+			// a spec compliant encoder should never generate sequences: ==, =\n and =\r, but we'll handle them to be spec compliant
+			// the yEnc specification requires any character following = to be unescaped, not skipped over, so we'll deal with that
+			// firstly, check for invalid sequences of = (we assume that these are rare, as a spec compliant yEnc encoder should not generate these)
+			if(LIKELIHOOD(0.0001, (mask & ((maskEq << 1) | escFirst)) != 0)) {
+				uint8_t tmp = eqFixLUT[(maskEq&0xff) & ~escFirst];
+				uint32_t maskEq2 = tmp;
+				for(int j=8; j<32; j+=8) {
+					tmp = eqFixLUT[((maskEq>>j)&0xff) & ~(tmp>>7)];
+					maskEq2 |= tmp<<j;
+				}
+				maskEq = maskEq2;
+				
+				// next, eliminate anything following a `=` from the special char mask; this eliminates cases of `=\r` so that they aren't removed
+				maskEq = (maskEq<<1) | escFirst;
+				mask &= ~maskEq;
+				escFirst = tmp>>7;
+				
+				// unescape chars following `=`
+				uint8x8_t maskEqTemp = vreinterpret_u8_u32(vmov_n_u32(maskEq));
+				cmpPacked = vbic_u8(cmpPacked, maskEqTemp); // `mask &= ~maskEq` in vector form
+#ifdef __aarch64__
+				uint8x16_t vMaskEqA = vqtbl1q_u8(
+					vcombine_u8(maskEqTemp, vdup_n_u8(0)),
+					vmakeq_u8(0,0,0,0,0,0,0,0, 1,1,1,1,1,1,1,1)
+				);
+				uint8x16_t vMaskEqB = vqtbl1q_u8(
+					vcombine_u8(maskEqTemp, vdup_n_u8(0)),
+					vmakeq_u8(2,2,2,2,2,2,2,2, 3,3,3,3,3,3,3,3)
+				);
+#else
+				uint8x16_t vMaskEqA = vcombine_u8(
+					vdup_lane_u8(maskEqTemp, 0),
+					vdup_lane_u8(maskEqTemp, 1)
+				);
+				uint8x16_t vMaskEqB = vcombine_u8(
+					vdup_lane_u8(maskEqTemp, 2),
+					vdup_lane_u8(maskEqTemp, 3)
+				);
+#endif
+				vMaskEqA = vtstq_u8(vMaskEqA, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128));
+				vMaskEqB = vtstq_u8(vMaskEqB, vmakeq_u8(1,2,4,8,16,32,64,128, 1,2,4,8,16,32,64,128));
+				
+				dataA = vsubq_u8(
+					dataA,
+					vbslq_u8(vMaskEqA, vdupq_n_u8(64+42), vdupq_n_u8(42))
+				);
+				dataB = vsubq_u8(
+					dataB,
+					vbslq_u8(vMaskEqB, vdupq_n_u8(64+42), vdupq_n_u8(42))
+				);
+			} else {
+				// no invalid = sequences found - we can cut out some things from above
+				// this code path is a shortened version of above; it's here because it's faster, and what we'll be dealing with most of the time
+				escFirst = (maskEq >> 31);
+				
+				dataA = vsubq_u8(
+					dataA,
+					vbslq_u8(
+						vextq_u8(vdupq_n_u8(42), cmpEqA, 15),
+						vdupq_n_u8(64+42),
+						yencOffset
+					)
+				);
+				dataB = vsubq_u8(
+					dataB,
+					vbslq_u8(
+						vextq_u8(cmpEqA, cmpEqB, 15),
+						vdupq_n_u8(64+42),
+						vdupq_n_u8(42)
+					)
+				);
+			}
+			yencOffset = vsetq_lane_u8((escFirst << 6) | 42, yencOffset, 0);
+			
+			// all that's left is to 'compress' the data (skip over masked chars)
+			uint32_t counts = 0x08080808 - vget_lane_u32(vreinterpret_u32_u8(vcnt_u8(cmpPacked)), 0);
+#ifdef __aarch64__
+			counts += counts >> 8;
+			vst1q_u8(p, vqtbl1q_u8(
+				dataA,
+				vld1q_u8_align((uint8_t*)(compactLUT + (mask&0x7fff)), 16)
+			));
+			p += counts & 0xff;
+			mask >>= 16;
+			vst1q_u8(p, vqtbl1q_u8(
+				dataB,
+				vld1q_u8_align((uint8_t*)(compactLUT + (mask&0x7fff)), 16)
+			));
+			p += (counts>>16) & 0xff;
+#else
+			// lookup compress masks and shuffle
+			vst1_u8(p, vtbl1_u8(
+				vget_low_u8(dataA),
+				vld1_u8_align((uint8_t*)(compactLUT + (mask&0xff)), 8)
+			));
+			p += counts & 0xff;
+			mask >>= 8;
+			vst1_u8(p, vtbl1_u8(
+				vget_high_u8(dataA),
+				vld1_u8_align((uint8_t*)(compactLUT + (mask&0xff)), 8)
+			));
+			p += (counts>>8) & 0xff;
+			mask >>= 8;
+			vst1_u8(p, vtbl1_u8(
+				vget_low_u8(dataB),
+				vld1_u8_align((uint8_t*)(compactLUT + (mask&0xff)), 8)
+			));
+			p += (counts>>16) & 0xff;
+			mask >>= 8;
+			vst1_u8(p, vtbl1_u8(
+				vget_high_u8(dataB),
+				vld1_u8_align((uint8_t*)(compactLUT + (mask&0xff)), 8)
+			));
+			p += (counts>>24) & 0xff;
+			
+#endif
+			
+		} else {
+			dataA = vsubq_u8(dataA, yencOffset);
+			dataB = vsubq_u8(dataB, vdupq_n_u8(42));
+			vst1q_u8_x2_unaligned(p, vcreate2_u8(dataA, dataB));
+			p += sizeof(uint8x16_t)*2;
+			escFirst = 0;
+#ifdef __aarch64__
+			yencOffset = vdupq_n_u8(42);
+#else
+			yencOffset = vcombine_u8(vdup_n_u8(42), vget_high_u8(yencOffset));
+#endif
+		}
+	}
+	
+	if(isRaw) {
+		if(len != 0) { // have to gone through at least one loop cycle
+			if(src[i-2] == '\r' && src[i-1] == '\n' && src[i] == '.')
+				nextMask = 1;
+			else if(src[i-1] == '\r' && src[i] == '\n' && src[i+1] == '.')
+				nextMask = 2;
+			else
+				nextMask = 0;
+		}
+	} else
+		nextMask = 0;
+}
+
+void decoder_set_neon_funcs() {
+	decoder_init_lut(eqFixLUT, compactLUT);
+	_do_decode = &do_decode_simd<false, false, sizeof(uint8x16_t)*2, do_decode_neon<false, false> >;
+	_do_decode_raw = &do_decode_simd<true, false, sizeof(uint8x16_t)*2, do_decode_neon<true, false> >;
+	_do_decode_end_raw = &do_decode_simd<true, true, sizeof(uint8x16_t)*2, do_decode_neon<true, true> >;
+}
+#else
+void decoder_set_neon_funcs() {}
+#endif