カテゴリー: 日記

　単純反復法でもやってみたが，Testsetの固い問題群は全然ダメ。簡易Newton法とNewton法はこれからの予定。明日にはできたらうれしいなっと。

　にしても，ステップ幅制御は相当真面目にやらんと使い物にならんなぁ・・・。

　あれこれ見直してたら下らんミスと思い込みが見つかって陽的Runge-Kutta法のプログラムを小修正。とりあえずLorenz方程式とPleiades問題はキッチリ解けていることを確認。

　が，そうなると他のTestset問題も扱いたくなるんだが，他は全部stiffなので陰的Runge-Kutta法は避けられず。来週ちと頑張ってみるとするか。Reductionに拘ると収拾がつかなくなるので，紀要に書いたシンプルな奴を土台にする予定。ただこれ遅いんだよなぁ～。どこまで速くできるかな。

　ちなみにLAPACKの解説書なのに最後がODEってのは前例があるのでそれを踏襲しているのである・・・ってプログラミングしててはっと気が付いた。まぁあれよりは必然性のある構成ではあるんだが。

実際のモデルから生成された常微分方程式を集めたTest Setに掲載されているODE問題群ををほぼそのままリンクしてODE solverにかけられるようにした。FortranソースをCからリンクすること自体はたやすいが，character*(*)で定義されている文字列をCにちゃんと受けた渡す方法が分からず，固定長＋char(0)追加という小修整を施すことに。なかなかうまくいかんもんじゃのう。

・・・ってワシが書いているのはLAPACKの本のハズなのだが・・・。

とりあえず動いて誤差も出ているようなので良しとする。古典的4段4次，Butcher 7段６次ときて，やっぱり埋め込み型も入れたいってことでDOPRI5も実装。

しかし陽的解法だとBLASだけで実装できちゃうからツマランな～。

リンクリストはチマチマ修正しております。

 

　とりあえずできるところからということで，CBLASだけを用いて倍精度陽的Euler法を実装。このまま古典的Runge-Kutta，埋め込み型公式と行きたいところだが今週中に上がるかしら？

　CBLASのリストも作らんといかんと気が付く。BLASはこの表(PDF)があるんだが，CBLASとなるとまた変数だの定数だのが違うもんで，結局cblas.hを見ながらプログラムを書くことになる。どういう表にすべきかはまた考えなきゃなぁ。でもオリジナルのBLAS表が使いやすいから，ほぼ同じ形のものにしようかなぁ。

　あと良く参照するリンクも書いておかないとな。どこかのエントリでまとめておいて，チマチマ更新するのが良いな。

　ということで「LAPACK入門（仮）参照URL一覧」を書きました。