派手なメッシュで周期的な境界を定義する方法を見つけようとしています。
サイズが 1x1x1 のボックスを定義し、その中に半径 0.25 の球を入れたとします。この球は中心ではありませんが、球の一部がボックスの反対側から出てくる必要があるように、境界に十分近いです。
たとえば、次のコードの場合
import numpy as np
x_ = np.linspace(0,1,100)
y_ = np.linspace(0,1,100)
z_ = np.linspace(0,1,100)
X,Y,Z = np.meshgrid(x_,y_,z_,indexing='ij')
I = (X-particle['x'])**2 + (Y-particle['y'])**2 + (Z-particle['z'])**2 < particle['r']**2
ブール値の 3D 配列を取得します。ここで、True 値は球の内側にあるメッシュ ポイントであり、False 値は球の内側にあるメッシュ ポイントです。ただし、これは私が望む周期的な境界を保証するものではありません。
各グリッドポイントをループすることなく、これにエレガントな方法はありますか