模拟晶体的螺旋生长中,用ndarray数组迭代,在计算到几百步后,开始报错:
D:/PycharmProjects/yyl1.py:91: RuntimeWarning: overflow encountered in double_scalars
+dt*(math.sin(pi*(phi[i,j]-phi_s(midx-i*h,midy-j*h)))+l*U[i,j]*(1+math.cos(pi*(phi[i,j]-phi_s(midx-i*h,midy-j*h)))))
D:/PycharmProjects/yyl1.py:92: RuntimeWarning: overflow encountered in double_scalars
sum_U[i,j]=10*(dt/(h**2))*(U[i+1,j]-2*U[i,j]+U[i-1,j]+U[i-1,j+1]-2*U[i,j]+U[i+1,j-1]+U[i+1,j+1]-2*U[i,j]+U[i-1,j-1]+U[i-1,j+1]-2*U[i,j]+U[i+1,j-1])+F*dt
Traceback (most recent call last):
File "D:/PycharmProjects/yyl1.py", line 91, in <module>
+dt*(math.sin(pi*(phi[i,j]-phi_s(midx-i*h,midy-j*h)))+l*U[i,j]*(1+math.cos(pi*(phi[i,j]-phi_s(midx-i*h,midy-j*h)))))
ValueError: math domain error
Process finished with exit code 1
求大神解决