27 int nSeg = 16/
sizeof(
T);
29 T * pp1 = (
T*)memalign( 16,
MATRIX_DIM*(queryLen+16)*
sizeof(
T));
31 profile->
len = queryLen;
33 profile->
min = DBL_MAX;
34 profile->
max = DBL_MIN;
50 int segLen = currentBlockSize / nSeg;
55 T * currentProfile = currentBlock+i*currentBlockSize;
56 for( j=0; nSeg*j < currentBlockSize; j++ ){
57 T * tmp = currentProfile + j*nSeg;
58 for( k=0; k<nSeg; k++ )
60 tmp[k] = (
T)matrix[ query[j + k*segLen +
blockStart] * MATRIX_DIM + i ];
61 if(profile->
min > tmp[k]) profile->
min = tmp[k];
62 if(profile->
max < tmp[k]) profile->
max = tmp[k];