为了比较Julia和C / C ++的性能,我编写了一个非常简单的程序。基本上,它返回的总质数小于给定数。问题是我的平均时间比c / c ++的平均时间差2.5倍(AFAIK,它应该稍微接近1.0倍)。我正在使用Julia 0.5.0(但我也在0.6.0上进行了测试,得到了相同的结果)。我发现更改%
为mod
进行模运算会有所帮助。另外,在任何地方设置类型都会给我带来一些好处。
有人可以帮我弄清楚在此示例中我还缺少哪些东西来改善Julia的表现吗?提前致谢。
module modIsPrime
export isprime, Test_isprime
@inline function isprime( nNumber :: Int64 ) :: Int64
s = Int( trunc( sqrt( nNumber ) ) )
if mod( nNumber, 2 ) == 0
return 0
else
for i = 3 : 2 : s
if mod( nNumber, i ) == 0
return 0
end
end
return 1
end
end
function Test_isprime( )
nPrimeNumbers :: Int64 = 0
n :: Int64 = 0
for n = 1 : 2 : 16000000
nPrimeNumbers += isprime( n )
end
display( [ "Prime numbers: " nPrimeNumbers ] )
end
end
编辑-正如罗恩提到的,这是我用来进行比较的C ++代码:
#include "stdafx.h"
#include <chrono>
#include <iostream>
#include <omp.h>
#include <math.h>
using namespace std;
using namespace chrono;
#define DEFAULT_MAX_TESTS 16'000'000
inline int isprime( unsigned long nNumber ) {
unsigned long
i
, s = static_cast< unsigned long >( sqrt( static_cast< double >( nNumber ) ) )
;
if ( nNumber % 2 == 0 ) {
return 0;
} else {
for( i = 3; i <= s; i += 2 ) {
if ( nNumber % i == 0 ) return 0;
} // for.
return 1;
} // else.
} // isprime.
int main( ) {
const unsigned long cstnMaxTests = DEFAULT_MAX_TESTS;
unsigned long nPrimeNumbers = 0;
steady_clock::time_point start = steady_clock::now( );
for( long n = 1; n < cstnMaxTests; n += 2 ) {
nPrimeNumbers += isprime( n );
} // for.
steady_clock::time_point end = steady_clock::now( );
auto time = duration_cast< milliseconds >( end - start ).count( );
cout
<< "milliseconds: " << time
<< endl << "Number of Primes: " << nPrimeNumbers
<< endl
;
return 0;
} // main.
这里有些代码的运行速度比您快4倍:
module modIsPrime
export isprime, Test_isprime
@inline function isprime(nNumber)
nNumber & 1 == 0 && return 0
for i = 3:2:isqrt(nNumber)
i*unsafe_trunc(Int,nNumber/i) == nNumber && return 0
end
return 1
end
function Test_isprime()
nPrimeNumbers = 1
for n = 3:2:16000000
nPrimeNumbers += isprime(n)
end
println("Prime numbers: ", nPrimeNumbers)
end
end
能够看到同一台机器上的时序比较会很好(我手边没有C ++编译器)。Test_isprime
正确计数有一个小的修正(结果相同,但是在您的代码中计数1
为素数和2
复合数)。
编辑:我已经移动了isqrt
测试后的计算是数字甚至节省了一些时间(虽然不是主要的事情),如果您想与C ++进行公平的比较,它也应该移到那里。
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句