だからここに私がやろうとしていることがあります、
N=1000
x=np.arange(0,1,1./float(len(N)))
XX,YY=np.meshgrid(x,x)
l=len(XX)
grid=np.array([ ([XX[i,i],YY[j,j],0. ]) for i in xrange(l) for j in xrange(l) ])
numpy ルーチンはかなり高速ですが、グリッドを別の形式にする必要があり、これにはかなり時間がかかります (numpy 配列のインデックス付けのためだと思います)。
提案をありがとう:)
乾杯