データ解析ツールocta ..
53:Pen4 2.4GHz
04/06/10 01:25
octave:1> A=rand(2000,2000);b=rand(2000,1);
octave:2> tic;A\b;toc
ans = 2.7709
最近のマシンは速いねぇ…
54:名無しさん＠お腹いっぱい。
04/06/10 01:33
>>53
octave:1> A=rand(2000,2000);b=rand(2000,1);
octave:2> tic;A\b;toc
ans = 193.02
(´･ω･`)ｼｮﾎﾞｰﾝ
Status of processor 1 as of: 06/10/04 01:33:21
Processor has been on-line since 05/25/04 18:25:10.
The sparcv9 processor operates at 296 MHz,
and has a sparcv9 floating point processor.
Status of processor 3 as of: 06/10/04 01:33:21
Processor has been on-line since 05/25/04 18:25:11.
The sparcv9 processor operates at 296 MHz,
and has a sparcv9 floating point processor.
55:名無しさん＠お腹いっぱい。
04/06/10 01:52
>>52
多分ならなかったと思うけど...
threadで組んでたかな?
56:名無しさん＠お腹いっぱい。
04/06/10 10:59
>>54
PentiumIII(katmai) 600MHz
octave-2.1.57:1> A=rand(2000,2000);b=rand(2000,1);
octave-2.1.57:2> tic;A\b;toc
ans = 17.742
sparcのv9って比較的最近のCPUじゃないの?
大昔(Coppermineですらない)のPenIIIよりも遅いのですか...??
57:名無しさん＠お腹いっぱい。
04/06/10 11:00
>>56
書くの忘れたけど，
OSは，FreeBSD 5-current で
CPUは 600MHz x 2のSMPです．
58:名無しさん＠お腹いっぱい。
04/06/10 11:10
ちなみに FreeBSD-4.9R で
ちょっと古いCPUですが，
Athlon-MP(多分Thunderbirdコア) 1.2GHzx2 のSMPで
octave-2.1.50:1> A=rand(2000,2000);b=rand(2000,1);
octave-2.1.50:2> tic;A\b;toc
ans = 5.1248
と5秒ちょっとのタイムをマークしました...
Octaveは多分 SMP は効いてないと思うので，速度は
>>53の Pen4 と比較してもクロック分の速度は出てるかな...
59:名無しさん＠お腹いっぱい。
04/06/10 23:19
>>49
それってかなり古いOctaveじゃない?
60:ThinkPad X31
04/06/11 01:05
octave:1> A=rand(2000,2000);b=rand(2000,1);
octave:2> tic;A\b;toc
ans = 5.6739
最近のB５サブノートは速いねぇ…
61:名無しさん＠お腹いっぱい。
04/06/11 02:05
Octaveって疎行列の計算はできますか？
周期境界条件を課した疎行列の計算パッケージがなかなか見つからなくて困っています。
FFTはNlogNのオーダーなんで遅いしなーって思うとNのオーダーで計算できる
疎行列パッケージが一番よいのですが。
62:SS S8
04/06/11 09:20
>>60
cygwinでコンパイルしたやつだけど、
ans = 34.755
だった orz
コンパイルオプションばりばりつけた方がいいのかなぁ
63:名無しさん＠お腹いっぱい。
04/06/11 14:52
>>61
自分で書くと感謝されるよ。
でもLAPACKとかで疎行列扱えないの？
最近はOctaveも外部ライブラリに頼っているからね。
64:名無しさん＠お腹いっぱい。
04/06/11 14:57
うちのThinkPad X31(PentiumM-1.4GHz)も，
FreeBSD 4.10-stable 上で
octave-2.1.57:1> A=rand(2000,2000);b=rand(2000,1);
octave-2.1.57:2> tic;A\b;toc
ans = 5.5506
という値を出しました．
65:名無しさん＠お腹いっぱい。
04/06/12 02:04
漏れも試そうとportsでoctaveをmakeしてるんだけど
atlasのビルドが一日たっても終わんないよう。
66:名無しさん＠お腹いっぱい。
04/06/12 02:31
atlasのインストール時ってコンパイルが終わらないというより、
数値計算を行っている時間が長いようですね。
ずっと数値計算プログラムが動いていますが、
計算精度のチェックでもしているのでしょうか？
ちなみにうちのThinkPad X31だと、とりあえず１日ってことはないな
３時間くらいで終わっていたように思えるけど...
67:名無しさん＠お腹いっぱい。
04/06/12 02:45
>>65
メモリが不安定だと計算処理が途中で異常終了して
計算を繰り返して丸一日以上かかることがありました．
もしも速いCPUを使っていて丸一日以上かかるようだったら、
メモリをチェックしてみた方が良いかもしれません．
普通に使っていて動いていても memtest とかが
通らない場合がありますよ．
私は atlas のコンパイル作業でメモリ異常に
気が付きました．
ちなみに Athlon等の1GHz以上のCPUなら数時間で計算処理が
終わります．
68:45
04/06/12 09:45
なにせ古いマシン、PentiumIII 550E ×2 です。
40GB HDDを買ってきて Fedora CORE 2を再インストールしました。
ご指摘のありました -g 付きの ./configure; make 版です。
$ uname -a
Linux localhost.localdomain 2.6.6-1.427smp #1 SMP Thu Jun 10 10:01:11 EDT 2004 i
686 i686 i386 GNU/Linux
$ octave
GNU Octave, version 2.1.57 (i686-pc-linux-gnu).
octave:1> A=rand(2000,2000);b=rand(2000,1);
octave:2> tic;A\b;toc
ans = 110.47
xcpustate で各プロセッサーの状態を監視しながら実行させました。
なるほど、片側のプロセッサーのみを占有していました。悲しい。
窓を２枚開いて同時に octave を起動して、上記処理を行なうと、
143.50 と 145.39 でした。二人で同時に動かしてもそれなりということか。
少し救われるかも。
47さんのアドバイスに従って再make 中です。
どんなプラットフォームでもちゃんと動作するマルチスレッドの
プログラムを作るのは難しいですからね。ＳＭＰを活かそうと思ったら
やはりガリガリ自分で作るしかなさそう。嬉しいような悲しいような。。。
69:名無しさん＠お腹いっぱい。
04/06/12 22:56
>>68
110秒はちょっと遅いですよね？
これはコンパイルオプションにかなり左右されるのでしょうか？
SSEとかどうなっているのでしょうか？
>>56 は 同じようスペックで17秒ちょいで終わっているようですが...
70:名無しさん＠お腹いっぱい。
04/06/13 03:03
研究環境を完全OpenSource化使用して、いろいろ画策中です。
特にMATLABからOctaveへ関数をいくつか移植しているのだけど、
けっこう大変だよね。
Octaveは関数が一つしか書けないし、微妙なところで関数の仕様が
違っていたりするし。
でもMATLABへのライセンス料を払わなくて済むなら
かなりの出費が押さえられる。
71:名無しさん＠お腹いっぱい。
04/06/13 03:04
＿＿＿＿
.／, - ､, - ､ ￣ ヽ
./-┤ ｡|｡ |―-､ ヽ
| ヽ｀- ○- ´ ／ ヽ |
| - | ― | |
| ´ | ｀ヽ . | |ヽ
∩ 人｀､ _ | ＿.- ´ | .| ＼
| ⌒ヽ ／ ＼ ￣￣￣ ノノ ＼
| |´ |￣―--― ´ヽ ＿ ／⌒＼
＼_ _/-――.| （ Ｔ ） ｀l Τ（ ）
￣ | ￣ } | ＼＿／
| ､--― ￣| /
72:名無しさん＠お腹いっぱい。
04/06/13 03:06
>>71
もしかしてキリ番狙っていたの?
73:名無しさん＠お腹いっぱい。
04/06/13 03:46
||
Λ||Λ
（ / ⌒ヽ
| | | 逝け＞７２
∪ 亅|
| | |
∪∪
：
：
‐ニ三ニ‐
74:＋＋＋
04/06/13 11:00
||
Λ||Λ
（ / ⌒ヽ
| | | 逝け>>７２
∪ 亅|
| | |
∪∪
：
：
‐ニ三ニ‐
75:名無しさん＠お腹いっぱい。
04/06/13 18:03
>>71, >>73, >>74
せっかくのAAなんだけどうちの環境，FreeBSD + mozilla デフォルト設定では
文字がずれて，きちんと見えません...
等幅フォントにすりゃいいんだけどね．
76:名無しさん＠お腹いっぱい。
04/06/13 18:30
このスレせっかくいい調子で話が進んでいたのに、
変な香具師が紛れ込むようになってきたね。
77:名無しさん＠お腹いっぱい。
04/06/13 18:33
必要もないのにageちゃうからだろ。
78:名無しさん＠お腹いっぱい。
04/06/13 19:15
さりげなくサイエンスな香りがするこのスレ、気に入っていたのに
79:名無しさん＠お腹いっぱい。
04/06/13 23:20
WindowsでOctaveって動きますか？
ちょっと面白そうなので立ち寄ってみました。
80:名無しさん＠お腹いっぱい。
04/06/13 23:24
>>79
動くよ。native 版もあるし、cygwin で動かす手もある
特に数値計算なら、こしょこしょコード書くよりまかせてしまったほうがいいと思う
81:名無しさん＠お腹いっぱい。
04/06/13 23:40
レスありがとうございます。
Windowsのnative版もあるんだ。
探してみます。
でもcygwinはややこしいからパスします。
82:45
04/06/14 22:01
途中報告。
-g を抜いたバイナリでSMPカーネルで同じ手順を行ないました。
バイナリサイズはナント 7MB でした。
実効結果は
ans = 105 くらいでした。大差ないです。
原因を引き続き調査中です。
そうそう、SPARC とかは float より double が早いそうですね。
ハードウェアがそう設計されているそうです。
83:名無しさん＠お腹いっぱい。
04/06/14 23:01
>>85
floatよりdoubleが速い処理系は多いよ．特に昔のRISC系CPUの場合は...
それと表示は，桁を省略しているけど，Octaveは基本的に double で
計算しているはずです．
それで105秒だと本当にそれくらいの速度かと...
atlasをコンパイルするとき -g を抜くだけではなくて，
gcc で -O3 とか -msse2 といった最適化オプションは
指定してありますか？
84:名無しさん＠お腹いっぱい。
04/06/15 00:48
>>83
FreeBSDのportsではSSE2対応のCPUでは，自動的に
SSE2が有効になるようにコンパイルオプションが
指定されているようです．
85:SS S8
04/06/15 01:12
gcc-3.3.1/cygwin/PenM 1G ですが、
-O2のみでは
ans = 34.755
だったのが、
FFLAGS, CFLAGS, CXXFLAGS全部に
-O3 -fomit-frame-pointer -mcpu=pentium3 -march=pentium4 -msse2 -funroll-loops
を指定したら、
ans = 14.412
になりました。ただし、一番効いてるのは -mcpu っぽい。
86:名無しさん＠お腹いっぱい。
04/06/15 02:36
よく覚えていないのですが、atlasってアセンブラで
SSE2のコード書いてませんでしたっけ？
たしかプリプロセッサでSSE2を有効にする変数を定義したら
SSE2が有効になったような...
87:名無しさん＠お腹いっぱい。
04/06/15 08:00
Opteron242だと
ans=11.8
やけにおそいな。
OS SuSE 9.1 X86_64
88:名無しさん＠お腹いっぱい。
04/06/16 04:15
>>87
B5サブノートのPentiumMに負けている...
Opteronって通常の数値計算が速いって聞いてますが、
実際のところどうなんでしょうか？
Octaveの使用を含む数値計算用に一台Dual CPUなのを購入予定なのですが...
89:名無しさん＠お腹いっぱい。
04/06/16 17:14
opteron 242*2(1.6GHz)
姫野ベンチ 2CPU
1313MFlops
Pentium200MHz 40.7倍
Xeon 2GHzしかないが
819MFlops
90:名無しさん＠お腹いっぱい。
04/06/19 03:03
今一つ速くないね．Opteron...
まだ最適化が十分ではないのかな？
91:名無しさん＠お腹いっぱい。
04/06/22 11:44
あげ
92:名無しさん＠お腹いっぱい。
04/06/22 11:48
しかしながらベンチマークスレに成り下がっているのは情けなくもある。
93:名無しさん＠お腹いっぱい。
04/06/22 12:14
演算の精度の制御はどうやればできますか？
たぶんデフォルトではdoubleというか64bitっぽいですが、
もっと精度が欲しいのです。
94:名無しさん＠お腹いっぱい。
04/06/22 22:44
scilab
95:名無しさん＠お腹いっぱい。
04/06/23 07:09
>>92
スーパーコンピューターの最大ユーザーはベンチマーク屋である。
という話を聞いたことあるけど。
もっともコンパイルオプションの違いで数倍のパフォーマンスが
違えばショックですね。
MATLABを買ったと思ってＰＣ買い換えれば幸せになれるかも。。。
96:名無しさん＠お腹いっぱい。
04/06/23 08:20
>>93
演算の精度制御は、できないみたい。
「GNU Octave Manual」によれば(p35)
numeric data is currently stored as double precision numbers.
use the IEEE floating format,
the relative precision is approximately 2.2e-16.
とある。
無論、表示の制御は、できる。C likeのフォーマットが使える。
97:名無しさん＠お腹いっぱい。
04/06/23 09:03
>>94,>>96
どうもありがとう。octaveじゃ無理なんですね。
scilabをインストールしてみたら、例のペンチマークがoctaveより速くて、
のりかえを検討していますが、こっちはどうやれば精度の制御ができるんでしょう？
98:名無しさん＠お腹いっぱい。
04/06/23 09:54
octaveは全く触ったことがない厨房ですが、octaveとRは機能的に
かなりかぶりますか？ それとも相互補完するものですか？
99:名無しさん＠お腹いっぱい。
04/06/24 12:56
若干かぶるけど，特定分野だとRの方が強いと思うよ...
Octaveは線形代数系の計算全般といったところか．
100:名無しさん＠お腹いっぱい。
04/06/24 12:59
>>95
MATLABのことを考えると，毎年PCを買い換えれますね．
MATLABは基本的に毎年ライセンス料を払わないといけないプレッシャーが
あります．
まぁ，払わなくてもいいんだろうけど，修正パッチとか
バージョンアップのことを考えるとライセンス料を払い続けるしか
ないシステムになっています．
101:名無しさん＠お腹いっぱい。
04/06/24 13:03
>>95
ベンチマークだいすきですぅ。
102:名無しさん＠お腹いっぱい。
04/06/25 14:28
MATLAB個人で買いましたよ。
octaveじゃできないこともおおいし。
ちなみに「データ解析ツール」じゃなくて「行列言語」な。
間違えるなよ。
103:名無しさん＠お腹いっぱい。
04/06/25 22:18
行列言語ってのもなんだな...
104:名無しさん＠お腹いっぱい。
04/06/25 22:28
Octave使っているのだけど，マニュアルとかが英語なのがちょっとね．
何かお奨めな参考書とかWebページとかりますか？
105:名無しさん＠お腹いっぱい。
04/06/25 22:28
とかありますか？
106:名無しさん＠お腹いっぱい。
04/06/26 09:24
あまり役に立たない本
「Octaveを用いた数値計算入門」
北本 卓也著
ピアソンエヂュケーション
ISBN4-89471-448-5
107:名無しさん＠お腹いっぱい。
04/07/01 07:59
卒研で構造物の強度計算をすることになりました。
それで手頃なツールをさがしているのですが、
Octaveでできますか？
研究室のOSはNetBSDというものだそうです。
108:名無しさん＠お腹いっぱい。
04/07/01 08:02
強度計算がどんな計算をするものなのかわからないので、わからない。
行列計算とか微分方程式とかに帰着させられるんならできるだろうけど。
OSは無関係。
109:名無しさん＠お腹いっぱい。
04/07/04 23:00
初めて立ち寄ってみました。
lsodeの中身は4段のRunge-Kuttaなのでしょうか？
何を見たら分かるのでしょう？どなたか教えてください。
110:名無しさん＠お腹いっぱい。
04/07/04 23:34
微分方程式は Runge-Kuttaを使っていたと思います．
今のOctaveは外部の数値計算ライブラリに依存しているから
そのソースコードを確認するのがいいのではないでしょうか？
でもFORTRANとかで書かれているのでちょっと困りますが...
111:名無しさん＠お腹いっぱい。
04/07/04 23:43
>>109
lsodeはRunge-Kuttaではなく、stiffな問題にはBDF法、nonstiffな問題にはAdams(predictor-corrector)と使い分けています。
いずれもRunge-KuttaのようなSingle step methodでは無く、Multi step methodと呼ばれる方法です。
URLﾘﾝｸ(www.llnl.gov)
もともとODEPACKというFORTRANのパッケージに含まれていて、ソースはwww.netlib.orgからダウンロード出来ていたと思います。
lsodeでググってみたら出てきますよ。
112:名無しさん＠お腹いっぱい。
04/07/05 07:45
二つの異なる情報ですが、ありがとうござました。
linspaceで刻み幅を固定して計算するのだからRKだろうか
と思っていました。ソースを読んで見ます。
BDFかあ。懐かしいなあ。ずっと以前でしたが、理解するのに
一苦労したっけ。小生は、現在は、刻み幅自動調整の
RKFの愛用者です。
113:名無しさん＠お腹いっぱい。
04/07/05 07:48
異なる内容でしたが、ありがとうございます。
linspaceで刻み幅を固定するので、RKかと思っていました。
小生は、現在は、刻み幅自動調整のRKFを愛用しているので
つい、そう思っていました。ソースを読んで見ます。
114:名称未設定
04/07/07 22:57
Matlabのtoolboxが必要なら、octeva-forgeをインストすれば
115:名無しさん＠お腹いっぱい。
04/07/14 03:05
>>114
こんなんあったんだね。
これはちょっと使えるかも。
116:名無しさん＠お腹いっぱい。
04/07/21 20:37
さてMATLAB R14が出たわけだが。
117:名無しさん＠お腹いっぱい。
04/07/21 21:18
MATLABは新しくなればなるほど重く不安定になってくる。
R12以前くらいかな良かったのは。
118:名無しさん＠お腹いっぱい。
04/09/28 22:30:37
久々にベンチマークネタ age
octave-2.1.57:1> A=rand(2000,2000);b=rand(2000,1);
octave-2.1.57:2> tic;A\b;toc
ans = 4.5008
CPU: AMD Athlon(tm) XP 1700+ (1472.93-MHz 686-class CPU)
Origin = "AuthenticAMD" Id = 0x681 Stepping = 1
そんなに速くないな...motherboardも古いし...
119:名無しさん＠お腹いっぱい。
04/09/29 04:41:37
octave-forgeってどーですか
120:名無しさん＠お腹いっぱい。
04/10/02 21:42:10
>>119
けっこう使えるよ。
公開されてるMatlab用のプログラムが変更無しで動く
121:名無しさん＠お腹いっぱい。
04/10/02 21:56:19
octave-forge ってそもそも何ですか?
122:名無しさん＠お腹いっぱい。
04/10/03 00:20:31
>>121
cat /usr/ports/math/octave-forge/pkg-descr
123:名無しさん＠お腹いっぱい。
04/11/13 20:43:16
あげ
124:名無しさん＠お腹いっぱい。
04/11/13 20:45:19
あげ
125:名無しさん＠お腹いっぱい。
04/11/14 03:58:42
日本語化パッチはまだでつか？
126:名無しさん＠お腹いっぱい。
04/11/14 04:15:32
Octave に日本語って必要なんですか?
127:名無しさん＠お腹いっぱい。
04/11/14 04:23:04
必要かどうかは解かりませんけど、
道具が人を選んではいけないと思います。
128:名無しさん＠お腹いっぱい。
04/11/14 04:35:04
>>127
意味分かりません
129:名無しさん＠お腹いっぱい。
04/11/14 05:48:38
英語読むのめどい
130:名無しさん＠お腹いっぱい。
04/11/19 11:26:10
馬鹿とハサミは使い方次第ですね
馬鹿には道具は使えないんでしょうかねぇ
131:名無しさん＠お腹いっぱい。
04/11/19 11:28:12
んなこといっちゃってはずかしいくせに
132:名無しさん＠お腹いっぱい。
04/11/19 16:55:41
Octaveって線形代数電卓って感じだなぁ．
普段何気なく使っていて，それほど必要だとは思わないんだけど
自分の使っているマシンに入ってないと気づいたときのショックは
計り知れない．
無いとかなり不便なんだよね．
133:名無しさん＠お腹いっぱい。
04/11/23 06:37:06
FreeBSDの ports の Octave 新しくなりましたね．
何が変わったのかな?
134:FreeBSD5.2.1-R
04/11/23 20:41:54
2.1.63にアップグレードしたら次のメッセージが出るようになった
warning: lo_ieee_init: unrecognized floating point format!
何でしょうねこれ
135:134
04/11/24 01:33:20
自己レス
最適化がきつすぎたみたい
コンパイルオプション変えたら出なくなった
136:名無しさん＠お腹いっぱい。
04/11/24 05:47:50
スレ違いでしたらすいません。
Octaveで以下の式を解きたいのですが、どの様な手順で進めるたらいいか、どなたか救いの手をよろしくお願いします。
dx(i)/dt = P(i) - k(i)x + a(i)x(i) - b(SUMx(j))
ここでxは2つあり、x＝１の時は以下の式で解を求めます。
When D <> 0
x１(t) = (Ps/2S + Pd/2D) + (x１(0) - Ps/2S)e^(-St) - (Pd/2D)e^(e-Dt)
When D = 0
x１(t) = xd(0) + Pdt/2 + Ps/2S + (xs(0) - Ps/2S)e^(-St),
where
S = k + B;
D = k - B;
xs=x１+x2;
xd=x１- x2;
Ps = P１+ P2;
Pd = P１- P2;
x=2の場合、Pd -> -Pd, x1 -> x2, xd -> -xdとし、同じ様に解を求めます。
137:名無しさん＠お腹いっぱい。
04/12/13 00:21:01
CPUとマザーボードを買い換えました．
AthlonからPen4に変えたのですが
octave って再コンパイルが必要だったんですね．
138:名無しさん＠お腹いっぱい。
04/12/14 20:04:09
-march=athlonとかやってたんでないのか?
139:名無しさん＠お腹いっぱい。
04/12/21 23:47:43
Octave は関係ないんだけど atlas がコンパイル時にSSEや3D Now を
自動的に認識して最適なコードを出力します．
Athlonだと3D Now が入ってしまうから Pen4だと動かなくなってしまいます．
140:名無しさん＠お腹いっぱい。
04/12/22 05:03:05
FreeBSD 5.3-stable + プレスコットなPen4で ATLASを
コンパイルしたら以下のようなエラーが出るようになりました．
何かエラーを解消する方法があれば教えてください．
同じOSで PentiumM なら問題なくコンパイルできたのに...
pre=s, loopO=JIK, ta=N tb=T, mb=0, nb=0, kb=0, lda=0, ldb=0, ldc=0
mu=4, nu=1, ku=8, muladd=0, lat=1 ====> mflop=1573.458410
*** Error code 255
Stop in /usr/ports/math/atlas/work/ATLAS/tune/blas/gemm/THREADED.
*** Error code 1
Stop in /usr/ports/math/atlas/work/ATLAS/bin/THREADED.
Assertion failed: (fp), function GetMMRES, file ../atlas_install.c, line 327.
141:140
04/12/23 09:52:22
PentiumM でコンパイルしたATLASを pkg_create してプレスコットな
Pen4 2.8GHz に入れました。
思ったより速くないですね。
PentiumM でコンパイルしているので SSE2 とかは有効になっている
はずです。
octave-2.1.63:1> A=rand(2000,2000);b=rand(2000,1);
octave-2.1.63:2> tic;A\b;toc
ans = 4.3339
octave-2.1.63:3> tic;A\b;toc
ans = 4.1330
142:名無しさん＠お腹いっぱい。
04/12/23 17:52:46
portsから入れるならWITH_OPTIMIZED_FLAGSを定義しておく
その他CFLAGSやFFLAGS等も適切に設定しておく
143:名無しさん＠お腹いっぱい。
04/12/23 18:21:14
>>141
SSE2の速度はPentiuMとPen4でぜんぜん違いますね。
>>142
WITH_OPTMIZIZED_FLAGSは定義しないでください。gcc-3.4.2の
バグで無限ループに陥ることアリ。CFLAGS/FFLAGSは意味なしです。
144:名無しさん＠お腹いっぱい。
04/12/23 19:31:54
意味無いこと無いだろ
145:名無しさん＠お腹いっぱい。
04/12/23 20:05:09
mahoたん降臨とみた。
146:名無しさん＠お腹いっぱい。
04/12/23 22:58:33
>>143
Pen4 ってSSE2 速いんですか?
それほどでもないと思うのですが...
147:143
04/12/23 23:06:00
>>64 で X31 のベンチがありましたが
こちらも同じマシンですがベンチの結果はほぼ同じです。
同じ環境で PenM 1.4GHz 5.5秒 Pen4 2.8GHz 4.3秒...
これをどうみるか...
148:名無しさん＠お腹いっぱい。
04/12/25 13:29:55
>>145
するどい。143です。
>>147
octaveではないですが、ATLASだとdouble(REAL*8)の
FLOPSがPen4, PenMで全然違います。厳密には覚えてないですが、
Pen4 2.4Gでdgeemが1600MFlopsでPenMが500MFlops程度
だったかと。dualでもtheoretical peak近いところまででます。
(宣伝したい訳ではないが) FreeBSD Expert 2005にATLASの
ベンチマークを行った結果を詳しく書いています。ただPentiumMに
関しては書いていません。REAL*8はあまりに低速だったからです。
single precisionということはないですか?
PenMのほうが同クロックだと相当高速ですよ。
149:名無しさん＠お腹いっぱい。
04/12/26 01:01:24
>>148
プレスコットでも北森でも同じなのでしょうか?
研究室で数値計算用に一台PCを組もうと思っているんですが
世間のPen4の不人気具合いからためらっています．
とりあえず FreeBSD Expert 買ってみます．
150:名無しさん＠お腹いっぱい。
04/12/26 17:44:12
>>149
Expert2005のXeon dualはプレスコットだろうか? すまぬ。厳密な
チェックは怠った...しかし、FreeBSD/amd64を入れてみたがブートはしなかった。
SSE2は理論性能は実クロックx2Flops
でます。簡単にはXeon dual (3.4G)ならば3.4x2(dual)x2(theoretical)
=12.8GFlops Opteronも同様。まぁそこまでは出ません。
数値には出にくいのですが、複雑な処理をさせるとXeonは不利で、
Opteronがキビキビしてるなというのが良く解ります。
FreeBSD ports+HPCというのに僕は期待しています。
ノウハウをみなで共有できたら良いと思いますし、よい記事執筆者
が出ることを期待もしています。
151:143
04/12/27 22:56:58
少しスレ違いな話で申し訳ないのですが
うちのプレスコなPen4とThinkPad X31 で FreeBSD 5.3-stable の
kernelコンパイル時間を計測してみたのですが
make buildkernel KERNCONF=MYPC
とやると ThinkPadX31(PenM 1.4GHz)，自作PC(プレスコPen4 2.8GHz)
X31は 10分35秒, 自作PCは8分42秒でした．
またuser time と CPU使用率は
X31 545.5秒 97.9%
自作PC 795秒 181%
でした．
正直いってプレスコ使うのが嫌になってきています．
152:140
04/12/27 22:57:38
>>151
ごめん名前間違えた 143 ではなく 140です．
153:名無しさん＠お腹いっぱい。
05/01/17 19:45:19 ,
すいません。
お聞きしていいですか？？
現在、octaveで１次元有限要素法のプログラムを書いているのですが、量子井戸の波動関数を求めようとした際に、どうしても解が振動してしまいます。
波動関数の２乗の値はまともなのに、２乗しないとだめです。
どうしてでしょうか？？
154:名無しさん＠お腹いっぱい。
05/01/18 17:46:38
>>153
ネタ？
155:名無しさん＠お腹いっぱい。
05/01/18 19:00:34
>>153
> 波動関数の２乗の値はまともなのに、２乗しないとだめです。
ここの意味がわかりません．
156:名無しさん＠お腹いっぱい。
05/01/18 21:58:47
>>153
とりあえずソースを
固有エネルギー・固有状態がほしいのか
波動関数の時間発展をみたいのか
何をしたいのか明確に説明プリーズ
157:153
05/01/19 04:45:16
すいません。
説明不足でした。
波動関数が振動するとは、１要素ごとにプラスとマイナスが入れ替わってるような、ギザギザした形状になってしまうということです。
以下にソースを載せますので、暇であればお試しください。
欲しいのは時間に依存しない方程式を解いて、固有関数と固有値です。
URLﾘﾝｸ(cyobe36.k-server.org)
158:153
05/01/19 04:46:57
URLﾘﾝｸ(cyobe36.k-server.org)
159:名無しさん＠お腹いっぱい。
05/01/19 21:31:54
>>157
ちゃんと解読してないけど
固有値が一番大きな方の波動関数をプロットしてない？
だったら１サイトごとに振動してても問題ないかと。
固有値が小さいやつをプロットしてみたらどうでしょう。
ただ、plot(x,v(:199))をしてみたけど
ポテンシャル井戸に局在してるようなモードでは無かったです。
対角化したときに固有値・固有状態がどういう順番で出るのか良く知らないんですけど。
160:名無しさん＠お腹いっぱい。
05/01/29 17:13:41
ATLASってiccでビルドするとどれくらい速くなるんでしょうか
161:名無しさん＠お腹いっぱい。
05/02/11 19:46:36
Octave GTK+
URLﾘﾝｸ(octave-gtk.sourceforge.net)
162:名無しさん＠お腹いっぱい。
05/02/11 22:01:01
すいません、宣伝になっちゃいますけど
liboctaveの使い方をまとめたLiboctaveWikiを立ち上げました
URLﾘﾝｸ(liboctave.halfmoon.jp)
ちょっとliboctave使ってみたいなっていう初心者の方から、
liboctaveなら俺に任せろ！な猛者まで、
Wikiの編集に参加していただけるとうれしいです。
なんでOctaveじゃなくてliboctaveの方なんだっていうのは無しの方向で・・(汗
163:名無しさん＠お腹いっぱい。
05/04/29 23:16:47
age
164:名無しさん＠お腹いっぱい。
05/04/30 21:37:06
sage
165:名無しさん＠お腹いっぱい。
05/05/02 03:31:08
% grep -1 BROKEN /usr/ports/math/octave/Makefile
.if ${OSVERSION} < 500000
BROKEN="Does not compile on FreeBSD 4.x"
.else
BROKEN="Does not compile on FreeBSD >= 5.x"
.endif
これはMakefile読めってことですか。
166:名無しさん＠お腹いっぱい。
05/05/02 18:59:44
FreeBSDではもうoctaveは使えません
167:名無しさん＠お腹いっぱい。
05/05/02 19:32:52
>>166
(´･ω･`)ｼｮﾎﾞｰﾝ
168:名無しさん＠お腹いっぱい。
05/05/03 20:28:50
>>166
マジ？
169:名無しさん＠お腹いっぱい。
05/05/08 07:09:16
だめじゃん
170:名無しさん＠お腹いっぱい。
05/05/08 18:22:22
焦って/usr/ports/distfilesを見たけどoctaveのファイルは残ってなかった
誰かFreeBSD5.4R用のi686バイナリパッケージくれよ
171:名無しさん＠お腹いっぱい。
05/05/08 19:39:29
こういう時こそportdowngradeの出番では
172:名無しさん＠お腹いっぱい。
05/05/10 12:26:18
FreeBSDのportsってたまに Octave が Broken になりますよね．
またしばらく待っていると復活しますよ．
173:名無しさん＠お腹いっぱい。
05/05/10 13:21:18
>171
joke かと思ったら本当にあった... sysutils/portdowngrade
174:名無しさん＠お腹いっぱい。
05/05/13 00:07:36
mahoﾀｿがその内直してくれるんじゃね
175:名無しさん＠お腹いっぱい。
05/05/30 22:52:46
なおった？
176:名無しさん＠お腹いっぱい。
05/06/01 10:16:48
少なくとも一昨日の段階ではまだ直ってなかったけど...
mahoタソたのむ
177:名無しさん＠お腹いっぱい。
05/06/01 10:51:19
5.x以降なら直ってるみたい。
4.xはBrokenだけど。
178:名無しさん＠お腹いっぱい。
05/06/01 12:51:40
漏れも使っているので、メンテ頼む mahoさん
179:名無しさん＠お腹いっぱい。
05/06/01 14:04:37
matlab だと、
x{1} = 1; x{2} = 2;
x{1} = [];
で、x の最初の cell を消去できますが、octave ではどうするんでしょうか?
同様に structure array の要素の消去の仕方も、[] を代入では駄目でした。
180:179
05/06/01 14:11:47
すみません。嘘付きました。
cell array の要素の消去は、[] を代入でできるんです。
しかし、structure array の要素の消去が分かりません。
181:179
05/06/01 14:36:59
解決しました。
web では見つけられませんでしたが、よく考えれば分かりました。
for i=1:10; x(i).f=i; end
とし、
x(7) の削除は、
x = x([1:6 8:10]);
ですね。
182:名無しさん＠お腹いっぱい。
05/06/03 07:45:41
今、FreeBSD で portsから octave を入れたのだけど
微妙に良くなっていっているよね octave
メンテや改良ってまだどんどん行われているの?
183:名無しさん＠お腹いっぱい。
05/06/03 09:17:17
FreeBSD 5.3でportsからoctaveを入れようとしたけれど、
===> octave-2.1.63 is marked as broken: Does not compile on FreeBSD >=5.x
って出て、インスコできなかったよ。
これってFreeBSD 5.x以上では使えないってことなの？
184:名無しさん＠お腹いっぱい。
05/06/03 10:01:00
>>183
portsを新しくしてください。
185:183
05/06/03 14:38:31
>>184
portsは新しくしたはずなんだが…ｵｶｼｲﾅ
186:名無しさん＠お腹いっぱい。
05/06/07 19:33:20
URLﾘﾝｸ(adlib.rsch.tuis.ac.jp)
に2.1.Xから多次元配列がサポートされたって書いてあるけど，具体的にどう使えるのでしょうか？
ChangeLogとかNEWSとかにも multi とか Arrayとか検索しても引っかからないし，
zeros(3,3,3)とかも使えないし…．
どなたか教えてください．
187:名無しさん＠お腹いっぱい。
05/06/07 20:35:35
octaveじゃないけど
math/sdpaも壊れてるような
188:名無しさん＠お腹いっぱい。
05/06/09 00:35:11
>>186
バージョンは？
2.1.50ではだめだけど、2.1.60以降ではzeros(3,3,3)でいける。
189:名無しさん＠お腹いっぱい。
05/06/09 22:07:53
>>188
2.1.71にしたら出来ました！ありがとうございました！
190:名無しさん＠お腹いっぱい。
05/06/21 11:04:05
imageコマンドで起動される外部画像ビューアーが、インストールされた計算機によって、
imagemagicだったり、xvだったりするのですが、この違いはなぜ生じるのでしょうか。
また、明示的にどちらかを使いたい場合、どうすればよいでしょうか。
191:名無しさん＠お腹いっぱい。
05/07/06 23:18:05
printf('kita---(.v.)-----!'); sleep(10);
これが１０秒後まで表示されないんですよね
どうしようもないですか？
192:名無しさん＠お腹いっぱい。
05/07/16 22:05:40
>>191
printf('kita---(.v.)-----!'); fflush(stdout);sleep(10);
193:名無しさん＠お腹いっぱい。
05/07/27 01:39:55
ちょっと板違いなんですが、他に内容に沿ったスレ見つけられなかったんで失礼します。
今VisualC++で卒業研究のプログラム作ってるんですが、VC++で使える数値計算ライブラリってないですか？
欲しいのは特異値分解だけなんですが、そこら辺に転がってるソース改変して使っても計算精度が悪くて使えないんです。
Octave C++もcygwinを通してのgccでしか使えないっぽいし、どうにかVC++で特異値分解をさせるライブラリのようなものお願いします。
194:名無しさん＠お腹いっぱい。
05/07/27 08:15:02
>>193
URLﾘﾝｸ(vision.kuee.kyoto-u.ac.jp)
これじゃだめ?
195:名無しさん＠お腹いっぱい。
05/07/27 14:59:39
(.v.)ｵﾚ
196:193
05/07/27 17:12:02
>>194
若干板違いだったのに、わざわざ有用な情報をありがとうございます。
結構色々調べたつもりだったんですけど、Octave C++をVC++で使えるとは知らなかったです。
Octaveを使うのあきらめてclapack使おうかと思っててたんですが、こっちの方が良さそうですね。
では、卒業研究頑張ってきますﾉｼ
197:193
05/07/28 02:03:21
なんかうまくいかねーorz
198:名無しさん＠お腹いっぱい。
05/10/08 07:15:57
Octave で補間や補外ってできますか?
199:名無しさん＠お腹いっぱい。
05/10/21 23:36:59
octave-2.9.3 に UMFPACK をリンクさせたいんですけど、どーやっても configure が libumfpack.a を見つけてくれません。
誰か助けてー。
200:名無しさん＠お腹いっぱい。
05/10/25 23:34:33
octaveの可視化機能は、どこまでMatlabに近いですか？
URLﾘﾝｸ(www.cybernet.co.jp)
こんな画像を出したいのですが、octave-forgeを入れていても、関数が無いようです。
201:199
05/10/26 00:15:58
>199
自己レス。
UMFPACKv4.4 の共有ライブラリを自分で作って CVS 版の octave で configure したら認識しました。
でも make がコンパイルエラーで止まる…今度はデバッグか…orz
202:199
05/10/27 06:28:25
>201
さらに自己レス。解決しますた。静的ライブラリでも大丈夫だった。
足りないライブラリを持ってきて Makefile をやっつけハックでいけた。
2048 * 2048 三重対角行列の LU 分解 - Pentium III Coppermine 1GHz です。
octave:1> a=eye(2048).*rand(2048)+diag(ones(2047,1),1)+diag(ones(2047,1),1)';
octave:2> t=time();[L U]=lu(a);time()-t
ans = 46.189
octave:1> b=sparse(eye(2048).*rand(2048)+diag(ones(2047,1),1)+diag(ones(2047,1),1)');
octave:2> t=time();[L U]=splu(b);time()-t
ans = 0.30429
(*´д｀*)
203:名無しさん＠お腹いっぱい。
05/12/01 11:28:50
portsからoctaveをインストールしようと思ったんだけど、
atlasのコンパイルに無茶苦茶時間かかりますね。
CPU毎に最適化したatlasのpackageってどこかに落ちてない?
204:203
05/12/01 11:31:15
>>203
ここのスレ的にはいちいち言わなくても問題ないと思いますが、
使っているOSはFreeBSD6.0Rです。
ここの板にいると、ついついみんなFreeBSDを使っていると思ってしまう。
205:名無しさん＠お腹いっぱい。
05/12/01 22:07:16
FreeBSD 6.0Stable, pentium-m 1.3GHzでAtlasのチューニングに36時間かかった
パッケージは保存したよ
デュアルブートでBSD使ってるからもうあんなのやりたくねえ
206:名無しさん＠お腹いっぱい。
05/12/02 00:03:18
今日はじめて使いました。
自分の作ったFFTのプログラムがあまりに遅かったので…。
ところで、起動時のフォルダって変更できるんですか？
微妙に不便
207:名無しさん＠お腹いっぱい。
05/12/03 12:19:24
？/.octaverc に好きなコマンドを突っ込んでおけば起動時に実行される
スタートアップファイルには他にもいろいろあって順番が決まっている
208:名無しさん＠お腹いっぱい。
05/12/03 20:42:54
>>207
なるほど。
ありがとうございました。
209:名無しさん＠お腹いっぱい。
05/12/07 14:24:02
しかし、Opteron dual な鯖で、昨夜の8時からatlasのインストールしているのに
まだ終わらない。
210:名無しさん＠お腹いっぱい。
05/12/26 17:28:06
liboctaveでodeをする方法だれかわかりませんか？
微分方程式がどけない、、、、、
211:名無しさん＠お腹いっぱい。
05/12/28 21:59:59
非線形だからMatlabとかという理屈がよくわからん
そんな非線形解なんてニュートン法でいいじゃん
ニュートン法自体はすげー簡単で単純
学部の頃はQuickBASICで書いてたぞ
奥村晴彦のコンピュータアルゴリズム事典にも載ってるし
ニュートン法でポテンシャルのパラメータを出して
インパクトパラメータを計算してた
212:名無しさん＠お腹いっぱい。
05/12/29 18:19:06
非線形っていうけど、Matlabって線形計算用なんだよね
もともとは...
213:名無しさん＠お腹いっぱい。
06/01/02 07:04:26
FreeBSDでoctaveを使いたいのだけど、atlasのコンパイルに時間がかかりすぎます。
それでpackagesを探しているのですが、FreeBSD6.0 or 6-stable用の
atlasのpackagesが見付かりません。
どこかにないでしょうか?
214:名無しさん＠お腹いっぱい。
06/01/02 10:22:57
ports にこう書いてあるから、自分で頑張る方が良いんじゃないかな...
NO_PACKAGE= Optimizes for the local machine, so a package is not generally useful. Set FORCE_PACKAGE to override
215:名無しさん＠お腹いっぱい。
06/01/02 13:35:50
>>213
オレ、年末にoctaveのportsをmakeしたけど
PentiumIII900MHzで30時間かかったよ
216:名無しさん＠お腹いっぱい。
06/01/04 10:13:40
LMItoolboxみたいのありますか？
217:名無しさん＠お腹いっぱい。
06/01/04 22:50:27
＞２１３
ATLASは実際にCPUにいろいろ計算させて、その結果を見て最適化していくんだよ。
だから、お前の希望するFreeBSD6.0 or 6-stable用、かつ、お前の使ってるCPU用のバイナリじゃなきゃだめじゃないの？
218:名無しさん＠お腹いっぱい。
06/01/04 22:55:46
まあatlasが必要ない用途も多いな。
219:名無しさん＠お腹いっぱい。
06/01/05 02:28:44
>>217
Atlas使うといっても何時間もかかるような数値計算をするわけでも
ないので、それほどカリカリに最適化してもらわなくてもいいんだけど...
220:名無しさん＠お腹いっぱい。
06/01/05 02:43:06
atlasってそんなにインストールに時間掛かるの？
openoffice位か？
221:名無しさん＠お腹いっぱい。
06/01/05 06:06:19
演算の最適化部分を省けばたいしたことないと思うけどね。
簡易インストールみたいなのないの?
222:名無しさん＠お腹いっぱい。
06/01/05 22:51:42
>221
最適化をあきらめる選択肢も必要だよね
チューニングに時間がかかりすぎるし
ATLASを使わないで、BLASだけでつかえたらいいのにね
223:名無しさん＠お腹いっぱい。
06/01/05 22:54:08
ところで、どうしてFreeBSDのPortsにあるATLASには-ffast-mathオプションが付いてるんだろう
はずさないとだめなんじゃないの？
224:名無しさん＠お腹いっぱい。
06/01/05 23:59:39
カリカリに最適化してくれるから俺はports
ひいてはFreeBSDを使っているわけだが
225:名無しさん＠お腹いっぱい。
06/01/06 00:03:13
>222
そういう make option の盛り込み案を maho たんに提案しる
226:名無しさん＠お腹いっぱい。
06/01/11 17:41:06 BE:81633582-
質問なのですが、save のファイル名の引数に文字列を保持した変数を
与えることはできるのでしょうか？
a=[1,2; 3,4];
filename = "test.dat";
save filename a;
とやると,"test.dat"ではなく"filename"というファイルに出力されてしまいます。
変数filenameの中身をみてくれる方法はないのでしょうか？
227:名無しさん＠お腹いっぱい。
06/01/13 23:14:55
>226
■俺も初心者。
とりあえず、save関数を使ってやってみてだめだった。つーか、マニュアルを見るとこうある。
たぶんだめなんじゃないのかな？代わりにfopenコマンドを見つけたから試してみたらどうかな？
１．fopenでファイルオープン ２．fputで書き込み ３．fcloseでファイルクローズ
perlのファイルハンドルみたいだし、""無かったら変数も展開してくれそうな雰囲気はあるよな。
俺も日頃からクソつまらない作業ばかりの社会人で試す前に寝てしまった。暇ができたらやってみるけど。
■saveコマンド
15.1.3 Simple File I/O
save options file v1 v2 . . . Command
Save the named variables v1, v2, . . . in the file file.
The special filename '-' can beused to write the output to your terminal.
■15.2.1 Opening and Closing Files
[fid, msg] = fopen (name, mode, arch) Built-in Function
fid list = fopen ("all") Built-in Function
file = fopen (fid) Built-in Function
For example,
myfile = fopen ("splat.dat", "r", "ieee-le");
fclose (fid) Built-in Function
■15.2.2 Simple Output
fputs (fid, string) Built-in Function
Write a string to a file with no formatting.
■これを参照した。
GNU Octave A high-level interactive language for numerical computations
Edition 3 for Octave version 2.1.x
228:名無しさん＠お腹いっぱい。
06/01/14 00:44:12
どうしてもsaveにこだわるのなら
eval(sprintf("save %s a",filename))
で可能。
こんな方法を使うメリットがあるかは疑問だが。
229:名無しさん＠お腹いっぱい。
06/01/14 01:29:15
>228
おおう
なるほど
230:名無しさん＠お腹いっぱい。
06/01/15 00:02:40
MatlabではM関数をエデットすればそれが次ぎの実行で反映されるのですが
Octaveでは一度読み込んだ関数はそのMファイルが書き変わっても反映されない様です
それでしょうがなく関数を書き換えたらexitしてoctaveを立ち上げ直しているのですが
関数を初期化する方法は無いでしょうか？
231:名無しさん＠お腹いっぱい。
06/01/15 08:34:26
clear all っぽい
でも何かエラーが出る
232:名無しさん＠お腹いっぱい。
06/01/15 13:00:18
x="output.txt"
y=[1,2,3,4]
save("x","y");
ってなかったけ
233:名無しさん＠お腹いっぱい。
06/01/15 20:07:34
>>230
マニュアルには更新日時を見て
再読込するよな事が書いてあるけど、
実際の動作は違うのか？
help -i reload
>>231
clear allは普通にできないか？
>>232
できない。
234:名無しさん＠お腹いっぱい。
06/01/21 18:48:32
Octaveで偏微分方程式は解けますか？
235:名無しさん＠お腹いっぱい。
06/02/07 13:45:49
octaveとscilabってどう違いますか？
初心者にお勧めはあります？
236:名無しさん＠お腹いっぱい。
06/03/05 16:18:36
cygwinでscilabをコンパイルできた人居る？
cygwinのコンソールでscilabを使いたいよ
237:名無しさん＠お腹いっぱい。
06/03/12 21:25:19
Windows用 Octave 統合環境
Octave Workshop
URLﾘﾝｸ(www.math.mcgill.ca)
インストールしたらちゃんと動いた！
でもグラフは gnuplot for win
エディタはメニューの File で open すれば立ち上がる
これの凄い所は全部 MinGW で作ってある所！
Cygwin 無しなので、UNIXになじみの無い職場でも使えるかも
238:名無しさん＠お腹いっぱい。
06/03/13 16:40:29
そういえば、Octave Workshopって、FreeBSDのportsに入っていないね。
239:名無しさん＠お腹いっぱい。
06/03/16 14:32:49
floating ｐｏｉｎｔの形式はｉｅｅｅじゃなきゃだめなのね。
故に、-mfancy-math-387 -ffast-ｍａｔｈなんかつけるとこけてしまう。
240:名無しさん＠お腹いっぱい。
06/03/17 23:32:02
FreeBSDのatlasのportsのfiles/{blas.h,atlas.h}って
+#ifdef __cplusplus
extern "C"
{
+#endif
--略--
+#ifdef __cplusplus
}
+#endif
なpatchをあてないとc言語の時に困るので
誰かsend-prしてください m(_ _)m
241:名無しさん＠お腹いっぱい。
06/03/17 23:52:31
cygwinでoctaveとatlasをインストールしようとしてたやついる？
gcc3.4でoctaveをコンパイルすると遅くなるって言ってたけど
gcc4.0とか普通にコンパイルできるぞ？
今からこれでoctaveをコンパイルしてみる
242:名無しさん＠お腹いっぱい。
06/03/18 06:12:02
>>241
ａｔｌａｓしらんが、ｏｃｔａｖｅは出来そうだな。
liboctave/lo-ieee.ｃｃかな？にCygwin_ｉｎｉｔとかいう関数がある。
243:名無しさん＠お腹いっぱい。
06/03/20 21:09:03
SciLabは4.0が出たみたいだな
使いにくいから使わないけど
244:名無しさん＠お腹いっぱい。
06/03/21 01:11:46
octaveの方がいろいろルーズでいいな
245:名無しさん＠お腹いっぱい。
06/03/22 12:30:18
>>244
確かにな
mファイルをコマンドとして使えないとか
関数mファイルを明示的に読み込まないといけないとか、、etc..
Octaveみたいに操作方法もMatLab互換にすればみんな使うのに
246:名無しさん＠お腹いっぱい。
06/04/03 00:21:50
俺はLMI tool(MATLAB非互換)が使いたいからSciLabも使うけどね。
慣れればSciLabもいいよ。
247:名無しさん＠お腹いっぱい。
06/04/05 23:06:34
octave-forgeのfillは頂点が4つまでしか描けないの？
有限要素法の研究をoctaveでやらにゃならんのに
8頂点の要素がかけないのはつらい・・・。
情報をください。ｴﾛい人！！
248:名無しさん＠お腹いっぱい。
06/05/14 07:20:32
>241
ATLASはgcc-3.3.3でmakeしました。3.4でもいいはずなのですが
ファイルの一部がmakeできなかった？（よく覚えていない）
octave以外はgcc-3.2でmake
(3.4はbroken octaveになるそう。3.3.3は例外処理が遅いため遅い。
cygwin1.dllの一月のverupより遅さが緩和されているような気がしますが)
gcc-3.2は普通にはgetできないので
www.nuce.nagoya-u.ac.jp/e8/Matsuoka/OctCpp.html
を参考にもってくる。
249:名無しさん＠お腹いっぱい。
06/05/14 07:23:26
>248 自己フォロー
>octave以外はgcc-3.2でmake
ごめん。
octaveはgcc-3.2でmake
の間違い
250:名無しさん＠お腹いっぱい。
06/05/20 12:12:17
Octave Workshop (>>237) を使っています。
起動直後に私のホームフォルダに自動的に移動しておいて欲しいのですが、
どのファイルにコマンドを書き込めばいいのでしょうか。
あちらこちらにある設定ファイルらしきものに書き込んでみたのですが
まったく変化がありませんでした。ヒントをもらえると嬉しいです。
251:名無しさん＠お腹いっぱい。
06/05/20 14:06:11
>250 自己解決しました。
普通にリンクファイルでフォルダを設定すればＯＫでした。
失礼しました。
252:名無しさん＠お腹いっぱい。
06/06/25 16:16:16
>>40のleastsqr.m を動かそうとしているんだけど、なかなかうまくいかない。
なにをやったか書くので、なにがいけないのか、分かる人がいたら教えてください。
x=[1 2 3 4 5];
y=[1 3 6 7 8];
function ret=F(x,p)
ret = p(1)*x*x + p(2)*x + p(3);
endfunction
pin=[1 1 1];
[f,p]=leasqr(x,y,pin,"F")
といった形で実行させると、
error: input(x)/output(y) data must have same number of rows
といったエラーが出る。leastsqr.mのソースを見ると、
% check data vectors- same length?
m=length(y); n=length(pin); p=pin;[m1,m2]=size(x);
if m1〜=m ,error('input(x)/output(y) data must have same number of rows ') ,end;
といった箇所があって、m=length(y)と[m1,m2]=size(x)のm1を比較している。
つまり、上の例だとmが5となり、m1が1となるので、両者がマッチしない。
だから、エラーになったということのようだった。
このあたり、データの与え方が間違っているのかと思って転置行列 x' で
実行させてみたりしたけれど、なかなかうまくいかない。Octave のバージョンは、
GNU Octave, version 2.1.53 (powerpc-apple-darwin7.5.0).
GNU Octave, version 2.1.69 (i386-pc-linux-gnu).
両者とも同じ結果だった。
次ページ最新レス表示スレッドの検索類似スレ一覧話題のニュースおまかせリスト▼オプションを表示暇つぶし2ch
5162日前に更新/221 KB
担当:undef