c++ - 使用 SIMD 命令的程序中的分段错误

Question

这里有什么问题？当我运行程序时，它说，Segmentation Fault (Core Dumped). 我使用了一些 SIMD 命令。

float function ( Point p1, Point p2, int dim )
{
      int k;
      float result=0.0;
      float *p3;
      p3 = (float*) malloc (16);
      k=dim%4;

      __m128 *v_p1 = (__m128*)p1.coord;
      __m128 *v_p2 = (__m128*)p2.coord;
      __m128 *v_p3 = (__m128*)p3;

      for (int i=0; i<dim/4; i++){
             *v_p3= _mm_sub_ps(*v_p1,*v_p2);
      }
      for(int i=0; i<dim; i++){
             result+=p3[i];
      }
      return(result);
}

Answer 1

正如评论所说，在使用 SIMD 内部函数时，内存中的数据必须对齐（在这种特殊情况下，16 字节对齐），以防您在 UNIX 系统中尝试使用以下方式分配数据posix_memalign()：

http://pubs.opengroup.org/onlinepubs/009695399/functions/posix_memalign.html

Answer 2

任何SIMD _ps指令都需要16字节对齐的数据。据我所知，至少没有正确对齐，所以如果你不使用正确对齐的数据，你p3肯定会得到一个。seg fault我自己无法运行此代码，但如果您__m128按值分配给变量，您应该没问题，因为它们应该正确对齐：

  __m128 v_p1 = _mm_set_ps( ... ); // not sure of the argument 
  __m128 v_p2 = _mm_set_ps( ... ); // not sure of the argument 
  __m128 v_p3 = _mm_set_ps1(p3) ;