Fortran与C:Mandelbrot基准

在该基准站点上获胜的C++版本是 手动矢量化 对于x86,使用simd intrinsics(SSE、AVX或AVX512) ,例如使用 _mm256_movemask_pd(v1 <= v2); 为了得到一个比较结果的整个向量的位掩码,让它检查4个像素是否超出界限。以及用于simd multiple等的gnu c本机向量语法 r2 + i2 将SIMD向量与正常C/C++算子相乘或相乘。

C++版本有一个为SIMD优化的循环条件:

 // Do 50 iterations of mandelbrot calculation for a vector of eight
 // complex values.  Check occasionally to see if the iterated results
 // have wandered beyond the point of no return (> 4.0).

Fortran只是使用OpenMP来实现自动并行化,而编译器的自动矢量化并不会创造出比手动调节循环条件更好的效果,因为手动调节循环条件可以继续执行源代码没有的冗余工作(因为这比更频繁地检查要便宜)。

有很多C和C++版本的程序和FORTRAN版本有类似的速度。 它们甚至对于C/C++源并不是手动矢量化的。

我不确定英特尔Fortran或其他编译器是否支持手动矢量化扩展。