(57)【要約】【課題】粒子法シミュレーションなどに用いる粒子探索方法において、対象粒子の影響半径の内部と共通部分を持たないバケットの粒子探索をしないようにすることで粒子探索を高速化する。【解決手段】本発明に関わる粒子探索方法は、解析空間をバケットで分割するバケット生成ステップと、前記バケットから対象粒子の近傍粒子探索の対象とする探索対象バケットを選択する探索対象バケット選択ステップと、前記探索対象バケットの前記対象粒子に対する前記近傍粒子探索を行う粒子探索ステップとを有することを特徴とする。【選択図】      図２
2014211832.tif 000002
【特許請求の範囲】
【請求項１】
解析空間をバケットで分割するバケット生成ステップと、
前記バケットから対象粒子の近傍粒子探索の対象とする探索対象バケットを選択する探索対象バケット選択ステップと、
前記探索対象バケットの前記対象粒子に対する前記近傍粒子探索を行う粒子探索ステップと、を有する粒子探索方法。
【請求項２】
前記バケット生成ステップは、前記バケットごとにバケット番号を付与する、請求項１記載の粒子探索方法。
【請求項３】
前記バケット生成ステップは、前記対象粒子の粒子番号と前記対象粒子が含まれる前記バケットの前記バケット番号とを関連づける、請求項１または２記載の粒子探索方法。
【請求項４】
前記探索対象バケット選択ステップは、前記対象粒子が含まれる前記バケットとその周辺バケットとを選択し、前記対象粒子を中心とする影響半径内の領域と前記周辺バケットの領域とが共通部分を持つか否かを判定することで前記探索対象バケットを選択する、請求項１から３の内の１項記載の粒子探索方法。
【請求項５】
前記探索対象バケット選択ステップは、前記対象粒子を中心とする前記影響半径の境界と前記周辺バケットの境界とが交差部を持つか否かを判定することで前記探索対象バケットを選択する、請求項４記載の粒子探索方法。
【請求項６】
前記粒子探索ステップは、前記探索対象バケット以外の前記バケットでの前記近傍粒子探索を行わない、請求項１から５の内の１項記載の粒子探索方法。
【請求項７】
解析空間をバケットで分割するバケット生成処理と、
前記バケットから対象粒子の近傍粒子探索の対象とする探索対象バケットを選択する探索対象バケット選択処理と、
前記探索対象バケットの前記対象粒子に対する前記近傍粒子探索を行う粒子探索処理と、を実行させる粒子探索プログラム。
【請求項８】
前記探索対象バケット選択処理は、前記対象粒子が含まれる前記バケットとその周辺バケットとを選択し、前記対象粒子を中心とする影響半径内の領域と前記周辺バケットの領域とが共通部分を持つか否かを判定することで前記探索対象バケットを選択する、請求項７記載の粒子探索プログラム。
【請求項９】
解析空間をバケットで分割するバケット生成部と、
前記バケットから対象粒子の近傍粒子探索の対象とする探索対象バケットを選択する探索対象バケット選択部と、
前記探索対象バケットの前記対象粒子に対する前記近傍粒子探索を行う粒子探索部と、を備える粒子探索装置。
【請求項１０】
前記探索対象バケット選択部は、前記対象粒子が含まれる前記バケットとその周辺バケットとを選択し、前記対象粒子を中心とする影響半径内の領域と前記周辺バケットの領域とが共通部分を持つか否かを判定することで前記探索対象バケットを選択する、請求項９記載の粒子探索装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、粒子法シミュレーションなどに用いられる粒子探索を行う粒子探索方法とこれを実行するプログラム、およびこれを実施する装置に関する。
【背景技術】
【０００２】
近年、流体や構造物の挙動を解析する手法として粒子法シミュレーションが注目されている。粒子法シミュレーションとは、対象とする物体を多数の粒子の集まりとして扱い、粒子の個々の挙動に注目して計算を行うシミュレーション手法である。粒子法シミュレーションの手法としては、ＳＰＨ（Ｓｍｏｏｔｈｅｄ  Ｐａｒｔｉｃｌｅ  Ｈｙｄｒｏｄｙｎａｍｉｃｓ）やＭＰＳ（Ｍｏｖｉｎｇ  Ｐａｒｔｉｃｌｅ  Ｓｉｍｕｌａｔｉｏｎ）が広く用いられている。
【０００３】
粒子法シミュレーションの長所の一つは、同じ大きさの複数の粒子でモデル化するため、要素でモデル作成しなければならない有限要素法などの格子法に比べて、複雑な形状のモデルでも短時間で容易にモデル作成ができることである。また、モデルが大きく変形する場合、格子法では要素が潰れるため解析するのに何度も要素作成をやり直すことが必要になるのに対し、粒子法ではその必要がないという長所もある。
【０００４】
その一方で、粒子法では、対象粒子の物理量を求めるのに周辺粒子の物理量を重み付け平均するが、周辺粒子を探索する処理に計算時間を要するという問題がある。ここで、周辺粒子とは、対象粒子を中心とした影響半径の内部、すなわち、２次元空間の場合は円内部、３次元空間の場合は球内部、に位置する粒子のことである。周辺粒子の探索処理に際しては、粒子１個に対して全粒子ｎ個から周辺粒子を探索しこれを全粒子にわたって行った場合、ｎ２回の粒子探索を繰り返すことになり効率的ではない。そのため、これまで様々な効率的な探索手法が提案されてきた。例えば、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）を用いた粒子探索に適合するツリー法やＧＰＵ（Ｇｒａｐｈｉｃａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）を用いた粒子探索に適合するバケット法などがある。
【０００５】
ツリー法の場合、１セルに１粒子が紐づくが、バケット法の場合は、１バケット（あるいはセル、グリッド、ボクセル、ブロックとも言う）に複数粒子を紐づかせることが可能である。非特許文献１によると、バケットの一辺を影響半径とした場合、対象粒子が属するバケットと近接したバケットについて、すなわち、２次元の場合は３２個のバケット、３次元の場合は３３個のバケットを粒子探索すればよい。また、特許文献１によると、バケットに対応するメモリを節約するため、粒子の中心が位置するバケットのみメモリを確保し、そのバケットに対し粒子探索すればよい。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特願２００９−６９９３０号公報
【非特許文献】
【０００７】
【非特許文献１】越塚誠一ら、“粒子法シミュレーション”、培風館、ｐｐ．４８−５０．
【発明の概要】
【発明が解決しようとする課題】
【０００８】
バケット法により粒子探索する場合、上記のように粒子探索する空間を狭めていくことで粒子探索を効率的にしてきている。しかしながら、非特許文献１や特許文献１の何れの場合でも、対象粒子の影響半径の内部と共通部分を持たない、本来、探索を必要としないバケットをも粒子探索するため、不必要な処理に計算時間を要している。
【０００９】
本発明は上記の問題に鑑みてなされたものであり、その目的は、対象粒子の影響半径の内部と共通部分を持たないバケットの粒子探索をしないようにすることで粒子探索を高速化する、粒子法シミュレーションなどに用いられる粒子探索方法とプログラムおよび装置を提供することにある。
【課題を解決するための手段】
【００１０】
本発明に関わる粒子法シミュレーションなどに用いられる粒子探索方法は、解析空間をバケットで分割するバケット生成ステップと、前記バケットから対象粒子の近傍粒子探索の対象とする探索対象バケットを選択する探索対象バケット選択ステップと、前記探索対象バケットの前記対象粒子に対する前記近傍粒子探索を行う粒子探索ステップとを有する。
【００１１】
本発明に関わる粒子探索プログラムは、解析空間をバケットで分割するバケット生成処理と、前記バケットから対象粒子の近傍粒子探索の対象とする探索対象バケットを選択する探索対象バケット選択処理と、前記探索対象バケットの前記対象粒子に対する前記近傍粒子探索を行う粒子探索処理と、を実行させる。
【００１２】
本発明に関わる粒子探索装置は、解析空間をバケットで分割するバケット生成部と、前記バケットから対象粒子の近傍粒子探索の対象とする探索対象バケットを選択する探索対象バケット選択部と、前記探索対象バケットの前記対象粒子に対する前記近傍粒子探索を行う粒子探索部とを備える。
【発明の効果】
【００１３】
本発明の粒子探索方法とプログラムおよび装置によれば、対象粒子の影響半径の内部とバケットとが共通部分を持つか否かを判定し、共通部分を持つバケットだけを粒子探索するよう構成されるので、粒子探索を効率化しその計算時間を短縮することが可能となる。これにより、粒子法シミュレーションなどの計算時間の短縮に寄与することができる。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施形態に関する粒子探索装置の構成を示すブロック図である。
【図２】本発明の実施形態に関する粒子探索方法の手順を示すフローチャートである。
【図３】本発明の実施形態に関するバケットとバケット番号の対応を示す模式図である。
【図４】本発明の実施形態に関する粒子探索方法における探索対象のバケットを選択する手順を示すフローチャートである。
【図５】本発明の実施形態に関する影響半径の境界とバケットとの交点の対応を示す模式図である。
【図６】本発明の実施例のバケットに含まれる粒子数を示す表である。
【図７】本発明の実施例の既存の粒子探索方法と比べた演算削減率を示す図である。
【発明を実施するための形態】
【００１５】
以下、図を参照しながら、本発明の実施形態を詳細に説明する。但し、以下に述べる実施形態には、本発明を実施するために技術的に好ましい限定がされているが、発明の範囲を以下に限定するものではない。
【００１６】
図１は、本発明の実施形態の粒子探索方法を用いて粒子探索を行う装置の構成を示すブロック図である。粒子探索装置１は、コンピュータ装置としての構成を備える。すなわち、プロセッサ２、主記憶手段３、固定記憶手段４、入力手段５、および出力手段６を備える。
【００１７】
プロセッサ２は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、ＧＰＵ（Ｇｒａｐｈｉｃａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）などであり、コンピュータプログラムを動作させることができる。主記憶手段３は、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）などであり、プロセッサ２が処理中のデータを一時的に保存することができる。固定記憶手段４は、ハードディスクドライブ、ソリッドステートドライブなどであり、コンピュータプログラムやデータを保存することができる。入力手段５は、キーボード、マウスなどである。出力手段６は、ディスプレイ、プリンタなどである。
【００１８】
プロセッサ２では、粒子探索実行プログラム４２が動作することにより、後述するバケット生成部２１、探索対象バケット選択部２２、粒子探索部２３、が順次機能する。主記憶手段３には、バケット情報３１が一時的に保存される。固定記憶手段４には、粒子探索に使用されるモデルデータ４１、粒子探索プログラム４２が予め保存されている。
【００１９】
バケット情報３１は、各粒子番号とバケット番号を関連づけたデータ、各バケットの境界と影響半径の境界との交点数、などのデータを含む。モデルデータ４１は、粒子の粒子番号、粒子の座標、粒子の速度、粒子の加速度、粒子の材料番号、材料特性値、影響半径、解析空間の範囲、バケットによる解析空間の分割数などのデータを含む。粒子探索実行プログラム４２には、バケット生成部２１、探索対象バケット選択部２２、粒子探索部２３が順次機能するプログラムが組み込まれている。
【００２０】
図２は、本発明の実施形態の粒子探索方法の手順を示すフローチャートである。図２を参照して、粒子探索装置１の動作を説明する。
【００２１】
粒子探索装置１のプロセッサ２で機能するバケット生成部２１は、ステップＳ１０１、Ｓ１０２、Ｓ１０３の処理を有する。探索対象バケット選択部２２は、ステップＳ１０４、Ｓ１０５の処理を有する。粒子探索部２３は、Ｓ１０６、Ｓ１０７、Ｓ１０８の処理を有する。
【００２２】
粒子探索装置１の動作としては、まず、バケット生成部２１のステップＳ１０１の処理へ進む。ステップＳ１０１では、入力手段５からプロセッサ２へ粒子探索実行プログラム４２の起動が命令され、プロセッサ２が粒子探索実行プログラム４２を起動すると、固定記憶手段４に保存されているモデルデータ４１を主記憶手段３へ一時的に保存する。
【００２３】
ステップＳ１０２では、プロセッサ２が、解析空間をバケットに分割し、各バケットにバケット番号を付与する処理を行う。この処理では、まず、モデルデータ４１を受けて、影響半径の大きさで解析空間の各軸方向を分割数ｄで等分割することで、バケットを生成することができる。次に各バケットにバケット番号を付与する。バケット番号は、各軸の正方向に位置する順番から算出し、その順番と関連づけて主記憶手段３へ一時的に保存する。例えば、２次元の場合、Ｘ軸の正方向にｉ番目、Ｙ軸の正方向にｊ番目の位置にあるバケット番号は、ｉ＋（ｊ−１）×ｄとして表す。
【００２４】
ステップＳ１０３では、プロセッサ２が、対象とする粒子の粒子番号と、対象粒子が含まれるバケットのバケット番号とを関連付ける処理を行う。この処理では、対象粒子の座標がどのバケットに含まれるかを特定し、対象粒子の粒子番号と対象粒子が含まれるバケットのバケット番号とを関連付けて、主記憶手段３へ一時的に保存する。
【００２５】
ステップＳ１０４では、プロセッサ２が、対象粒子が含まれるバケットとその周辺のバケットとを選択する処理を行う。例えば、２次元の場合、図３に示すバケットとバケット番号を対応させた模式図のように、対象粒子が含まれるバケット番号をｉ＋（ｊ−１）×ｄとすると、周辺バケットのバケット番号は、
ｉ−１＋（ｊ−１）×ｄ、ｉ＋１＋（ｊ−１）×ｄ、
ｉ−１＋（ｊ−２）×ｄ、ｉ＋１＋（ｊ−２）×ｄ、ｉ＋（ｊ−２）×ｄ、
ｉ−１＋ｊ×ｄ、ｉ＋１＋ｊ×ｄ、ｉ＋ｊ×ｄ、
の８つとなる。これらの周辺バケットのバケット番号と対象粒子が含まれるバケットのバケット番号とを関連付けて主記憶手段３へ一時的に保存する。
【００２６】
ステップＳ１０５では、プロセッサ２が、対象粒子を中心ｘ０とする影響半径Ｒの領域内部Ｐ°（２次元の場合は円内部、３次元の場合は球内部）と対象バケットＣとの共通部分があるか否かを判定する処理を行う。ここで、Ｐ、Ｐ°は、
【数１】
2014211832.tif 000003
（太字はベクトルを表す）
であり、それぞれ、Ｐは対象粒子を中心ｘ０とする影響半径Ｒの領域、Ｐ°はＰの境界を除いた内部の領域を表す。
【００２７】
このステップＳ１０５の詳細を、図４に示すフローチャートによって説明する。図４は、本実施形態の粒子探索方法における探索対象のバケットを選択する手順を示すフローチャートである。
【００２８】
ステップＳ２０１では、対象粒子を中心（ｘ０，ｙ０）とする影響半径の境界∂Ｐとバケットの境界∂Ｃとが共通部分、すなわち、交点を持つか否かを判定する処理を行う。ここで、∂Ｐは∂Ｐ＝Ｐ−Ｐ°であり、影響半径の境界を表す。∂Ｃは対象バケットＣの境界を表す。
【００２９】
例えば、２次元の場合、バケットの１辺の長さは影響半径Ｒと同じとして、対象粒子が含まれるバケットの番号をｉ＋（ｊ−１）×ｄとすると、交点を持つか否かを判定するバケットの境界は、
ｘ＝ｘｉ＝Ｒ×（ｉ−１）、ｘ＝ｘｉ＋１＝Ｒ×ｉ
ｙ＝ｙｊ＝Ｒ×（ｊ−１）、ｙ＝ｙｊ＋１＝Ｒ×ｊ
である。境界ｘ＝ｘｉとの交点は、（ｘｉ，ｙ０±｛Ｒ２−（ｘｉ−ｘ０）２｝１／２）である。これがどのｙの区間［ｙｊ−１，ｙｊ］、［ｙｊ，ｙｊ＋１］、［ｙｊ＋１，ｙｊ＋２］に含まれるかを調べることで、交点を持つバケットが特定される。
【００３０】
例えば、図５に示すように、ｙ＝ｙ０−｛Ｒ２−（ｘｉ−ｘ０）２｝１／２は［ｙｊ−１，ｙｊ］の区間に含まれるので、交点Ｂを持つバケットの番号は、［ｙｊ−１，ｙｊ］の区間において境界ｘ＝ｘｉを隔てて隣接する、ｉ−１＋（ｊ−２）×ｄ、ｉ＋（ｊ−２）×ｄとなる。また、ｙ＝ｙ０＋｛Ｒ２−（ｘｉ−ｘ０）２｝１／２は区間［ｙｊ＋１，ｙｊ＋２］に含まれるので、交点Ａを持つバケットの番号は、区間［ｙｊ＋１，ｙｊ＋２］において境界ｘ＝ｘｉを隔てて隣接する、ｉ−１＋ｊ×ｄ、ｉ＋ｊ×ｄとなる。他の３つの境界、すなわち、ｘ＝ｘｉ＋１、ｙ＝ｙｊ、ｙ＝ｙｊ＋１についてもこれと同様の処理を行う。
【００３１】
ステップＳ２０２では、交点が含まれたバケットのフラグに１を加算する処理を行う。
【００３２】
ステップＳ２０３では、対象バケットの全境界についてＳ２０１、Ｓ２０２の処理が終了したか判定し、終了していればステップＳ２０５に進み、終了していなければステップＳ２０４に進む。
【００３３】
ステップＳ２０４では、次の境界に進む処理を行い、Ｓ２０１に進む。
【００３４】
ステップＳ２０５では、バケットのフラグが２以上であるかを判定する処理を行う。フラグが２以上のバケットは、影響半径の境界との交点を２つ以上持っているため、対象粒子を中心とする影響半径の領域内部と共通部分を持つことになる。よって、このバケットは粒子探索対象となり、ステップＳ１０６に進む。バケットのフラグが１であれば、影響半径の境界とバケットの境界が１点で接しており、バケットは影響半径の領域内部と共通部分を持たないため、粒子探索対象ではなく、ステップＳ１０７に進む。フラグが０であれば、バケットは影響半径の領域内部と共通部分を持たないため、粒子探索対象ではなく、ステップＳ１０７に進む。
【００３５】
ステップＳ１０６では、プロセッサ２が、粒子探索対象のバケットの粒子探索をする処理を行う。対象粒子の中心の座標から粒子間距離を計算し、粒子間距離が影響半径未満であれば周辺粒子とし、対象粒子と周辺粒子の粒子番号を関連付けて主記憶手段３へ一時的に保存する。
【００３６】
ステップＳ１０７では、プロセッサ２が、探索対象のバケットが全て探索されたか否かを判定する処理を行う。探索対象のバケットが全て探索されていれば、対象粒子の周辺粒子が探索されているため、処理を終了する。探索対象のバケットが全て探索されていなければ、ステップＳ１０８に進む。
【００３７】
ステップＳ１０８では、プロセッサ２が、次のバケットに進む処理を行い、Ｓ１０５に進む。
【００３８】
本実施形態の粒子法シミュレーションなどに用いられる粒子探索プログラムは、上記の本実施形態の粒子探索方法を実行するプログラムであり、粒子探索実行プログラム４２として固定記憶装置４に格納されている。
【００３９】
本実施形態の粒子探索方法、粒子探索プログラムおよび粒子探索装置によれば、対象粒子の影響半径の内部とバケットとが共通部分を持つか否かを判断し、共通部分を持つバケットだけを粒子探索するよう構成されるので、粒子探索を効率化しその計算時間を短縮することが可能となる。これにより、粒子法シミュレーションなどの計算時間の短縮に寄与することができる。
【００４０】
本実施形態により粒子探索が効率化されることを、以下の実施例により説明する。本実施例では、２次元空間において、対象粒子が存在するバケットとその周辺バケットの合計９つのバケットに対し、バケットの１辺の長さは影響半径Ｒと同じとして、Ｒが粒子間距離ｒの４倍、８倍、１６倍、３２倍、６４倍の場合を各々挙げる。
【００４１】
図６は、１つのバケットに含まれる粒子数、および９つのバケットに含まれる粒子数を表す。対象粒子の位置は、図５に示すように、９つのバケットのうち中央に位置するバケット（ｉ＋（ｊ−１）×ｄ）の中で最も左下とする。この場合、本発明の実施形態を適用すると、最も右上に位置するバケット（ｉ＋１＋ｊ×ｄ）に対する粒子探索が削減できることとなる。
【００４２】
対象粒子を中心とする影響半径Ｒの内部に位置する周辺粒子を探索し、Ｓ１０５〜Ｓ１０８の処理における加減演算、乗除演算、平方根演算、ＩＦ分岐演算について演算回数を数え、９つのバケットに含まれる粒子数毎に、既存の粒子探索方法と比べた演算削減率の粒子数依存性を図７に示す。ここで、演算削減率とは、（１−本発明の実施例の演算回数÷既存の方法の演算回数）×１００（％）を表す。
【００４３】
Ｒ＝４ｒのとき、加減演算、ＩＦ分岐演算が削減されない。これは、Ｓ１０５のバケットと影響半径の境界が交点を持つか判定する演算が、１バケット４粒子に対する粒子探索の演算より多くなってしまうためと考えられる。
【００４４】
一方、Ｒが８ｒより大きい場合は、加減演算、乗除演算、平方根演算、ＩＦ分岐演算とも削減される。これは、影響半径の境界が交点を持つか否かを判定する演算が、１バケットに対する粒子探索の演算より少なくなるためと考えられる。さらに、Ｒが大きくなるにつれて、削減率は増加していく。これは、１バケットに含まれる粒子数が多くなるためと考えられる。
【００４５】
Ｒ＝６４ｒのとき、加減演算の削減率が１０％程度に達し、乗除演算、平方根演算、ＩＦ分岐演算の削減率が各々１１％程度に達する。従って、１バケットに含まれる粒子数が十分多くなると、影響半径の境界とバケットの境界が交点を持つか否かを判定する演算を実施しても、既存の粒子探索方法に比べて最大１１％程度の演算削減が可能となる。
【００４６】
以上、本発明の粒子法シミュレーションなどに用いられる粒子探索方法とプログラムおよび装置は、上記の実施形態や実施例に限定されることなく、特許請求の範囲に記載した発明の範囲内で、種々の変形が可能であり、それらも本発明の範囲内に含まれるものであることはいうまでもない。
【００４７】
また、上記の実施形態の一部又は全部は、以下の付記のようにも記載され得るが、以下には限られない。
【００４８】
付記
（付記１）
解析空間をバケットで分割するバケット生成ステップと、
前記バケットから対象粒子の近傍粒子探索の対象とする探索対象バケットを選択する探索対象バケット選択ステップと、
前記探索対象バケットの前記対象粒子に対する前記近傍粒子探索を行う粒子探索ステップと、を有する粒子探索方法。
（付記２）
前記バケット生成ステップは、前記バケットごとにバケット番号を付与する、
付記１記載の粒子探索方法。
（付記３）
前記バケット生成ステップは、前記対象粒子の粒子番号と前記対象粒子が含まれる前記バケットの前記バケット番号とを関連づける、付記１または２記載の粒子探索方法。
（付記４）
前記探索対象バケット選択ステップは、前記対象粒子が含まれる前記バケットとその周辺バケットとを選択し、前記対象粒子を中心とする影響半径内の領域と前記周辺バケットの領域とが共通部分を持つか否かを判定することで前記探索対象バケットを選択する、付記１から３の内の１項記載の粒子探索方法。
（付記５）
前記探索対象バケット選択ステップは、前記対象粒子を中心とする前記影響半径の境界と前記周辺バケットの境界とが交差部を持つか否かを判定することで前記探索対象バケットを選択する、付記４記載の粒子探索方法。
（付記６）
前記影響半径の境界は、２次元空間では円、３次元空間では球である、付記５記載の粒子探索方法。
（付記７）
前記粒子探索ステップは、前記探索対象バケット以外の前記バケットでの前記近傍粒子探索を行わない、付記１から６の内の１項記載の粒子探索方法。
（付記８）
前記バケット生成ステップは、前記バケットを前記影響半径の大きさで分割することにより生成する、付記１から７の内の１項記載の粒子探索方法。
（付記９）
解析空間をバケットで分割するバケット生成処理と、
前記バケットから対象粒子の近傍粒子探索の対象とする探索対象バケットを選択する探索対象バケット選択処理と、
前記探索対象バケットの前記対象粒子に対する前記近傍粒子探索を行う粒子探索処理と、を実行させる粒子探索プログラム。
（付記１０）
前記探索対象バケット選択処理は、前記対象粒子が含まれる前記バケットとその周辺バケットとを選択し、前記対象粒子を中心とする影響半径内の領域と前記周辺バケットの領域とが共通部分を持つか否かを判定することで前記探索対象バケットを選択する、付記９記載の粒子探索プログラム。
（付記１１）
解析空間をバケットで分割するバケット生成部と、
前記バケットから対象粒子の近傍粒子探索の対象とする探索対象バケットを選択する探索対象バケット選択部と、
前記探索対象バケットの前記対象粒子に対する前記近傍粒子探索を行う粒子探索部と、を備える粒子探索装置。
（付記１２）
前記探索対象バケット選択部は、前記対象粒子が含まれる前記バケットとその周辺バケットとを選択し、前記対象粒子を中心とする影響半径内の領域と前記周辺バケットの領域とが共通部分を持つか否かを判定することで前記探索対象バケットを選択する、付記１１記載の粒子探索装置。
【符号の説明】
【００４９】
１    粒子探索装置
２    プロセッサ
３    主記憶手段
４    固定記憶手段
５    入力手段
６    出力手段
２１    バケット生成部
２２    探索対象バケット選択部
２３    粒子探索部
３１    バケット情報
４１    モデルデータ
４２    粒子探索実行プログラム
【図１】
2014211832.tif 000004
【図２】
2014211832.tif 000005
【図３】
2014211832.tif 000006
【図４】
2014211832.tif 000007
【図５】
2014211832.tif 000008
【図６】
2014211832.tif 000009
【図７】
2014211832.tif 000010