ruby-list

佐藤です。

Vine2.5において、これらをインストールして、

blas        http://l5web.laser5.co.jp/rpm2html/7.1/blas-3.0-9.i386.html
lapack        http://l5web.laser5.co.jp/rpm2html/7.1/lapack-3.0-9.i386.html
g2c        ftp://ftp.ring.gr.jp/pub/linux/Vine/VinePlus/2.5/RPMS/i386


/sbin/ldconfigをやってからインストールしたらうまく行きました。

西浦さんのzheevtes.rbもエラーなく実行できました。

ただ、申し分けないのですが、使い方について教えてください。
練習として

 1 0 -1
 1 2  1
 2 2  3

この行列の固有値、固有ベクトルを求めようとして

require 'fortranarray'

original_H = FortranArray.new(2, 3, 3, 4, 'D', 'GE','lower').set_all(0.0)

original_H[1,1] = 1.0
original_H[1,2] = 0.0
original_H[1,3] = -1.0

original_H[2,1] = 1.0
original_H[2,2] = 2.0
original_H[2,3] = 1.0

original_H[3,1] = 2.0
original_H[3,2] = 2.0
original_H[3,3] = 3.0

p original_H

ary_H = original_H.clone

w = ary_H.zheev('V') # w -> eigenvalues, ary_H -> eigenvectors

p w

p ary_H

というスクリプトをつくってみました。
結果は

# 2, 1..3, 1..3, 1..4, 1..3, 'D', 'GE', 'LOWER', 'NON-UNIT', 0, 0
1
1
2
0
2
2
-1
1
3

# 1, 1..3, 1..3, 'D', 'GE', 'UPPER', 'NON-UNIT', 0, 0
-1.909515966155951
1.390699272883653
4.518816693272298

# 2, 1..3, 1..3, 1..4, 1..3, 'D', 'GE', 'LOWER', 'NON-UNIT', 0, 0
-0.6647535743902899
0
0.2708020827196853
0
0.6962534864042075
2
-0.4656102811879915
0
0.57864919010172

となりました。

固有値は1,2,3となるはずなのですが、
そうなっていないようです。使い方が間違っているからだと
思いますが、どうスクリプトを書いていいのかわかりません。

もうしわけないですが教えていただけないでしょうか。

Thread

Prev Next

In This Thread

Prev Next

[#39102] pstoreのバックアップファイルについて — 堀川 久 <vzw00011@...>

[#39104] Internet Explorer で印刷 — Shin-ichiro HARA <sinara@...>

[#39106] What processor do you run ruby on? — "K.Sasada" <ko1@...>

[#39111] *a=aの挙動の違い — Tamiyasu Murakami <tamiyasu_m@...>

[#39118] YAML without Syck — IWATSUKI Hiroyuki <don@...>

[#39129] InternetExplorer ってインターフェースとして使える？ — Shin-ichiro HARA <sinara@...>

[#39142] WORD Graphのデータシートにデータを貼り付けようとしています — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39148] [BUG] ruby-mode indentation bug — Kenta MURATA <muraken2@...>

[#39154] Re: WORD Graphのデータシートにデータを貼り付けようとしています — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39158] setup.rb on mswin32 — KUBO Takehiro <kubo@...>

[#39170] [ANN] Kwartz -- a template system, available in Ruby, PHP and Java — "maco@..." <kwa@...>

[#39188] debugger on RubyUnit — URA Takefumi <hage@...>

[#39189] MinGW — EMOTO Masahiko <emo@...>

[#39192] Re: WORD Graph のデータシートにデータを貼り付けようとしています — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39197] [ANN] Nora 0.0.20040216 released — MoonWolf <moonwolf@...>

[#39198] RHG 全章公開 — Minero Aoki <aamine@...>

[#39201] undefined reference to `__imp__ruby_safe_level' — EMOTO Masahiko <emo@...>

[#39202] Re: WORD Graph のデータシートにデータを貼り付けようとしています — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39204] Re: WORD Graph のデータシートにデータを貼り付けようとしています — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39212] 特異メソッドでのクラス変数 — "eclipse" <eclipse@...>

[#39213] IconvでUTF-16を変換 — Takashi & Kayoko Sano <tksano@...3.kcn.ne.jp>

[#39232] クラス変数をモジュールのメソッドから参照したい — ruby@...

[#39233] MinGW / socket.so — EMOTO Masahiko <emo@...>

[#39235] Zlib memory leak — sheepman <sheepman@...>

[#39236] uriってどう使うんでしょうか？ — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39237] erb コマンドが compile error になる — Kouhei Yanagita <sugi@...>

[#39239] Re: uriってどう使うんでしょうか？ — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39243] How about 'irb -I path'? — IKEGAMI Daisuke <ikegami@...>

[#39245] Re: uriってどう使うんでしょうか — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39248] Re: uriってどう使うんでしょうか — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39250] 気象庁の．．． Re: uriってどう使うんでしょうか — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39251] ipv6 patch for native win32 — "U.Nakamura" <usa@...>

[#39252] ruby-mode.el のインデントの問題 — Kouhei Yanagita <sugi@...>

[#39254] [ANN] Nora 0.0.20040222 released — MoonWolf <moonwolf@...>

[#39255] 気象庁の．．． Re: uriってどう使うんでしょうか — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39257] Re: 気象庁の．．． Re: uriってどう使うんでしょうか — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39260] File#path.tainted? — 堀川 久 <vzw00011@...>

[#39265] 1.8.1 の File.basename — FUJIMOTO Hisakuni <hisa@...>

[#39269] Ruby 0.95 — Masayoshi Takahashi <maki@...>

[#39270] OSX + REXML — Mitsuru Koyama <mit@...>

[#39273] Re: 気象庁の．．． Re: uriってどう使うんでしょうか — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39275] DnD on win32 — Shinichiro HIDA <shinichiro@...>

[#39286] Re: [ruby-dev:16312] Re: ライブラリ拡大計画 — "NAKAMURA, Hiroshi" <nakahiro@...>

[#39288] 固有値、固有ベクトルの計算 — Masahiro Sato <msato@...>

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

Thread

In This Thread

[#39102] pstoreのバックアップファイルについて — 堀川久 <vzw00011@...>

[#39142] WORD　Graphのデータシートにデータを貼り付けようとしています — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39154] Re: WORD　Graphのデータシートにデータを貼り付けようとしています — Hideyuki Oya <h_oya@...7.dion.ne.jp>

[#39201] undefined reference to `impruby_safe_level' — EMOTO Masahiko <emo@...>

[#39260] File#path.tainted? — 堀川久 <vzw00011@...>