pythonの話
indexで回すことが前提のプログラムになってます.多分癖ですね.なかなか抜けないです.多分すべてをindexループなしで書こうと思うとイテレータでしょうね.可読性の問題がありますが.参考:Pythonのイテレータとジェネレータ - Qiita
indexのfor文ですが,以下のサイトを見ているとそんなに遅くないのかと思ってます.mapで書けるケースは限られているので,リスト内包で書くくらいが無難です.それにしてもNumbaとCythonはやばすぎです. 参考:どこまで速くできる? 達人に学ぶPython超高速データ分析~PyData.Tokyo Meetup #4イベントレポート (1/3):CodeZine(コードジン)
また以下のサイトより.私もリスト内包の可読性には疑問をもっております.シンプルなものなら良いですが山ほどif文があると... 参考:Pythonのリスト内包表記をdisる - shkh's blog
よってindexループに対する個人的なルールは,
- map,lambdaが可能なら使う
- リスト内包で綺麗に書けるなら使う.
- 諦めてfor文.
ご存知かと思いますが,index付きで回すenumerateというものもありますが,rangeより遅いようです. 参考:python - range(len(list)) or enumerate(list)? - Stack Overflow
Section Quasi 2D Uniform Flow model
- 今日中に終わると思ったのですがもう一息です.かなりいい感じですが,壁面付近がQuasi 2Dの限界かと.
- は水深平均渦動粘性係数です.いつものです.本来は水深方向に放物線分布を持ちますが断面内で一定値にしてます.