CLAPACKの動作について

このエントリーをはてなブックマークに追加
93山崎渉:03/05/21 23:23
━―━―━―━―━―━―━―━―━[JR山崎駅(^^)]━―━―━―━―━―━―━―━―━―
94山崎渉:03/05/28 14:32
     ∧_∧
ピュ.ー (  ^^ ) <これからも僕を応援して下さいね(^^)。
  =〔~∪ ̄ ̄〕
  = ◎――◎                      山崎渉
95山崎 渉:03/07/12 12:33

 __∧_∧_
 |(  ^^ )| <寝るぽ(^^)
 |\⌒⌒⌒\
 \ |⌒⌒⌒~|         山崎渉
   ~ ̄ ̄ ̄ ̄
96名無しさん@3周年:03/07/18 04:05
ATLAS使うとどれくらい速くなるんだ?
97age:03/07/22 04:35
(ノ゜ー゜)ノage!
98名無しさん@3周年:03/07/22 15:53
サンプルムービーだ
http://www.k-514.com/sample/sample.html
  _、_
( ,_ノ` )y━・~~~ 見たらどうだ・・・
99織田:03/07/23 10:57
[簡単に稼げる]頑張れば、頑張るほど稼げるビジネスであります。
http://www.c-gmf.com/index1.htm
「クーポン屋」「アクセスが1日2000人ほどある人なら、1ヶ月274万円
ほど稼げるビジネスであります。」1日110円で宣伝できるサイトであります。
契約者は年間広告費40,000円を当社に振込んでいただければ、OKで
あります。そこで、紹介者として貴方様専用のCGIフォームを用意しますので
料金回収無し、ノルマ無し、HPにバナーを貼ればOKであります。
パートナー様の目標は1ヶ月300件ほど契約を成立していただければ
貴方様の口座に274万円を振込しますので、頑張って営業をしてください。
どんな人でも稼げますので、1度でも良いので見てくださいね。
http://www.c-gmf.com/index1.htm
100山崎 渉:03/08/02 02:24
   ∧_∧
  (  ^^ )< ぬるぽ(^^)
101名無しさん@3周年:03/08/02 02:30
>>96
だいたい >>94 の3倍程度と考えればいいかな。
     ∧_∧  ∧_∧
ピュ.ー (  ・3・) (  ^^ ) <これからも僕たちを応援して下さいね(^^)。
  =〔~∪ ̄ ̄ ̄∪ ̄ ̄〕
  = ◎――――――◎                      山崎渉&ぼるじょあ
103山崎 渉:03/08/15 18:51
    (⌒V⌒)
   │ ^ ^ │<これからも僕を応援して下さいね(^^)。
  ⊂|    |つ
   (_)(_)                      山崎パン
104名無しさん@3周年:04/02/11 22:46
age
105名無しさん@3周年:04/07/19 00:43
Ax=bを計算したときの、xとLとUがほしいのですが、
CLAPACKのdgbsv_に入れた係数行列って、なにになって返されるのでしょうか?
コメントを見ると計算後のAにはP*L*Uが入っていると書いてあるのですが、
それだとAを置換した行列になるはずが、実際に計算してみるとそうはなっていないようです。

http://www.a.mei.titech.ac.jp/~kabe/calcsoft/library/lapack/clapack/usage/dgbsv/
↑のページにも紹介されているのですが、零の項の位置が異なっています。
ソースも見てみたんですが解読できませんでした。
どなたか詳しい方教えてもらえませんでしょうか。
ちなみにCLAPACKではなくてCPPLAPACKを使っています。
106名無しさん@3周年:04/07/19 00:44
実行結果はこんな感じになってます
実行前のA
1 2 1 x x
2 9 8 9 x
1 8 4 5 2
x 9 5 1 6
x x 2 6 7

実行後のA
2 9 8 9 0
0.5 9 5 1 6
0.5 0.388889 2 6 7
x -0.277778 -0.805556 5.94444 6.47222
x x -0.972222 0.102804 6.64019

別の方法で実行したL
1 0 0 x x
2 1 0 0 x
1 1.2 1 0 0
x 1.8 1.38095 1 0
x x -0.47619 -0.450331 1

別の方法で実行したU
1 2 1 x x
0 5 6 9 x
0 0 -4.2 -5.8 2
x 0 0 -7.19048 3.2381
x x 0 0 9.4106
107名無しさん@3周年:04/07/19 02:55
だんだん分かってきました。
IPIVで部分枢軸選択をおこなっているんですね。
でもIPIVの数列を置換行列を指定する仕組みがわかりません(汗
5x5の行列のばあい、IPIV[5]の例としては(2,2,3,4,5)とか(3,4,5,5,5)だったりします。
i 列を IPIV[i] 列と置換するというのなら同じ数字は表れないんじゃないかと思うんですが。
108名無しさん@3周年:04/07/19 11:44
A.dgesv(b)で、AはPAをLU分解した後のLとUの成分行列になり、
対角項で切り離してL,Uを求めてL*Uを求めるとPAと一致しました。

ところが、同じことをバンドマトリックス dgbmtrix Bを用いて、
B.dgbsv(b)を行うと、PBとLUが一致しません。
dgbsvも行置換されたBのLとUの成分行列を返すのですが、Bがバンドマトリクスでも
PBはバンドマトリクスにはなりませんよね。そう考えると
PBのL,U成分行列をdgbmatrixに入れて返すのは無理があると思うのですが・・・

そもそもdgbmatrixクラスのヘッダファイルを見てみるとコンストラクタでは
dgbmatrix(m,n,kl,ky)としたときに、(kl+ku+1)*nの領域しか確保していません。
対して、Cのdgbsrv_の係数行列は(2kl+ku+1)*nの大きさが必要です。
上の質問を書きながら気がつきました。

CPPLAPCKじゃなくてCLAPACK使ったほうがいいのかな。
109名無しさん@3周年:04/07/20 01:46
CPPLAPACKってそんなにできの良い物合ったかな、CLAPACKのほうがよさげ
110名無しさん@3周年:04/07/20 19:00
CLAPACKにはLU分解をするだけの関数、LUを用いて解くだけの関数、があるんですね。
エルミート分解もできるようなので、CPPLAPACKはやめてCLAPACKを使うことにします。
どうもお騒がせしました。
111名無しさん@3周年:04/07/20 19:07
エルミートじゃなくてコレスキーでした
112名無しさん@3周年:04/07/21 01:46
>>108
> dgbmatrix(m,n,kl,ky)としたときに、(kl+ku+1)*nの領域しか確保していません。
> 対して、Cのdgbsrv_の係数行列は(2kl+ku+1)*nの大きさが必要です。
dgbmatrix.dgbsvの中身を見ると,関数の中で一度(M,N,KL,KU+KL)の大きさを
持つバンドマトリックスを作って,データをそこにコピーして,解いて,それ
を元のマトリックスに浅くコピーしている模様. 試しに,dgbsvした後の
マトリックスのサイズを見てみると,勝手に変更されてるのが分かります.
大きさも(2kl+ku+1)*nになってます. 結局,CLAPACKでdgesvを使った時に
得られるものと同じものが得られるハズでは?
113名無しさん@3周年:04/07/21 10:20
>>112
本当だ!再確保してる。
見落としてました。
だから>>106の実行後のAは右上のxがなくなっているんですね。
レスありがとうございます。

でもどうしてL*U計算させると合わないんだろう・・・
もう少しバンドマトリクスの格納方法の絵とにらめっこしてみます。
114名無しさん@3周年:04/09/08 12:48
cblasのルーチン使うときの引数の説明などが日本語化されているページって
ありませんか?小国さんの本にはcblasルーチンの引数の説明が全く無いので
全く分かりません。英語ページでもいいので、分かりやすいところはありま
せんか?
115名無しさん@3周年:04/09/30 09:08:23
Macでこのライブラリーを使いたいのですが、どうしてもコンパイルが上手くいきません。。
CodeWarrior用のものもダウンロードしましたが、ファイルが壊れている
ようなんです。。

MacでClapackなどを使われている方、教えてください。
116名無しさん@3周年:04/09/30 12:34:35
>>115
XCodeがインストール環境であれば、vecLib.frameworkを読み込めばOKです。
117名無しさん@3周年:04/09/30 13:05:03
もはや投稿されたことの無いカキコなどというものは存在しない。
ならば、自分の頭で考え出したものよりも、2年3年の2ちゃんねらーの
蓄積が宿ったコピペやその改変を投稿した方がずっと良い。
118115:04/09/30 14:19:05
>>116
ありがとうございます。

別件なのですが、今、XCodeってダウンロード出来なくなっていませんか?
例のアップルのサイトに行き、会員になってダウンロードのページまで
行ったのですが、Xcodeの古いのから新しいバージョンすべて、
ダウンロード出来なかったです。。
(T-T)

折角、トライしようと思ったんですが。。。
119名無しさん@3周年:04/10/12 17:13:44
zgeevを使うとうまく固有ベクトル出してくれないんすけど.....。

なんか初歩的なミスでもしてるんでしょうか?

サンプルプログラムが落ちている場所など御存知ではないでしょうか?
120名無しさん@3周年:05/01/22 20:06:41
       /                    .\
     /                       .\
 \ /                          \ /
   \      ___________     . /
    \   ./                \   /
     | \/                   \/...|
     |  |                       |  |
     |  |                       |  |
  .._ |_|                       |_|_..
     |  |     >>1さあ‥ 逝こうか‥‥     |  |
     |  |                       |  |
     |  |                       |  |
     |  |                       |  |
     |  |/                   ..\|  |
     |  |                       |  |
   .._|/|      _ _ _ _ _ _ _ _     ...|\|_
  / ..|  |   _−_−_−_−_− _   |  |  \
     |  |/_/━ ━ ━ ━ ━ ━\_\|  |
   .._| /_/━ ━ ━ ━ ━ ━ ━ \_\.|_
     /_/━  ━  ━  ━  ━  ━  ━ \_\
   /  /━  ━  ━  ━  ━  ━  ━\  \
121山.崎 渉:05/02/22 20:57:05
...これからも僕を応援して下さいね(^^)。   
  
━―━―━―━―━―━―━―━―━[JR山崎駅(^^)]━―━―━―━―━―━―━―━―━―
         
     ∧_∧
ピュ.ー (  ^^ ) <これからも僕を応援して下さいね(^^)。                         
  =〔~∪ ̄ ̄〕                                            
  = ◎――◎                      山崎渉                       
                                
 __∧_∧_                                                 
 |(  ^^ )| <寝るぽ(^^)      
 |\⌒⌒⌒\                                
 \ |⌒⌒⌒~|         山崎渉             
   ~ ̄ ̄ ̄ ̄                            
                            
   ∧_∧                                       
  (  ^^ )< ぬるぽ(^^)      
                                                       
    (⌒V⌒)                    
   │ ^ ^ │<これからも僕を応援して下さいね(^^)。   
  ⊂|    |つ                                
   (_)(_)                      山崎パン 
                                         
     ∧_∧  ∧_∧
ピュ.ー (  ・3・) (  ^^ ) <これからも僕たちを応援して下さいね(^^)。
  =〔~∪ ̄ ̄ ̄∪ ̄ ̄〕                          
  = ◎――――――◎                      山崎渉&ぼるじょあ
122名無しさん@5周年:2006/03/14(火) 12:50:23
最近LAPACKてどうなんですか?
123porn:2006/04/26(水) 20:01:54
124CBLAS???:2006/05/20(土) 12:44:48
125名無しさん@5周年:2006/05/20(土) 12:50:40
CBLASのAPIについては、ここの blast-forum/cinterface.ps を見てください。
http://www.netlib.org/blas/
126名無しさん@5周年:2006/06/26(月) 23:16:51

【数学】フェルマーの最終定理に反例が・・・ワイルズ会見
http://news18.2ch.net/test/read.cgi/news7/1149592020/

127バネ足ジョップリン:2007/06/15(金) 18:15:30
バネ足ジョップリンです
CLAMPかと思って釣られました
128名無しさん@5周年:2007/06/17(日) 00:08:20
CLAMP
129名無しさん@5周年:2007/07/19(木) 01:59:12
lapack本家が最近微妙にバージョンアップを始めているな
更新履歴を読んでもよく分からぬのでどなたか解説頼むよ
130馬場:2008/04/24(木) 00:08:01
あのーバンドマトリックスを一番精度よく解く方法をおしえていただきたいんですけど,,
131名無しさん@5周年:2008/12/30(火) 14:42:10
テスト
132名無しさん@5周年:2008/12/31(水) 23:08:25
池■田大■作の本名はソン・テチャク。小泉純一郎、小沢一朗は朝鮮人。
911では小型の水爆が使用されている。
http://ri■ch■ardkosh■im■izu.at.webry.info/
創価の保険金殺人事件。
オウム事件は、統一・創価.北朝鮮の共同犯行である。CIAが監修している。
http://www15.ocn.ne.jp/~oy■ako■don/kok_web■site/ir■iguc■hi.htm
与党も野党もメディアも全部朝鮮人だった。
http://jb■bs.li■vedo■or.jp/b■bs/read.cgi/news/20■92/115794■1306/

2ちゃんねるは統一■教会が個人情報を集めるための道具。運営には統一■教会がいる。
駅前で「手相を見せてください」と勧誘してるのが統一教■会。(カルト宗教)  
カルト宗教の下にいる人と上にいる奴を分けて考えないといけない。下にいる利用される人は上がどんなことをしてるか知らない。

ユダヤ権力の子分→2ちゃん運営=統一協会上層部=層化上層部=自民党清和会=野党の朝鮮人ハーフの政治家=
与党の朝鮮人ハーフの政治家=金 正日(キム・ジョンイル )=読売サンケイ=小沢十朗
ユダヤ人=ロックフェラー=ブッシュ=クリントン=ヒラリー=アドルフ・ヒトラー=オサマ・ビンラディン

毎日新聞スレを荒らしてる奴らも統■一教会の可能性が高い。
荒らしは洗脳するために「ネトウヨ」などのレッテル付けレスを何千回もする。
現実には「ネトウヨ」などは存在しない。
http://changi.2ch.net/test/read.cgi/ms/1230363385/
133名無しさん@5周年:2009/01/14(水) 16:30:36
lapackの導入について調べてたら良スレハケーン
134名無しさん@5周年:2009/01/14(水) 17:05:58
CLAPACKをwindowsに導入したいんですが、うまくできません。
http://d.hatena.ne.jp/blono/20080919/1221783259
これみたけど、よくわからん
135名無しさん@5周年:2009/01/14(水) 21:48:38
>>134
コンパイラは何じゃ。

VC++なら
Netlibから CLAPACK-3.1.1-VisualStudio.zip を落として解凍して
INCLUDE と LIB\Win32 または LIB\x64 にパスを通したり通った所にコピーするなりして
マルチスレッド用スタティックライブラリ使用(/MT or /MTd)でコンパイルするだけだ。
136134:2009/01/15(木) 12:08:21
>>135
Microsoft Visual Studio 2005 Version 8.0 の C++のヤツです
下の2行の意味がわからない場合はどこから出直せばいい?
137名無しさん@5周年:2009/01/16(金) 01:40:36
>>136
VSのIDE使用を前提でいくと

>INCLUDE と LIB\Win32 または LIB\x64 にパスを通したり通った所にコピーするなりして
ttp://www.ne.jp/asahi/hishidama/home/tech/vcpp/setting.html
この辺の
「インクルードファイルの検索パス」、「ライブラリの検索パス」、「ライブラリの追加」
辺りを参照のこと。
追加するライブラリは構成が「Debug」の場合は「clapackd.lib blasd.lib tmglibd.lib libf2cd.lib」
「Release」の場合は「clapack.lib blas.lib tmglib.lib libf2c.lib」

>マルチスレッド用スタティックライブラリ使用(/MT or /MTd)でコンパイルするだけだ。
プロジェクトのプロパティの「構成プロパティ」→「C/C++」→「コード生成」の「ランタイム ライブラリ」を
構成が「Debug」の場合は「マルチスレッド デバッグ(/MTd)」
「Release」の場合は「マルチスレッド(/MT)」
に変更する。

>下の2行の意味がわからない場合はどこから出直せばいい?
Visual Studioの使い方がわからなければ「Visual Studio ライブラリ 追加」でググるところから。
「ディレクトリって何?食い物?」ならばWindowsのファイル操作入門から。
138134:2009/01/16(金) 22:23:02
>>137
まじでありがとうございます
あまりの優しさに胸が熱くなった
Visual Studioの使い方がよくわかってないんでその辺から出直してきます
うまくいったら報告します
139134:2009/01/16(金) 22:42:34
普段は「Debug」でバグを確認してて、
計算時間を計るときは「Release」にするんですけど、
そういうときは「マルチスレッド デバッグ(/MTd)」と「マルチスレッド(/MT)」を
毎回切り替えるということですか?
140名無しさん@5周年:2009/01/17(土) 19:15:56
>>139
プロジェクトの「構成プロパティ」は構成(Debug, Release等)毎に設定できるので
最初に各構成毎に使うランタイムライブラリの種類や追加の依存ファイルを設定しておけば
後はリストボックスでDebug/Releaseを切り替えるだけでOK
141名無しさん@5周年:2009/01/17(土) 20:37:21
>>140
そうだったのかw
VCって便利ですね
142名無しさん@5周年
cpplapack-2010.03.27-1-setup.exe あげ