同じ構造を持つ2つの異なるVTKファイルに保存された2つのフィールド間の数値の差を計算するにはどうすればよいですか？
構造化されたグリッド形式の2つのVTKファイルがあるとします。構造化されたグリッドは同じであり（同じ順序で同じポイントのリストを持っています）、各VTKファイルに「Phi」と呼ばれるフィールドがあります。同じ構造のグリッドを使用して、3番目のVTKファイルを作成し、最初のVTKファイルのPhiと2番目のVTKファイルのPhiの差であるフィールドをプロットします。 これを手動で行う方法を知っています。2つのVTKファイルの未加工テキストを解析し、データを配列にコピーし、一方の配列を他方から減算し、正しい形式のデータを新しいファイルにダンプできます。この差を計算してVTKにエクスポートするより良い方法はありますか？Python、またはVisItやParaviewのような視覚化ソフトウェアのソリューションは、C ++のようなコンパイルされた言語を使用するよりも望ましいでしょう。 この差を計算する目的は、PDEの解を計算するためのさまざまな数値手法を比較することです。同じソフトウェアを使用してソリューションを生成しているため、生成する各ファイルでフィールドPhiを除くすべてのデータが同じになることを保証できます。