@@ -180,10 +180,8 @@ static __m128i* qP_byte (const int8_t* read_num,
180
180
j = i ;
181
181
for (segNum = 0 ; LIKELY (segNum < 16 ) ; segNum ++ ) {
182
182
* t ++ = j >= readLen ? bias : mat [nt * n + read_num [j ]] + bias ;
183
- fprintf (stderr , "%d\t" , * (t - 1 )); // Mengyao
184
183
j += segLen ;
185
184
}
186
- fprintf (stderr , "\n" ); // Mengyao
187
185
}
188
186
}
189
187
return vProfile ;
@@ -209,7 +207,6 @@ static alignment_end* sw_sse2_byte (const int8_t* ref,
209
207
is set to 0, it will not be used */
210
208
uint8_t bias , /* Shift 0 point to a positive value. */
211
209
int32_t maskLen ) {
212
- fprintf (stderr , "in sw_sse2_byte function\n" ); // Mengyao
213
210
214
211
// Put the largest number of the 16 numbers in vm into m.
215
212
#define max16 (m , vm ) (vm) = _mm_max_epu8((vm), _mm_srli_si128((vm), 8)); \
@@ -259,7 +256,6 @@ static alignment_end* sw_sse2_byte (const int8_t* ref,
259
256
step = -1 ;
260
257
}
261
258
for (i = begin ; LIKELY (i != end ); i += step ) {
262
- uint16_t val [8 ]; // Mengyao
263
259
int32_t cmp ;
264
260
__m128i e , vF = vZero , vMaxColumn = vZero ; /* Initialize F value to 0.
265
261
Any errors to vH values will be corrected in the Lazy_F loop.
@@ -284,10 +280,6 @@ static alignment_end* sw_sse2_byte (const int8_t* ref,
284
280
vH = _mm_max_epu8 (vH , e );
285
281
vH = _mm_max_epu8 (vH , vF );
286
282
vMaxColumn = _mm_max_epu8 (vMaxColumn , vH );
287
- memcpy (val , & vMaxColumn , sizeof (val )); // Mengyao
288
- printf ("j: %d\tvMaxColumn1: %i %i %i %i %i %i %i %i \n" , // Mengyao
289
- j , val [0 ], val [1 ], val [2 ], val [3 ], val [4 ], val [5 ],
290
- val [6 ], val [7 ]);
291
283
292
284
/* Save vH values. */
293
285
_mm_store_si128 (pvHStore + j , vH );
@@ -321,11 +313,6 @@ static alignment_end* sw_sse2_byte (const int8_t* ref,
321
313
}
322
314
323
315
end :
324
- memcpy (val , & vMaxColumn , sizeof (val )); // Mengyao
325
- printf ("vMaxColumn2: %i %i %i %i %i %i %i %i \n" ,
326
- val [0 ], val [1 ], val [2 ], val [3 ], val [4 ], val [5 ],
327
- val [6 ], val [7 ]);
328
-
329
316
vMaxScore = _mm_max_epu8 (vMaxScore , vMaxColumn );
330
317
vTemp = _mm_cmpeq_epi8 (vMaxMark , vMaxScore );
331
318
cmp = _mm_movemask_epi8 (vTemp );
@@ -335,7 +322,6 @@ static alignment_end* sw_sse2_byte (const int8_t* ref,
335
322
max16 (temp , vMaxScore );
336
323
vMaxScore = vMaxMark ;
337
324
338
- fprintf (stderr , "temp: %d\n" , temp ); // Mengyao
339
325
if (LIKELY (temp > max )) {
340
326
max = temp ;
341
327
if (max + bias >= 255 ) break ; //overflow
@@ -855,7 +841,6 @@ s_align* ssw_align (const s_profile* prof,
855
841
// Find the alignment scores and ending positions
856
842
if (prof -> profile_byte ) {
857
843
bests = sw_sse2_byte (ref , 0 , refLen , readLen , weight_gapO , weight_gapE , prof -> profile_byte , -1 , prof -> bias , maskLen );
858
- fprintf (stderr , "bests[0].score: %d\n" , bests [0 ].score ); // Mengyao
859
844
if (prof -> profile_word && bests [0 ].score == 255 ) {
860
845
free (bests );
861
846
bests = sw_sse2_word (ref , 0 , refLen , readLen , weight_gapO , weight_gapE , prof -> profile_word , -1 , maskLen );
0 commit comments