Почему load_ps () работает на одном ПК, а не на другом? - PullRequest
Новая
       68

Почему load_ps () работает на одном ПК, а не на другом?

0 голосов
/ 22 мая 2018

Я написал следующий код для масштабирования набора чисел: 

 #include <stdio.h>
 #include <stdlib.h>
 #include <math.h>
 #include "immintrin.h"

 void scale(struct problem_param prob_param, float* features)
 {
    int i,j,k;
    for (j = 0; j < prob_param.Nr_ft; j++)
    {
        __m256      range_vec ,low_up_vec , low_vec,tmp_vec;
        __m256      feat_min_vec, feat_vec;
        unsigned    count           = prob_param.Size;
        unsigned    offset          = j * prob_param.Size;
        float       feature_max     = features[offset];
        float       feature_min     = features[offset];
        /*
         * Look for min and max of each feature.
         */
        for ( i = 1; i < prob_param.Size ; i++)
        {
            if (features[i + offset] > feature_max )        feature_max = features[i + offset];

            if (features[i + offset] < feature_min )    feature_min = features[i + offset];
        }
        printf("feature : %u \t min = %f \t max = %f \n",j,feature_min,feature_max);
        /*
         * Set the range.
         * Set constant vectors for the vector instructions.
         */
        float       range   = feature_max - feature_min;
        feat_min_vec        = _mm256_set1_ps (feature_min);
        range_vec           = _mm256_set1_ps (range);
        low_up_vec          = _mm256_set1_ps (prob_param.upper_limit - prob_param.lower_limit);
        low_vec             = _mm256_set1_ps (prob_param.lower_limit);
        /*
         * Normalising
         * -----------
         * Head
         */
        for ( i = 0; i < prob_param.Size && count >= 7 ; i+=8)
        {
            feat_vec    = _mm256_load_ps(&features[i + offset]);
            tmp_vec     = _mm256_sub_ps(feat_vec,feat_min_vec);
            tmp_vec     = _mm256_mul_ps(tmp_vec,low_up_vec);
            tmp_vec     = _mm256_div_ps(tmp_vec,range_vec);
            feat_vec    = _mm256_add_ps(tmp_vec,low_vec);

            _mm256_store_ps (&features[i + offset], feat_vec);

            count -=8;
        }
        /*
         * Normalising
         * -----------
         * Tail
         */
        for ( k = i; k < prob_param.Size ; k++)
        {
            features[k + offset] = prob_param.lower_limit + (prob_param.upper_limit - prob_param.lower_limit) * (features[k + offset] - feature_min) / range;
        }
    }

Это функция, которая заботится о масштабировании, и я называю это так: 

#include <stdio.h>
#include <stdlib.h>

#include "data.h" 
#include "common.h"

#define     training_size       3089
#define     number_features     4
#define     low                 -1.0
#define     up                  1.0
float* feature_array;

int main()
{
  struct problem_param pp;

  pp.Size           =       training_size;
  pp.Nr_ft          =       number_features;
  pp.lower_limit    =       low;
  pp.upper_limit    =       up;

posix_memalign((void **) &feature_array, 32, (size_t) training_size * number_features *sizeof(float));

scale(pp,feature_array);

return EXIT_SUCCESS;
}

Я тестировал этот код на своем MacBook Pro Core i5 Haswell, и он работает, но когда я тестирую его на ASUS Core I7 Haswell, он показывает ошибку сегментации для загрузки.Я что-то упустил?

1 Ответ

0 голосов
/ 22 мая 2018

Значение offset (и, следовательно, i + offset) не всегда будет кратным 8 (в приведенном выше примере оно равно 0, 3089, 6178, 9267), поэтому ваши встроенные функции загрузки и хранения будутв общем случае не выровняйте.

Самое простое решение - использовать _mm256_loadu_ps вместо _mm256_load_ps и _mm256_storeu_ps вместо _mm256_store_ps.

.работать на вашем MacBook Pro, я предполагаю, что clang генерирует невыровненные инструкции по загрузке / хранению за вашей спиной, тем самым скрывая проблему, пока вы не попытаетесь запустить код в системе с другим компилятором.

Обновление: Iтолько что подтвердил приведенную выше гипотезу, скомпилировав и разобрав сгенерированный код (на Haswell MacBook Pro с macOS 10.13.4 и Xcode 9.3.1): 

>>> vmovups (%r14,%r13,4), %ymm0
    vsubps  192(%rsp), %ymm0, %ymm0 ## 32-byte Folded Reload
    vmulps  448(%rsp), %ymm0, %ymm0 ## 32-byte Folded Reload
    vdivps  384(%rsp), %ymm0, %ymm0 ## 32-byte Folded Reload
    vaddps  416(%rsp), %ymm0, %ymm0 ## 32-byte Folded Reload
>>> vmovups %ymm0, (%r14,%r13,4)

Обратите внимание на использование vmovups вместо vmovaps.

...