[ruby-list:39306] Re: 固有値、固有ベクトルの計算

西松です.

まだまだひどく未完成のFortranArray-0.01.24.tar.gzを
Webに載せておきました.
RubyでFortranの配列を扱うことができます.
Moduleはやめたみたいです（まるで他人事.）.
プログラムインターフェースが大幅に変わりました.
ごめんなさい.

>>> [ruby-list:39303] 佐藤さん
> ruby: error while loading shared libraries:
> /usr/local/lib/ruby/site_ruby/1.6/i686-linux/FortranArrayModule/FAM.so:
> undefined symbol: zheev_

% man zheev
などがうまくいかない場合は, おそらくBLAS, LAPACKがインストール
されていません. まずBLAS, LAPACKをソースからもしくはパッケージ
などからman pagesを含めてインストールして下さい. libblas.so,
liblapack.so などが必要です. Linuxでしたらldconfig(8), また
MANPATHを通すことも忘れずに.

$ ldd /usr/local/lib/ruby/site_ruby/1.8/i686-linux/fortranarray/fa.so 
        liblapack.so => /usr/local/lib/liblapack.so (0x40016000)
        libdl.so.2 => /lib/libdl.so.2 (0x40459000)
        libcrypt.so.1 => /lib/libcrypt.so.1 (0x4045c000)
        libm.so.6 => /lib/libm.so.6 (0x40489000)
        libc.so.6 => /lib/libc.so.6 (0x404a7000)
        /lib/ld-linux.so.2 => /lib/ld-linux.so.2 (0x80000000)
$ cat zheevtes.rb
# zheevtes.rb
require 'fortranarray'
original_H = FortranArray.new(2, 1..4, 1..4, 1..5, 'Z', 'HE', 'lower').set_all(0.0)
original_H[1,1] = 1.1
original_H[2,2] = 2.2
original_H[3,3] = 3.3
original_H[4,4] = 4.4

original_H[2,1] = Complex(0.2, 0.1)
original_H[3,1] = Complex(0.3, 0.1)
original_H[4,1] = Complex(0.4, 0.1)

original_H[3,2] = Complex(0.3, 0.2)
original_H[4,2] = Complex(0.4, 0.2)

original_H[4,3] = Complex(0.4, 0.3)

ary_H = original_H.clone
w = ary_H.zheev('V') # w -> eigenvalues, ary_H -> eigenvectors
#print ary_H[2,4], "\n"
#print ary_H[4].to_s
#print w[4], "\n"
y = FortranArray.new(1,4,'Z')
y.mv!(original_H, ary_H[4])
y.mv!(original_H, ary_H[3])

ary_H = original_H.clone
w = ary_H.zheev('V', 3) # w -> eigenvalues, ary_H -> eigenvectors
y = FortranArray.new(1,3,'Z')
y.mv!(original_H, ary_H[3], Complex(1.0,0.0), Complex(0.0,0.0), 'LO', 1, 1, 3, 1234567, 5)
   # y -> H * eigenvectors[3] = eigenvalues[3] * eigenvectors[3]
$ ruby --version zheevtest.rb
ruby 1.8.0 (2003-06-23) [i686-linux]
$ 

-- 
 love && peace && free_software
 西松 毅