有各种说法认为,在某些情况下,Fortran 可以比 C 更快,例如在别名方面,我经常听说它比 C 更好地进行自动矢量化(请参阅此处进行一些很好的讨论)。
但是,对于像计算斐波那契数和 Mandelbrot 这样的简单函数,使用直接的解决方案,没有任何技巧和编译器的额外提示/关键字,我希望它们确实执行相同的操作。
C实现:
int fib(int n) {
return n < 2 ? n : fib(n-1) + fib(n-2);
}
int mandel(double complex z) {
int maxiter = 80;
double complex c = z;
for (int n=0; n<maxiter; ++n) {
if (cabs(z) > 2.0) {
return n;
}
z = z*z+c;
}
return maxiter;
}
Fortran 实现:
integer, parameter :: dp=kind(0.d0) ! double precision
integer recursive function fib(n) result(r)
integer, intent(in) :: n
if (n < 2) then
r = n
else
r = fib(n-1) + fib(n-2)
end if
end function
integer function mandel(z0) result(r)
complex(dp), intent(in) :: z0
complex(dp) :: c, z
integer :: n, maxiter
maxiter = 80
z = z0
c = z0
do n = 1, maxiter
if (abs(z) > 2) then
r = n-1
return
end if
z = z**2 + c
end do
r = maxiter
end function
朱莉娅实现:
fib(n) = n < 2 ? n : fib(n-1) + fib(n-2)
function mandel(z)
c = z
maxiter = 80
for n = 1:maxiter
if abs(z) > 2
return n-1
end
z = z^2 + c
end
return maxiter
end
(可以在这里找到包括其他基准函数在内的完整代码。)
根据Julia 主页,Julia 和 Fortran (with ) 在这两个函数上的-O3
性能优于 C (with )。-O3
这个怎么可能?