(57)【要約】【課題】ＦＤＴＤ法を用いた電磁界解析において、計算時のメモリアクセスの頻度を低減し、計算時間を短縮すること。【解決手段】離散化された解析空間の複数モデルに関する計算に使用されるパラメータを読み込む読込部と、パラメータを複数モデル間で共通化する共通化部と、複数モデルに対して、読込部によって読み込まれたパラメータを用いてＦＤＴＤ法により計算を行う計算部と、を備える。共通化部は、複数モデルの全てのモデルにおいて共通するパラメータをグループ化する。読込部は、複数モデルの計算に必要なパラメータを、共通化部によってグループ化された単位でまとめて読み込みを行う。【選択図】図１
2015125492.tif 000002
【特許請求の範囲】
【請求項１】
離散化された解析空間の複数モデルに関する計算に使用されるパラメータを読み込む読込部と、
前記パラメータを前記複数モデル間で共通化する共通化部と、
前記複数モデルに対して、前記読込部によって読み込まれた前記パラメータを用いてＦＤＴＤ法により計算を行う計算部と、
を備え、
前記共通化部は、前記複数モデルの全てのモデルにおいて共通するパラメータをグループ化し、
前記読込部は、前記複数モデルの計算に必要なパラメータを、前記共通化部によってグループ化された単位でまとめて読み込みを行う、解析装置。
【請求項２】
離散化された解析空間の複数モデルに関する計算に使用されるパラメータを読み込む読込ステップと、
前記パラメータを前記複数モデル間で共通化する共通化ステップと、
離散化された解析空間の複数モデルに対して、前記読込ステップによって読み込まれた前記パラメータを用いてＦＤＴＤ法により計算を行う計算ステップと、
を有し、
前記共通化ステップでは、前記複数モデルの全てのモデルにおいて共通するパラメータをグループ化し、
前記読込ステップでは、前記複数モデルの計算に必要なパラメータを、前記共通化ステップにおいてグループ化された単位でまとめて読み込みを行う、解析方法。
【請求項３】
請求項１に記載された解析装置としてコンピュータを機能させるためのコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＦＤＴＤ（Finite Difference Time Domain method）法を用いた電磁界解析の解析技術に関する。
【背景技術】
【０００２】
近年、時間領域の有限差分法であるＦＤＴＤ法を用いて電磁界を解析する手法が提案されている。ＦＤＴＤ法は、セルに分割された解析空間において適用される。ＦＤＴＤ法は、計算対象となる解析空間の全体について、計算を実行する周波数に応じたサイズのボクセルに離散化して実行される。ボクセルに離散化された解析空間のモデルをボクセルモデルという。
【０００３】
ＦＤＴＤ法を用いて電波の伝搬特性解析等の電磁界シミュレーションを実施する場合、解析空間内に配する電波の波源（例えばアンテナ）や周囲環境を直方体の微小ブロック（Ｙｅｅセル）で再現したボクセルモデルを構築する必要がある。ＦＤＴＤ法は、ボクセルモデルを構成する全てのボクセルに対して、マックスウェル方程式を空間・時間領域での差分方程式に展開して逐次計算することで、電場・磁場を決定する。ＦＤＴＤ法を用いた電磁界解析は、一般的に計算量が多く時間を要する。このような問題を解決するために、ＦＤＴＤ法を用いた電磁界解析に要する時間を削減する技術が提案されている（例えば、特許文献１参照。）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００８−２５０９５６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、ＦＤＴＤ法では、計算に必要なボクセル毎のパラメータがメモリに保持される。そのため、計算時に大量のメモリアクセスが発生し、計算時間が長くなるという問題があった。
上記事情に鑑み、本発明は、ＦＤＴＤ法を用いた電磁界解析において、計算時のメモリアクセスの頻度を低減し、計算時間を短縮する技術を提供することを目的としている。
【課題を解決するための手段】
【０００６】
本発明の一態様は、離散化された解析空間の複数モデルに関する計算に使用されるパラメータを読み込む読込部と、前記パラメータを前記複数モデル間で共通化する共通化部と、前記複数モデルに対して、前記読込部によって読み込まれた前記パラメータを用いてＦＤＴＤ法により計算を行う計算部と、を備え、前記共通化部は、前記複数モデルの全てのモデルにおいて共通するパラメータをグループ化し、前記読込部は、前記複数モデルの計算に必要なパラメータを、前記共通化部によってグループ化された単位でまとめて読み込みを行う、解析装置である。
【０００７】
本発明の一態様は、離散化された解析空間の複数モデルに関する計算に使用されるパラメータを読み込む読込ステップと、前記パラメータを前記複数モデル間で共通化する共通化ステップと、離散化された解析空間の複数モデルに対して、前記読込ステップによって読み込まれた前記パラメータを用いてＦＤＴＤ法により計算を行う計算ステップと、を有し、前記共通化ステップでは、前記複数モデルの全てのモデルにおいて共通するパラメータをグループ化し、前記読込ステップでは、前記複数モデルの計算に必要なパラメータを、前記共通化ステップにおいてグループ化された単位でまとめて読み込みを行う、解析方法である。
【０００８】
本発明の一態様は、上記の解析装置としてコンピュータを機能させるためのコンピュータプログラムである。
【発明の効果】
【０００９】
本発明により、ＦＤＴＤ法を用いた電磁界解析において、計算時のメモリアクセスの頻度を低減し、計算時間を短縮することが可能となる。
【図面の簡単な説明】
【００１０】
【図１】ＦＤＴＤ法を用いて解析を行う対象のボクセルモデルの具体例を表す図である。
【図２】ＦＤＴＤ法を用いて解析を行う対象のボクセルモデルの具体例を表す図である。
【図３】モデル１のボクセルモデルを簡易的に表した模式図である。
【図４】モデル２のボクセルモデルを簡易的に表した模式図である。
【図５】従来手法を用いた解析装置１００の機能構成を表した機能ブロック図である。
【図６】ＦＤＴＤ法を用いた計算手法に従来手法１を用いた場合の、媒質定数の読み込みのタイミングを表した図である。
【図７】ＦＤＴＤ法を用いた電磁界解析の計算に従来手法１を用いた場合の処理の流れを表すフローチャートである。
【図８】ＦＤＴＤ法を用いた計算手法に従来手法２を用いた場合の、媒質定数の読み込みのタイミングを表した図である。
【図９】ＦＤＴＤ法による電磁界解析の計算に従来手法２を用いた場合の処理の流れを表すフローチャートである。
【図１０】実施形態の解析装置１００ａの機能構成を表す機能ブロック図である。
【図１１】実施形態の解析装置１００ａにおける媒質定数の読み込みを表した図である。
【図１２】実施形態の解析装置１００ａが電磁界解析の計算を行う流れを表すフローチャートである。
【図１３】媒質定数の読み込み処理の流れを表すフローチャートである。
【発明を実施するための形態】
【００１１】
［概要］
ＦＤＴＤ法を用いた電磁界解析では、ボクセルモデルを構成する全てのボクセルについて、マックスウェル方程式が空間・時間領域での差分方程式に展開される。展開された式に基づいて逐次計算が行われることにより、電場・磁場が決定される。ＦＤＴＤ法では、解析対象の時間が複数の単位時間（以下、「時間ステップ」という。）に分割される。ボクセルモデルを構成する全てのボクセルについて電場・磁場の計算が時間ステップ毎に行われる。ある時間ステップｎ（ｎは１以上の整数）における電場・磁場が算出されると、時間ステップｎにおいて決定された電場・磁場に基づいて、次の時間ステップ（ｎ＋１）における電場・磁場が決定される。そのため、ＦＤＴＤ法では、次の時間ステップに処理が移行する際に、解析対象のボクセルモデルを構成する全てのボクセルの媒質定数の読み込みが行われる。媒質定数とは、波動が伝播する場を構成する物質・物体における、波動の伝わり易さを表す定数である。例えば、媒質定数は、誘電率、透磁率である。ＦＤＴＤ法では、ボクセル毎の媒質定数が必要である。そのため、解析の対象となる空間によっては膨大な回数のメモリアクセスが発生する。
実施形態の解析装置は、計算に必要な媒質定数の読み込みおいて発生するメモリアクセスの回数を低減することによって、ＦＤＴＤ法を用いた電磁界解析に要する計算時間を短縮する。実施形態の手法によって得られる効果を明確にするために、まず、従来手法における媒質定数の読み込みの方法について説明する。
【００１２】
図１及び図２は、ＦＤＴＤ法を用いて解析を行う対象のボクセルモデルの具体例を表す図である。図１及び図２は、同じ列車の車両内部の空間を解析の対象としたボクセルモデルである。図１は、車両内部に乗客がいない状態のボクセルモデル（以下、「モデル１」という。）である。図２は、車両内部の全座席に乗客が着席した状態のボクセルモデル（以下、「モデル２」という。）である。
図３は、モデル１のボクセルモデルを簡易的に表した模式図である。図３は、モデル１のボクセルモデルが車両を構成するモデル（車両モデル）のみで構成されていることを表している。
図４は、モデル２のボクセルモデルを簡易的に表した模式図である。図４は、モデル２のボクセルモデルが、車両モデルと、車両内部の座席に着席した人体部分を構成するモデル（人体モデル）とで構成されていることを表している。
【００１３】
ＦＤＴＤ法を用いた複数モデルの計算手法には、２つの従来手法がある。一方は、各モデルの計算を直列で行う手法（以下、「従来手法１」という。）である。他方は、各モデルの計算を時間ステップ毎に並列で行う手法（以下、「従来手法２」という。）である。
図５は、従来手法を用いた解析装置１００の機能構成を表した機能ブロック図である。解析装置１００は、制御部１０、ＦＤＴＤ算出部１１、解析空間情報メモリ１２、ボクセルモデル保存メモリ１３及び媒質定数読込部１４を備える。
【００１４】
制御部１０は、解析装置１００の全体の動作を制御する。
ＦＤＴＤ算出部１１は、ＦＤＴＤ法を用いた電磁界解析の各種計算処理を実行する。ＦＤＴＤ算出部１１は、ＦＤＴＤ算出部１１ａと、電界成分吸収境界条件計算部１１ｂと、磁界成分計算部１１ｃと、磁界成分吸収境界条件計算部１１ｄとを備える。
【００１５】
解析空間情報メモリ１２は、磁気ハードディスク装置や半導体記憶装置などの記憶装置を用いて構成される。解析空間情報メモリ１２は、解析空間の構成を座標系の情報とともに記憶する。解析空間情報メモリ１２には、解析空間上の各点に対応する媒質定数の値も記憶される。
ボクセルモデル保存メモリ１３は、磁気ハードディスク装置や半導体記憶装置などの記憶装置を用いて構成される。ボクセルモデル保存メモリ１３は、ボクセルモデルの情報を記憶する。
媒質定数読込部１４は、計算対象のボクセルに対応する媒質定数を解析空間情報メモリ１２から読み出す。
【００１６】
図６は、ＦＤＴＤ法を用いた計算手法に従来手法１を用いた場合の、媒質定数の読み込みのタイミングを表した図である。従来手法１では、モデル１及びモデル２の計算が直列で行われる。また、従来手法１では、モデル毎に全ての時間ステップについて計算が行われる。
図７は、ＦＤＴＤ法を用いた電磁界解析の計算に従来手法１を用いた場合の処理の流れを表すフローチャートである。まず、制御部１０は、ボクセルモデル保存メモリ１３から、計算を行う対象のボクセルモデルを読み込む（ステップＳ１０１）。
【００１７】
制御部１０は、解析空間情報メモリ１２から、計算対象のボクセルモデルに適用する３次元空間の座標の絶対値の情報を読み込む（ステップＳ１０２）。３次元空間の座標の絶対値の情報とは、例えば、解析空間の構造を規定する３次元ＣＡＤデータ（３Ｄ−ＣＡＤデータ）などである。制御部１０は、入力された３次元空間の座標の絶対値をボクセルモデルに適用し、ボクセルの座標に対応する３次元配列を生成する。制御部１０は、ボクセルモデルに対応する配列を生成すると、各座標値に対応する媒質定数を割り当てる（ステップＳ１０３）。
【００１８】
制御部１０は、計算を行う対象の時刻に時間ステップを更新する（ステップＳ１０４）。制御部１０は、現在の時間ステップに所定の単位時間を加算した時刻を新しい時間ステップとする。媒質定数読込部１４は、計算対象のボクセルモデルの計算対象のボクセルについて、媒質定数の読み込みを行う（ステップＳ１０５）。
【００１９】
ＦＤＴＤ算出部１１は、現在の時間ステップの電磁界成分の計算のために、計算対象のボクセルの直前の時間ステップの電磁界成分を読み込む（ステップＳ１０６）。ＦＤＴＤ算出部１１は、直前の時間ステップの電磁界成分に基づいて、計算対象のボクセルの現在の時間ステップにおける電磁界成分を計算する（ステップＳ１０７）。
【００２０】
制御部１０は、計算対象のボクセルにおける電磁界成分の計算が終了すると、現在の時間ステップにおいて、計算対象のボクセルモデルの全てのボクセルについて電磁界成分の計算が終了したか否かを判定する（ステップＳ１０８）。現在の時間ステップにおいて計算が終了していないボクセルが存在している場合（ステップＳ１０８−ＮＯ）、処理がステップＳ１０５に戻り、計算が終了していないボクセルについて電磁界成分の計算が継続される。
【００２１】
一方、現在の時間ステップにおいて全てのボクセルの計算が終了した場合（ステップＳ１０８−ＹＥＳ）、制御部１０は、電磁界成分の計算が収束したか否かを判定する（ステップＳ１０９）。電磁界成分の計算が収束していない場合（ステップＳ１０９−ＮＯ）、処理がステップＳ１０４に戻り、電磁界成分の計算が収束するまで計算が継続される。
【００２２】
一方、電磁界成分の計算が収束した場合（ステップＳ１０９−ＹＥＳ）、制御部１０は、全てのボクセルモデルについて計算が終了したか否かを判定する（ステップＳ１１０）。計算が終了していないボクセルモデルが存在する場合（ステップＳ１１０−ＮＯ）、処理がステップＳ１０１に戻り、計算が終了していないボクセルモデルについての計算が継続される。一方、全てのボクセルモデルについて計算が終了した場合（ステップＳ１１０−ＹＥＳ）、解析装置１００は、解析処理を終了し結果を出力する。
【００２３】
図８は、ＦＤＴＤ法を用いた計算手法に従来手法２を用いた場合の、媒質定数の読み込みのタイミングを表した図である。従来手法２では、モデル１及びモデル２の計算が並列で行われる。また、従来手法２では、時間ステップ毎に全てのボクセルモデルについて計算が行われる。
図９は、ＦＤＴＤ法による電磁界解析の計算に従来手法２を用いた場合の処理の流れを表すフローチャートである。まず、制御部１０は、ボクセルモデル保存メモリ１３から、計算を行う対象の全てのボクセルモデルを読み込む（ステップＳ２０１）。
【００２４】
制御部１０は、解析空間情報メモリ１２から、各ボクセルモデルに適用する３次元空間の座標の絶対値の情報を読み込む（ステップＳ２０２）。制御部１０は、入力された３次元空間の座標の絶対値を各ボクセルモデルに適用し、ボクセルの座標に対応する３次元配列を生成する。制御部１０は、ボクセルモデルに対応する配列を生成すると、各モデルの各座標値に対応する媒質定数を割り当てる（ステップＳ２０３）。
【００２５】
制御部１０は、計算を行う対象の時刻に時間ステップを更新する（ステップＳ２０４）。制御部１０は、現在の時間ステップに所定の単位時間を加算した時刻を新しい時間ステップとする。媒質定数読込部１４は、計算対象の各ボクセルモデルの計算対象のボクセルについて、媒質定数の読み込みを行う（ステップＳ２０５）。
【００２６】
ＦＤＴＤ算出部１１は、現在の時間ステップの電磁界成分の計算のために、計算対象のボクセルの直前の時間ステップの電磁界成分を読み込む（ステップＳ２０６）。ＦＤＴＤ算出部１１は、直前の時間ステップの電磁界成分に基づいて、計算対象のボクセルの現在の時間ステップにおける電磁界成分を計算する（ステップＳ２０７）。
【００２７】
制御部１０は、計算対象のボクセルにおける電磁界成分の計算が終了すると、現在の時間ステップにおいて、全てのボクセルモデルの全てのボクセルについて電磁界成分の計算が終了したか否かを判定する（ステップＳ２０８）。現在の時間ステップにおいて計算が終了していないボクセルが存在している場合（ステップＳ２０８−ＮＯ）、処理がステップＳ２０５に戻り、計算が終了していないボクセルについて電磁界成分の計算が継続される。
【００２８】
一方、現在の時間ステップにおいて全てのボクセルの計算が終了した場合（ステップＳ２０８−ＹＥＳ）、制御部１０は、電磁界成分の計算が収束したか否かを判定する（ステップＳ２０９）。電磁界成分の計算が収束していない場合（ステップＳ２０９−ＮＯ）、処理がステップＳ２０４に戻り、電磁界成分の計算が継続される。一方、電磁界成分の計算が収束した場合（ステップＳ２０９−ＹＥＳ）、解析装置１００は、解析処理を終了し結果を出力する。
【００２９】
以上のように、従来手法によるＦＤＴＤ法を用いた電磁界解析では、従来手法１及び従来手法２のいずれにおいても、同じ回数の媒質定数の読み込みが発生する。すなわち、一つのモデルの計算に要する時間をＴ、計算対象のボクセルモデルのボクセル数をＮ、計算対象のモデルの数をＭとした場合、一つのモデルの計算で発生する媒質定数の読み込みによるメモリアクセスの回数は、Ｍ×Ｔ×Ｎ回となる。
実施形態の解析装置は、計算に必要な媒質定数の読み込みおいて発生するメモリアクセスの回数を低減することによって、ＦＤＴＤ法を用いた電磁界解析に要する計算時間を短縮する。
以下、実施形態の解析装置について説明する。
【００３０】
［実施形態］
図１０は、実施形態の解析装置１００ａの機能構成を表す機能ブロック図である。解析装置１００ａは、媒質定数共通化部１５を備える点で、図５に示した従来手法を用いた解析装置１００と異なる。
【００３１】
媒質定数共通化部１５は、全てのボクセルモデル間で共通する媒質定数を持つボクセルを共通部分としてグループ化する。また、媒質定数共通化部１５は、共通部分以外のボクセルを非共通部分としてグループ化する。例えば、図３及び図４の例の場合、媒質定数共通化部１５は、人体部分（人体モデル）を非共通部分とし、人体部分（人体モデル）以外の車両部分（車両モデル）を共通部分としてグループ化する。媒質定数共通化部１５は、このようにグループ化された共通部分及び非共通部分の情報をグループ化情報として、解析空間情報メモリ１２に保存する。
【００３２】
図１１は、実施形態の解析装置１００ａにおける媒質定数の読み込みを表した図である。実施形態の解析装置１００ａでは、従来手法２と同様にモデル１及びモデル２の計算を並列で行い、時間ステップ毎に全てのボクセルモデルについて計算を行う。
【００３３】
図１２は、実施形態の解析装置１００ａが電磁界解析の計算を行う流れを表すフローチャートである。まず、制御部１０は、ボクセルモデル保存メモリ１３から、計算を行う対象の全てのボクセルモデルを読み込む（ステップＳ３０１）。
【００３４】
制御部１０は、解析空間情報メモリ１２から、各ボクセルモデルに適用する３次元空間の座標の絶対値の情報を読み込む（ステップＳ３０２）。制御部１０は、入力された３次元空間の座標の絶対値を各ボクセルモデルに適用し、ボクセルの座標に対応する３次元配列を生成する。制御部１０は、ボクセルモデルに対応する配列を生成すると、各モデルの各座標値に対応する媒質定数を割り当てる（ステップＳ３０３）。
【００３５】
媒質定数共通化部１５は、ステップＳ３０３においてボクセルに割り当てられた媒質定数を、各ボクセルモデル間で共通化する（ステップＳ３０４）。具体的には、媒質定数共通化部１５は、各ボクセルモデル間で同じ座標に位置するボクセルに割り当てられた媒質定数を比較する。媒質定数共通化部１５は、比較の結果、全てのボクセルモデル間で同一の媒質定数を有するボクセルを共通部分として定義する。また、媒質定数共通化部１５は、共通部分として定義されなかったボクセルを非共通部分として、ボクセルモデル毎に定義する。媒質定数共通化部１５は、共通部分及び非共通部分のボクセルの定義と、共通部分及び非共通部分に対応する媒質定数とを組み合わせた情報を、媒質定数の共通化情報として解析空間情報メモリ１２に保存する。
【００３６】
制御部１０は、計算を行う対象の時刻に時間ステップを更新する（ステップＳ３０５）。制御部１０は、現在の時間ステップに所定の単位時間を加算した時刻を新しい時間ステップとする。
媒質定数読込部１４は、計算対象の各ボクセルモデルの計算対象のボクセルについて、媒質定数の読み込み処理を行う（ステップＳ３０６）。
【００３７】
図１３は、媒質定数の読み込み処理の流れを表すフローチャートである。
まず、媒質定数読込部１４は、計算対象のボクセルモデルについて、計算対象のボクセルに対応する媒質定数が読み込まれているか否かを判定する（ステップＳ４０１）。既に媒質定数が読み込まれている場合（ステップＳ４０１−ＹＥＳ）、媒質定数読込部１４は、媒質定数の読み込みを行うことなく処理を終了する。一方、媒質定数が読み込まれていない場合（ステップＳ４０１−ＮＯ）、媒質定数読込部１４は、解析空間情報メモリ１２に記憶されている媒質定数の共通化情報を参照する（ステップＳ４０２）。
【００３８】
媒質定数読込部１４は、媒質定数の共通化情報から、対象のボクセルが共通部分であるか否かを判定する（ステップＳ４０３）。対象のボクセルが共通部分である場合（ステップＳ４０３−ＹＥＳ）、媒質定数読込部１４は、対象のボクセルに対応する媒質定数を読み込む（ステップＳ４０４）。そして、媒質定数読込部１４は、読み込んだ媒質定数を、共通する全てのボクセルモデルの対応するボクセルの媒質定数に設定する（ステップＳ２０４）。
一方、対象のボクセルが共通部分でない場合（ステップＳ４０３−ＮＯ）、媒質定数読込部１４は、対象のボクセルが属するボクセルモデルの非共通部分の媒質定数を読み込む（ステップＳ４０５）。
【００３９】
図１２の説明に戻る。ＦＤＴＤ算出部１１は、現在の時間ステップの電磁界成分の計算のために、計算対象のボクセルの直前の時間ステップの電磁界成分を読み込む（ステップＳ３０７）。
ＦＤＴＤ算出部１１は、直前の時間ステップの電磁界成分に基づいて、計算対象のボクセルの現在の時間ステップにおける電磁界成分を計算する（ステップＳ３０８）。
【００４０】
制御部１０は、計算対象のボクセルにおける電磁界成分の計算が終了すると、現在の時間ステップにおいて、全てのボクセルモデルの全てのボクセルについて電磁界成分の計算が終了したか否かを判定する（ステップＳ３０９）。現在の時間ステップにおいて計算が終了していないボクセルが存在している場合（ステップＳ３０９−ＮＯ）、処理がステップＳ３０６に戻り、計算が終了していないボクセルについて電磁界成分の計算が継続される。
【００４１】
一方、現在の時間ステップにおいて全てのボクセルの計算が終了した場合（ステップＳ３０９−ＹＥＳ）、制御部１０は、電磁界成分の計算が収束したか否かを判定する（ステップＳ３１０）。電磁界成分の計算が収束していない場合（ステップＳ３１０−ＮＯ）、処理がステップＳ３０５に戻り、電磁界成分の計算が継続される。
一方、電磁界成分の計算が収束した場合（ステップＳ３１０−ＹＥＳ）、解析装置１００ａは、解析処理を終了し結果を出力する。
【００４２】
このように構成された実施形態の解析装置１００ａでは、電磁界解析の計算に必要なボクセル毎の媒質定数の読み込みをボクセルモデル間で共通化することによって、並列して計算されるボクセルモデルの媒質定数の読み込みの重複が排除される。そのため、媒質定数の読み込みに要するメモリアクセスの回数が削減される。その結果、ＦＤＴＤ法による電磁界解析の計算速度を向上させることが可能となる。
【００４３】
例えば、図３及び図４の例において、ボクセルモデルの数をＭ、ボクセルモデルを構成するボクセルの数をＮ、ボクセルモデルを構成するボクセルのうち非共通部分のボクセルの数をＰ、ＦＤＴＤ法による計算を終了するまでに必要な時間ステップ数をＴとした場合を考える。
従来のＦＤＴＤ法では、計算時の媒質定数の更新によって発生するメモリアクセスの回数は、Ｍ×Ｎ×Ｔ回であったのに対し、本実施形態では、｛（Ｎ−Ｐ）＋Ｐ×Ｍ｝×Ｔ回となる。
Ｍ×Ｎ×Ｔ−｛（Ｎ−Ｐ）＋Ｐ×Ｍ｝×Ｔ＝（Ｍ−１）×Ｔ×Ｐ
すなわち本実施形態のＦＤＴＤ法では、従来のＦＤＴＤ法と比較して（Ｍ−１）×Ｔ×Ｐ回のメモリアクセスを削減することができ、並列して計算を行うボクセルモデルが多いほど高い削減効果を得ることができる。また、本実施形態のＦＤＴＤ法では、解析の初期段階のみなど限定的なタイミングで効果が得られるのではなく、処理の全時間において効果が得られる。
【００４４】
上述した実施形態における解析装置１００ａをコンピュータで実現するようにしてもよい。その場合、この機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することによって実現してもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含んでもよい。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよく、ＦＰＧＡ（Field Programmable Gate Array）等のプログラマブルロジックデバイスを用いて実現されるものであってもよい。
以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
【符号の説明】
【００４５】
１００、１００ａ…解析装置，  １０…制御部，  １１…ＦＤＴＤ算出部，  １１ａ…電界成分計算部，  １１ｂ…電界成分吸収境界条件計算部，  １１ｃ…磁界成分計算部，  １１ｄ…磁界成分吸収境界条件計算部，  １２…解析空間情報メモリ，  １３…ボクセルモデル保存メモリ，  １４…媒質定数読込部，  １５…媒質定数共通化部
【図１】
2015125492.tif 000003
【図２】
2015125492.tif 000004
【図３】
2015125492.tif 000005
【図４】
2015125492.tif 000006
【図５】
2015125492.tif 000007
【図６】
2015125492.tif 000008
【図７】
2015125492.tif 000009
【図８】
2015125492.tif 000010
【図９】
2015125492.tif 000011
【図１０】
2015125492.tif 000012
【図１１】
2015125492.tif 000013
【図１２】
2015125492.tif 000014
【図１３】
2015125492.tif 000015