非構造化メッシュでのローカル適応メッシュ細分割（ローカルAMR）の戦略は何ですか？

非構造化メッシュのローカルAMRに興味があります。現在、OpenFOAMライブラリを使用しています-完全に構造化されていないローカルAMRをサポートしています。

• セル絞り込み条件は、切り取られるセルのリストを決定します
• 選択したセルがリファインされます：メッシュ全体が再構築されます
• 古いメッシュから新しいメッシュにマップが作成されます
• 接続性が再計算されます（フェースセル、エッジフェースなど）
• フィールドは新しいメッシュにマッピングされます

関係するデータ構造は基本的にC ++ベクトルであるため、メッシュは膨張してコピーされます。

静的データ構造を使用するメッシュ上に構築できる代替アプローチについて学ぶ必要があります。それらの一つは、中に存在する並列オクツリーフォレストローカルAMR、あるp4estとデンドロ。

誰かが非構造化メッシュのローカル適応AMR戦略に関する最近のレビューペーパーを私に指摘できますか？

経験に基づくアドバイスはさらに優れています。どのローカルAMRエンジンが、固定データ構造ベースの非構造化メッシュに最適ですか？

論文の最初のページでツリー間の通信のバランスについて読む前に、概要が必要です。:)

上記のコメントから、セルを追加するときにベクターをコピーしないようにすることを理解しています。最も簡単な方法は、必要なセルの最大数のためにスペースを予約することです。

std::vector<YourCellType> myVectorOfCells;
vectorOfCells.reserve(maxNoCells);

ベクトルはmaxNoCellsセルを作成するためのスペースを割り当てましたが、まだセルは作成されていません。これでO(1)、ベクター自体をコピーすることなく、maxNoCellsをベクターに追加できます。ただし、C ++標準では、push_backオペレーションを償却 する必要がありO(1)ます。maxNoCellsを超える数を追加すると、ベクターはそれ自体をコピーし、以前のk倍の数のセルのスペースを予約します（通常の実装では、1.4から2の間でkを選択します）。これにより、ベクターをセルに追加し続けることができますO(1)。このサイズ変更操作はそうで はありませんO(1)。

O(1)$n/2$$n/2$O(1)時間...最初にメモリを予約する限り、O（1）時間で要素を追加することもできます。バンガース教授が前述したように、ツリーのような階層的なデータ構造も、データを格納するために内部的にベクトルを使用します。

ただし、シミュレーションの最初にメモリを割り当てることをお勧めします。十分なメモリがあるかどうかを確認するには、必要なセルの数を把握する必要があります。200.000プロセッサでのシミュレーションでは、データ構造を再割り当てしたり、メモリを使い果たしたり、ディスクにスワップしたりする必要はありません。それが起こった場合、私の意見は、ユーザーの入力エラーのためにプログラムが大声で失敗するはずだということです。