四分木および八分木グリッドの視覚化

いわゆるクワッドツリーおよびオクツリーグリッドは、適応メッシュリファインメントを必要とするアプリケーションにとって非常に魅力的です。たとえば、GerrisやParameshで使用されています。そのようなグリッドの適切なファイル形式と、サポートする視覚化ソフトウェアを知っている人はいますか？Gerris FAQのこの回答も参照してください。私が現在気付いている唯一の潜在的な候補はVTK HyperOctreeクラスですが、ドキュメントはほとんどないようです。

メッシュの具体的な例として、それぞれがセルを含む個のボックスで構成される、まばらに洗練された八分木グリッドを考えてみましょう。Visit / Paraviewでの私の現在の戦略は次のとおりです。$10^4$$8 \times 8 \times 8=512$

1. グリッドを非構造化VTKファイルとして書き込みます。これは、大規模なデータセットではコストがかかり、クワッド/オクツリーの「構造化」を利用しません。
2. ツリー内の隣接する「ボックス」をより大きなブロックに集め、ブロック構造のグリッドを記述します。これにはかなりの追加コードが含まれ、クワッド/オクツリーに埋め込まれた接続情報が失われます。

更新現在、octreeグリッドを直接書き込んで視覚化するための実行可能な形式がない場合は、octreeをブロック構造のグリッドとして書き込むために使用できるファイル形式について提案をいただければ幸いです。理想的には、このフォーマット：

1. 定義/実装が比較的シンプルで、理想的にはオープンソースで、C / Fortranと互換性があります。
2. 多数のグリッドブロック（1000以上など）を効率的に処理できます。
3. 特に細分割境界でグリッドブロックの接続を簡単に指定できます。

私は現在、この目的でSiloを使用していますが、上記の3つの点でより良いスコア（より小さく、より速く、より単純）を探しています。

私の場合、マルチレベルの改良を行っていたとき、以前は異なるレベルのグリッドファイルを個別に作成していました。通常、レベル0は固定されたままのベースの粗いメッシュです。レベル1以上では、パッチごとに異なるファイルを使用しました。レベル1に5つの細かいパッチがある場合、私の戦略は、プログラムの最後にLevel_1_Patch_1、Level_1_Patch_2 ...のように5つの異なるファイルを書き込むことでした（Cでは動的に実行できます（Fortranでは機能しません））。

        char filename1[64];
        sprintf(filename1, "Patch%d.dat", patch_number);

        FILE *file;
        file= fopen(filename1, "wb");

ここでは、有限差分構造化グリッドシステムを使用したことを述べておきます。

したがって、すべてのレベルとそれらの.PLTまたは.DATファイルを取得したことがありました。TecPLOTにすべてインポートします。たとえば、レベル1にいくつかの細かいパッチを使用して2レベルのコードを記述していた場合、メッシュはTecplotでは次のようになります。

私は構造化されていないグリッドの扱いにはあまり詳しくありませんが、構造化されたグリッドはTecplotでうまく処理されます