python加速循环嵌套的方法_我如何加速python嵌套循环？

我试图通过计算身体两侧的影响来计算埋藏物体的重力效应,然后总结贡献以在一个站点获得一次测量,重复一些站点.代码如下(正文是一个正方形,代码围绕它顺时针计算,这就是为什么它从-x回到-x坐标)

grav = []

x=si.arange(-30.0,30.0,0.5)

#-9.79742526 9.78716693 22.32153704 27.07382349 2138.27146193

xcorn = (-9.79742526,9.78716693 ,9.78716693 ,-9.79742526,-9.79742526)

zcorn = (22.32153704,22.32153704,27.07382349,27.07382349,22.32153704)

gamma = (6.672*(10**-11))#'N m^2 / Kg^2'

rho = 2138.27146193#'Kg / m^3'

grav = []

iter_time=[]

def procedure():

for i in si.arange(len(x)):# cycles position

t0=time.clock()

sum_lines = 0.0

for n in si.arange(len(xcorn)-1):#cycles corners

x1 = xcorn[n]-x[i]

x2 = xcorn[n+1]-x[i]

z1 = zcorn[n]-0.0 #just depth to corner since all observations are on the surface.

z2 = zcorn[n+1]-0.0

r1 = ((z1**2) + (x1**2))**0.5

r2 = ((z2**2) + (x2**2))**0.5

O1 = si.arctan2(z1,x1)

O2 = si.arctan2(z2,x2)

denom = z2-z1

if denom == 0.0:

denom = 1.0e-6

alpha = (x2-x1)/denom

beta = ((x1*z2)-(x2*z1))/denom

factor = (beta/(1.0+(alpha**2)))

term1 = si.log(r2/r1)#log base 10

term2 = alpha*(O2-O1)

sum_lines = sum_lines + (factor*(term1-term2))

sum_lines = sum_lines*2*gamma*rho

grav.append(sum_lines)

t1 = time.clock()

dt = t1-t0

iter_time.append(dt)

任何帮助加速这个循环将不胜感激谢谢.