Не выровненный доступ по выровненному указателю?

Question

У меня есть этот код функции:

bool interpolate(const Mat &im, float ofsx, float ofsy, float a11, float a12, float a21, float a22, Mat &res)
{         
   bool ret = false;
   // input size (-1 for the safe bilinear interpolation)
   const int width = im.cols-1;
   const int height = im.rows-1;
   // output size
   const int halfWidth  = res.cols >> 1;
   const int halfHeight = res.rows >> 1;
   float *out = res.ptr<float>(0);
   const float *imptr  = im.ptr<float>(0);
   for (int j=-halfHeight; j<=halfHeight; ++j)
   {
      const float rx = ofsx + j * a12;
      const float ry = ofsy + j * a22;
      #pragma omp simd
      for(int i=-halfWidth; i<=halfWidth; ++i, out++)
      {
         float wx = rx + i * a11;
         float wy = ry + i * a21;
         const int x = (int) floor(wx);
         const int y = (int) floor(wy);
         if (x >= 0 && y >= 0 && x < width && y < height)
         {
            // compute weights
            wx -= x; wy -= y;
            int rowOffset = y*im.cols;
            int rowOffset1 = (y+1)*im.cols;
            // bilinear interpolation
            *out =
                (1.0f - wy) * ((1.0f - wx) * imptr[rowOffset+x]   + wx * imptr[rowOffset+x+1]) +
                (       wy) * ((1.0f - wx) * imptr[rowOffset1+x] + wx * imptr[rowOffset1+x+1]);
         } else {
            *out = 0;
            ret =  true; // touching boundary of the input            
         }
      }
   }
   return ret;
}

Который является частью кода, скомпилированного с icpc и следующие варианты:

INTEL_OPT=-O3 -simd -xCORE-AVX2 -parallel -qopenmp -fargument-noalias -ansi-alias -no-prec-div -fp-model fast=2 -fma -align -finline-functions
INTEL_PROFILE=-g -qopt-report=5 -Bdynamic -shared-intel -debug inline-debug-info -qopenmp-link dynamic -parallel-source-info=2 -ldl

В файле отчета об оптимизации (.optr) есть такие строки:

  remark #15389: vectorization support: reference *out has unaligned access   [ /home/luca/Dropbox/HKUST/CloudCache/cloudcache/CloudCache/Descriptors/hesaff/helpers.cpp(259,14) ]
  remark #15389: vectorization support: reference *out has unaligned access   [ /home/luca/Dropbox/HKUST/CloudCache/cloudcache/CloudCache/Descriptors/hesaff/helpers.cpp(263,14) ]
  remark #15389: vectorization support: reference *out has unaligned access   [ /home/luca/Dropbox/HKUST/CloudCache/cloudcache/CloudCache/Descriptors/hesaff/helpers.cpp(263,14) ]

Линии 259 и 263 являются те, где, очевидно, out значение изменено.

Насколько мне известно, это означает, что out не выровнен (больше информации здесь [ https://software.intel.com/en-us/articles/fdiag15126] и здесь).

Тем не менее, это уже странно, потому что, как я уже узнал здесь, cv::Mat объекты уже выровнены. Но это может быть законно, так как компилятор не может знать о выравнивании...Однако сообщение появляется, если я специально выделяю out как выровненный указатель, если я добавлю следующую строку:

out = (float*) _mm_malloc(res.cols*res.rows*sizeof(float), 32);

Этого должно быть достаточно, чтобы сообщить компилятору, что данные выровнены. Я ошибся? В противном случае, как я могу сказать, что out уже выровнен?

Answer 1

Компилятор видит out++, Это не может оставаться на одной линии.

(То есть он остается выровненным для float* доступ, но сообщение относится к векторизации SIMD, которая требует более строгого выравнивания)