１、粒子の初期位置を原点として、 ２、適当なシェルの半径（吸収境界条件の位置a）を決めてdrawTimeする ３、もしdrawTimeした時間と現在時刻の和がシミュレーション終了時刻より以前ならば ４、平面状での角度（0-2pi）を乱数で決定 ５、シェルの半径と角度から現在位置からのdisplacementを得て、位置を更新する ６、時間をdrawTimeで得た時刻だけ進める ７、2に戻る (edited) ３、もしdrawTimeした時間と現在時刻の和がシミュレーション終了時刻よりも後ならば ４、シミュレーション終了時刻引く現在時刻で得た時間幅（drawTimeで得た時刻より手前なはず）でdrawRする ５、平面状での角度を乱数で決定 ６、drawRで得た半径と角度から現在位置からのdisplacementを得て、位置を更新する ７、時間をシミュレーション終了時刻まで進める ８、end (edited)
先の処理に加えて、 １、今自分が直方体のどの平面状にいるのかを位置に加えてに保持 ２、シェルが平面のどの辺も含まない場合は普通に位置更新。ただし、displacementをRと角度から計算するときには所属する平面の法線ベクトルをもとに行う。 ３、辺をひとつだけ含む場合はその辺を軸として折り返した点を新たな位置として、所属する平面の番号を隣り合った平面に移動する。 ４、点を含む場合（辺をふたつ含む場合）は、直方体の場合は３と同様に行う→必要であれば2度折り返す。 ５、シェルの半径を選ぶときは点を含まない最大距離を目指す。そのシェルの半径があるthresholdよりも小さい場合はあきらめて点を高々ひとつ含むようにとる。 ６、点を高々ひとつ取るようにとったときの最大半径がthresholdよりもなお小さい場合はエラーで落とす。 基本的にはこのような処理になります。 おそらく必要になるポリゴンがらみの処理は、 ＊２の法線ベクトルを利用したdisplacementの計算 ＊３，４における辺をdisplacementが横切った際の折り返しの処理 ＊最も近い点または辺とそこまでの距離計算（シェルの半径を決定するのに使う）