Test example MSA 14

Author: nicodr97

scripts/generate_trimmed_msas.sh

@@ -6,6 +6,10 @@ mkdir test_msas/automated1
 mkdir test_msas/nogaps
 mkdir test_msas/noallgaps
 
+
+./bin/trimal -in dataset/example.014.AA.EggNOG.COG0591.fasta -strict -out test_msas/strict/example.014.AA.EggNOG.COG0591.fasta
+
+'
 ./bin/trimal -in dataset/example.004.AA.fasta -gappyout -out test_msas/gappyout/example.004.AA.fasta
 ./bin/trimal -in dataset/example.004.AA.fasta -strict -out test_msas/strict/example.004.AA.fasta
 ./bin/trimal -in dataset/example.004.AA.fasta -strictplus -out test_msas/strictplus/example.004.AA.fasta
@@ -510,3 +514,4 @@ mkdir test_msas/noallgaps
 ./bin/trimal -in dataset/example.091.AA.strNOG.ENOG411BWBU.fasta -automated1 -out test_msas/automated1/example.091.AA.strNOG.ENOG411BWBU.fasta
 ./bin/trimal -in dataset/example.091.AA.strNOG.ENOG411BWBU.fasta -nogaps -out test_msas/nogaps/example.091.AA.strNOG.ENOG411BWBU.fasta
 ./bin/trimal -in dataset/example.091.AA.strNOG.ENOG411BWBU.fasta -noallgaps -out test_msas/noallgaps/example.091.AA.strNOG.ENOG411BWBU.fasta
+'

source/Platform/x86/AVX2.cpp

@@ -109,12 +109,22 @@ class AVX2Vector {
   inline uint16_t sum() const {
     __m256i sum256 = _mm256_sad_epu8(vector, _mm256_setzero_si256());
     __m128i sum128 = _mm_add_epi64(_mm256_extractf128_si256(sum256, 1), _mm256_castsi256_si128(sum256));
+    uint16_t res32 = _mm_extract_epi32(sum128, 0) + _mm_extract_epi32(sum128, 2);
+    uint16_t res64 = _mm_extract_epi64(sum128, 0) + _mm_extract_epi64(sum128, 1);
+    if (res32 != res64) {
+      std::cout << "uint16_t res32 = " << res32 << "\n";
+      std::cout << "uint16_t res64 = " << res64 << "\n";
+      exit (1);
+    }
+
+
     return _mm_extract_epi64(sum128, 0) + _mm_extract_epi64(sum128, 1);
   }
 #else
   inline uint16_t sum() const {
     __m256i sum256 = _mm256_sad_epu8(vector, _mm256_setzero_si256());
     __m128i sum128 = _mm_add_epi64(_mm256_extractf128_si256(sum256, 1), _mm256_castsi256_si128(sum256));
+     std::cout << "NOT CPU_FEATURES_ARCH_X86_64 \n";
 
     return _mm_extract_epi32(sum128, 0) + _mm_extract_epi32(sum128, 2);
   }