小编usr*_*r_g的帖子

我开始编写一个程序进行非线性光束计算.我之所以选择Python,是因为它的Matlab就像代码一样,我正在进行速度测试(确保python是正确的语言来进行快速数值计算)并尝试熟悉python3.我尝试了一种算法,计算从t = 1到n的1/t ^ 2之和(来自书Julia High Perfromance),以比较python3与julia的速度.现在我有一些问题:

1)在我的计算中,朱莉娅没有预期的那么快？Julia ist JIT编译.它应该非常快.为什么python更快？

2)为什么python中的循环如此缓慢？

3)为什么python sum方法比numpy.sum方法慢

4)为什么python geting的sum函数与numpy.sum函数的解决方案略有不同？

我希望你能帮助我.

代码:

# Benchmark Python vs Julia
# (from Julia High Performance page 7 from Avik Sengupta)

import scipy as sp
import time

# Sum of 1/t^2 from t = 1 to n by using loops:
# --------------------------------------------
def pisum(w,n):
    u_sum = 0
    for vi in range(w):
        u_sum = 0
        for vj in range(1,n+1,1):
            u_sum += 1.0/(vj*vj)
    return u_sum

# Sum of 1/t^2 from t …