computational-sediment-hyd.hatenablog.jp
の続きです.
数値計算は,特に変わったことはせずに普通の平面二次元計算です.
ドライベッドの処理
ちょっとサボって連続式だけに以下のような条件を付け加えました.
if == 0
if
else :
else if == 0
if
else :
ぎりぎり質量が保存される感じです.
今後修正します.
計算速度
computational-sediment-hyd.hatenablog.jp
numba.jitでまたまた衝撃的な速度です. 時間がかかり過ぎるので,計算時間を120分の1にしています.
なんと600倍です.意味がわかりません.
コードはこちら
Wuさんの本のDry Bedの処理のところ(6.1.4 Wetting and drying technique)が面白いです.
- 作者: Weiming Wu
- 出版社/メーカー: Routledge
- 発売日: 2007/11/15
- メディア: ペーパーバック
- この商品を含むブログを見る