【手続補正書】
【提出日】令和1年6月4日(2019.6.4)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００５０
【補正方法】変更
【補正の内容】
【００５０】
本開示は、その特定の実施態様を参照して詳細に説明されたが、当業者であれば、これらの実施態様の精神及び範囲から逸脱することなく様々な変更及び改良を行うことができることは明らかであろう。従って、本開示は、添付の特許請求の範囲及びその等価物の範囲に含まれるのであれば、この開示の改良及び変更を包含するものとする。
本件出願は、以下の構成の発明を提供する。
（構成１）
物理的システムを構築するプロセッサ実行方法であって：
物理的システムに関連した特性データを受け取るステップ；
該特性データに基づいて該物理的システムのモデルを構築するステップ；
該モデルを使用して該物理的システムの非線形過渡シミュレーションを行うステップであって、該シミュレーションが、複数の時間間隔に亘って行われ、該ステップが：
該複数の時間間隔を、それぞれが複数の時間間隔を含む時間間隔の群に分割するステップ；
時間間隔の各群に関連したデータを異なる解ユニットに送信するステップ；及び
該時間間隔の群のそれぞれを、該異なる解ユニットを用いて並列に解くステップ、を含む、該ステップ、
該シミュレーションの結果に基づいて物理的特性を決定するステップ、を含み、
該物理的システムを、該シミュレーションで決定された物理的特性に基づいて構築又は修正する、前記方法。
（構成２）
前記物理的システムが、周期的システムであるか又は非周期的システムであるかを決定するステップをさらに含み、該物理的システムが周期的システムである場合に第1のモデルが構築され、該物理的システムが非周期的システムである場合に第2のモデルが構築される、構成1記載の方法。
（構成３）
前記決定するステップが、コンピュータシステムによって自動的に行われる、構成2記載の方法。
（構成４）
前記物理的システムが周期的システムである場合は、前記モデルが：
（数１）
2017107539.tif 000001
に従って構築され、
式中、Ki及びMiは、ヤコビ行列であり、Δxiは、非線形反復の増分であり、かつbiは、非線形反復の残差である、構成1記載の方法。
（構成５）
Mn行列の階数を閾値と比較するステップをさらに含み；
該Mn行列の階数が、該閾値未満である場合は、直接ソルバーを使用して前記モデルの時間間隔を解き；
該Mn行列の階数が、該閾値以上である場合は、反復ソルバーを使用して該モデルの時間間隔を解く、構成4記載の方法。
（構成６）
前記反復ソルバーが、収束基準が満たされるまで前記モデルの時間間隔を処理するループを含む、構成5記載の方法。
（構成７）
並列で解くための前記モデルの時間間隔の分割の前に、前処理行列を用いてA行列を調整するステップをさらに含む、構成6記載の方法。
（構成８）
元のシステム行列のMnブロックの階数が閾値未満である場合は、該Mnブロックを0に設定し、残りの下ブロック三角行列を、直接ソルバーによって並列にLU分解し、該元の行列の解がウッドベリー式によって得られる、構成5記載の方法。
（構成９）
前記物理的システムが非周期的システムである場合は、前記モデルが：
（数２）
2017107539.tif 000002
に従って構築され、
式中、Ki及びMiは、ヤコビ行列であり、Δxiは、非線形反復の増分であり、かつbiは、非線形反復の残差である、構成1記載の方法。
（構成１０）
ブロック直接ソルバーを使用して、前記非周期的システムモデルの時間間隔を解く、構成9記載の方法。
（構成１１）
前記ブロック直接ソルバーが、複数のA行列のブロックのそれぞれに対してLU分解を行い、続いて複数のブロックの前進代入演算を行う、構成10記載の方法。
（構成１２）
前記時間間隔の群が並列に解かれたら、収束チェックを行い；
収束が達成されていない場合は、各群に関連したヤコビ行列を、並列に更新して再び解き；
収束が達成されている場合は、後処理を行って前記物理的特性を決定する、構成1記載の方法。
（構成１３）
前記群を、ニュートン・ラフソン法を使用して調整する、構成12記載の方法。
（構成１４）
前記モデルに関連したデータの第1の部分を、どの特定の群の時間間隔が第1の解ユニットに割り当てられたかに基づいて該第1の解ユニットに送信し、該モデルに関連したデータの第2の異なる部分を、どの特定の群の時間間隔が第2の解ユニットに割り当てられたかに基づいて該第２の解ユニットに送信する、構成1記載の方法。
（構成１５）
各解ユニットが、別個のデータ処理ユニットである、構成1記載の方法。
（構成１６）
前記シミュレーションが、入力に対する前記物理的システムの物理的反応を解析する、構成1記載の方法。
（構成１７）
前記物理的反応が、媒体中の波動伝播、熱伝導、又は電気機械装置の過渡挙動である、構成16記載の方法。
（構成１８）
物理的システムを構築するためのコンピュータ実行システムであって：
1つ以上のデータプロセッサ；
複数のステップを実行するために1つ以上のデータプロセッサに指令するための命令で符号化されたコンピュータ可読媒体であって、該複数のステップが：
物理的システムに関連した特性データを受け取るステップ；
該特性データに基づいて該物理的システムのモデルを構築するステップ；
該モデルを使用して該物理的システムの非線形過渡シミュレーションを行うステップであって、該シミュレーションが、複数の時間間隔に亘って行われ、該ステップが：
該複数の時間間隔を、それぞれが複数の時間間隔を含む時間間隔の群に分割するステップ；
時間間隔の各群に関連したデータを異なる解ユニットに送信するステップ；及び
該時間間隔の群のそれぞれを、該異なる解ユニットを用いて並列に解くステップ、を含む、該ステップ、
該シミュレーションの結果に基づいて物理的特性を決定するステップを含む、該コンピュータ可読媒体を含み、
該物理的システムを、該シミュレーションで決定された物理的特性に基づいて構築又は修正する、前記システム。
（構成１９）
前記複数のステップが：
前記物理的システムが、周期的システムであるか又は非周期的システムであるかを決定するステップをさらに含み、該物理的システムが周期的システムである場合に第1のモデルが構築され、該物理的システムが非周期的システムである場合に第2のモデルが構築される、構成18記載のシステム。
（構成２０）
前記物理的システムが周期的システムである場合は、前記モデルが：
（数３）
2017107539.tif 000003
に従って構築され、
式中、Ki及びMiは、ヤコビ行列であり、Δxiは、非線形反復の増分であり、かつbiは、非線形反復の残差である、構成18記載のシステム。
（構成２１）
前記複数のステップが：
M行列の階数を閾値と比較するステップをさらに含み；
該M行列の階数が、該閾値未満である場合は、直接ソルバーを使用して前記モデルの時間間隔を解き；
該M行列の階数が、該閾値以上である場合は、反復ソルバーを使用して該モデルの時間間隔を解く、構成20記載のシステム。
（構成２２）
物理的システムを構築するための方法のステップを実行するために1つ以上のデータプロセッサに指令するための命令で符号化されたコンピュータ可読媒体であって、該方法が：
物理的システムに関連した特性データを受け取るステップ；
該特性データに基づいて該物理的システムのモデルを構築するステップ；
該モデルを使用して該物理的システムの非線形過渡シミュレーションを行うステップであって、該シミュレーションが、複数の時間間隔に亘って行われ、該ステップが：
該複数の時間間隔を、それぞれが複数の時間間隔を含む時間間隔の群に分割するステップ；
時間間隔の各群に関連したデータを異なる解ユニットに送信するステップ；及び
該時間間隔の群のそれぞれを、該異なる解ユニットを用いて並列に解くステップ、を含む、該ステップ、
該シミュレーションの結果に基づいて物理的特性を決定するステップ、を含み、
該物理的システムを、該シミュレーションで決定された物理的特性に基づいて構築又は修正する、前記コンピュータ可読媒体。