72 num_points = (num_points > USHRT_MAX) ? USHRT_MAX : num_points;
73 const uint32_t eighthPoints = num_points / 8;
75 float* inputPtr = (
float*)source;
77 __m256 indexIncrementValues = _mm256_set1_ps(8);
78 __m256 currentIndexes = _mm256_set_ps(-1, -2, -3, -4, -5, -6, -7, -8);
80 float min = source[0];
82 __m256 minValues = _mm256_set1_ps(min);
83 __m256 minValuesIndex = _mm256_setzero_ps();
84 __m256 compareResults;
90 for (uint32_t number = 0; number < eighthPoints; number++) {
92 currentValues = _mm256_load_ps(inputPtr);
94 currentIndexes = _mm256_add_ps(currentIndexes, indexIncrementValues);
96 compareResults = _mm256_cmp_ps(currentValues, minValues, _CMP_LT_OS);
98 minValuesIndex = _mm256_blendv_ps(minValuesIndex, currentIndexes, compareResults);
99 minValues = _mm256_blendv_ps(minValues, currentValues, compareResults);
103 _mm256_store_ps(minValuesBuffer, minValues);
104 _mm256_store_ps(minIndexesBuffer, minValuesIndex);
106 for (uint32_t number = 0; number < 8; number++) {
107 if (minValuesBuffer[number] < min) {
108 index = minIndexesBuffer[number];
109 min = minValuesBuffer[number];
110 }
else if (minValuesBuffer[number] == min) {
111 if (index > minIndexesBuffer[number])
112 index = minIndexesBuffer[number];
116 for (uint32_t number = eighthPoints * 8; number < num_points; number++) {
117 if (source[number] < min) {
119 min = source[number];
122 target[0] = (uint16_t)index;
140 __m128 currentIndexes = _mm_set_ps(-1, -2, -3, -4);
197 num_points = (num_points > USHRT_MAX) ? USHRT_MAX : num_points;
198 const uint32_t quarterPoints = num_points / 4;
200 float* inputPtr = (
float*)source;
202 __m128 indexIncrementValues = _mm_set1_ps(4);
203 __m128 currentIndexes = _mm_set_ps(-1, -2, -3, -4);
205 float min = source[0];
207 __m128 minValues = _mm_set1_ps(min);
208 __m128 minValuesIndex = _mm_setzero_ps();
209 __m128 compareResults;
210 __m128 currentValues;
215 for (uint32_t number = 0; number < quarterPoints; number++) {
217 currentValues = _mm_load_ps(inputPtr);
219 currentIndexes = _mm_add_ps(currentIndexes, indexIncrementValues);
221 compareResults = _mm_cmplt_ps(currentValues, minValues);
223 minValuesIndex = _mm_or_ps(_mm_and_ps(compareResults, currentIndexes),
224 _mm_andnot_ps(compareResults, minValuesIndex));
225 minValues = _mm_or_ps(_mm_and_ps(compareResults, currentValues),
226 _mm_andnot_ps(compareResults, minValues));
230 _mm_store_ps(minValuesBuffer, minValues);
231 _mm_store_ps(minIndexesBuffer, minValuesIndex);
233 for (uint32_t number = 0; number < 4; number++) {
234 if (minValuesBuffer[number] < min) {
235 index = minIndexesBuffer[number];
236 min = minValuesBuffer[number];
237 }
else if (minValuesBuffer[number] == min) {
238 if (index > minIndexesBuffer[number])
239 index = minIndexesBuffer[number];
243 for (uint32_t number = quarterPoints * 4; number < num_points; number++) {
244 if (source[number] < min) {
246 min = source[number];
249 target[0] = (uint16_t)index;
294 num_points = (num_points > USHRT_MAX) ? USHRT_MAX : num_points;
295 const uint32_t eighthPoints = num_points / 8;
297 float* inputPtr = (
float*)source;
299 __m256 indexIncrementValues = _mm256_set1_ps(8);
300 __m256 currentIndexes = _mm256_set_ps(-1, -2, -3, -4, -5, -6, -7, -8);
302 float min = source[0];
304 __m256 minValues = _mm256_set1_ps(min);
305 __m256 minValuesIndex = _mm256_setzero_ps();
306 __m256 compareResults;
307 __m256 currentValues;
312 for (uint32_t number = 0; number < eighthPoints; number++) {
314 currentValues = _mm256_loadu_ps(inputPtr);
316 currentIndexes = _mm256_add_ps(currentIndexes, indexIncrementValues);
318 compareResults = _mm256_cmp_ps(currentValues, minValues, _CMP_LT_OS);
320 minValuesIndex = _mm256_blendv_ps(minValuesIndex, currentIndexes, compareResults);
321 minValues = _mm256_blendv_ps(minValues, currentValues, compareResults);
325 _mm256_storeu_ps(minValuesBuffer, minValues);
326 _mm256_storeu_ps(minIndexesBuffer, minValuesIndex);
328 for (uint32_t number = 0; number < 8; number++) {
329 if (minValuesBuffer[number] < min) {
330 index = minIndexesBuffer[number];
331 min = minValuesBuffer[number];
332 }
else if (minValuesBuffer[number] == min) {
333 if (index > minIndexesBuffer[number])
334 index = minIndexesBuffer[number];
338 for (uint32_t number = eighthPoints * 8; number < num_points; number++) {
339 if (source[number] < min) {
341 min = source[number];
344 target[0] = (uint16_t)index;