前回のプログラムを要素数492ほどで試してみたら、ネットブック(Acer Aspire One HAPPY2）では、1ステップで2.55秒もかかりました。お、おう。ネットブックで解析を回すつもりは無いとはいえ、最低でも1万ステップくらい回すことを踏まえるとさすがに遅い。ノートPC（Dell Inspiron14 core i3-4010U）でも0.63秒で4倍くらい速いものの、全然遅い。勉強用とはいえ、5000要素で1万ステップくらいは１分以内でやってほしい。1ステップで0.005秒くらいとなると、100倍以上の高速化が必要だ。
Pythonが他のプログラムに比べても演算の処理速度が遅いことは知っていたので、最終段階で高速化しようと思っていたけども、現時点で、高速化しないと色々遊んでみることもままならなさそう・・・。っということで、Cython化してみました。
意味も分からず線とか追加。厳密な当たり判定すらできず、ズルしてみたけど、うまく動いとる。ゲームとかならまだしも、解析でやっちゃいかんでしょー。
FEMやったりDEMやったり節操がないですが、手っ取り早くできそうなことには手をつけていくスタンスでさ。
個別要素法に関して、コードを公開してくれている講義資料を見つけましたので、とりあえずPython化。ありがたい。
大学の講義資料なので、一部、意図的に書かれていないこととかあって、あってることやら。
http://www.fml.t.u-tokyo.ac.jp/lecture/lecture_6/text_DEM2014.pdf
粒子間の動作だけだけど、とりあえず動いた。
おー
FEM、FEMと言うけど、これから使うのは差分法を使った解析ソフト。構造分野に絞らず、数値解法として共通する部分も学ぶべきかな？と思い、下記の本を購入しました。
すると、まー、びっくりするほど、まったく分からない。宇宙語のごとしで、まるで頭に入ってこない。「このように～」って言われても何をどうしてどうなったんだー。
前回のVBAによる有限要素法入門を呼んだときも難しいなーっと思っていたけど、内容はなんとなく分かったので、実は分かりやすい本だったんだなぁと今になって実感。最近はネット通販に頼りがちですが、やはり本屋で中身を見ることは大事だなぁ。
さりとて、買ってしまったものは活用せねば。とりあえず・・・、例のごとくPythonで書き換えてみよっと。