逆行列を求めたかっただけなんだ。何度も何度も求めなければならないのでできるだけ速い方がよかろうと、ATLAS と CLAPACK をインストールした。色々苦労してインストールしたのち、このようなライブラリでは共通の悩みとなる、必要な関数探しの旅にでた。
で、LAPACK のマニュアルをみつけて、matrix inversion というのがあったのでクリっと。すると、
These driver routines cover all the functionality of the computational routines for linear systems, except matrix inversion. It is seldom necessary to compute the inverse of a matrix explicitly, and it is certainly not recommended as a means of solving linear systems.
と書いてある。逆行列求めるルーチンないじゃん。
がはーん!となっかりしていたら、A と C を与えて A・B = C を解くルーチンで C を単位行列にすればいいんじゃないの?と言われた。なるほど。これはしかし、余計な計算をしちゃったりなんかしてないのかしらん...
前のコンタクトレンズがいまいちフィットしてなかったので交換してもらったやつの検査。度は合ってるけど、左のコンタクトがなんかずれるなぁ〜と思っていたらやっぱり形が合ってないとのこと。
つーことで再び交換。三度目の正直なるか。