1
下面是(非常)shorterned我的程序的vesrion理解困难numpy的阵列行为
的normvecs矢量没有覆盖掉了,你woud能够看到自己只是一旦你运行该代码段中,参数是否正确 GEOM函数接受N + 1×2阵列和n作为输入,我想我做一些非常愚蠢的(我认为我是),或者我不明白这种行为
import numpy as np
#geometry calculations
def geom(pts,n):
r = np.zeros(n)
normvecs = np.zeros((n,2))
tgtvecs = np.zeros((n,2))
alphap = np.zeros(n)
cpts = np.zeros((n,2))
#collocation points
cpts[:,0] = (pts[0:n,0]+pts[1:n+1,0])/2
cpts[:,1] = (pts[0:n,1]+pts[1:n+1,1])/2
#length of panels
r[:] = np.sqrt((pts[0:n,0]-pts[1:n+1,0])**2 + (pts[0:n,1]-pts[1:n+1,1])**2)
#angle of each panel with the horizantal(chord)/refernce axis
alphap[:] = np.arctan2(pts[1:n+1,1]-pts[0:n,1],pts[1:n+1,0]-pts[0:n,0])
#normal vectors
normvecs[:,0] = -np.sin(alphap[:])
normvecs[:,1] = np.cos(alphap[:])
xx = np.cos(alphap[:])
print(normvecs[:,0],normvecs[:,1],np.cos(alphap[:]),xx)
return
geom(np.random.rand(31,2),30)
normvecs [:,1]仅显示nomvecs [:,0]的负数 我不能覆盖像normvecs [:,1] = xx(没有错误..但只是打印相同的东西)
我可能听起来不清楚,一旦你运行你可能会看到我在尝试告诉。
什么问题?(如果是这样的!)
是没有问题,这是一个看不见的逻辑错误我跳过代码中的一行(考虑它无关紧要)张贴在我有这个(错误,我的坏!)赋值`tgtvecs [:,0] = np。 cos(alphap [:]); normvecs [:,] = np.sin(alphap [:])` – fedvasu 2011-08-17 13:25:05