[論文レビュー] Sea-Rise Flooding on Massive Dynamic Terrains

Herman Haverkort, Jeffrey Janssen|arXiv (Cornell University)|Nov 8, 2012

Computer Graphics and Visualization Techniques参考文献 8被引用数 7

ひとこと要約

本稿では、Z順序曲線走査とキャッシュに配慮したセパレータ手法を活用することで、O(Scan(N))のI/O操作を達成する、I/Oに配慮した新規なアルゴリズムを提示している。これは、時間順処理手法に比べて実際の性能が最大10倍向上する。空間局所性を活用することで、コストの高いソート処理や座標の保存を回避し、単一フロー方向モデルの正しさや一般性を損なわずに、大規模な地形データセットにおける性能を顕著に向上させている。

ABSTRACT

Predicting floods caused by storm surges is a crucial task. Since the rise of ocean water can create floods that extend far onto land, the flood damage can be severe. By developing efficient flood prediction algorithms that use very detailed terrain models and accurate sea-level forecasts, users can plan mitigations such as flood walls and gates to minimize the damage from storm surge flooding. In this paper we present a data structure for predicting floods from dynamic sea-level forecast data on dynamic massive terrains. The forecast data is dynamic in the sense that new forecasts are released several times per day; the terrain is dynamic in the sense that the terrain model may be updated to plan flood mitigations. Since accurate flood risk computations require using very detailed terrain models, and such terrain models can easily exceed the size of the main memory in a regular computer, our data structure is I/O-efficient, that is, it minimizes the number of I/Os (i.e. block transfers) between main memory and disk. For a terrain represented as a raster of N cells, it can be constructed using O(N/B log_M/B N/B) I/Os, it can compute the flood risk in a given small region using O(log_B N) I/Os, and it can handle updating the terrain elevation in a given small region using O(log²_B N) I/Os, where B is the block size and M is the capacity of main memory.

研究の動機と目的

• 大規模なグリッド地形における水文的解析の性能ボトルネックを解消するため、フローアキュムレーションにおけるI/O操作を削減すること。
• 実際の地形データセット（主記憶容量を超えるもの）に対して実用的であり、かつI/Oに配慮したアルゴリズムを設計すること。
• 空間的データレイアウト（例：Z順序）とキャッシュに配慮した設計が、従来の時間順処理法に比べてI/O複雑度と実行時間において優れているかどうかを検討すること。
• ソートに依存するパイプラインを、正しさとスケーラビリティを維持したまま、より単純でキャッシュ最適化された代替手法に置き換える可能性を評価すること。

提案手法

• 行単位の処理中に空間局所性を向上させるためにZ順序走査を用い、キャッシュミスとI/O操作を削減する。
• キャッシュに配慮したセパレータベースのアルゴリズムを採用し、グリッドをΘ(√M) × Θ(√M)のサブグリッドに分割し、サイズΘ(√M)の境界保持型代替グラフを計算し、効率的に統合する。
• M や B の知識を必要としないキャッシュオブliviousなセパレータアルゴリズムの変種を適用し、一般化されたキャッシュポリシーに依存する。
• トポロジカル順序の要件を排除することで、入力をソートする必要がなくなり、座標の保存・復元の必要がなくなり、I/Oオーバーヘッドが削減される。
• 行単位またはZ順序形式でデータを処理し、標準パラメータ（M = 2^30、B = 2^14 または 2^16）を仮定したI/Oモデル下での性能分析を実施する。
• Arge らが提唱した時間順処理法（O(Sort(N))のI/Oを要し、再順序付けに必要な座標を保存）と、すべてのアルゴリズムを比較する。

実験結果

リサーチクエスチョン

• RQ1 ソートや座標の保存に依存せずに、グリッド地形におけるフローアキュムレーションをO(Scan(N))のI/Oで行うことは可能か？
• RQ2 Z順序データレイアウトは、行優先順序に比べて、フローアキュムレーションにおけるI/O効率をどの程度向上させるか？
• RQ3 セパレータに基づくキャッシュに配慮した、またはキャッシュオブリビアスなアルゴリズムは、現実の地形データにおいて時間順処理法に比べてより優れたI/O性能を達成できるか？
• RQ4 集水域仮定が、I/OナーディブなアルゴリズムとZ順序走査アルゴリズムの実際の性能にどのような影響を与えるか？
• RQ5 グリッド構造を活用することで、フロールーティングや階層的ウォータ shed ラベル付けなどの、水文パイプラインの他の段階に対しても同等の効率性を持つアルゴリズムを設計できるか？

主な発見

• 提案されたキャッシュに配慮したセパレータベースのアルゴリズムは、O(Scan(N))のI/O複雑度を達成し、35億セルのグリッドで2.0～6.6分で実行され、時間順処理法を著しく上回る。
• I/OナーディブなアルゴリズムのZ順序バージョンは、同じデータセットで111分の行優先処理に比べて41分に短縮された。
• キャッシュオブリビアスなセパレータアルゴリズムは、O(Scan(N))のI/O複雑度を達成し、118分で処理を完了し、M や B のチューニングなしで高い頑健性を示した。
• TerraStream が採用した時間順処理法は、同じタスクに455分を要し、実際の性能差が最大10倍に達する可能性があることが示された。
• 座標の保存とソートを回避することでI/Oオーバーヘッドが削減され、新規アルゴリズムは漸近的および実用的に高速であることが実証された。
• 予備的な結果から、Z順序レイアウトを用いることで、単純な走査戦略でも洗練されたキャッシュに配慮した設計と同等の性能向上が得られる可能性がある。