前提・実現したいこと
問題の解決では無く、pythonに関する質問です。
発生している問題・エラーメッセージ
別な人から頂いたソースコードの一部です。
python
1indgas = np.where(dist2(posgas[:,0]-haloPos[j][0],posgas[:,1]-haloPos[j][1],posgas[:,2]-haloPos[j][2],boxSize) < r200[j]**2) 2allgasindex[j] = indgas 3inddm, = np.where(dist2(posdm[:,0]-haloPos[j][0],posdm[:,1]-haloPos[j][1],posdm[:,2]-haloPos[j][2],boxSize) < r200[j]**2) 4 alldmindex[j] = inddm 5 ret = np.hstack([posdm[inddm], idxdm[inddm].reshape(np.size(inddm),1)]) 6 sys.stdout.flush() 7 return ret 8if __name__ == "__main__": 9 p = Pool(NTask) 10 tmp = p.map( store_ID, over300idx ) 11# print(tmp) 12 for i in range(np.size(over300idx)): 13 if str(species) == 'H': 14 alldmpos_H [over300idx[i]] = tmp[i][:,0:3] 15 alldmidx_H [over300idx[i]] = tmp[i][:, 3]
質問は2点ございます。
まず、haloPos[j][1]とありますが、こちらはndarrayの表記でいうと何という表記の仕方でしょうか。(定義によるとhaloPos[j][1]はj番目のhaloの位置(2行目にかかれてある)私はndarrayでj-1行目2行目の表記はhaloPos[j, 1]だと思っていました…
さらに、 tmp[i][:,0:3]に関する質問です。これはそもそも2次元ndarrayになっているのでしょうか? この表記法を何と呼ぶのか分からず、ネットでもうまく検索できなかった為、こちらにて質問をさせていただきました。
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。