いわて駐在研究日誌2。

NEVER STAND BEHIND ME

iadj計算の見直し

Gridgenで作成したメッシュデータには、隣接セル情報がないため、ある意味力技で計算している。
したがって、もとの計算ループは、(セル数×セルフェース数)×セル数のループになっている。

  • 一番内側ループでの不要なテンポラリ配列へのコピーを削減 1396.056s → 1221.330s (微妙)
  • index比較をサブルーチンを呼ぶのをやめて直接展開して記述 → 604.64s