それと Mac OS X で BLAS の関数、例えば、dgemm を呼び出すのにちょっとてこずったので、やり方を書いておく。 まず、ライブラリ含まれているかどうかを調べる。/usr/lib/libatlas.dylib に含まれているはず。 $ otool -v -t /usr/lib/libatlas.dylib|grep …

BLAS, LAPACK の subroutine がややこしいので wrapper を書いた。 しかし、BLAS や LAPACK は破壊的な subroutine なので、 operator だと内部でコピーするので余計な生成がおきてしまう。 *= += -= とかは、メンバに一時的な配列を加えておいて、copy して…