データ解析ツールoctaveを語ろう

知る人ぞ知る強力なデータ解析ツール octave ...

unixを科学研究で使っている香具師には非常に心強い味方のはずだ。

matlab互換を主張しながら、微妙にバージョンアップを続ける
octaveについて語ろう。

リンクを貼るとか、

せめてものネタふりとして自分の使用感を晒すとか、

　今度、卒研で使わざるを得なくなりました
　インストールのしかたから教えてください

とか、

書くことがあるんじゃネーノ？ｗ

>>2
どうも申し訳ない。

http://www.octave.org/

がリンクです。

実は、研究でmatlabを使っていたんですが、
サポート料金を払う金がおしくなって octave にしました。

octave 意外と使えるので、というか軽いし gnuplot との融合性が
高くて本家のmatlabよりも使いやすいのではないかと
思っています。

もちろん、パッケージなんかを考えるとmatlabの方が良いのだけど
私の使い方だと、高価なパッケージ使うこともあまりないので
それほどデメリットにはなりませんでした。

あくまで数値解析ソフトなわけで、数式解析ソフトな matlab
と比較するのはどうかと思う

>>4
mathematica と間違えていませんか?

matlab で数式解析するという話は聞いたことがありませんが?

私もFreeBSD上で octave 使っています．
最近では数学関数も増えているようでけっこう使いやすく
なっています．

非線型の最小自乗法が扱える関数があればベストなんだけどね．

matlabの影響力は日本は少ないからね

octaveいいね。卒論で他人がFFTやらのプログラム書いてるときに、
俺はoctaveで楽させてもらった。

>>5
その通り。失礼した
あんまり使ってないけど、行列計算の検算とか変形に便利ではあった

octaveって使ってみたのですが，
行列計算ってそんなに必要なんですか？

>>9
失礼、上げてしまった。

必要なところには必要、不要なところには不要。

scilabというのもあるので、どっちがいいか比較キボン。
http://www-rocq.inria.fr/scilab/

直接LAPACKを使う場合と比べてどの位遅くなるのでしょうか

>>12
scilabの方がmatlab互換という点では良いかも知れないけど，
純粋に線形解析に使うのだったらoctaveの方が使いやすいと思います．

私的感想だと数学関数はoctaveの方が良いと思う．

>>13
速度に関してはインタプリタの分だけ遅いという感じですね．
最近のoctaveは中間コードに変換しているような感じなので
それほど速度の低下は感じられません．

octaveと数学関係の外部ライブラリの間をいったり来たりするような
感じの処理だと遅くなるようですが...

Perl とか awk とかでやるにはややめんどい数値処理
(FFTとか統計処理とか) にいい。liboctave というのも
あるようだが未見。

科学関係のソフトがわりと丁寧に整理・紹介されている。octave もあるよ
http://www.mlb.co.jp/linux/science/
http://www.mlb.co.jp/linux/science/octave/index.html

>>16
けっこう参考になるリンクですね。

最近は、貧乏研究室に配属になったので、Open source系の
科学技術計算ソフトの情報はありがたい。

最近は、FreeBSDにoctaveをインストールしてデータ解析の毎日です。
しかし、portsから入れたらoctaveのインストールに時間のかかること．．．

昔学生だったときに三田村さんのページをよく見た。

…と思ったらまだあった。:-)

>>18
だれそれ?

オープンソースなツールを使って科学技術計算やっている人々の話が聞きたい．

ほとんどが学術系だと思うけど企業とかの研究者でもoctaveとか使っている人
いますか？

ところで画像解析にmatlabを使っている人多いと思うけど，
octaveで画像解析している人いますか？

どうしても画像になると octave　だと若干貧弱になると思うのだけど...

特に表示するときとかインデックスドカラーなので
カラーパレットの設定などが大変なんですが，
うまくやっている人いたらコツを教えてください．

>>4-5
matlabにも数式解析機能はある(toolboxとして)。
その目的のために買う人は少ないだろうけど。

>>20
技術屋としての数値計算に僕が主に使うのはemacsとgcc。
matlabとoctaveも使っているけど、自分で使う道具としてではなくて
プログラミングをあまり知らない人とのやりとりのために使っている。

Cの読み書きができない人(客)でもMatlabやOctaveのコードなら自分でいじれ
ることが多い(オレのまわりでは)。

>>22

数値計算アルゴリズムを開発する場合，octaveみたいに試作が直ぐにできる
ツールを使った方が効率が良くないですか？

最終的にはC言語で実装するにしても，アルゴリズムの開発で
試行錯誤を繰り返している段階の場合は，やはり最初は
octaveのようなツールが有効だと思いますが...

>>23
＞試行錯誤を繰り返している段階の場合は，
このあたりは主に紙の上でやります。

>>24
それは凄い．

私はコンピュータの上でしか試行錯誤ができなくなってしまっています．
コンピュータに頼りきるのも考え方が雑になってしまうという
弊害もあるのかも知れないけれども便利ですよ．

私のところのoctaveの使用目的というのは，
(1) 特にソフトウェアは必要なくて，
　　対象（物理系の実験データ）の解析結果データのみが
　　必要なとき

(2) 新しいアルゴリズムを試作するとき

の２つの時ですね．

>>25
アルゴリズムとして定式化する以前の試行錯誤のところじゃないですかねぇ
もちろん、分野によって異なる計算によるわけですが。

私もoctave+gnuplotはかなり便利に使わせてもらってます。

Octaveで非線形計算やるパッケージをご存じないですか?

誰か作ってそうなんだけど...

非線形方程式を解くコマンドはfsolveだったはづだけど、計算..？

みんながOctaveをどのように使っているか、聞きたい。

私は、某工学分野の研究員だが、
Octaveは、アルゴリズム設計の前段階で使うことが多い。

ちょっと変な使い方かも知れないけど、Octaveでデータを統計解析して、
その傾向に合わせてC言語でプログラムを組んでいる。
つまり、直接アルゴリズム開発に使っているのではなくて、
アルゴリズムを開発するときの参考というか理論的な根拠にしている。

こういう使い方にはけっこう向いていると思うのだが。

>>28

非線形の最小自乗法が欲しい。
fsolveは、微分方程式を解く関数じゃなかったっけ？

>>30
gnuplot

>>30
そのへんはRとかどうなん？

>>30
でもあれはデータフィッティングしてグラフをプロットするだけでしょう？
ちゃんとパラメータとかを求めたいってことだと思うのだが...

スマソ，リンク間違えた

>>31 ね...

>>32
Rにしちゃうと Octaveの他の恩恵が受けられなくなる...
という問題が...

Octaveをインストールするとき，FreeBSD 4.x だと
readlineのバージョンが古いっていって libreadlineが使えません．
コマンドライン編集や履歴を使いたいのでこれでは
不便なのですが...

それで手動でlibreadlineの新しいやつをインストールして
使っているのですが，octaveのバージョンが変ると
パッチ当てとかがめんどうです．

FreeBSD 5.2くらいなら libreadlineのバージョンが
新しいので問題ないのだけど...

RはしらんけどSはやたらdirect mailを送ってきた印象があるな

非線形一般均衡系のシミュレーションに使おうとしてうまくいきまへんでした > octave
で､MATLAB買っちゃいました｡toolboxマンセｰ

>>37
こういう分野では MATLABというか サードパーティ製のtoolboxが
かなり強力だよね．

でも価格が大変だよね．
貧乏な大学の研究室だと変えないところが多そうだ．

うちも非線形解析用のパッケージ買っていたけど，
本体のライセンス代金も含めて大変なことになっていた．

octaveはプログラム組むのに、１ファイル１関数っていうのが
不便ですね。
最近のバージョンだとこの制限なくなっているのかな？

MATLABも昔はそうだったけど、最近のはそういう制限無くなっているし

>>30
FAQだとよ。
http://www.octave.org/octave-lists/archive/help-octave.1999/msg00770.html

>>40
これは良いかも知れない．
丁度Levenberg-Marquardt法のアルゴリズム探してました．

あとは，これをパラメータの範囲指定ができるように改造して
使わないと解がlocal minimum に陥ってしまって
収束しないんですよね．

もしかしてそれもあったりして...

cygwin で liboctave 使えますか?
下のようなテストプログラムでエラーになってしまいます.

$ cat test.cpp
#include <iostream>
#include <octave/config.h>
#include <octave/Matrix.h>

using namespace std;

int main()
{
Matrix m(2, 2, 1.0);
cout << m << endl;
return 0;
}

$ g++ -I/usr/include/octave-2.1.49 -L/usr/lib/octave-2.1.49 -loctave -lcruft -lreadline -ltermcap -llapack -lblas -lg2c -lfftw -lm test.cpp -o test
/cygdrive/c/DOCUME~1/user/LOCALS~1/Temp/ccdhMv8V.o(.text+0x8e):test.cpp: undefined reference to `operator<<(std::basic_ostream<char, std::char_traits<char> >&, Matrix const&)'
/cygdrive/c/DOCUME~1/user/LOCALS~1/Temp/ccdhMv8V.o(.text$_ZN6Array2IdED2Ev+0xd):test.cpp: undefined reference to `Array<double>::~Array [not-in-charge]()'
/cygdrive/c/DOCUME~1/user/LOCALS~1/Temp/ccdhMv8V.o(.text$_ZN6Array2IdEC2EiiRKd+0x1b):test.cpp: undefined reference to `Array2<double>::get_size(int, int) const'
/cygdrive/c/DOCUME~1/user/LOCALS~1/Temp/ccdhMv8V.o(.text$_ZN6Array2IdEC2EiiRKd+0x33):test.cpp: undefined reference to `Array<double>::Array[not-in-charge](int, double const&)'
collect2: ld returned 1 exit status
make: *** [test] Error 1

>>42

octave は知らないが、与える引数間違ってるんじゃないの？

>>43
どう間違っているのでしょうか?
Linux と同じじゃダメなのでしょうか?

素晴らしーいソフトウェアです。
Fedora CORE 2 でSMPマシンにインストールしました。
コンパイルに１時間強、出来上がったバイナリ− octave が60MBを越えたのには
驚きましたけど。
紹介していただいたみなさんに感謝感激です。
しばらくマイブームとなりそうです。

>>42

俺が試した限りでは、

octave-2.1.57 を ./configure; make install した状態で、

$ g++ -I/usr/local/include/octave-2.1.57 test.cpp -L/usr/local/lib/octave-2.1.57 -loctave -lcruft -lg2c -o test

でコンパイルできたが。

>>45
60MBは、たぶん -g つけたままコンパイルしたんじゃないか？
liboctaveも巨大なままでは？
make clean; make CFLAGS=-O2 CXXFLAGS=-O2 をおすすめする。
あと、strip octave。

>>46
上ので通りました. -llapack とかは不要なんですね.

$ g++ -I/usr/local/include/octave-2.1.57 -L/usr/local/lib/octave-2.1.57 -loctave -lcruft -lg2c -o test test.cpp

こう書くと undefined reference がいっぱい出ますが,
順番が大事だとわかりました.

どうも, ありがとうございました.

octave ってスパース行列扱うの弱くない？

最近はCPUもメモリも強力だけど、古いマシンで
計算させたら死にそうになった。

Scilab使へば

がびーん。
./configure; make しただけだったのですが、そういえばご丁寧に -g
が付いていたような。。。
気を取り直してもう一度 build　してみます。
その前に HDD 増設かも。

連続書き込みすみません。
octave って SMP マシンの恩恵を得られるのでしょうか？
HT対応のＣＰＵを２個挿せば、体感？４つになります。
個人でも 4 way も持てるなんて、想像しただけで幸せです。

octave:1> A=rand(2000,2000);b=rand(2000,1);
octave:2> tic;A\b;toc
ans = 2.7709

最近のマシンは速いねぇ…

>>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.

>>52

多分ならなかったと思うけど...
threadで組んでたかな?

>>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よりも遅いのですか...??

>>56
書くの忘れたけど，
OSは，FreeBSD 5-current で
CPUは 600MHz x 2のSMPです．

ちなみに 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 と比較してもクロック分の速度は出てるかな...

>>49

それってかなり古いOctaveじゃない?

octave:1> A=rand(2000,2000);b=rand(2000,1);
octave:2> tic;A\b;toc
ans = 5.6739

最近のB５サブノートは速いねぇ…

Octaveって疎行列の計算はできますか？
周期境界条件を課した疎行列の計算パッケージがなかなか見つからなくて困っています。
FFTはNlogNのオーダーなんで遅いしなーって思うとNのオーダーで計算できる
疎行列パッケージが一番よいのですが。

>>60
cygwinでコンパイルしたやつだけど、
ans = 34.755
だった orz
コンパイルオプションばりばりつけた方がいいのかなぁ

>>61
自分で書くと感謝されるよ。

でもLAPACKとかで疎行列扱えないの？
最近はOctaveも外部ライブラリに頼っているからね。

うちの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

という値を出しました．

漏れも試そうとportsでoctaveをmakeしてるんだけど
atlasのビルドが一日たっても終わんないよう。

atlasのインストール時ってコンパイルが終わらないというより、
数値計算を行っている時間が長いようですね。

ずっと数値計算プログラムが動いていますが、
計算精度のチェックでもしているのでしょうか？

ちなみにうちのThinkPad X31だと、とりあえず１日ってことはないな
３時間くらいで終わっていたように思えるけど...

>>65
メモリが不安定だと計算処理が途中で異常終了して
計算を繰り返して丸一日以上かかることがありました．

もしも速いCPUを使っていて丸一日以上かかるようだったら、
メモリをチェックしてみた方が良いかもしれません．
普通に使っていて動いていても memtest とかが
通らない場合がありますよ．

私は atlas のコンパイル作業でメモリ異常に
気が付きました．

ちなみに Athlon等の1GHz以上のCPUなら数時間で計算処理が
終わります．

なにせ古いマシン、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 中です。

どんなプラットフォームでもちゃんと動作するマルチスレッドの
プログラムを作るのは難しいですからね。ＳＭＰを活かそうと思ったら
やはりガリガリ自分で作るしかなさそう。嬉しいような悲しいような。。。

>>68

110秒はちょっと遅いですよね？
これはコンパイルオプションにかなり左右されるのでしょうか？
SSEとかどうなっているのでしょうか？

>>56 は　同じようスペックで17秒ちょいで終わっているようですが...

研究環境を完全OpenSource化使用して、いろいろ画策中です。

特にMATLABからOctaveへ関数をいくつか移植しているのだけど、
けっこう大変だよね。

Octaveは関数が一つしか書けないし、微妙なところで関数の仕様が
違っていたりするし。

でもMATLABへのライセンス料を払わなくて済むなら
かなりの出費が押さえられる。

　　　　　　　　　　　 　　 ＿＿＿＿
　　　　　　　　　　　.／, - ､, - ､　　￣ ヽ
　　　　　　　　　　./-┤　｡|｡　 |――-､　ヽ
　　　　　　 　 　　|　ヽ｀- ○- ´ ／　　ヽ　 |
　　　　　　 　 　　|　-　　 |　 　 　―　　 |　 |
　　　　　　 　 　　|　 ´　　|　 　　｀ヽ　 .　|　 |ヽ
　　∩　　　　　　人｀､ _ 　|　 　　＿.- ´ |　.|　　＼
　　|　　⌒ヽ　／　 ＼　￣￣￣　 　　 ノノ 　 　 　 ＼
　　|　 　 　 |´　 　 　　|￣―--―― ´ヽ 　 　 ＿　　／⌒＼
　　＼_ 　 _/-―――.|　（ Ｔ ）　　　　　　｀l 　 　 Τ（　 　 　 ）
　 　 　 ￣　　　　　　 |　　￣　　　　　　 　}　 　 　 |　＼＿／
　　　　　　　　　　　　 | ､--―　￣|　 　　/

>>71
もしかしてキリ番狙っていたの?

　　 ||
　Λ||Λ
（ / ⌒ヽ
　| | 　　| 逝け＞７２
　∪ 亅|
　 |　|　|
　 ∪∪
　　 ：
　　 ：

　‐ニ三ニ‐

　　 ||
　Λ||Λ
（ / ⌒ヽ
　| | 　　| 逝け>>７２
　∪ 亅|
　 |　|　|
　 ∪∪
　　 ：
　　 ：

　‐ニ三ニ‐

>>71, >>73, >>74

せっかくのAAなんだけどうちの環境，FreeBSD + mozilla デフォルト設定では
文字がずれて，きちんと見えません...

等幅フォントにすりゃいいんだけどね．

このスレせっかくいい調子で話が進んでいたのに、
変な香具師が紛れ込むようになってきたね。

必要もないのにageちゃうからだろ。

さりげなくサイエンスな香りがするこのスレ、気に入っていたのに

WindowsでOctaveって動きますか？

ちょっと面白そうなので立ち寄ってみました。

>>79
動くよ。native 版もあるし、cygwin で動かす手もある
特に数値計算なら、こしょこしょコード書くよりまかせてしまったほうがいいと思う

レスありがとうございます。

Windowsのnative版もあるんだ。
探してみます。

でもcygwinはややこしいからパスします。

途中報告。
-g を抜いたバイナリでSMPカーネルで同じ手順を行ないました。
バイナリサイズはナント 7MB でした。
実効結果は
ans = 105 くらいでした。大差ないです。
原因を引き続き調査中です。

そうそう、SPARC とかは float より double が早いそうですね。
ハードウェアがそう設計されているそうです。

>>85
floatよりdoubleが速い処理系は多いよ．特に昔のRISC系CPUの場合は...

それと表示は，桁を省略しているけど，Octaveは基本的に double で
計算しているはずです．
それで105秒だと本当にそれくらいの速度かと...

atlasをコンパイルするとき -g を抜くだけではなくて，
gcc で　-O3 とか　-msse2 といった最適化オプションは
指定してありますか？

>>83

FreeBSDのportsではSSE2対応のCPUでは，自動的に
SSE2が有効になるようにコンパイルオプションが
指定されているようです．

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 っぽい。

よく覚えていないのですが、atlasってアセンブラで
SSE2のコード書いてませんでしたっけ？
たしかプリプロセッサでSSE2を有効にする変数を定義したら
SSE2が有効になったような...

Opteron242だと
ans=11.8
やけにおそいな。
OS SuSE 9.1 X86_64

>>87
B5サブノートのPentiumMに負けている...

Opteronって通常の数値計算が速いって聞いてますが、
実際のところどうなんでしょうか？

Octaveの使用を含む数値計算用に一台Dual CPUなのを購入予定なのですが...

opteron 242*2(1.6GHz)
姫野ベンチ　2CPU
1313MFlops
Pentium200MHz 40.7倍

Xeon　2GHzしかないが
819MFlops

今一つ速くないね．Opteron...

まだ最適化が十分ではないのかな？

あげ

しかしながらベンチマークスレに成り下がっているのは情けなくもある。

演算の精度の制御はどうやればできますか？
たぶんデフォルトではdoubleというか64bitっぽいですが、
もっと精度が欲しいのです。

scilab

>>92
スーパーコンピューターの最大ユーザーはベンチマーク屋である。
という話を聞いたことあるけど。
もっともコンパイルオプションの違いで数倍のパフォーマンスが
違えばショックですね。
MATLABを買ったと思ってＰＣ買い換えれば幸せになれるかも。。。

>>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のフォーマットが使える。

>>94,>>96
どうもありがとう。octaveじゃ無理なんですね。

scilabをインストールしてみたら、例のペンチマークがoctaveより速くて、
のりかえを検討していますが、こっちはどうやれば精度の制御ができるんでしょう？

octaveは全く触ったことがない厨房ですが、octaveとRは機能的に
かなりかぶりますか？　それとも相互補完するものですか？

若干かぶるけど，特定分野だとRの方が強いと思うよ...

Octaveは線形代数系の計算全般といったところか．

>>95
MATLABのことを考えると，毎年PCを買い換えれますね．

MATLABは基本的に毎年ライセンス料を払わないといけないプレッシャーが
あります．
まぁ，払わなくてもいいんだろうけど，修正パッチとか
バージョンアップのことを考えるとライセンス料を払い続けるしか
ないシステムになっています．

>>95

ベンチマークだいすきですぅ。

MATLAB個人で買いましたよ。
octaveじゃできないこともおおいし。

ちなみに「データ解析ツール」じゃなくて「行列言語」な。
間違えるなよ。

行列言語ってのもなんだな...

Octave使っているのだけど，マニュアルとかが英語なのがちょっとね．
何かお奨めな参考書とかWebページとかりますか？

とかありますか？

あまり役に立たない本
「Octaveを用いた数値計算入門」
北本　卓也著
ピアソンエヂュケーション
ISBN4-89471-448-5

卒研で構造物の強度計算をすることになりました。
それで手頃なツールをさがしているのですが、
Octaveでできますか？

研究室のOSはNetBSDというものだそうです。

強度計算がどんな計算をするものなのかわからないので、わからない。
行列計算とか微分方程式とかに帰着させられるんならできるだろうけど。
OSは無関係。

初めて立ち寄ってみました。
lsodeの中身は4段のRunge-Kuttaなのでしょうか？
何を見たら分かるのでしょう？どなたか教えてください。

微分方程式は Runge-Kuttaを使っていたと思います．

今のOctaveは外部の数値計算ライブラリに依存しているから
そのソースコードを確認するのがいいのではないでしょうか？

でもFORTRANとかで書かれているのでちょっと困りますが...

>>109
lsodeはRunge-Kuttaではなく、stiffな問題にはBDF法、nonstiffな問題にはAdams(predictor-corrector)と使い分けています。
いずれもRunge-KuttaのようなSingle step methodでは無く、Multi step methodと呼ばれる方法です。
http://www.llnl.gov/CASC/odepack/

もともとODEPACKというFORTRANのパッケージに含まれていて、ソースはwww.netlib.orgからダウンロード出来ていたと思います。
lsodeでググってみたら出てきますよ。

二つの異なる情報ですが、ありがとうござました。
linspaceで刻み幅を固定して計算するのだからRKだろうか
と思っていました。ソースを読んで見ます。
BDFかあ。懐かしいなあ。ずっと以前でしたが、理解するのに
一苦労したっけ。小生は、現在は、刻み幅自動調整の
RKFの愛用者です。

異なる内容でしたが、ありがとうございます。
linspaceで刻み幅を固定するので、RKかと思っていました。
小生は、現在は、刻み幅自動調整のRKFを愛用しているので
つい、そう思っていました。ソースを読んで見ます。

Matlabのtoolboxが必要なら、octeva-forgeをインストすれば

>>114

こんなんあったんだね。
これはちょっと使えるかも。

さてMATLAB R14が出たわけだが。

MATLABは新しくなればなるほど重く不安定になってくる。
R12以前くらいかな良かったのは。

久々にベンチマークネタ 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も古いし...

octave-forgeってどーですか

>>119
けっこう使えるよ。
公開されてるMatlab用のプログラムが変更無しで動く

octave-forge ってそもそも何ですか?

>>121
cat /usr/ports/math/octave-forge/pkg-descr

あげ

あげ

日本語化パッチはまだでつか？

Octave に日本語って必要なんですか?

必要かどうかは解かりませんけど、
道具が人を選んではいけないと思います。

>>127
意味分かりません

英語読むのめどい

馬鹿とハサミは使い方次第ですね
馬鹿には道具は使えないんでしょうかねぇ

んなこといっちゃってはずかしいくせに

Octaveって線形代数電卓って感じだなぁ．

普段何気なく使っていて，それほど必要だとは思わないんだけど
自分の使っているマシンに入ってないと気づいたときのショックは
計り知れない．
無いとかなり不便なんだよね．

FreeBSDの ports の Octave 新しくなりましたね．
何が変わったのかな?

2.1.63にアップグレードしたら次のメッセージが出るようになった
warning: lo_ieee_init: unrecognized floating point format!
何でしょうねこれ

自己レス
最適化がきつすぎたみたい
コンパイルオプション変えたら出なくなった

スレ違いでしたらすいません。
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とし、同じ様に解を求めます。

CPUとマザーボードを買い換えました．
AthlonからPen4に変えたのですが
octave って再コンパイルが必要だったんですね．

-march=athlonとかやってたんでないのか?

Octave は関係ないんだけど atlas がコンパイル時にSSEや3D Now を
自動的に認識して最適なコードを出力します．

Athlonだと3D Now が入ってしまうから Pen4だと動かなくなってしまいます．

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.

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

portsから入れるならWITH_OPTIMIZED_FLAGSを定義しておく
その他CFLAGSやFFLAGS等も適切に設定しておく

>>141
SSE2の速度はPentiuMとPen4でぜんぜん違いますね。
>>142
WITH_OPTMIZIZED_FLAGSは定義しないでください。gcc-3.4.2の
バグで無限ループに陥ることアリ。CFLAGS/FFLAGSは意味なしです。

意味無いこと無いだろ

mahoたん降臨とみた。

>>143

Pen4 ってSSE2 速いんですか?

それほどでもないと思うのですが...

>>64 で X31 のベンチがありましたが
こちらも同じマシンですがベンチの結果はほぼ同じです。

同じ環境で PenM 1.4GHz 5.5秒 Pen4 2.8GHz 4.3秒...
これをどうみるか...

>>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のほうが同クロックだと相当高速ですよ。

>>148

プレスコットでも北森でも同じなのでしょうか?

研究室で数値計算用に一台PCを組もうと思っているんですが
世間のPen4の不人気具合いからためらっています．

とりあえず FreeBSD Expert 買ってみます．

>>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というのに僕は期待しています。
ノウハウをみなで共有できたら良いと思いますし、よい記事執筆者
が出ることを期待もしています。

少しスレ違いな話で申し訳ないのですが

うちのプレスコな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%
でした．

正直いってプレスコ使うのが嫌になってきています．

>>151

ごめん名前間違えた 143 ではなく 140です．

すいません。
お聞きしていいですか？？
現在、octaveで１次元有限要素法のプログラムを書いているのですが、量子井戸の波動関数を求めようとした際に、どうしても解が振動してしまいます。
波動関数の２乗の値はまともなのに、２乗しないとだめです。
どうしてでしょうか？？

>>153
ネタ？

>>153

> 波動関数の２乗の値はまともなのに、２乗しないとだめです。

ここの意味がわかりません．

>>153
とりあえずソースを

固有エネルギー・固有状態がほしいのか
波動関数の時間発展をみたいのか
何をしたいのか明確に説明プリーズ

すいません。
説明不足でした。
波動関数が振動するとは、１要素ごとにプラスとマイナスが入れ替わってるような、ギザギザした形状になってしまうということです。
以下にソースを載せますので、暇であればお試しください。
欲しいのは時間に依存しない方程式を解いて、固有関数と固有値です。

http://cyobe36.k-server.org/newpage1.htm

ttp://cyobe36.k-server.org/newpage1.htm

>>157
ちゃんと解読してないけど
固有値が一番大きな方の波動関数をプロットしてない？
だったら１サイトごとに振動してても問題ないかと。
固有値が小さいやつをプロットしてみたらどうでしょう。
ただ、plot(x,v(:199))をしてみたけど
ポテンシャル井戸に局在してるようなモードでは無かったです。
対角化したときに固有値・固有状態がどういう順番で出るのか良く知らないんですけど。

ATLASってiccでビルドするとどれくらい速くなるんでしょうか

Octave GTK+
http://octave-gtk.sourceforge.net/

すいません、宣伝になっちゃいますけど

liboctaveの使い方をまとめたLiboctaveWikiを立ち上げました
http://liboctave.halfmoon.jp/
ちょっとliboctave使ってみたいなっていう初心者の方から、
liboctaveなら俺に任せろ！な猛者まで、
Wikiの編集に参加していただけるとうれしいです。
なんでOctaveじゃなくてliboctaveの方なんだっていうのは無しの方向で・・(汗

age

sage

% 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読めってことですか。

FreeBSDではもうoctaveは使えません

>>166
(´･ω･`)ｼｮﾎﾞｰﾝ

>>166
マジ？

だめじゃん

焦って/usr/ports/distfilesを見たけどoctaveのファイルは残ってなかった
誰かFreeBSD5.4R用のi686バイナリパッケージくれよ

こういう時こそportdowngradeの出番では

FreeBSDのportsってたまに Octave が Broken になりますよね．
またしばらく待っていると復活しますよ．

>171
joke かと思ったら本当にあった... sysutils/portdowngrade

mahoﾀｿがその内直してくれるんじゃね

なおった？

少なくとも一昨日の段階ではまだ直ってなかったけど...

mahoタソたのむ

5.x以降なら直ってるみたい。
4.xはBrokenだけど。

漏れも使っているので、メンテ頼む mahoさん

matlab だと、
x{1} = 1; x{2} = 2;
x{1} = [];
で、x の最初の cell を消去できますが、octave ではどうするんでしょうか?
同様に structure array の要素の消去の仕方も、[] を代入では駄目でした。

すみません。嘘付きました。
cell array の要素の消去は、[] を代入でできるんです。
しかし、structure array の要素の消去が分かりません。

解決しました。
web では見つけられませんでしたが、よく考えれば分かりました。

for i=1:10; x(i).f=i; end
とし、
x(7) の削除は、
x = x([1:6 8:10]);
ですね。

今、FreeBSD で portsから octave を入れたのだけど
微妙に良くなっていっているよね octave

メンテや改良ってまだどんどん行われているの?

FreeBSD 5.3でportsからoctaveを入れようとしたけれど、
===> octave-2.1.63 is marked as broken: Does not compile on FreeBSD >=5.x
って出て、インスコできなかったよ。

これってFreeBSD 5.x以上では使えないってことなの？

>>183
portsを新しくしてください。

>>184
portsは新しくしたはずなんだが…ｵｶｼｲﾅ

ttp://adlib.rsch.tuis.ac.jp/~akira/unix/octave/update.html
に2.1.Xから多次元配列がサポートされたって書いてあるけど，具体的にどう使えるのでしょうか？
ChangeLogとかNEWSとかにも multi とか　Arrayとか検索しても引っかからないし，
zeros(3,3,3)とかも使えないし…．
どなたか教えてください．

octaveじゃないけど
math/sdpaも壊れてるような

>>186
バージョンは？
2.1.50ではだめだけど、2.1.60以降ではzeros(3,3,3)でいける。

>>188
2.1.71にしたら出来ました！ありがとうございました！

imageコマンドで起動される外部画像ビューアーが、インストールされた計算機によって、
imagemagicだったり、xvだったりするのですが、この違いはなぜ生じるのでしょうか。
また、明示的にどちらかを使いたい場合、どうすればよいでしょうか。

printf('kita---(.v.)-----!'); sleep(10);
これが１０秒後まで表示されないんですよね
どうしようもないですか？

>>191
printf('kita---(.v.)-----!'); fflush(stdout);sleep(10);

ちょっと板違いなんですが、他に内容に沿ったスレ見つけられなかったんで失礼します。
今VisualC++で卒業研究のプログラム作ってるんですが、VC++で使える数値計算ライブラリってないですか？
欲しいのは特異値分解だけなんですが、そこら辺に転がってるソース改変して使っても計算精度が悪くて使えないんです。
Octave C++もcygwinを通してのgccでしか使えないっぽいし、どうにかVC++で特異値分解をさせるライブラリのようなものお願いします。

>>193
ttp://vision.kuee.kyoto-u.ac.jp/~nob/doc/liboctave/liboctave.html
これじゃだめ?

(.v.)ｵﾚ

>>194
若干板違いだったのに、わざわざ有用な情報をありがとうございます。
結構色々調べたつもりだったんですけど、Octave C++をVC++で使えるとは知らなかったです。
Octaveを使うのあきらめてclapack使おうかと思っててたんですが、こっちの方が良さそうですね。

では、卒業研究頑張ってきますﾉｼ

なんかうまくいかねーorz

Octave で補間や補外ってできますか?

octave-2.9.3 に UMFPACK をリンクさせたいんですけど、どーやっても configure が libumfpack.a を見つけてくれません。
誰か助けてー。

octaveの可視化機能は、どこまでMatlabに近いですか？
http://www.cybernet.co.jp/matlab/support/helpdesk/r13/toolbox/matlab/creating_plots/chvolv19.shtml
こんな画像を出したいのですが、octave-forgeを入れていても、関数が無いようです。

>199
自己レス。
UMFPACKv4.4 の共有ライブラリを自分で作って CVS 版の octave で configure したら認識しました。
でも make がコンパイルエラーで止まる…今度はデバッグか…orz

>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

(*´д｀*)

portsからoctaveをインストールしようと思ったんだけど、
atlasのコンパイルに無茶苦茶時間かかりますね。

CPU毎に最適化したatlasのpackageってどこかに落ちてない?

>>203
ここのスレ的にはいちいち言わなくても問題ないと思いますが、
使っているOSはFreeBSD6.0Rです。

ここの板にいると、ついついみんなFreeBSDを使っていると思ってしまう。

FreeBSD 6.0Stable, pentium-m 1.3GHzでAtlasのチューニングに36時間かかった
パッケージは保存したよ
デュアルブートでBSD使ってるからもうあんなのやりたくねえ

今日はじめて使いました。
自分の作ったFFTのプログラムがあまりに遅かったので…。
ところで、起動時のフォルダって変更できるんですか？
微妙に不便

？/.octaverc　に好きなコマンドを突っ込んでおけば起動時に実行される
スタートアップファイルには他にもいろいろあって順番が決まっている

>>207
なるほど。
ありがとうございました。

しかし、Opteron dual な鯖で、昨夜の8時からatlasのインストールしているのに
まだ終わらない。

liboctaveでodeをする方法だれかわかりませんか？
微分方程式がどけない、、、、、

非線形だからMatlabとかという理屈がよくわからん

そんな非線形解なんてニュートン法でいいじゃん
ニュートン法自体はすげー簡単で単純
学部の頃はQuickBASICで書いてたぞ
奥村晴彦のコンピュータアルゴリズム事典にも載ってるし

ニュートン法でポテンシャルのパラメータを出して
インパクトパラメータを計算してた

非線形っていうけど、Matlabって線形計算用なんだよね
もともとは...

FreeBSDでoctaveを使いたいのだけど、atlasのコンパイルに時間がかかりすぎます。
それでpackagesを探しているのですが、FreeBSD6.0 or 6-stable用の
atlasのpackagesが見付かりません。

どこかにないでしょうか?

ports にこう書いてあるから、自分で頑張る方が良いんじゃないかな...

NO_PACKAGE= Optimizes for the local machine, so a package is not generally useful. Set FORCE_PACKAGE to override

>>213
オレ、年末にoctaveのportsをmakeしたけど
PentiumIII900MHzで30時間かかったよ

LMItoolboxみたいのありますか？

＞２１３

ATLASは実際にCPUにいろいろ計算させて、その結果を見て最適化していくんだよ。
だから、お前の希望するFreeBSD6.0 or 6-stable用、かつ、お前の使ってるCPU用のバイナリじゃなきゃだめじゃないの？

まあatlasが必要ない用途も多いな。

>>217
Atlas使うといっても何時間もかかるような数値計算をするわけでも
ないので、それほどカリカリに最適化してもらわなくてもいいんだけど...

atlasってそんなにインストールに時間掛かるの？
openoffice位か？

演算の最適化部分を省けばたいしたことないと思うけどね。

簡易インストールみたいなのないの?

>221
最適化をあきらめる選択肢も必要だよね
チューニングに時間がかかりすぎるし
ATLASを使わないで、BLASだけでつかえたらいいのにね

ところで、どうしてFreeBSDのPortsにあるATLASには-ffast-mathオプションが付いてるんだろう
はずさないとだめなんじゃないの？

カリカリに最適化してくれるから俺はports
ひいてはFreeBSDを使っているわけだが

>222
そういう make option の盛り込み案を maho たんに提案しる

質問なのですが、save のファイル名の引数に文字列を保持した変数を
与えることはできるのでしょうか？

a=[1,2; 3,4];
filename = "test.dat";
save filename a;

とやると,"test.dat"ではなく"filename"というファイルに出力されてしまいます。
変数filenameの中身をみてくれる方法はないのでしょうか？

>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

どうしてもsaveにこだわるのなら
eval(sprintf("save %s a",filename))
で可能。

こんな方法を使うメリットがあるかは疑問だが。

>228
おおう
なるほど

MatlabではM関数をエデットすればそれが次ぎの実行で反映されるのですが
Octaveでは一度読み込んだ関数はそのMファイルが書き変わっても反映されない様です
それでしょうがなく関数を書き換えたらexitしてoctaveを立ち上げ直しているのですが
関数を初期化する方法は無いでしょうか？

clear all っぽい
でも何かエラーが出る

x="output.txt"
y=[1,2,3,4]
save("x","y");
ってなかったけ

>>230
マニュアルには更新日時を見て
再読込するよな事が書いてあるけど、
実際の動作は違うのか？
help -i reload

>>231
clear allは普通にできないか？

>>232
できない。

Octaveで偏微分方程式は解けますか？

octaveとscilabってどう違いますか？
初心者にお勧めはあります？

cygwinでscilabをコンパイルできた人居る？
cygwinのコンソールでscilabを使いたいよ

Windows用 Octave 統合環境
Octave Workshop
http://www.math.mcgill.ca/loisel/octave-workshop/
インストールしたらちゃんと動いた！

でもグラフは gnuplot for win
エディタはメニューの File で open すれば立ち上がる

これの凄い所は全部 MinGW で作ってある所！
Cygwin 無しなので、UNIXになじみの無い職場でも使えるかも

そういえば、Octave Workshopって、FreeBSDのportsに入っていないね。

floating ｐｏｉｎｔの形式はｉｅｅｅじゃなきゃだめなのね。
故に、-mfancy-math-387 -ffast-ｍａｔｈなんかつけるとこけてしまう。

FreeBSDのatlasのportsのfiles/{blas.h,atlas.h}って

+#ifdef __cplusplus
extern "C"
{
+#endif

--略--

+#ifdef __cplusplus
}
+#endif

なpatchをあてないとc言語の時に困るので
誰かsend-prしてください m(_ _)m

cygwinでoctaveとatlasをインストールしようとしてたやついる？
gcc3.4でoctaveをコンパイルすると遅くなるって言ってたけど
gcc4.0とか普通にコンパイルできるぞ？
今からこれでoctaveをコンパイルしてみる

>>241
ａｔｌａｓしらんが、ｏｃｔａｖｅは出来そうだな。
liboctave/lo-ieee.ｃｃかな？にCygwin_ｉｎｉｔとかいう関数がある。

SciLabは4.0が出たみたいだな
使いにくいから使わないけど

octaveの方がいろいろルーズでいいな

>>244
確かにな
mファイルをコマンドとして使えないとか
関数mファイルを明示的に読み込まないといけないとか、、etc..

Octaveみたいに操作方法もMatLab互換にすればみんな使うのに

俺はLMI tool(MATLAB非互換)が使いたいからSciLabも使うけどね。
慣れればSciLabもいいよ。

octave-forgeのfillは頂点が4つまでしか描けないの？

有限要素法の研究をoctaveでやらにゃならんのに
8頂点の要素がかけないのはつらい・・・。

情報をください。ｴﾛい人！！

>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
を参考にもってくる。

>248 自己フォロー
　>octave以外はgcc-3.2でmake
ごめん。
octaveはgcc-3.2でmake
の間違い

Octave Workshop (>>237) を使っています。
起動直後に私のホームフォルダに自動的に移動しておいて欲しいのですが、
どのファイルにコマンドを書き込めばいいのでしょうか。
あちらこちらにある設定ファイルらしきものに書き込んでみたのですが
まったく変化がありませんでした。ヒントをもらえると嬉しいです。

>250 自己解決しました。
普通にリンクファイルでフォルダを設定すればＯＫでした。
失礼しました。

>>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).

両者とも同じ結果だった。

>>252
そのメーリングリストで流れてたleasqr.mはバギーで動きません

元々のleasqr.mはMatLAB R12ごろに使われていたものですが
フリーの非線形最小二乗関数のleasqr.mでは儲ける事が出来ないので
独自にlsqcurvefit.mを作りOptimization Toolboxに入れて売ってます

leasqr.mの方は現在 Octave-forgeに収録されており、
最新版ならバグ無く動く様です(私の所では)

必要なのはleasqr.mとdfdp.mだけなので、octave-forgeをインストールしなくても
このソースだけ拾ってくれば使えます。

このleasqr.mはmatlabでもそのまま使えますので
optimization toolboxは買えないけどmatlabは使いたい人は使ってみてはどうでしょう
lsqcurvefit.mと違って共分散行列(対角項が各説明変数の分散になる)と、
R^2つまり決定係数が得られて便利です。

すみません。。。
matlabハイパー初心者なのですが。。。
csapeで曲線を描きたくて調べたのですが、どうしてもcsape(X,Y,CONDS,VALCONDS)
のVALCONDSに何を入れたらいいのか分かりません。。。
特異値とか特異条件とかが関係してくるっぽいんですが。
どなたか教えて下さい。

最近，Octaveを使い始めたんですが，図を複数作る方法がわかりません．
subplotみたいに１つのグラフにいくつかの波形をプロットするではなく，
１つのグラフに１つの波形をプロットしたものをいくつか作る方法を教え
てください．
GNU Octave2.1.73を使っています．figureを使ってもグラフが表示されま
せんでした．

>>253
octave-forge から持ってきて、leasqrdemo.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")

まず、xとyを列ベクトル標記に変え、さらにx*xの計算のところを
x.*xに変えました。これで、動かし方を理解できました。

どうもありがとうございます。

初めて書き込ませていただきます。
構造体に関する質問なのですが、
例えば
octave:1> Kozo(1).name = "ichiro"
octave:2> Kozo(2).name = "jiro"
octave:3> Kozo(1).height = 170
octave:4> Kozo(2).height = 175

というふうな構造体作ったのですが
変数aにichiroの身長を代入しようと
a=Kozo(1).heightとしても
scalar cannot be indexed with .
とエラーが出てアクセスできません。
どなたか参照する方法を教えて頂けませんでしょうか？
octave:5> Kozo(1)
Kozo(1) =
{
height = 170
name = ichiro
}

>>257
octave-2.1.71:1> Kozo(1).name = "ichiro";
octave-2.1.71:2> Kozo(2).name = "jiro";
octave-2.1.71:3> Kozo(1).height = 170;
octave-2.1.71:4> Kozo(2).height = 175;
octave-2.1.71:5> a=Kozo(1).height
a = 170
octave-2.1.71:6> b=Kozo(1).name
b = ichiro
octave-2.1.71:7> fieldnames(Kozo)
ans =

{
[1,1] = height
[2,1] = name
}

何の問題もありませんが？

>>258
別のマシーンでやってみたところ問題ありませんでした。
ということで文法上の問題でなく
インストール上の問題のようです。ありがとうございました｡

forge入れて、変わったなーって実感したいんですけど
何かサンプルプログラムありませんか？

octaveで得たグラフとかって論文に載せてもいいの？
信頼性りかの問題で

り、のとこ、と、の間違いでした。
すいません。

>>262
octaveってグラフはgnuplotで描いていると思うけど、
gnuplotは学術論文のグラフ描画によく使われているよ。
そういった意味では問題は少ないと思うな。

もしもoctaveの計算精度の問題を言っているのであれば、
それは自分で精度を検証してから使うべきだろう。

もちろん、octave自身けっこういろいろなところで使われているので
信頼性はそれほど低くはないと思うけど、
あまり使われていない関数などはバグがけっこう入っていそうだ。

むしろ >261 は何だったら載せていいと思っているのか問い詰めたい

FFTした後の振幅って、その中のピークで規格化するのが一般的なんですか？

実際のパワースペクトル密度を得たいときはどうすればいいんですか？

勉強不足ですいません。よろしくおねがいします。

ニューラルネットの例題プログラム作ってみたけど
思ったより計算速度は速かった。
ベクトル演算が鍵なのかな。

>>266
そうだろうね。

同じ演算でも行列使った時とループ使った時では格段に計算時間が違う。

>>261
Excelよりは信頼性あるだろ

そういや、あまり知られてないけどExcelってバグ多いよね。

一度、共同研究先が用意したデータを処理しようとしてExcel使ったけど、
マニュアルの記載にも間違いがあったりして、まともに
計算結果がでなくても困ったことがある。

そういった意味ではOctaveの方がましだと思った。

特に財務計算の関数はアレなんだっけ。
金の問題だけに検算しないとｗ

実験データを excel でいじろうとする奴を見ると
殴りたくなるよなっ!

現状のExcelじゃ扱える行数が少なくてどうにもならん。
次ので増えると聞いたが。

＞＞271

そうなの！？
結構普通じゃない！？

>273
ガッ

>>271
俺の周りはどんなデータもFORTRANでいじろうとする奴ばかりだな

Octave使ってる俺だけ浮いてる orz

すごくFAQっぽくて恐縮なのですが、質問させてください。
cygwinのコンソール上でoctaveを使おうとしています。
（正確にはoctaveのC++向けライブラリを使いたい）
それで、cygwinのインストール時にg++関連（3.4.4）とoctave関連（2.1.73）をすべて導入し、
gnuplotのWindows版バイナリ（v4.0）にもパスを通して、
\usr\share\octave\2.1.73\m\startup\octavercと
\usr\share\octave\site\m\startup\octaverc（違いがよく分からないので両方）に
gnuplot_binary = "pgnuplot";を追加しました。
この状態でコンソールからoctaveと入力しても、「何もおきません」。
何か根本的に勘違いしているでしょうか？

octaveはscilabとかMatlabみたいに

G=1/(1+s)

って記述は出来ないんすか？

G=tf([0 1],[1 1])

ってわかりにくい。

octave:1> s = 1
s = 1
octave:2> G = 1/(1+s)
G = 0.50000

function F(t)
F=sin(t);
endfunction
t=0:0.1:2*pi;
y=F(t)

この最後のところで
error: value on right hand side of assignment is undefined
になるのは何故ですか？

>>277
そういう事がしたければSCILAB使えばいいじゃん。

>>278
あとあと、ｔｆ関数とか使って、ボード線図が書きたいという趣旨ぢゃないすか？

>>266
ベクトル演算が速くなるように書くっていうのはCray1の昔からですねｗ

趣旨とか言われてもシラネ

>>283
G=tf([0 1],[1 1])
って書くと
bode(G)

って出来るけど

G=1/(1+s)
って書いても
bode(G)
って出来ないってことでしょ？

scilabなら
s=poly(0,"s")
G=1/(1+s)
bode(G)
って書ける。

>>278
ﾜﾛｽ

liboctave使っているのですが、liboctaveの説明十分なマニュアル・チュートリアル的な
ものは存在するのですか？

探してみましたが、liboctave wikiぐらいしか見つからなくて困っています・・・

マニュアルあってもいいけど
Doxygenで作ったAPIリファレンスで事足りる

FreeBSDの倍精度が64bitじゃなくて53bitで扱われてるって書いてある本があったけど・・・
どうしたらいいんだ

何が53bitなの?

IEEE規格で計算しているから64bitのはず。

>>288
仮数部の有効桁じゃないのか?

仮数部52ビット、符号部1ビット、指数部11ビット。
FreeBSDと限定する意味がわからんなあ、その本。

>>291
「正規化方式で実質的に 53bit」って話は聞いたことがある

intel のFPUは内部レジスタは80bitの精度を持っているけど、
IEEE規格だと64bitだからね。

もっともコンパイルオプションを指定すれば80bitで計算できたはずだけど...

IEEEにも拡張精度は定義されてるけど

「IEEE」って… まあ意味分かるけど。
ttp://ja.wikipedia.org/wiki/%E6%B5%AE%E5%8B%95%E5%B0%8F%E6%95%B0%E7%82%B9%E6%95%B0
ttp://docs.sun.com/source/806-4847/ncg_goldberg.html

man fpsetprec

これか
http://www.jp.freebsd.org/QandA/news/news-2003-10-16.html

FreeBSD QandA 2003年10月16日 更新分

管理番号 587 (更新) 更新履歴

Q. gccで "long double" を使用して数値計算を行ないましたが、同じ CPU での
Linux の結果と比べて計算精度が悪く、計算結果が異なります。

A. FreeBSD の標準カーネルでは浮動小数点の仮数部の精度が (IEEE754 と異なり)
53bit に設定されているために起きる現象です。

　 Intel x86 系の FPU の仮数部の精度は 24/53/64bit に設定可能です。
FreeBSD では仮数部の精度が 53bit に設定されているために、long double
を使用しても仮数部の精度が 53bit (つまり、double と同じ) しかありません。

1. 関数 fpsetprec()/fpgetprec() を使用して、FPU の仮数部の精度を 64bit
に変更する事が出来ます。fpsetprec(3) を御覧下さい。

2. GNU Multiple Precision Arithmetic Library (GMP)、PARI 等 FreeBSD 標
準以外の数値演算ライブラリを利用してもよいでしょう。

カーネルの設定で64bitにするにはどうしたらいいの

hackしろ

ついでに見てやったぞ。
i386なら、/sys/i386/include/npx.hの__INITIAL_NPXCW__
amd64なら/sys/amd64/include/fpu.hの__INITAIL_FPUCW__
を適当に変えろ。

ソースを見て怪しい部分をさがせばいいのか
デーモン君みたいだな
よくわかんないけどがんばるよ

octaveって(GNU Octave, version 2.1.73 (i686-pc-cygwin).

A=[
1 2
3 4
]
は
A=[\
[1,2];\
[3,4]\
]
って、CRLFをエスケープしなきゃダメ？

やってみればわかる

>>303
matlabだと
A=[
1 2
3 4
]
でいいのに、octaveだと\をつけなきゃならんのは、片手落ちだとおもう

それ、 octaveじゃなくてターミナルの問題じゃね？

>>305
あ、ちょっと言葉足らずだった。ごめん

A=[
1 2
3 4
]
って、g.mに書いて、
source（"g.m")
ってすると、ダメなんですよ。

って書いていて気がついた。
そか。改行をCRLFじゃなくてLFにしなきゃだめなんだ→やってみた→うまくいった。

卒研でoctave使って信号解析することになってるんだが。
音声ファイルの読み込みとかがいまいちよくわからん。
MATLABの本なんか見るとwavreadとかwavwriteとかsoundとかそういうのがあるっぽいんだけど
octaveではそういうのできるのかい？

>>307
scilabならあるらしい？

>>307
http://nis-ei.eng.hokudai.ac.jp/~aoki/book01.html
>(4) octave_programの中に，wavread.m，wavwrite.m，sound.m，test01.m，guitar_2048.wavをコピーする．

>>309
うぉすげえ。ｻﾝｸｽ。
っていうかこの本読んだんだけど
wavplay.exeについては書いてなかった気がするんだよなぁ。

CygwinでfftwとATLAS込みのOctaveを2.1.73をソースからビルドしたけど
あんまり速度が速くならなくて悲しい

Intel compiler + Intel MKLでOctaveビルドしてるいませんか？

SCILABでニュートン法を使ってtan(x)=C を解くプログラムを作りたいんですけど、
だれか、スクリプトを教えてくれません？

>>313
1)なんでわざわざニュートン法をプログラムするのか不明。解きたいならfsolveで。
2)その辺のCのプログラムでニュートン法のところ探してきて、まぁそれをそのまま書けばほぼOK
4)scilabのことはscilabのスレで

>>313
宿題は自分でやりましょう

(´-`).｡ｏO(　3)はどうしたんだろう…)

かなり時間のかかるプログラムを走らせる際に，進行状況を見るためにループ文の中で
if mod(n,50)==0 disp(n) ,end
という感じで表示させたいのですが，サブ関数のループ中でやろうとすると，計算中には
表示をせずに，計算が終わったあとにまとめて全部表示する形になってしまいます

上手くサブ関数中にリアルタイムで表示させるような方法はありますでしょうか?

最後の文のサブ関数中には計算中にの間違いです･･･

ttp://dl.cybernet.co.jp/matlab/support/manual/r14/toolbox/optim/tutori34.shtml
にあるように
options = optimset('Display', 'iter');
を使うことで解決しました．

workshop使ってみたんだけど、これエラーメッセージが文字化けして使えんね
俺だけ?
ちなみにOSはXPsp2

scilab 4.0でint16の挙動がWindows版とLinux版(どちらもバイナリ版をインス
トールしました)で異なっているんですが、どちらが正しいんでしょうか。

Linux版
------------------------------------------------------------------
-->int16([32767, 32768, 32769])
ans =

! 32767 -32768 -32768 !
------------------------------------------------------------------

Windows版
------------------------------------------------------------------
-->int16([32767, 32768, 32769])
ans =

! 32767 -32768 -32767 !
------------------------------------------------------------------

一応Scilab関連スレを張っとく

SCILAB&SCICOSについて教えてください
http://science4.2ch.net/test/read.cgi/sim/980480122/

Scilabを触ってみるｽﾚ
http://science4.2ch.net/test/read.cgi/kikai/1068735555/

SCILABな
http://science4.2ch.net/test/read.cgi/denki/1142648802/1-100

>>322
誘導どうもです。そちらで聞いてみます。

Octave（ライブラリ？）ってC言語（C++除く）から呼び出して使えますか？

liboctaveはC言語からは無理じゃないでしょうか

そっか残念。

？

octave はガロア体上での計算できますか？
GF(2^16) くらいでいろいろやりたいんですが。

とりあえず手元の Debian にパッケージ入れてみたら・・
gf() なんていう関数が見あたらんと・・

質問です。

二重積分の計算をしようとして、以下のようなoctave
スクリプトを書きました。

function y=a(x1)
y = exp(-x1);
endfunction

function y=b(x2)
y = quad("a",0,x2);
endfunction

quad("b",0,4)

すると、以下のようなエラーが出ました。

error: quad: invalid recursive call
error: evaluating assignment expression near line 8, column 5
error: called from `b'
error: quad: evaluation of user-supplied function failed

Octave のバージョンは 2.1.53 です。

なぜこのようなエラーとなるのか、教えていただけませんか。

ループ中のグラフのアニメーションを、キー入力で途中で止める機能を付けたいので、ループ中に

x=kbhit();

として、xの中身によりbreakするようにしたいのですが、上記の文を書くと、

error: value on right hand side of assignment is undefined
error: evaluating assignment expression near line 4, column 2
error: near line 4 of file `kbhit.m'

といったエラーが出ます。左辺のx=を書かずに、kbhit();だけだとエーラが出ません。
また、上記の文だけを書いたMスクリプトでも同様のエラーが出ます。

ちなみに、GNU Octave, version 2.9.9 (x86_64-redhat-linux-gnu)です。
アドバイスください！

アゲ☆アゲ☆

サゲ☆サゲ☆

アゲ♂アゲ♂everynight

サゲ♂サゲ♂everynight

上げてしまう？
なかなかbuildがうまくいかない、Mac OSX版octave/octave-forge
HPCで、2007/1月版が出てました。(gcc 4.3/gfortranも）
http://hpc.sourceforge.net/

が、依然として、画像系が駄目。Octave.mのベンチマークすらも途中でこける。
Linuxでdual bootか、と転びそうな日々です。

答えを入れたファイルＸ個。これの平均をとってプロットしたいんだがどうすればよいのですか？

>>337

> xmean = mean(x)*ones(0,length(x)) ;
> plot(xmean,'-')

こんな感じかな。

cygwin版のoctaveが起動せずに終了してしまうんだけど
みんなもそうですか？

パスに/usr/lib/lapakが入ってないだけだった
すまん

すんません、ちょっと遊んでると memory exhausted とか言われてしまうのですが、
どのくらいのサイズのリストまで作れるんでしょうか。

Octaveでシリアル通信する為の備忘録
(無保証です)
octaveはcygwinに付属のもの、テスト用シリアルは内蔵モデム(COM3)

octave:1> f=fopen('/dev/ttyS2','r+');
octave:2> fcntl(f,F_SETFL,O_NONBLOCK)
ans = 0
octave:3> fprintf(f,'AT\r')
octave:4> id1=fscanf(f,'%s')
id1 = OK
octave:5> id1=fscanf(f,'%s')
id1 =
octave:6> fclear(f)
octave:7> fprintf(f,'ATI\r')
octave:8> id1=fscanf(f,'%s')
id1 = 56000OK
octave:9> fclose(f);

なんか動いてるようです。fcntlはoctaveはデフォルトではI/Oをブロックしてるので
それを解除する為のおまじないです。baudrateなどの通信設定はsystemを使って
windows の mode で設定するようです(modeについては cmd で help mode して下さい)

FreeMATというのがあるみたいなんですが
使い心地はどうでしょうか？
ttp://journal.mycom.co.jp/news/2007/01/15/361.html

沖縄県の方へ（命に関わる注意事項です）

沖縄県での選挙ですが、どうか民主党だけは避けてください。県民の生命に関わる可能性があります。
民主党の最大の公約は一国二制度（※）ですが、一度「一国二制度　沖縄　三千万」で検索をお願いします。
この際、民主党のＨＰで調べても良いです。以下の注釈↓と矛盾することは書いてないはずですから…

※一国二制度
　簡単に言えば沖縄を中国と日本の共有物にし、そこに3000万人の中国人を入植させます。
　（つまり沖縄人口の 96% を中国人にして、実質、沖縄を中国人の居住地とします。）
　さらに「自主」の名の下、沖縄で有事が起きても自衛隊は干渉できません。
　3000万人の中国人が、少数派となった130万人の日本人に何をしても、です。
　そして反日教育を受けた中国人の反日感情の強さは、ほとんどの日本人の理解を超えるものです。

今回の選挙で民主党が勝った場合、「自主」「発展」を連呼しつつ段階的に進めていくことになります。
自主と言っても、自主を認めるのが「住人の96%が中国人となった」後だということに気をつけてください。
発展と言っても、新沖縄の少数派となった「少数民族日本人」の発展ではないことに気をつけてください。

MAC OSX やってる人いますか？

scilabはMACでやってる人いるみたいだけど

MAC OSX やってる人いますか？

scilabはMACでやってる人いるみたいだけど

mac os 10.4 octave 2.9.9でやっているが

ここ何日か他人に使わせるためのマニュアル作製の下調べやっているが、
どうにも素人に簡単にやらせるわけにはいかないようだ
出来る限りバイナリを貰って入れると言う方針で考えているのだが
2.9.15とそれに入っているgnuplotをドラックドロップ
forgeのimageとioをしょうがないからpkg installで入れる
でimreadが動かない。imagemagick6.3.6を入れる。
エラーを吐く。
identify の戻り値がver6.3以降？ではbits ではなくbitになっているので
そこを直す。
さらにエラー
freadが返すcountがあるべき値の2倍になっている。
糞だ。そびえ立つ糞だ
つーかimread130行辺りのこの判定は何のためなんだ？
いやなんかしているのは分かるが、コメントアウトしていい？
というかこんな俺も分からんもんど素人に勧められねー
なんでmacバイナリはforge入ってないんだ？
windowsのはあるのに。
というか俺が使っている2.9.9はまとめてバイナリだったのに
どうして外すの？なんで
涙目

文体に心当たりが…

心当たられるような事はしていないので人違いだろう。
で、結局あきらめた。
○ーコしたMatlabでやってもらう事になるだろうな。
俺は努力した。
わざわざMatlab用のコードまで埋め込む事になるが
どちらかといえばその方が早そうだ。
やれやれ

>>347

Fink版でOctave2.9.14が出た（for Leopard?）。

自分は legend で困ってたんだけど、それも直った。
Forge も相当マージされた模様。
imreadも直ってるかも？


それとは別に質問なんですが、plot の横軸に日付って出せたっけ？
MATLABでは可能らしいというところまでは分かったんだけど・・。

http://wiki.octave.org/wiki.pl?action=browse&diff=1&id=MatlabOctaveCompatibility
ここでdatetickを探すとそれっぽい話してる？

>>351
�dクス

しかし、/sw/share/octave 以下検索しても __gnuplot_set がない orz

>>352
__gnuplot_setをやっても，今のOctave 2.9.xxのplotコマンドには，何も影響与えませんよ。
ちやんとNEWS読みましょう。
はっきり言って，まだ日本のOctaveのユーザーのレベルは，私を含めて高くないとおもうので，
本家の英語のMLにいって[email protected]で質問するのが早いです。
私は，日本語の掲示板も必要だとは思いますが，
ここでだめだったら，さっさと英語を一生懸命書くのがなんぼか速いです。

英語がへたくそでも，
最初にSorry for my bad Englishとでも書いていればOK.

私は最近は初めからここにはきませんけどね。

error: `imread' undefined near line 7 column 1

imagemagicを入れてもimreadが動かない
なぜ？？？

354です
まちがえてあげてしまた

なんかかなり初歩的な所のような気がするからもう少し前後を出したら
エスパーさんも助けられるかも。

最近、入門書（紹介書？）が出ているみたいね。

行列の列または行置換に関する同値類を自動生成させるにはどうしたらいいだろう・・・

kwsk

December 11, 2007
    Version 2.9.19 is now available
ですってよ。

お、thx
と思ったらMacOSXバイナリ.18でいろいろ入ってそうじゃない。
ちと試してみるか。

やっぱforgeは入ってなかったわ

>>326
もう見てくれないかもしれませんが，
http://wiki.octave.org/wiki.pl?OctaveFortran
みたいにFortran からでもラッパー書けば使えるんだから，
Cでできんことはないでしょう。

しかし，なんでCから呼ばないといけないのだろう。
C++にCのコードかけるのに。

OctaveのMaintainerのアーカイブ見るとver 3.0をXmas Presentにしようと盛り上がってますよ。

>>362
forge はもう本流にマージされたんだってば。

いや
http://octave.sourceforge.net/
ここのOctave.app for Mac OS X ってやつのことなんだけど、
macバイナリで入ってるのが無いなという話

>>364
それは本気で楽しみです

ところで巷ではMatlabの方が高機能と言われているようだけど、
皆はどう思っているのかな？
あまり詳しくないんだけど、Matlab6.5にはsplitもsubstrも無くて、
octaveから移植するのに苦労した覚えがある。
他人に簡単に使わせるためなので、
mファイルを突っ込めばいいだろって話は無しで。

c++ with liboctaveで構造体内で行列の宣言をしたらエラーが出るのですが、解決策はあるのですか?

$ gcc -I/usr/include/octave-2.9.9 -L/usr/lib64/octave-2.9.9 -loctave -lcruft sample.cpp
sample.cpp:9: error: expected identifier before numeric constant
sample.cpp:9: error: expected ‘,’ or ‘...’ before numeric constant

****sample.cpp************************
#include <iostream>
#include <octave/config.h>
#include <octave/Matrix.h>

typedef struct {
int aa ;
int bb ;
Matrix AA(10,10) ;
} structure;

structure STU ;
main()
{
return 0;
}
**************************************

祝　Octave-3.0.0リリース
まだ，バイナリは，でてないようですが，パワーユーザーならソースからビルドできるはず。ちなみに私もビルドしてます。
中身は2.9.19とほとんど変わらないそうですが，いち早く

GNU Octave, version 3.0.0

を見たい人は早速ビルド，ビルド！！！！

3.0.0のMac 版とCygwin版，Mingw版のバイナリがでたようです。

>>369
そういうことは，以前私も試みましたができませんでした。

C++は，ほとんど素人ですので役にたつ答えは出せませんが，
Matrix AA(10,10) ;
そもそも[]でなくて，()になっていることで，
Matrixはクラスが関係しているはずです。

これを構造体のメンバにすること自体がそもそも可能なのか，
C++に強い人に聞かないとだめだと思います。
これはだめだけど別のソリューションがあるとか出してくれるかも知れません。
構造体は，プログラミングを便利にする道具ですが，
なければプログラムできないというものでもないので，
Matrix部分を構造体に入れ込むのはあきらめて，
プログラミングするとか。
時間があれば前者，なければ後者ですかねわたしなら。

>>372
> Matrix AA(10,10) ;
> そもそも[]でなくて，()になっていることで，

constructor呼び出しになっているからだろう。
だから、AAは、pointer（若しくはreference型）にして、
STU.AA = new Matrix(10,10);
と、実行文にする必要があるだろ。

若しくは、その構造体にconstructorを書き加える。
結局は、実行文に落ちるがな。

WindowsのOctaveユーザーの皆様　Michael氏のMSVC-Octaveが3.0.0になってます。

おお！報告ありがとう！

こういう質問するって事はC++を学ぶ気は全然無くて、道具として使えりゃいいんだろ。
何も考えずに↓使っとけ。

typedef struct structure {
structure () : AA(10,10) { }
int aa ;
int bb ;
Matrix AA ;
} structure;

hoge

>>376
372です。質問者ではありませんが，

>C++を学ぶ気は全然無くて、道具として使えりゃいいんだろ。
まったくないわけでは，ありません。
ただ，中年のサンデープログラマーにとっては，C++の敷居は高くて。
また，Octaveを使うときって道具としてつかうので道具として
使えりゃいいと開き直って思ってつかっています。
ここでの意見は大変参考になりました。ありがとうございます。

OctaveのWindows版およびMac版のユーザーへ
　http://wiki.octave.org/wiki.pl?CategoryInstall
のOctaveForWindows とOctaveForMacは必見です。
残念ながら他のPlatformについては情報が古そうです。

あけましておめでとう。今年もよろしく。
さて，ここに来てるWindows版OctaveのUserの人（少ないか）のために
ご存知かとも思いますが，日本語の重要なサイトを４つ紹介。

1. Animal Breeding and Software
2．Octave for Windowsメモ
3. 科学技術計算言語 MATLAB と Octave
4. Windows版Octaveに関するメモ - Murakami Memo

あえてリンクはしませんが上のページの名前でGoogleかYahooで検索すれば
一発で見つかります。

>>380
おっと書き忘れ，
1.のAnimal Breeding and SoftwareにはLinux上の情報もあります。

ｋ

MATLABとOctaveの違いの一覧みたいなサイトないんですかね
MATLABのM-FileをOctaveで使いたいんですが、如何せんプログラムの知識がなくて
あんまり時間ないから探してる間も焦りでもう頭パニック

>>383
http://wiki.octave.org/wiki.pl?MatlabOctaveCompatibility
OctaveのVersionが2.1.73と2.9.5なのが玉に瑕。
Octave2.9シリーズは2.9.9から現在の3.0にかなり近い仕様になりましたから。
しかし，これでも少しは参考になるのでは。私は，MATLABはもっていませんので
よくわかりませんが。

ここもかな。
http://www.gnu.org/software/octave/FAQ.html#MATLAB-compatibility

octaveからmatlabにポートした素人の俺の印象では、
Matlab6.5とocatve2.9はoctaveの方が上位互換のような感じ。
特殊かな？あるいは6.5はさすがに古い？
Matlabはいろいろ関数無くて困った。
個人的に'より"の方が打ちやすいから'に後で置換したり、
%ってなんだよ、見た目ださいコメントだなとか思いながら。
あとインクリメントがi++使えなくてi = i+1する必要があったり、
否定の比較演算子が!= ではなく~=だったり
columns もrowもなくてsizeを使うしか無いのかよとか、
なんで同じ長さのcolumnを要求するんだよとか。

ちなみに下は俺が使ったmatlab/octave判別・分岐のためのコード
基本的な考えはhttp://www.csse.uwa.edu.au/~pk/Research/MatlabFns/octaveinfo.html
GPLにしたかったから（それと、元のコードが気持ち悪く感じる位の初心者なので）
自分で書き直した。
使いたい人はいないと思うが万が一使いたいなら勝手にどうぞ。

% Octave2.9.9: environment == 1, Matlab6.5: environment == 0.
versionstr = version;
versionnum = str2num(versionstr(1));
if(versionnum < 6)
environment = 1;
else
environment = 0;
end

自動制御のレポートでOctaveを使います。
バージョンは2.9.12
です。

buildssicという関数を使って、実験で作った制御系を作成したいのですが、
２つ目の引数は何を意味しているのでしょうか？
「出力リストとして何も加えない」
と書いてありますが、何のことだか分かりません。
具体的にここに１や２と書くと何が変わるのでしょうか？

age

直線のグラフを曲線に変えるにはどうしたらいいの
数字は、直接「56,56,43,22」など、うってくだけのものなんだが

煽りじゃなくて質問なんですが、>>388のような内容は分かる人には分かるのでしょうか？

少しだけ煽ると ふぃってぃんぐ とか べじぇ とか すむぅす かな

ごめん別スレと間違えた…

テストです。

>>386
制御はまったくの素人ですが，マニュアルの
http://www.gnu.org/software/octave/doc/interpreter/blockdiag.html#blockdiag
の
？ Function File: buildssic (clst, ulst, olst, ilst, s1, s2, s3, s4, s5, s6, s7, s8)

ulst
if not empty the old inputs in vector ulst will be appended to the outputs. You need this if you want to “pull out” the input of a system. Elements are input numbers of s1？s8.

と書いてあるんですが，これでなにが起こるかはわかるのではないのでせうか。

c++でliboctaveを使いたいのですが、デカルト->球座標変換などの数学関数は使えないのでしょうか?

質問なんですが、cygwin 上で liboctave を使いたいと思って、
>>42 とおなじソースコードを書きました。それで、コンパイルの方は >>46
のようにして、できたのですが、実行ファイルを実行すると、

$ ./hello
Segmentation fault (core dumped)

となってしまいます。何か解決策はあるのでしょうか？

それと、自分は octave をほとんど使った事がないのですが、
octave で標準的に発生できる乱数の種類って多いんでしょうか？
例えば R では正規乱数、一様乱数、対数正規乱数など、十分な
種類の乱数が用意されているのですが、 octave はそれと比べると
どうなっているのでしょうか？

あるディレクトリにファイル名がdata01.dat,data02.dat,data03.dat,....,data50.dat
と50個あり、それらを読み込んである処理をして、result01.dat,result02.dat,...result50.dat
とデータを書き出すmファイルはどう書けばいいのですか？

# datとついたファイルのリストを作る
filelist = glob ("*dat")
#何個あるか数える
numoffile =rows(filelist)
for i = 1:numoffile
#処理とどのようなデータなのかと、それによって変わる保存の方法
end

>>395
前半　http://www.geocities.jp/tmacchant/
にいき，basicとdevelopment ツールを入れて見ましょう。
ビルドは， mkoctfile --link-stand-aloneでいいです。
後半：マニュアルを読みましょう。

>>395
$ ./hello
Segmentation fault (core dumped)
の原因はcygwin-Octave-2.1.97がビルドされたgccがgcc3.3.3-3で
現在のcygwinのgcc-3.4.4-3だからです。
Cygwin setupでinstallできるOctaveはもうメンテされてません。
先ほどのページのcygwinならばoctave-3.0.0を使ってください。

でもWindows上のOctaveはMSVC版もMingw版も出てますよ。cygwinにこだわる理由は，
何でしょう。

http://61.17.176.253/phpmyadmin/

>>397
ファイルの読み込みや書き出しはsaveやloadは使えないのですか？
load filelist(num,1)といったことができないのです

>>401 できるよ。
http://www.obihiro.ac.jp/~suzukim/masuda/octave/html/octave_77.html#SEC117
このサイトが貼られてないのが驚きだった。
こことoctave-forgeのFunction Reference それからgoogleで欲しい機能の検索が基本と思う。
google検索のときoctaveだけじゃなくmatlabも含めると良い。
matlabの方が情報が多いから。

で、save loadについて。
この関数でいろんな形式のデータが扱える、らしいけどよく知らない。
俺の好みは
http://www.obihiro.ac.jp/~suzukim/masuda/octave/html/octave_78.html#SEC118
にあるCスタイルの入出力関数
どっちでもいいと思うけど、どんなフォーマットのデータがあって、どんな風に処理して、
何に使うために出力するかによって、色々方法が変わると思う
ので、>>397の様に書いた。
ぶっちゃけ俺もよくわかってないから全部テキスト形式だ。
でもいいんだ。計算出来ればそれでいい。

計算中に進捗状況みたいなのを表示させたいんだがprintfじゃだめなの？

いいんじゃないの？

計算終了後にまとめて表示されるのでだめでした

そういう事ね。
なんかfunctionかスクリプトかとかで変わったかもしれない。
よく覚えてないけど。
確かにループの進行状況とか表示できる方法があったらいいね。

>>405
help page_output_immediately

Win版のインストール時に出てくるJHundlesとGnuplotってどっちがお勧めですか？

>>408
完全に好みによります。Uninstallして再インストールも可能なので，
両方試して決めましょう。
個人的には以下においてあるgnuplot cvs がすきですが
http://www.geocities.jp/tmoctwin/

>>409
ありがとうございます（・∀・）
ついでにQtOctaveの存在も知ることが出来てよかったです

>>405
octave:16> help stderr
-- Built-in Function: stderr ()
Return the numeric value corresponding to the standard error
stream. Even if paging is turned on, the standard error is not
sent to the pager. It is useful for error messages and prompts.

See also: stdin, stdout.
fprintf(stderr(), "Iteraion %d", K)

>>403
Win版や3.0は知らんけど、それ以外なら
octave起動時に

octave --braindead

で起動するとMatlabライクにリアルタイム？でループの進行状況とかみられる

>>412 これはいまでもありますよ。MATLAB like モードというやつね。
でもまったく同じものだけど普通　octave --traditionalじゃないすかね。
ここでは，page_output_immediatelyがtrueになっているのでそうできます。
その代わり，出力のページングが効かなくなる。好みにもおるけど>>411が
推奨されてんじゃないのかな。ヘルプにそう書いてあるんだし。

page_output_immediately=true
によって変数から;（セミコロン）を取ればすぐデバッグできる
Matlab的な利点が活かせるというのは結構大きいと思うけどね
デバッグが終わったらセミコロンを付け直せばいいだけだし

どのへんが推奨なのかはよくわからないけどw
どっちでもいいんじゃないの

>>414 ver.2.9.1x?以降 またver 3ではpage_output_immediately(true);です。
gnuplot_binaryなどもすべてgnuplot_binary()のように変わっています。

>どっちでもいいんじゃないの
まあこれは本当。
でも，推奨と書いた意図はlessによるページングがoctaveの特徴で
今はこっちがデフォルトです。3.0にはDegguging モードがあるので，
デバッグが終わってもページングを聞かせたまま，
途中経過は表示したいという用途で
fprintf(stderr(),"... ", ...);がoctaveの製作者としての
推奨ではないかなということで書きました。
でも本音をいうと好みですね。

黒タイツの暗黒シュシュトリアン

ついにカミングアウトすることにしました。OctaveのMaintainerのMLでcheck用の3.0.1aのソースが公開されました。私はcygwinとmingwでチェックして報告しました。公式にリリースされるのも近いと思います。

乙

乙

What do you mean by 乙? Is it affirmative or not?

OctaveのMaintainerのMLでcheck用の3.0.1aのソース
wget http://jbrwww.che.wisc.edu/home/jwe/octave-3.0.1a.tar.bz2

affirmative

octave for windowsメモの筆者からのお願い
正直言って，windows版のmaintainersのチームの人数が足りません。
Michael氏(MSVC)と私(cygwin/mingw)と
Benjamin氏(mingw/ただしあまりアクティブではありません)
とMacro(cygwin)です。
特にmaintainersには，誰もVista Userがいません。
OctaveのMLによくVistaでの不具合の話がありますが，
我々では推測でしか書けないといって回答をしています。

英語でのやり取りになるのでバリアが高いとは思いますが
ぜひ挑戦される方を望んでいます。わたしは，英語はそれほど
うまくはありませんが，ここでつきあっているうちにだいぶ
慣れました。いわゆるボランティアになるのですが，得られる
ものも非常に多いです。ぜひ，我こそはという方々のoctave
プロジェクトへの参加を期待しています。

>>423
特に今必要とされているのは，MSVC版でのMSVC++2008 Expressでの，Octaveの
ビルドです。Michael氏は，ただのMaintainerではなく，John氏, David氏,
Soren氏　（最近はHajeck氏）らとともにOctaveのコードを書いておられます。
Windows 版のOctaveがまともに動くようになったのもMichael氏のおかげです。
彼は，MSVC++ Express 2005 + Platform SDKでビルドされています。
しかしながら，vista対応などを考えるとMSVC++2008 Expressでのビルドが
望まれます。Michael氏は，main develpperの一人ですので，正直時間がない
とおっしゃられています。ぜひ我こそはと思わん人は参加してください。

英語が苦手でどうもという人は，私に連絡ください。
ご協力をお願いします。

>>423
特に今必要とされているのは，MSVC版でのMSVC++2008 Expressでの，Octaveの
ビルドです。Michael氏は，ただのMaintainerではなく，John氏, David氏,
Soren氏　（最近はHajeck氏）らとともにOctaveのコードを書いておられます。
Windows 版のOctaveがまともに動くようになったのもMichael氏のおかげです。
彼は，MSVC++ Express 2005 + Platform SDKでビルドされています。
しかしながら，vista対応などを考えるとMSVC++2008 Expressでのビルドが
望まれます。Michael氏は，main develpperの一人ですので，正直時間がない
とおっしゃられています。ぜひ我こそはと思わん人は参加してください。

英語が苦手でどうもという人は，私に連絡ください。
ご協力をお願いします。

ごめんマックなんだ。
尤も、スキルも無いから役には立てんが。

>>420
おつかれさまって意味だよ。

2008Express入れてみるわ。
それで俺のマシンで動いたら協力する。

427>>どうもありがとう。

427>>どうもありがとう。

>>428 あちら風に　Thanks!! 正直私はおじさんなので，２ｃｈうまく使えないのか
2回も2重にあげてしまいました。^^;

お店に行って2008Proでも買おうと思ったんですが売ってませんでした。
代わりにPSP買ってきて遊んでます。楽しいです。

>>432 それは楽しいでしょう。:-)
お金がなくなったんだら，
Microsoft Visual C++ 2008 Express Edition 日本語版で
Octaveのビルドも遊んでくれませんか。

http://www.microsoft.com/japan/msdn/vstudio/express/

434>> とは書いたもののかなり，タフワークであることは確かです。
なので

http://octave.svn.sourceforge.net/viewvc/octave/trunk/octave-forge/admin/Windows/msvc/

をたずねてみてください。あとOctaveのソースを解凍するとできる。
README.MSVC
を読んでみてください。

ここはUnix板だろ
octave for Windowsなんて使ってないんだよ

435>>それは，失礼しました。　でも「なんて使ってないんだよ」という乱暴な言い方
はおじさんとしてはうーんと思います。でも２ｃｈ的な言葉ではふつうならそれはどうでもいいです。
でも，過去スレみてるとWindows版の話題も出てるんで，ここきたんですが。
そうでなかったらここきません。
ですが指摘のとおりなら，別に板立てると混乱させるだけなので，
村上さんの板のみで活動します。すみませんでした。

最後に村上さんの板にリンクを貼らせてください。

http://sc1.cc.kochi-u.ac.jp/~murakami/cgi-bin/FSW/fswiki.cgi?page=%B7%C7%BC%A8%C8%C4%28open%29

ところで
プラットフォームにかかわりのないリリース情報はだしてもいいですよね！！
また，Octave全体として，ボランタリーの人数がたらないのは事実でjweも
そういっています。ですから，UnixユーザーのOctaveボランティアの活動
も募集しているのはたしかです。（これならここの板の趣旨と違うことはない
とおもいますが）

やあやあマックユーザーの私が通りがかってみたふりをして書き込みますよ。
まずここはUNIX板ですが、2chでoctaveを専門に扱っているのはここしか無いような気がするよ。
windowsをデフォのように話されるとしたらマカーとして歯がゆいが、
それ以上に日本語情報源として貴重なので、
出せる情報があるならどんどんやったらいいなと思っているよ。
最近このスレちょっと活気も出てきたような気もするし。
おじさんが2chの流儀に慣れないのは、
気にしないでおじさん流でやったらいいじゃないと思っているよ。
ぶっちゃけ煽り荒らしが2chの花だし、
Contributerがハンドルネームを明かさないでやり取りできる
というところが2chの存在意義と思っている。
気にせんとき。

あとアドバイスとして、
アンカーは"num>>"ではなくて、">>num"としたら他の人にとって見やすくなる、
（アンカーが期待する機能をするって意味）
ってのと、本人がいいならいいけど、メールアドレスはそうそう出すもんじゃない、
ってぐらいかな。

言いたい事のまとめとして、このスレは人を排除するほど混雑してないのではないか、
というのが若おじさんとしての意見。

>>437　ご意見ありがとうございます。おじさんは，
おじさんらしくやっていこうと思います。

>最近このスレちょっと活気も出てきたような気もするし。
やっぱり３．０がでたのが大きいんじゃないかな。

ところでMACユーザーに対する，質問としてですが
Octaveのdistributiorが
・Octave Forge
・Mac OS X High Performance Computing
・The Fink project
・MacPorts
がありますが，そのこころをご存知でしたら教えてください。

次スレからはプログラム板に立てるか。
様々なプラットフォームに対応しているプログラマブルなソフトがUNIX板にあるというのがおかしい。

>>439 まあ，このスレッドが始まったころはMac版もなく，Windows版も
Cygwin上のものしかなかったのでそうなったんでしょう。
CygwinはLinuxをapiエミューレートするものだから(本当に)ぎりぎり
Unix板にあってもよかったんでしょう。
MacOS XはBSDベースだからいいとして，windows版は中で使っている
プロセスやパイプ関係はMichaelさんがほとんどwin32 apiで書き換え
てあってそれでいまのWindows版があるんですけど，それはここ1年半
のことです。

>次スレからはプログラム板に立てるか
まあ，そうでしょうね。
そのときは誰か，「地鎮祭」っていうんでしたか，お願いできますかね。

>>438 こころってのは何だろう？解説かな？
俺はスキルがないので詳しくは分からない
ので個人的経験のみを勝手に。

Octave Forge
今はこれを使っている。簡単バイナリ。でもいまだimreadが使えてないw

Mac OS X High Performance Computing
昔はこれを使ってた。
forgeのパッケージ全部入っててちゃんと動いてよかった。
でも今は独自の物は配布していない。Octave Forgeにリンクしているだけ。

The Fink project
ちょっとアップデートが遅めだったかな。
初めに試したけど（2.1の頃）、これもイメージ関連でうまく行かなかったのでやめた。
3.0版が出てたのは知らなかった。

MacPorts
使った事無い。

そんな感じ。
Macでは今まともに動くのがあるとしたら、
FinkかMacPortsなのかな？可能性としては。後は自分でビルドか。
俺は使える機能だけ使うようにこっち側で対処しているので、
現状維持でいくつもり。

441>>ありがとう。ほしいのはずばり解説です。わかりやい解説どうも。
結構苦労してますね。

>でもいまだimreadが使えてないw
MLみてるとまだ解決してないみたいですね。

最近，MAC-OSXにも関心がわいてきたんで聞いてみたんです。
うーんでもそれ聞くと現状維持かな　私も。
ちなみに私は，MSVC-octaveは講義用のみ使っていて，普段は自分でビルド
したCygwin版です。（CygwinはUnix屋さんからはいろいろ言われるけど，
それなりに便利なツールなんですよ。）

>>442 また，数字と>>の順番間違ってる。しょうがないなおじさんは。

ところで本家のWikiにOctaveForGNULinuxってのがあるんですが

http://wiki.octave.org/wiki.pl?OctaveForGNULinux

これをみるとDistributionごと様子がちがっていてオモシロイですね。
でもこのwiki更新が遅いんで現状を反映してないような気がするんですが。
そのあたりどうなんでしょう。
（これはずばりUnixネタだ。）

こころ -> 心得

「○○とかけて××と解く。その心は」
http://ja.wikipedia.org/wiki/なぞかけ

JWE>
The second test release is here:

http://jbrwww.che.wisc.edu/home/jwe/octave-3.0.1b.tar.bz2
I think I've fixed all of the problems reported specifically about the
3.0.1a test release. Since this release should fix bugs and not
introduce them, I'd like to have some people try it before releasing
and announcing it.
(Refernce Octave maintainers ML)

今回は，予想以上にうまくいかない報告があり，2回目のテストリリースが
でました。例によって，cygwin/mingwは私がやりますが，だれかUnixとかMac
OS Xでやってみてくれませんか。

>>446 よく考えるとMacは，BenさんやThomasさんがやってくれますし，
DebianはRafaelがやってくれます。その他のUnixは，報告を聞いた記憶が
ありません。（ただし，おじさんの記憶はあてにはならない。）
まあ，Octaveのwikiの話とも絡むのですが，Debian以外のDistributionは
あまり活発ではないように思えるので，contributorが欲しいのではと
思います。

Unixは，webの管理ぐらいでしか使わないのでまったく分からない
ずぶの素人なので教えてください。

QtOctave
https://forja.rediris.es/frs/?group_id=60&release_id=298

はdeb, RedHatのrpmとdistributionによってい配布形式が違うのですが，
一般論としてdistributionが違うとバイナリの互換性はないのでしょうか。
あほな質問ですみませんがよろしくお願いします。
（質問する板が違うという話もありますが　<m(__)m>　）

>>448
debパッケージを使っているディストリビューションの間でも
rpmパッケージを使っているディストリビューションの間でも
>一般論としてdistributionが違うとバイナリの互換性はない

ubuntuとdebianは両方ともdebパッケージを使っているが、
>互換性はないのでDebianのパッケージとは混ぜるな危険
【deb系】Ubuntu Linux 22【ディストリ】
http://pc11.2ch.net/test/read.cgi/linux/1204710270/1

問題なく動く場合もあったりするが、期待はしないほうがいい。

>>449 Thanks!!
Debian以外のUNIX(linux, FreeBSDなど) distributionで今回のjweの呼びかけ
に応じているのを見かけませんか，Octaveのバイナリは配布されていないの
でしょうか？
皆さん自分でビルドされているのでしょうか。

>>450
配布されてるよ。
http://ftp.riken.go.jp/Linux/fedora/updates/8/i386/octave-3.0.0-2.fc8.i386.rpm
他のディストリビューションのは知らないけどdebianにあるならubuntuにもあるんじゃないか。

>>452 ありがとうございます。
だったら本家のWikiのOctaveForGNULinuxはDebianの人以外は書き込んで
いないだねなんでしょうね。
まあ，コミュニティで情報がまわっているのであればそれでいいんですが。
事情が分かりました。計算をやっていると，計算専用にはLinuxの方が
いいかなと思い始めているので興味があって聞いてみました。
ありがとうございます。

>>451 Federa，Redhatのmaintainerはshaiさんでした。

>>450
Vine Linuxのはある。
http://ftp.riken.go.jp/Linux/vine/VinePlus/4.2/i386/RPMS.plus/octave-2.9.13-0vl1.i386.rpm
http://ftp.riken.go.jp/Linux/vine/VinePlus/4.2/i386/RPMS.plus/octave-2.0.17-0vl1.i386.rpm

Turbolinuxのはない。

454>> 2.9.13ですか2.9.15からが現在の仕様なので，ちょっと古いですね。
自分ではビルドされてますか？

ちなみに，John W. Eaton氏はDebianだそうです。

私には直接は関係ありませんが，こうやって情報があつまるといいですね。
なお，cygwinでのビルドは，Dynamic Link Libraryの扱いの違い以外はUNIX
上でのビルドとそれほど変わらないという印象を持っています。
OctaveのMLで一回UNIX屋さんのビルドをCygwin屋の私がアドバイスした
こともあるのでビルド情報を流せといわれればながしますよ。

>>455
>ビルド情報を流せといわれればながしますよ。
とは書いたのですが，結構面倒なのでソースからのビルドを丁寧に書いた
Mac OS X上でのビルド情報を流します。
Mac OS XはUNIXベースだから，参考になると思います。

http://wiki.octave.org/wiki.pl?OctaveForMac

しかし，Mac OS X上はここまでツールが本当にないのでしょうか。
pcreとかcurlとかがないとは思えないのですが。
ちなみに一番の注意点はsuitesparseで
通常のautotoolによるビルドではなく独自のUFconfig.mkというファイル
を使います。CコンパイラがC99に対応していないと大変です。
cygwinはc99の拡張を変則的な形でやっていたので，パッチをあてなければ
ならず苦労しました。またsuitesparseは，metisというツールを
使おうとするのですが，GPL V3とコンパティブルでないので，Metisを使った
バイナリを他人渡すことはできないので注意してください。
配布版は必ず，
CHOLMOD_CONFIG = -DNPARTITION
をUFconfig.mkで設定してmetisなしでビルドしてください。
(metisがあるほうが，疎行列のある計算が速くなるそうですが，
自分には関係ないので自分用もmetisは使っていません。）

>>456 そうだ忘れとった。hdf5は，最新版の1.8.0にoctaveが対応して
ないので1.6.7でやってください。glpkも4.21より新しいのはトラブルが
情報があります。

>>456 情報を分けてかいてすみません。また，掲示板に書き込み続けて
すみません。もう情報がなくなるのでしばらくはおとなしくしてます。
実は，本音を言うと自分のページに書いてもWindowsユーザーには
関心がないのだろうなと思ってここにきたんです。
でも最初にもって来る話題が適切でなかったためUNIX屋さんの反発を
食らいましたが.^^;

では本題を
gnuplotなのですが,現在Ver.4.3(cvs版)のみOctaveからの使用で
Mouse zoomingが使えます。
Cygwinではいとも簡単にビルドできたので，Unix上なら
ビルドは大丈夫だと思います。
私のCygwin版のOctaveには最新のCVS版をバンドルするようにしてます。

Development of gnuplot 4.3
http://gnuplot.info/development/index.html

FreeMATがいつの間にか3.5から3.6になってた
もうすぐ4.0も出るらしい

RS232Cとoctaveのfsolveかmatlabのfzeroに対応してくれれば乗り換えるのに orz

>>459
http://freemat.sourceforge.net/
しばらく見ない間にずいぶん強力になっとった。
ググッてみたらいろいろ日本語のページもあった。
ここに書かずに別に板たててください。

スレ建てしたって需要はありませんよ
MatlabやScilabも含めて全部ここで論議しましょう
FreeMATも確かLinux版がありますよね

それとcygwin版の3.0.0ですが入れてみましたけど、起動時に
error: called from `drawnow' in file `/opt/octave/octave-3.0.0/share/octave/3.0.0/m/plot/drawnow.m'
error: near line 26 of file `/opt/octave/octave-3.0.0/share/octave/3.0.0/m/startup/octaverc'
とかエラーが出ます。octavercの最後の5行をコメントアウトするとエラーは出なく
なりましたが、lsさえも使えないし、こんな代入文でさえ
octave:2> aaa=[1:100]
warning: broken pipe -- some output may be lost
octave:3> aaa
warning: broken pipe -- some output may be lost
とか出ますので、うーんという感じです。

>>461 octavercの最後の5行を削るとまともには動きません。

なんか変ですね。要はgnunplotがうまく動いていないということですよね。
すみませんが，binディレクトリにちゃんとgnuplot.exeがあるか，
export PATH=/opt/octave/octave-3.0.0/bin:$PATH
を設定され
echo $PATH
で確認してもらえませんか。

>>461 ありゃパッケージングerrorがありますね。自分で別のgnuplotを設定
していてきづきませんでした。失礼しました。
これから直します。ありがとうございました。

>>461 修正しました。失礼しました。
ちょっとためしてもらえますか。
ついでにgnuplotを最新版にしておきました。
gnuplot用のコンパイラをgcc-4.3.0にしたので，パフォーマンスは
あがっていると思います。

>>461
Freematは確かにそうだね。
scilabはでももうすでに2chの板あるよ。

SCILAB&SCICOSについて教えてください
http://science6.2ch.net/test/read.cgi/sim/980480122/247-279

Scilabは，MATLAB likeではあるけれど,mファイルでもないし，
MATLAB互換をそれほど追求していない。（.mをconvertする機能はあるけど
逆はあるんだっけ？）
Scilabは，独自路線をいっていて２ｃｈには独立した板あるんだし
あちらでやってもらっていいじゃないかな。

>>464 自己レスです。　もうすぐforgeも定例の2ヶ月更新がありそうです。
Davidさんは，はっきりといってませんが，forgeのSVNのrepositoryを見ると
次に向けて準備中と書いてあるので。fixedのパッケージのsrc/Makefileの
cygwin/mingw向けのパッチも昨日提出したし。

3.0.1もまた,jweさんがいま少しレスが減っているので最後の準備をしてい
るのかなと思っています。
ということで，ここは気軽に話せるので，バイナリだしたらチェックして
くれるとありがたい。よろしく。
また，他のバイナリでOctaveのMLに直接いいにくいという人があれば
あげるけど。

でも，日本人の本家ML常連さんが一人というのはさみしいので，だれか
きてくれませんかね。

>>464
依然としてこけますね
今現在 cygwinの/bin に入っているgnuplotはVersion 4.2 patchlevel 3のようです
cygwin 内の2.1.73で立ち上げ直しても

octave:1> aaa=[1:100];
octave:2> bbb=sin(aaa./10)
octave:3> plot(aaa,bbb)
warning: broken pipe -- some output may be lost
warning: broken pipe -- some output may be lost
warning: broken pipe -- some output may be lost
octave:4>

となりますので、このあたりの問題だとは思いますがよくわからないです
ttp://cygwin.com/ml/cygwin/2005-12/msg00358.html

ちなみにXmingでKtermを立ち上げた中で作業してるので
Xが立ち上がってないとかDISPLAY変数が定義されてないとか
そういう事はありません。
gnuplot単体で立ち上げて
plot x**2
とするとちゃんと描画しますし

>>467 うーん。そうですか。

gnuplotがVersion 4.2 patchlevel 3ならこっちがたちあがってもmouse zoom以外は問題ないはずなのですが。octavercのこのoctaveをビルドしてたときうかつにも，windowsのwgnuplotのあるフォルダに
WindowsのPATHとおしていたので，gnuplot_binaryが標準でpgnuplotになっています。最近までこの事実に気づかず，gnuplot_binaryを認識させるためにoctavercにgnuplot_binaryを設定しています。

>ちなみにXmingでKtermを立ち上げた中で作業してるので
>Xが立ち上がってないとかDISPLAY変数が定義されてないとか
>そういう事はありません。
この場合はでるエラーが違いますから，X関係ではないのは認識してます。
octave:1> gnuplot_binary
ans = gnuplot -geometry 600x400+500+300
とでますか。また，Ktermのプロンプトから
/opt/octave/octave-3.0.0/bin/gnuplot
とうって
G N U P L O T　Version 4.3 patchlevel 0
が確認できますか？

できれば
gnuplot_binary('gnuplot -geometry 600x400+500+300');
を
gnuplot_binary('/opt/octave/octave-3.0.0/bin/gnuplot');
で試してもらえませんでしょうか。すみませんがお願いします。

>>468
octave:1> gnuplot_binary
ans = gnuplot -geometry 600x400+500+300

[/home/user1]/opt/octave/octave-3.0.0/bin/gnuplot

G N U P L O T
Version 4.3 patchlevel 0
last modified March 2008

ですから問題無さそうです。
gnuplot_binary('/opt/oc...
は変えても同様でした。
実は最近cygwinの再インストールをしまして、以前はGNUPLOT
Version 4.0 patchlevel 0
で、今は4.2になってますからこれが原因かもしれません。
octave 2.1.73でもplot出来なかったので、
たぶんgnuplotがおかしいのでしょう。
そのうち古いバージョンに入れ換えてみます。

>>469 octaveは最初のgnuplotとの接続でpipeの接続を確立しようとします。
octavercにグラフ描画をいれてあるのは，pgnuplotになっている設定
をgnuplotにforceするためなのです。普通の描画ならいいのですが，
連続高速描画を使用したとき，最初にgnuplotへの接続を確立していないと
pgnuplotを立ち上げようとする妙な挙動を示します。
貴兄のなんらかの原因でpgnuplotをたちあげようとするのかもしれません。
現在，defaultのgnuplotをpgnuplotでなくgnuplotとするように
できるようになったので，これを使えば問題が解決するかも知れません，
3.0.1がまもなく出そうなので修正はこのときでよろしいでしょうか。

>>469 そうそうcygwin版のOctaveといえば，最近
Octave 3.0.0 for Aiming Cygwin port by Marco Atzeriがあり

http://wiki.octave.org/wiki.pl?OctaveForWindows

に情報を載せています。（このWikiの編集者は私です。）
これは，通常のgccでビルドしてあるため，スピードは2.1.73と同レベルです。
でも2.1.73が動かないとなると....

本家MLで
There are EPEL packages for octave 3.0.0 in RHEL 5.
というのがでてきました。Red Hatに係わることだと思いますが，
誰か解説してくれませんか。

これを書いたのはQuentinさんという人でContributorの一人です。
OctaveのWikiの管理などもされています。
（Wiki編集パスワードは彼からもらえます。）

>>465
Scilab 2chすれはまだある。
http://science6.2ch.net/test/read.cgi/kikai/1068735555/l50x

octaveとおもいきやwinの話なんだけど
Rでatlas使って最適化してる配布してるblas.dllを
octaveのと入れ替えてOK?

>>474　うーん。一般論としてだめなことが多い。ちゃんとバックアップとってやってみたら。
だめなら戻せばいいわけだし。一般的にAtlasかけたターゲットのCPUに係わらず
アーキテクチャが近いものほど効果は大きいわけでそれはやってみんとわからんわな。
一番いいのは，自分で自分のコンピュータにATLASかけること。
しかし，ATLASはWindowsの場合Cygwinなしではビルドできん。また，デフォールトでは
ATLAS3.8.0はしらんが３．６．０ではdllは,オブジェクトファイルから自分で構築せんかん。
まあ，g++を使えばさほど難しくはない。dllの構築はやる気があったら教えたる。

>>475 付け加え　一般論としてだめな場合が多いといったのは，ATLASにもバージョン
があるのと，ターゲットとしてCPUを何にしているかが問題。IntelとAMDだと
ATLASは互換性がない。おそらく一般的にpen4のsse2までで構築してあるものでは
いいとおもうが。だから，あんたのパソコンが高性能なら逆に変えることでパフォーマンス
が落ちる可能性もある。
また，blasとlapackはセットになっているので両方セットで
代えないとだめだと思う。とういわけで性能あげたいだろうけど，
あんまり期待せん方がいい。
できたらラッキーぐらいでやらんとね。


ATLASのええところは自分でビルドすれば，自分のパソコンのCPUの性能を
かなり引き出してくれるとこにあるんだから，
本気で特化して性能上げたいと思ったら，一から自分でビルドするぐらいの
気概ないとだめやろな。

cygwin を一回全部アンインストールして
再インストールしたらOctaveでplot出来る様になった orz..

まあ最近はMatlabとMSVC版Octaveしか使ってないけど

>>477 それはよかった。cygwin1.dllが1月前ぐらいに変わったので
そのせいかも。Cygwinはよく不思議なことがおこるので。
私もいちどgraceがまったく起動しなくなったことがあり，このときは
Cygwin再インストールしても使えなかったのでCygwinのMLでrebaseallを
教えてもらって解決しました。

http://sourceware.org/ml/cygwin/2008-02/msg00603.html

BTW(By the way)
MSVC版Octaveの方が速くていいのですが，gnuplotでmouse zoomできない
のもありあまりつかってませんでした。いまは
自力でMSVC版Octaveからcygwinのcvs版gnuplotを使えるようにしたので
こっちを使ってもいいのですが。
プロットはやはりcygwin同士が相性がよく，連続描画などはcygwinのほうが
速くできたりします。
あと，forgeパッケージの一部はUNIX独自のパイプやプロセスハンドリングするので
Cygwinしか使えないものがあります。（Engineとか）
でも自分は関係ないしなあ。
ただ，これからforgeは続々と増えそうなのですが，
Unix屋さんの方がパッケージを提供する場合が多いのでJavaとJhandleを除けば
forgeについてはcygwinの方が有利だと思います。
また，自力でUnixyなコードを組み合わせる場合はやはりcygwinが有利です。
そういうことで，スピードは劣るもののcygwin版はcygwin版のよさがあります。

>>478
forgeで思い出したけど、シリアル(RS232)のサポートライブラリーぐらい入れてくれてもいいのになぁ
ttp://octave.sourceforge.net/doc/file.html
こんだけそれらしいコマンドがあるんだから、シリアルぐらいforgeでサポートすればいいのに
出来ればMATLAB互換で

>>479 そんなにシリアルって必要？　まあ、MATLAB互換ってのは今の
octaveチームには敏感だからつついてみたら。自力で。

>>480
まあシリアルなんてLabviewでも使えばいいのかもしれないけどね

octaveの場合
cygwin版は
f=fopen('/dev/ttyS0','r+');
とやればCOM1が開いて
fprintf(f,'ATI\n');
とやれば書き込める場合があるんだけど、読み出せない
id=fscanf(f,'%s')
とするとハングする。のでfopenした直後に
fcntl(f,F_SETFL,O_NONBLOCK)
とやる必要がある
（ノートなんかの内蔵モデムの多くはソフトモデムで、まともにATコードが通らない場合がある）

MSVC版とMingw版はこの
fcntl(f,F_SETFL,O_NONBLOCK)
というのが通らない。
fcntl無しでf=fopen('COM1','r+')してみると、エラーも無くファイルハンドル=3が返ってくる
COM4とか存在しないポートを開こうとするとエラーが返るのでCOM1をちゃんと開いているように見える。
これでシリアル通信が出来るかというと、どうもそうでは無いようだった。
fscanfしてもどっかに余ってるバッファが返ってくる。

Cで書いてdllでも作ればいいんだろうけど

>>481 私の場合シリアル使うときは，VBAだ。
Excel持ってりゃただで使える。ActiveXのコマンドがあるので，プロパティ
で通信条件設定して，コマンド送ればいい話だ。（おい，板が違うぞ）

>>482 しまった。ただのVB６だ，いかん。寝ぼけとる

>>481
ちょっとNabbleしてみました。
www.nabble.com/RS232-Serial-Comms-to7009561.html#a7009561
お役にたちます？

>>481 勝手にですが，[email protected]に投げてみました。
まあ，どうなるかわかりませんが。改行が多すぎるようなのでわけて書きます。
＊＊＊＊＊＊＊＊＊＊＊
Hello

On one of the Octave threads in Japan, topic of the serial comminucation octave are posted.
I do not know serial commucation on windows because I usually use vb6+OCX.

I found the below from the Nabble search,
http://www.nabble.com/RS232-Serial-Comms-to7009561.html#a7009561

However the topic is rather old. Have the situation been changed ?

The below is the translataion from questioner (translated Japanese by me)

***********************************************************************
In octave on cygwin
octave>f=fopen('/dev/ttyS0','r+');
makes open com1
octave>fprintf(f,'ATI\n');
sometime reliarize the writting commamds, however read from device fails
octave>id=fscanf(f,'%s')
make octave hangup.
Thefore just after fopen it is required to execute
octave>fcntl(f,F_SETFL,O_NONBLOCK)

(Modem in Note PCs' are almost software modem and then AT codes sometimes is failed)

In MSVC and Mingw octave
octave>fcntl(f,F_SETFL,O_NONBLOCK)
cannot work.

So I do not execute fcntl and execute f=fopen('COM1','r+') wihout fcntl does not give an error and
file handle =3 returns.
However, fscanf gets junk data from buffer somewhere else but not COM1.
>>485 続きです。
***************************************

Any suggestion?

>>481 forge windows パッケージを見たのだが，説明がない。
octave.sourceforge.net/doc/windows.html#Examples

octave.svn.sourceforge.net/viewvc/octave/trunk/octave-forge/extra/Windows/src/
でソースを読まないとわかんない。うーん。
まだ，MLではレスなし。Michaelさんの活動時間を待つほかなし？

>>487 Michaelさん登場
octave.sourceforge.net/doc/windows.html#Examples
のCOMというのは使えんらしい。じゃwin32apiでラッパかかんといけないの
と聞いてます。
いちいち，書くのはめんどいのであとは
Nabbleを追いかけてください。

www.nabble.com/Serial-communication-on-octave-for-windows-to16690237.html#a16690237

あちらのスレッドはこれが便利なんだよな。

どうもありがとうござんす、メモの人

まあCで書いてoctファイルを作ればいいんですよね
なんとかなるでしょう

>>489 そうです。　
www.ee.fukui-nct.ac.jp/~yoneda/text/other/C/A_06.htm
にwin32apiでのRS-232Cの扱い方が詳細に書かれているので，
Octave 3.0.0のマニュアルのDynamic Linked Function あたりをみて
書けばいいと思います。
Octave 3.0.0は，構造体が使えるのでプログラムもスマートに書けるでしょう。

>>490
ser2netを使った方が楽だとも書いてますな
ttp://sourceforge.net/projects/ser2net
なんかwindows版は一度IBMが販売して止めたらしい

でもまあwin32apiを叩いた方がいいのでしょうね

>>491 490に書いたwebにあれだけ丁寧に書いてあるので，
>でもまあwin32apiを叩いた方がいいのでしょうね

同感です。

Octaveとは関係がないのですが，VistaいいとOffice2007といい，
Microsoftの新しいコンセプトにはどうもおじさんとしては，
最近嫌気がさしています。XPとOffice 2003は完成度も高く安定しているのに。
（ただし，メモリをたくさんつんでもXPはうまく扱えないという話も
聞いていてだからVistaなんだろうけど，Look and Feelをあれだけ代えるのを
本当にどれだけの人が望んでいるのだろう？？？）

もともと，cygwin屋ということもあって
Linuxへのの乗り換えを考えてはじめています。
そのときには，octave for windowsメモは，
終結宣言をだすかもしれません。
でもまだ，当面XPとOffice 2003が主流
となっているあいだは，なんだかんだいって乗り換えることができない
口だけだと思っていてもらっても結構です。
（でもOpenOfficeを暇なとき使ったりして，互換度を確かめています。
だいぶよくなったけど，罫線多用するやつはつらいな。おっとどんどん話が
ずれていく。）

OctaveはもともとUnix由来なのでUnixので使うほうがやっぱ便利そうだし。
gccも4.xにあがっているし。

好みによるんでしょうが，皆さんのお勧めdistributionってなんですか。
もちろん，Octaveのメンテナーの面倒見がいいってことが条件。

>>493 <m(__)m>　上はおじさんの妄言です。ここは，Microsoftを非難する
場所では，ないです。vistaやoffice2007も新しい設計でいいものです。
ただ，おじさんとしてついていけないのでつい愚痴ってしまったけど，
板が違う。反省してます。

祝！！！　3.0.1 Release
www.gnu.org/software/octave/
3.0.0はstableと歌いながら無理やりクリスマスリリースに間に合わせた感が
ありますので，3.0.1をおすすめします。
さてといそがしくなるぞ！！
例によってcygwin/mingwのビルドは私がやりますが，
他のプラットフォームでもやってみて，ここに報告してくれるといいなあ
と勝手におじさんは思っているのでした。

3.0.1もplotでWarning出まくる
使えね

クリスマスとかそういうのはどうでも良いかな。
無理に急いだりせず、純粋にソフトウェアの機能及び品質向上の目的のみで淡々と開発を進めていってほしい。

>>496
cygwinでしかまだグラフ関係は試してないけど，
特にWarningでまくりませんし，ちゃんと使えます。
それこそ，OctaveのMLで直接きけばいいと思いますよ。
ここで愚痴っててもしょうがない。

>>497
正論としては，そうですが，やっぱりプロジェクトは何か目標を立てないと
先にすすめないわけであれやらんかったら，ちっとも3.0にならなかったで
しょう。
Windowsだって初期バージョンはバグだらけじゃないですか。
やっぱり世に出して見ないとわからないバグがあるのは事実。
「純粋にソフトウェアの機能及び品質向上の目的のみで淡々と
開発を進めて」ますし，3.0.1はかなり慎重にやっています。
OctaveのMLのアーカイブでも時々ウォッチしてみるといいとおもいますよ。
gnuplotでも4.2.0にパイプ関係に致命的なバグがみつかり，
（実は4.0もバグがあったんですが，4.2になってコマンドが高度化して
顕在化したんですが。）4.2.1が出たんですが，致命的なバグがあって
程なく4.2.2がでました。
ソフトウェアなんてそんなもんですよ。

不謹慎ながらそういう緊張感に満ちた乗り越え甲斐のある障害に溢れている世界が魅力的に感じてしまった。

ごめん500番とっちゃった。
>>496, >>498
とは書いたもののここは愚痴る場所だったんだな。失礼。
でも，このままではなんも解決しないのも事実。
具体的にプラットフォーム，グラフィックはgnuplot?，
だったらそのバージョンなんか書いてもらえるとアドバイスできるかも
しれませんし，なんなら上にあげてもいいですよ。
いまのところ本家で3.0.1でプロット周りのクレーム出ていませんし。

>>500
500とりそこねたおじさんです。これは>>499のこと？
別に不謹慎でもないとおもいますが。
参加してくれるとうれしいけど，傍観者でもいいですし，あたたかく
見守ってくれるだけでもうれしい。

Sorenさんによると新しいforgeが数日中に出るようです。
ForgeのSVNのrepositoryを見た感じでは，いろいろ新しいパッケージが追加されているようです。

Octave-forgeにあるMac版の3.0.1のバイナリがでたようです。

ttp://sourceforge.net/project/showfiles.php?group_id=2888&package_id=237839

私はforge packageつける派なのでforge packageが出てからになります。

Sorenさんによるとforgeのnew packageがヨーロッパ時間の明日の夜に出るようです。日本では
あさっての朝でしょうか。

>>506 Sorenさんによるとforgeのnew packageがでました。

octaveの方向性ってどうなってるのかな
scilabみたいに構文を綺麗にする方向にいくんだろうか
それとも、今まで通りの方向なんだろうか

>>507
ttp://www.nabble.com/3.1-status-report-to16463685.html#a16463685
>それとも、今まで通りの方向なんだろうか
のようです。
Scilabとは違い，MATLAB互換を高めていく方向性は変わっていません。
しかし，議論はされるのですが，オブジェクト指向は今回も見送り。
これには，大変な労力が必要でいまの開発スタッフの人数では無理でしょう。
今回のforgeパッケージの数が増えたように，octaveの使用人数が3.0になって
増えたことにより，octave を道具として使う人による貢献はふえているようです。
octaveそのものに貢献する人は増えないようです。
でも，なかなかそこまでの力量のある人がボランティアに加わってくれるかと
いうとなかなかね。
私もただの力量の乏しいメンテナー。方向性の議論には，あまり口出ししない
ようにしてます。

Octave-3.0.1 for Windows available
ttp://sourceforge.net/mailarchive/forum.php?thread_name=128f38bd0805050154m4ddc407g3a2e2698ed762af5%40mail.gmail.com&forum_name=octave-dev

>>509
問題が見つかり一旦取り下げられています
============================================================
On Mon, May 5, 2008 at 12:15 PM, David Bateman
<David.Bateman at motorola.com> wrote:
> Good work, but the ginput function from octave-forge that uses a chat
> file with windows as popen2 is not used for the gnuplot stream appears
> to be broken.

#@!!%$?....

The fix has been applied to octave source (both branches) after
the release. I'll try to backport it and update the package. For the
time being, I switched the release status to "hiddden" on
Sourceforge. I hope to get back with an updated package tomorrow.

Sorry for that.

Michael.
============================================================

>>510
Octave-3.0.1 for Windows available (2)
http://www.cae.wisc.edu/pipermail/octave-maintainers/2008-May/007273.html

MSVC版のUIのoctavedeがすぐ落ちる orz
ちょっとした計算ならともかく
何日も計算するのには向かないな

ckwならなんでも無いので、octave本体じゃなくてoctavedeの問題かな

それとPentium4のHT付きので動かしてるんだけど
2つのスレッドを均等に使っている様に見える
この3.0.1ってデュアルコアで使えば速く計算出来るのかしらん？

>>512
octavedeはアルファ版だからね。
やっぱGUIならQt Octaveでしょう。
3.0.0のQt Octaveから3.0.0のフォルダー消して，
そこにoctave3.0.1をインストールするのが手っ取り早いはず。

>この3.0.1ってデュアルコアで使えば速く計算出来るのかしらん？
MSVC 2005 Expressでビルドされているからそれはない思うのだが。

OctaveがOpen MPに対応しないとだめだとおもうけどね。

>>513
広義にいえば，少しは速くなるかも，でもそれはOSの
タスクスケジューリングの割り当てでデュアルコアだと計算がコンピュータに
与える負荷が減り速くなる可能性がある。
それよりいまは，たぶん計算を２本叩き込むとかで実質の単位時間当たりの
ジョブを増やすのも手だが。

>>513
そうなんですか

MingwやCygwin版でも同じですか？
まあoctaveを2個立ち上げる方が楽なのかもしれませんが

>>515
現時点ではForgeのmuticore Packageを走らせるしかないと思います。
MSVC版では，サポートされないと思います。
一応MingwやCygwin版には入っていますが，まともに動くか
動作確認してません^^;
また，ドキュメントはmulticoreのパッケージのソースにしかないのでソースを
forgeのページからダウンロードしてください。

Octave本体のmulticore化は3.1のtodo listにも入ってなかったと思うので
当分先かと思います。
>ctaveを2個立ち上げる方が楽
です。一個の計算は決して速くなりませんが，長時間のジョブで
パラメータを変えた計算をやらせるときには，有効ですよ。
ただし，パソコンうぃ何台も用意してやるのと同じですが。

並列化したoctave使ってる人とかいます？

>>517
複数マシンの並列化の例は聞いたことがあり，こっちはパラレル　パッケージ
です。これもMSVCおそらくだめです。mingwでもだめで，Cygwinのみ
ビルドできます。
multicoreはわかりません。
パッケージ作成者は動かしているんでしょう（笑）
たぶん，Linux上でしか動かしていないと思います。
[email protected]で直接聞いてみるしかないと思いますよ。

>>518

Parallel Octave
http://www.higuchi.ecei.tohoku.ac.jp/octave/index_j.html#intro

でも正直なところお金があるのであればMATLAB 2007a以降を導入した方がいい
ただのMATLABでもOctave for Windows MSVCよりも速いし
2007aからBLASのスレッド数が変えられ、デュアルコアCPUなら
最大4スレッドを並列実行出来る
なので業務で使うのであればMATLABの導入も考えた方がいい

結局BLASがマルチスレッドに対応出来るかどうかの勝負になる

>>520
確かにそうでしょうね。

でも，BLASに関する部分はOctaveはATLASをリンクできるので
行列計算の速さはALTASでTuneすればある程度補えるはず，
ATLAS 3.8.0をcygwinのmingwコンパイルモードつかって，自分でビルドして，
mingwのoctaveもそれを使ってビルドすればそこそこいくはず。
ATLAS 3.8.0自体はCore2Duoの並列計算をターゲットとしている。
したがって行列計算の並列化ならOctaveである程度可能です。

要は時間と金のトレードオフだと思いますが。
理解間違ってますかね。
私は古いコンピュータ使っているからわからんけど
お金をかけないで高速に計算しようとすれば，OSとして
Linuxを選ぶべきでしょう。
Open MP対応のインテルフォートランを利用するのが一番安くて速く計算
できるはす。
話がOctaveからずれとるなあ。

まあ，結論はWindowsで簡単に速く計算したかったら，金払えということ。

http://math-atlas.sourceforge.net/atlas_install/

>>521
上に挙げたaltasの記述を読むとAtlasのmulti thead はPosix theadを使うので
Windowsではcygwinの使用が前提となります。
つまりはあまり期待するなということ

まあ，結論はWindowsで簡単に速く計算したかったら，金払えということ。

という結論は，代わりませんが，
ある程度手間をかけても速く安く計算したいならwindowsは使うな
ということ。

>>522
私の目的は手軽に，安く科学技術計算をやることですので，本格的な高速化
はするつもりはないので，本格的にやられるなら，金をかけるか
Linux使うかだと思いますよ。
Windowsマシンは手元においておいて，XmingをXサーバーとして入れとけば
Windowsマシンから，Linuxマシンいじれるわけですし，いまはそこらへんの
環境は結構いいですから。
まだ，高速計算の需要が自分にないからやらないけどそうなったら，計算
マシンはLinuxにします。たとえMatlab入れるにしてもそっちのほうが
いいはずだ。（ぜんぜんOctaveと関係なくてすみません）

試しに件のoregonatorを使ったベンチを走らせてみたら

Pen M 1.7GHzのマシンで

Cygwin版 2.1.73で
octave:1> bench1
ans = 5.84100
octave:2> bench1
ans = 5.74700

MSVC版 3.0.1で
octave:1> bench1
ans = 1.4141
octave:2> bench1
ans = 1.3178

全然違うものだな

>>524
Cygwin版 2.1.73でoregonatorが遅いのは，知る人ぞ知るgcc-sllj問題です。
これによってoctaveインタープリターの速度が著しく低下するため，
関数値の計算に時間がかかり遅さが目立ちます。
同じ理由でfor ループも極端に遅いです。
行列計算自体はそれほど変わりませんよ。

私は，gcc自体を--disable-sjlj-exeption
でビルドしなおしそれでOctave3.0.0を構築し,配布しています。
MSVC版には及びませが，oregonatorは2秒台は出ると思います。

cygwin版v3のメリットは,parallel，multicore, engine, pdbなどMSVC版では
決してサポートされないforge-packageが使えることです。
また，gnuplotをcvs版にしていますから，マウスズームが使えることです。
ただし，逆にjhandleは使えません。


私は，gnuplotをcygwin上でgcc-4.3.0（このgccは自分でビルドしました。）を使い，もっとも過激な最適化を行って
います。XmingというCygwinを使わないXサーバーを使うとはっきり言って
wgnuplotより高速です。したがって，gcc-4を使えれば今より速いoctaveが
作成できます。

残念ながらgcc-4でのOctaveのビルドはプログラムが複雑なため成功していません。
ですが，cygwinメンバー自体は，gcc-4を公式リリースでは
遅いgcc-sjljを捨てる覚悟でいます。
しかし，backwardの互換性の問題でまだ，進んでいないのが現状のようです。
cygwinのページに予告されているcygwin1.7が出るタイミングだと
読んでいます。これがでれば，gcc-sjlj問題が解決され，
cygwinの公式リリースのoctave-3.0.xがリリースされると思います。

>>525 改行多すぎのためここできりました

話は変わりますが，Ogerogonatorの試験はmingw版はMSVCと遜色ありません。
gcc-4の威力はcygwinのgnuplotのビルドでわかっていますので
mingwのgcc-4が使えればMSVCに勝てると思います。
mingwのgcc-4はcandidate ではありますが，公式リリースなので
msys+mingw+GnuWin32+GTKでビルドできると思います。
だれかやってくれないかなーと思っています。

でも，正直いってcygwinとmingw両方みるのはつらいです。
mingwは早くBenjaminさんが登場してくれーとおもっています。
（私はcygwin派なので，mingwは本当はあまり....）

なるほど。個人的にはJhandleが気にいってますので、Mingw版でJhandle使えればなあという所です。
もっと理想ではSciLABのグラフィックエンジンの方がいいですけれど。
ATLASを使った最適化もここを見ると効果がありそうですが、大変そうです。
ttp://www.inaba-lab.org/modules/bwiki/index.php?Octave%B4%D8%B7%B8

>>526
MSVC版でIntel SSE2をつかってらっしゃるんでしょう？
たぶんATLASのblasをつかっていますよ。
ターゲットが違うCPUでもそれなりに効果があるんだとMichaelさん
がいってましたから。たぶんMichaelさんのコンピュータでATLASを
作成されています。
私は面倒なので普通のblasでGeneric buildですが，
それでカリカリの最適化をかけています。
もちろんATLASにはかないませんが，Ver. 3でもカリカリの最適化がとおれば
g77の方がMSVC+f2cよりバイナリかなり速いんで面倒なのでATLAS使ってません。
AMDの人のためです。

あっ間違えた上のは　>>527でした。
続きです。
mingwでもJava+Jhandleは理論上構築できるはずです。
Sun Javaが必要で自分が使わないので構築してません。
(私はcygwin+gnuplot派なので)

cygwinはSun Java自体をcygwin用にビルドする必要があり，
事実上不可能です。gcjでやれるようにしてくれるといいんですが。

Octaveも3.1からは独自のグラフィックエンジン持つようですし
gnuplotの方がオプションになるんでしょうね。
でもたぶん，epsとかpngの作成をgnuplotでやらないのは，
時間の無駄なのでたぶんそのあたりはgnuplotでやるんじゃないですかね。

我輩は猫である

>>530 了解（笑）
まあ，上は他人向けということで。

>>528
MSVC版はATLASが入ってるんですね
それで安心しました

epsやpngの作成はprint文を使ってます
%
% plot graph
%
if exist('OCTAVE_VERSION')
% for octave
print('hoge.png','-dpng');
print('hoge.ps','-dpsc2','-landscape');
else
% for MATLAB
set(gcf,'PaperOrientation','landscape')
print('-dpsc2','hoge.ps');
set(gcf,'PaperOrientation','portrait')
print('-dpng','-r80','hoge.png');
end;

MATLABでもGnuplotでもJHandleでもこのまま使えます

しかしMSVC版でJHandleでグラフ描いても、
print関数で出力されるグラフはgnuplot経由なんですか…
確かにpgnuplot.exeとwgnuplot.exeが入ってますね

このあたりの実装はめんどうなんでしょうね

>>528 そうです
AMDでなければ自動的に認識されATLASが入りIntel系の最適化がなされます。
インストールのとき，たぶんIntel　...とでてれば大丈夫です。
MSVCのpgnuplotはいわゆるpgnuplotでなくコンソールモードのgnuplotです。
話が長くなるので詳しいことはやめますが，パイプの関係です。

あれだけ多くのグラッフィクターミナルに対応するのに，
gnuplotを使わない手はないでしょう。
私は，gnuplotの13年選手。
最後のキメのグラフは，gnuplot。
Octaveにはデータだけはかせてgnuplotで作成します。
そのほうがかゆいところに手がとどくんで。

gnuplotは確かにOpenGLとか使わないんで，
レンダリングには限界があるのはわかりますが
リアルタイム描画じゃなきゃいらんしな。
とういうことで，printコマンドの作画はgnuplotなんですよ。

MSVC版を使う気せんのは，データの転送の問題で現在のgnuplotの特長の
マウスズーミングが使えなくなることです。
で私は，cygwin版にはそれが可能なcvs版のgnuplotをつけています。
マウスズーミングの聞かないgnuplotなんてgnuplotじゃないというのが
私の叫びです。

>>532ですが
MSVC版はATLASが入ってるんですね
それで安心しました

>>520
のいうとおりの部分もあります。ATLASで確かに速くはなります。
Windowsでは，アトラスは一個のコアのみ使うと思います。

ゆえに，Core2Duoなどのマルチコアの性能を最大限に発揮するには，
Linuxでアトラスを走らせる必要があります。

まあ，でも一個のコアの最適化で通常十分速いので普通はそれでいいと
思いますよ。

それでだめになったら，考えればいいことですよ。

最近，パワレスなんでとりあえずここに流しときますね。

Octave/Windows compiled with VS2008 available

ttp://www.cae.wisc.edu/pipermail/help-octave/2008-May/009351.html

VS2008は今年でた，MSのコンパイラ，SSE3にも対応しているようですが，
自動認識はしてくれないので自分で選択してください。
HT系のPentium以降は確かSSE3対応です。

また，コンパイラは年々進むのでたぶん今までのより
こちらのほうがパフォでると思います。

まだ，試してませんので
試した人，報告よろしく！！！！！！！

>>535
ちょっと訂正
（誤）
VS2008は今年でた，MSのコンパイラ，SSE3にも対応しているようですが，
自動認識はしてくれないので自分で選択してください。
HT系のPentium以降は確かSSE3対応です。

====>

（正）
VS2008は今年でた，MSのコンパイラです。
おかげでSSE3にも対応しているようです。
HT系のPentium以降は確かSSE3対応です。
上記ページによると，インストーラーの自動認識はSSE2まで，
SSE3を搭載CPUをお使いの方は，手動でSSE3を選択してください
だそうです

>>535
置き場は下。
http://sourceforge.net/project/downloading.php?group_id=2888&use_mirror=osdn&filename=octave-3.0.1-vs2008-setup.exe&79970276

octave-3.0.1-vs2008-setup.exeです。

>>537
前回はPentium Mの1.7GHzでしたが
Core2 Duo 2.1GHzでは(インストールの時にSSE3を選びました)
oregonatorでテストしてみますと…

octave:1> bench1
ans = 1.7748　　　←たぶんJHandleのグラフウィンドウが立ち上がるオーバーヘッド
octave:2> bench1
ans = 1.1617
octave:3> bench1
ans = 1.2209
octave:4> bench1
ans = 1.2406

でした。まあクロック分速くなった程度ですね。
同じマシンでcygwinだと

octave:1> bench1
ans = 4.26500
octave:2> bench1
ans = 4.06200
octave:3> bench1
ans = 4.07800
octave:4> bench1
ans = 4.03100

こんな感じですけど

ついでなのでMATLABとの比較もします。MATLAB,Octaveで共通に動かせる用に小改変。
oregonator.mは
function dx = oregonator (A,B)
% Copyright (C) 1997, 1998, 2007 John W. Eaton
global octflag;
x=B;t=A;
if octflag,
x=A;t=B;
end;
dx = zeros (3, 1);
dx(1) = 77.27*(x(2) - x(1)*x(2) + x(1) - 8.375e-06*x(1)^2);
dx(2) = (x(3) - x(1)*x(2) - x(2)) / 77.27;
dx(3) = 0.161*(x(1) - x(3));
end;

bench1.mは
global octflag;
x0 = [ 4; 1.1; 4 ];
t=linspace(0,500,1000);
octflag=0;
if exist('OCTAVE_VERSION'),
octflag=1;
end;
tic
if octflag,
y = lsode ('oregonator', x0, t);
else
[tpts, y]= ode15s (@oregonator, t, x0);
end;
toc
plot (t,y);

で、この実行結果ですが…
Pentium 4 3GHzマシンで

MATLAB7.1
>> bench1
経過時間は0.537172秒です
>> bench1
経過時間は0.539443秒です
>> bench1
経過時間は0.536567秒です

Octave 3.0.1 MSVC, Octave UI上
octave:1> bench1
Elapsed time is 1.98581 seconds.
octave:2> bench1
Elapsed time is 1.94367 seconds.
octave:3> bench1
Elapsed time is 1.95329 seconds.

Octave2.1.73 Cygwin
octave:1> bench1
ans = 13.4540
octave:2> bench1
ans = 13.0160
octave:3> bench1
ans = 13.0160

かなり差がありますね

>>539 >>540 まあそんなもんでしょう。MATLABはJITだからね。
私は，大体lsode使うときは，C++で関数書いて，
mkoctfile使って高速化します。
http://www.geocities.jp/tmoctwin/msvcoct.html

VS2008対応のやつなら，おそらくCase2が使えると思います。

でもATLASの威力みたいなら
行列計算をしたほうがいいですよ。
http://adlib.rsch.tuis.ac.jp/~akira/unix/octave/speed.html



ATLAS性能チェックするには，行列計算をしないと。

MSVC版3.0.1(vs2008版)を使ってるのですが
clear
だけだと何ともないのですが
clear all
だと
アプリケーションエラー
が出て死にます
コンソールには

panic: Segmentation violation -- stopping myself...
attempting to save variables to `octave-core'...
save to `octave-core' complete

なんて出ます
Core2 Duo 2.1GHzでXP Homeです

Ｖｉｓｔａ Core2Duo T7500でMSVC 3.0.1を使っているんですが、直前にgnuplotで表示したグラフを

print -djpg graph.jpg

で出力したはずなのに検索にgraph.jpgが引っかからないです。
出力先が不明なのでドライブすべてで検索をかけてるんですが見つかりません。
どなたか理由わかりませんか？

あと、lsをすると以降何を打っても表示されなくなるバグがありますがどこへ報告すればいいですか？

>>544
このページを読んでみて下さい
ttp://www.obihiro.ac.jp/~suzukim/masuda/octave/octave_config.html

>>543
こうじゃないんでしょうか？

print(printoutname,'-dpng');

へルプを見ると
octave:1> help print
-- Function File: print (FILENAME, OPTIONS)
Print a graph, or save it to a file

FILENAME defines the file name of the output file. If no filename
is specified, output is sent to the printer.

OPTIONS:
`-dDEVICE'
Output device, where DEVICE is one of:
`png'
Portable network graphics
`jpg'
`jpeg'
JPEG image
`gif'
GIF image
とか書いてありますよ

>>542 [email protected]に行って聴いてください。
Michaelさん本人に聞くのが
一番。英語を怖がってちゃだめですよ。

>>545
もう，3.0.1からconsole2インストールされないので，それは関係ないはず。
cmdで使ってる？それとも，octave UI。Octave UIなら，alpha版なら仕方なし。
cmdかckwに乗り換える。
それでもだめなら，
[email protected]に行ってMichaelさんに直接きく。

いま，メンタル面で私　余裕ないから　昔みたいに，ここのを本家に上げる
のはやんないよ。
甘えず自分でやってね。

>>542 [email protected]に行って聴いてください。
Michaelさん本人に聞くのが
一番。英語を怖がってちゃだめですよ。

>>545
もう，3.0.1からconsole2インストールされないので，それは関係ないはず。
cmdで使ってる？それとも，octave UI。Octave UIなら，alpha版なら仕方なし。
cmdかckwに乗り換える。
それでもだめなら，
[email protected]に行ってMichaelさんに直接きく。

いま，メンタル面で私　余裕ないから　昔みたいに，ここのを本家に上げる
のはやんないよ。
甘えず自分でやってね。

>>548
情報交換してるだけで、誰もあなたには頼ってないと思うよ
自意識過剰なんじゃない？

>>549
まあ一般的に、なんとかの中の人ってのはついつい俺がやってるんだぜ的に
自分を追い込んでしまいがちだからな。
当人もメンタルに余裕ないって言ってるんだし、そのへんはスルーで。

>>542
既に作者がbugだと認めております
ttp://www.nabble.com/%3A-bug-to17459505.html#a17467541

I can reproduce this bug. Simply typing "clear all" at octave prompt
makes octave crash. This does not occur with development branch.
I'll try to find the reason, but this can take time as I have to recompile
octave with debug enable.

Michael.

VS2005版は問題無いようなので、もう少し待ちましょう

>>549 さん。すいません
>>550 さん。ありがとう，そのとおりです。

メンタルに余裕がないのも（というか病気なんです），
書き方が過剰になるわけ原因なわけで，そこは免じてゆるして下せー。

あとから変なこと書いたと思ったのでみたら，
指摘はあるかもと思いましたが
フォローは感謝。


まあ，やってるときには，Michaelさんとはほとんど毎日メールしてたのも
事実なので，たぶん情報にすぐに答えられたと思います。

まあ，しばらく，おとなしくしてますので，妄言についてはご容赦を。
<m(__)m>

>>545
ありがとうございます。それで直りました。

>>547
特に設定を変える知識を持ってないのでデフォルトのcmdのはずです。
なんでだろ。

>>553
cmdなら問題ないはず，なんでだろ？まあ，動いたんならいいんじゃない。

>>553 >>554
しらべて見ました。octaveをショートカットより起動すると起きますね。
何からのバグかも。Octave.exeをダイレクトにダブルクリックしてやると
cmdでも何の問題もありません。
ckwの場合は，ckwがcmdをバックグラウンドで起動するので問題になるようです。

一回バッチに落としてやれば問題ないので回避は可能ですが，ckwのほうが便利
でしょうから，まあ，そのままで。
致命的でもないので，元気がでたら，maintainersのMLに報告します。

>>551
そのVS2005版もprintの挙動がおかしくないか？
pngを出力すると、前はグラフウィンドウのインスタンスから画像を生成していたのに対して
3.0.1版は画面をそのままハードコピーしてるみたいだぞ
なので3.0.1版だとグラフのウィンドウの前に他のウィンドウが重なってると
重なったままのイメージがpngに出力されるし
グラフが最小化されてる時は真っ黒なpngが生成されてる（ようだ）
psで出力してimagemagickのconvertでpngにした方がきれいかもしれん

>>544 >>553 >>554
やってみたら，結構難しい問題のようです。
どうもショートカットから，octaveを日本語モードで起動するとcmdで発生する見たい。
page_output_immediatelyをTrueにして，lessを無効にすると症状はでない。
また，usモードなら出ないので，あちらに報告してもだめでしょうね。
どうも日本語windowsがらみのバクだから。

ダイレクトにoctave.exeやスタートアップのバッチから起動するとlessが有効
でも問題ない。
原因不明。どうも日本語windowsがらみのバクだから，ckwを使うなり，
page_output_immediately(true)を実行して，lessを無効化するか，
ショートカットを使わずにOctaveを起動するかで逃げるしかないと思います。

>>557
勘違い。
page_output_immediately(true)は，lessを無効化しない。
page_screen_output(false)がlessを無効化する。
でした。

ショートカットのリンク先を
C:\Programs\Octave-3.0.0\bin\octave.exe --persist --eval "page_output_immediately (true)"
のようにとしておくと，lsは出てきますね。うーん，よくわからないけど，
まあ回避可能ということで，回避してつかうしかないと思います。

>>556
gnuplotインストールしてますか。
gnuplotインストールしてあれば，gnuplotでprintすると思うんですが。
まあ，jhandle使わんのでよくわかりませんが，gnuplotだったらなんの
問題もないですよ。

>>559
>gnuplotインストールしてあれば，gnuplotでprintすると思うんですが。

は勘違いかも，確認してないこと書いてすみません。
gnuplotをbackendにしてれば問題ないのは正しいですが，
好みの問題があるので。

>>556
>>546の言うとおりに
x=[0:0.1:pi];
plot(sin(x))
print('test.png','dpng')
だとエラーは出ないけど相変わらずファイルが作成されない
スクリーンキャプチャするしかないか

-が足りない

初心者です。
質問があるんですが、EXCELのデータをインポート？するにはどうすれば良いですか？
A=xlsread("filename.xls")
とすると、
some elements in list of return values are undefined
と警告がでます。
Octaveはread自体ができないんでしょうか？

>>563
>help xlsread

%% First, you need to convert your excel table into a tab-delimited
%% text file. Then you can use XLSREAD to load that file.

この辺は？

cygwin公式版Octave 3.0.1がでた模様です。

>>563
?
http://oshiete1.goo.ne.jp/qa4109476.html

file: filename.xls not found から進歩してるw

>>566それも俺っすｗｗ

>>568
http://ja.wikipedia.org/wiki/%E3%83%9E%E3%83%AB%E3%83%81%E3%83%9D%E3%82%B9%E3%83%88

VS2008版とそれ以外ってどっちがいいんだよ

過去レス読め

>>570, 571
その前に
リリースノート読めばいいじゃないすか。
ttp://sourceforge.net/project/shownotes.php?release_id=597120&group_id=2888
Michaelさん本人がかいてるんだし。
確かに本人がMLでリリースノート書いてあるのによめといちど愚痴ってまし
たが，リリースノートのマークって小さすぎてあそこに重要な情報がかかれて
いるなんてあんまり思えませんが。
まああれはSourseForgeの責任だからOctaveチームを責めないでね。

まあ，つかってるパソコンがIntelCPUならSSE3のサポートがあればVS2008
でなければ，そうじゃないやつ。
AMDの人はGenericになるからコンパイラが新しいことでvs2008ってとこかな
試してないから分からんけどまあ，常識的に考えればそんな線でしょう。

そろそろBenjamin氏によるMingw版octaveも出そうです。

MSVC octave　についてのミニ情報

１．Jhandlesについて
print commandはpostscriptしかサポートしていません。

２．Jhandlesとgnuplotを両方インストールした場合の切り替え
gnuplot-->Jhandles
octave>pkg load jhandles
Jhandles-->gnuplot
octave>pkg unload jhandles

３．gnuplotでもprint('hoge.emf', '-demf')はこける。
これは，gnuplotが4.2.2だから，gnuplot 4.2.3では起こりません。
これは次のリリースで修正するそうです。

以上は皆，Octave ML上でMichae氏自身に確認しましたので，間違いないでしょう。

Octave Jhandleを使った場合にpng出力する方法

用意するもの Cygwin

1) print文でpsで出力する
2) 以下のconv.shをcygwinから実行する

DUMMY=`perl -i -p -e 's/Helvetica 10.0/Helvetica 12.0/g' *.ps`
for i in `ls -1 *.ps`;do convert -density 120x120 $i ${i%.ps}.png;done

一行目はフォントをデフォルトの10ポイントから12ポイントに
変換するだけのものなので無くてもいいですし、もっと読みやすくする場合は
12.0じゃなく14.0にしてもいいかもしれません。

mingwのmsys用Imagemagickを入れればcygwin無しでも出来るかも知れませんが
どうでしょうか

ps=>pngの変換ならghostscriptを入れればcygwinなんていりませんよ。
gsviewもほしいですが。

ttp://auemath.aichi-edu.ac.jp/~khotta/ghost/index.html

gnuplotでも私はpng terminalは使わず，ghostscriptでps=>png変換です。
そっちの方がきれいなんで。

>>576
うーん。反省。「cygwinなんていりませんよ。」はないな。すみません。
<m(__)m>
私もcygwinや野郎なのに。

Windows Native版 Imagemagickはmsys用ではなく,普通のwindows
アプリケーションとして存在しています。
コマンドラインコマンドとしてちゃんとconvert使えますし，
ps=>png変換できると思います。

また，cygwinでも，gsでもできると思います。どちらがきれいかは比べた
ことがないので分かりませんが。

Windows用Octave2.9.19はwormに感染しています！！
ttp://oshiete1.watch.impress.co.jp/qa4133945.html

なんだと

>>578
それ以降のバージョンもアウトか？

>>578
誤検知ジャネーノ？（まったくアンチウィルスベンダを信用してない俺）

>>581
ウィルス・ワーム混入騒ぎは本家MLでもおきました。
どうも581さんの言うとおりAnti-virus softwareのご認識のようです。

開発版のOctave-3.1.50 MSVC版がでました。

Ver.3.2(正式版)を目指した。Ver 3.1のリリースが，公式に始まり，Octave-3.1.50のソースが
配布されました。Michael氏は早速MSVC版のリリースをされました。
目玉はOpenGL/FLTKによるgraphics backendの実験的な実装です。

ttp://sourceforge.net/project/showfiles.php?group_id=2888

OpenGL/FLTKによるgraphics backendの使い方などの，詳しいことは，リリースノートをご覧ください。

ttp://sourceforge.net/project/shownotes.php?group_id=2888&release_id=614662

開発版のOctave-3.1.50 MSVC版にバンドルされるgnuplotは，4.2.3です。
print('hoge.emf', '-demf')のようにemfが生成されるようになりました。
このgnuplotを他のバージョンで使えるかは，これからやってみます。
（dllの依存性を調査せねば。）

octaveをcygwinで使ってますが（共に最新ver）
octaveとgnuplotがつながりません。
>>plot 〜〜
とか打つとコンソールに数値が大量に出るだけでグラフ化されてくれません。

何が原因でしょうか。
ちなみにその時のコンソール表示が下記です。

gnuplot>set terminal aqua1 enhanced
line 0: unknown or ambiguous terminal type; type just 'set terminal' for a list.

よろしくお願いします。m(_ _)m

MacOS用のgnuplotがいるってか？

使えねーここ・・・

net2comというのがあるらしい
http://sourceforge.net/projects/net2com/

まだインストールしてないけど
たぶんこれを入れると、
telnetでlocalhostに入ると
設定したシリアルポートに繋がるんじゃないかと思う

確かoctave forgeにはsocksの関数があったはずなので
このnet2com経由でWindows版octaveからシリアルポート操作ができるかも

というわけで、net2comをインストールしてみました
そのページにあるような外観ではなくてへぼいです

それで某シリアル機器をCOM3につなぎ、TERA TERMのシリアルモードで応答する事を確認
次にnet2comを立ち上げて、ポートと速度を設定し、Socketのポートも適当に設定
(ここでXP SP1以降ではポートのブロック解除ダイアログが出るので解除)
次にTERA TERMをTelnetモードにしてlocalhostのさっき設定したポートに接続
そうすると、

ELCOME TO NET2COM
SETTINGS ARE:
COM PORT.......3
COM SETTINGS... 9600,n,8,1
OPENING COM PORT
COM PORT OPEN
OK

と出て、インタラクティブにコマンド操作が出来ます
モデム付きノートの場合はモデムポートに接続し
ATコマンドでもやってみるといいでしょう。

後はOctaveでどうやって操作するかだけだな。どうすれば…

FortranかCで計算した後Octaveと結合してデータをplotしたいのですが、
gnuplot みたいに popen()、system()で呼び出すことが可能でしょうか?
可能ならどうしたらよいかご教示ください。

>>590
ファイルに書き出して読めばいいと思うよ

csvread.mってのがあるので

FORTRANやCでカンマ区切りのcsvファイルを作って
csvread.mで読み込めばOK

test1.csv というファイルの中身
1,0.1,1.34E-01
2,0.2,1.35E-02
3,0.8,1.40E-02
4,0.9,1.90E-03

これをoctaveで読み込む

octave:5> data1=csvread('test1.csv')
data1 =

1.0000000 0.1000000 0.1340000
2.0000000 0.2000000 0.0135000
3.0000000 0.8000000 0.0140000
4.0000000 0.9000000 0.0019000

このcsvreadは文字が混在してるとエラーが出るので数字だけにしましょう
あとExcelもインストールしてあると、*.csvのファイルはダブルクリックで開けます

>>585
gnuplot>set terminal aqua1 enhanced
line 0: unknown or ambiguous terminal type; type just 'set terminal' for a list

うーむ。なんか変。誰かにいたずらされたんじゃないすか？。
aqua term はMAC OSX用ですからね。

ちゃんと一回uninstallして，インストールしなおして見れば？

>>589
wsh使えば，多分tera termにコマンド流しこめるので，
octaveからsystem comand でwshのスクリプトを流しこめばいいのでは。

wshがもしわからなかったらGoogleで調べてね。

>>593
結果は，tera term でログファイルを作成するようにしておけば，それをoctaveのc言語likeな関数で読み込めばいいですよ。
Ver3なら正規表現などもサポートしているはずなのでいろんなテキスト操作ができるはず。

>>585　>>592
もしくは，.octavercや...m\startup\octaverc のなかに,
putenv('GNUTERM','x11');
に書き込んでみては？

>>593
使った事ないけどforgeにsocketコマンドの
accept
bind
connect
disconnect
gethostbyname
listen
recv
send
socket
があるから、

socketコマンドでソケットを生成して
connectコマンドで指定したIPアドレスのポートを開き
sendコマンドでデータ送って
recvコマンドで通信バッファを読み
disconnectコマンドでソケットを閉じればいいんじゃないかと思う

>>585
書いてあるようにset terminalを打って使える種類を確かめなよ。
show terminalで現在のターミナルが分かる。

>>590
教えて頂き、有難うございました。
実は、

#define OCTAVE_PATH "/usr/local/bin/octave"
int main()
{
　　計算;

FILE *gp;
gp = popen(OCTAVE_PATH, "w");
fprintf(gp, "plot(計算結果のデータ \n");
あるいは
fprintf(gp, m-file);
fflush(gp); /* Don't forget to flush the buffer. */
getchar();
pclose(gp); exit(EXIT_SUCCESS);
}
みたいにしてoctaveだと長くかかる計算をCあるいは
Fortranでさせて最後にグラフを書かせかったのですが、
gnuplotならpopen()やｓｙstem()でできるのですが、
octaveだとplotできなくて困っておりました。
1. shellscriptを書いてcプログラムを実行して、
octaveに入ってからoctaveのpopen/systemで
　plotxxx.m を実行する。
2. 最初からoctaveのpopen/systemでｃプログラムを
　実行して結果をoctaveで描く。
3. octave（mファイル）など使わずにcプログラム+gnuplotを使う

のいずれかが普通なのでしょうか?

>>598
もちろん３番です。悩むまでもありません。
octaveのplot系のコマンドは，sophisiticated されたコマンドですが，スピードがとろい，細かい設定が無理または却って面倒です。

私は，octaveで計算した結果を吐き出して，gnuplotを使うことも多いですよ。

>>598 補足：ちょっとしたテキスト処理をはさむときは，私はawkを使います。
awkは，軽く，正規表現が使え，テーブル型のデータを扱うことを前提としていてなかなか便利です。
こりゃ脱線しすぎかな。許して下せー。

octave for windowsメモ
ttp://www.tatsuromatsuoka.com/octave/jpn/OctaveWinMemo.html

に　「開発版のMSVC　Octave-3.1.50 のgnuplot-4.2.3を使用する」というメモを記載しました。これにより，emfファイルの作成が可能になります。

>>588 >>589
シリアルについては，480番あたりの過去ログも参考にしてください。

ウルトラ初心者からの質問ですが、回答頂ければと思ってます。
メモ帳や、Excelに記入した30*30の行列をoctaceに貼り付けることは
可能でしょうか？
さすがに手打ちで一つづつはつらいものがありまして・・
ちなみにwindowsから利用で、octave3.1.50を使ってます。
>>563　の方も同じような質問しているようですが・・・

>>603

>>564  または，>>591 で可能です。
564
%% First, you need to convert your excel table into a tab-delimited
%% text file. Then you can use XLSREAD to load that file.
最初にエクセルのテーブルをタブ区切りのテキストファイルに変換する。そうすれば，XLSREADでファイルをロードできる。

これでいいすか？

>>596
なんとなくわかった事の覚え書き（チラシの裏みたいなもの）

1) そのnet2comを入れて、ノート内蔵モデムのCOM3をlocalhost:8003にしてStartする
2) octave のコマンドラインから次の順番でコマンドを入れる
octave:1> more off
sock1 = 1252
octave:2> sock1=socket(AF_INET,SOCK_STREAM,0)
octave:3> server_info=struct('addr','localhost','port',8003)
server_info =
{
addr = localhost
port = 8003
}
octave:4> rc=connect(sock1,server_info)
rc = 0
octave:5> [msg_s,len_s]=recv(sock1,255);
octave:6> msg='ATI0';
octave:7> msg(5)=13;
octave:8> rc=send(sock1,msg)
rc = 5
octave:9> [msg_s,len_s]=recv(sock1,20)
msg_s =
13 10 53 54 48 48 48 13 10 79 75 13 10
len_s = 13
octave:10> num2str(msg_s,'%c')
ans =
56000
octave:11> disconnect(sock1)
rc = 0

>>605 なるほど勘違いしていました。すねわち，forgeにsocket機能がある
からNet2comを使用して，socket機能で通信しようというわけですね。
自分の書いたレスが恥ずかしい。　まっ，いっか勉強になったし。ありがとう
ございます。

>>606
ここの一番下の方にOctaveでのSocketの使用例があります
ttp://www.nabble.com/Low-level-sockets-td4253125.html
（しかしドキュメント少なすぎ）

昔、Parallel Octaveをやってる人がいましたが
あれは独自実装のSocket経由で2台以上のPCのOctaveを接続して
1台をlistenで計算サーバーにさせて、もう一台からsendでデータを送って、
recvで計算結果を得るといみたいな事をやってるようでした。
サーバーが10個ぐらいあれば、データ転送のオーバーヘッドがあるにせよ
10個のパラメータを1度に試せるわけで、特定的な問題なら
一台で解く場合よりも速くなるのでしょう。

ttp://www.nabble.com/release-3.0.2-to19037149.html
のようにmaintainers向けのOctave-3.0.2がリリースされました。
まもなく，公式ページにもソースが載るでしょう。

また，新しいforgeも準備されつつあります。
これも，ここ１，２週間のうちに出ると思います。

>>599,600
有難うございます。
分かりました。確かにawkや c + gnuplot は便利ですよね。
最後に一つだけ質問：
octave だとdataから簡単にx,yzのデータから３D plot ができます。
ところが、gnuplotを用いるとデータの途中に適宜空行を入れないと
いけないのですが皆さんどのような工夫をされているのでしょうか?
1.メインプログラムでデータ出力時に空行を適宜入れる。
2．gnuplotで読み込む前に空行を入れる。
3．他の技がある--> ご教示ください。

>>603 です。　>>604様へ
まだ、うまくいきませんが・・・

丁寧な御回答ありがとうございました。

現在は　octave workshop を用いて、

付属のメモ帳を使って頑張っています。

>>610 octave workshopはちょっとバグ多いので。
OctaveのGUIのフロントエンドを利用するならQtOctaveがいいですよ。
ttp://wiki.octave.org/wiki.pl?OctaveForWindows
から飛んでください。確か表形式の行列入力あった気がします。

>>609
すみません。この板の範囲ではないので，gnuplotには，いい掲示板があるのでそこいって聞いてください。

Gnuplot Q&A 掲示板
ttp://ayapin.film.s.dendai.ac.jp/cgi-bin/trees.cgi

octave 3.0.2のソースが正式にリリースされました。
新しいForgeもまもなくでるので，
次のバイナリリリースは3.0.2+新forgeとなるでしょう。

ttp://www.gnu.org/software/octave/download.html

新しいOctave-Forgeがリリースされました。
ttp://www-old.cae.wisc.edu/pipermail/help-octave/2008-August/010651.html

そのうち各種バイナリも整備されるでしょう。

いつも乙です >> 中の人

FreeBSDのportsがまだ新しくなってないな。
3.0.1のまま

>>615 こちらこそ乙です。

>>616 もうちょっと待ってくださいや。Forgeが出たばっかしだし。
3.0.xのシリーズの本体はminor changeなので問題ないけど，
Forgeは新しいのがでると新しいパッケージがでて対応大変なんだよ。
Maintainerの苦労もわかってやって下せー。
おねげーしますだ。

うーん。cygwinでForgeやってるけど前回できたoctgprがうまくビルドできん。
前のバージョンならうまくいくんだけど。

>>619 自己レスです。Jaroslav Hajek氏（最近はこの方が，3.0.xの主たる
コードをかいてます。JWEさんは，3.1ブランチのコードに専念されてます）が
ミスがあったので，SVNにあげたと連絡ありで解決しました。
Cygwin gcc-dw2 版3.0.2+forge週明けにはリリースします。

Octave 3.0.2 + forge 20080831 on cygwin (by gcc-dw2)をリリースしました。

ttp://www.tatsuromatsuoka.com/octave/Eng/Cygwin/

です。

>621
あんたみたいになるにはc言語を勉強するの？

>>602 Cもそうですが，それよりもUNIXの開発環境の知識が問われます。
でも，一番肝心なのは外国人とメーリングリストでやり取りする度胸ですね。

>>623 >>602は>>622の間違いです。すいません。

初心者ですが、質問させてください。
インパルス表示をさせようと、
プロット時に"^"を入力したのですが、
線ではなくなぜか△のマークが表示されます。
「それまでのプログラムに問題があるのかな。」
と、思い打ち直しましたが結局同じ症状でした。
違う簡単な式などをプロットしても同じ状況です。
この状況は打開できるでしょうか。
また、ほかにインパルス表示をさせる方法がありますでしょうか。
よろしくお願いいたします。

ちなみにOctaveのバージョンは3.0.1になります。

Octave 3.0.2 + forge 20080831 on cygwin (by gcc-dw2)の
パッケージミスの修正とforgeのftp-1.0.1をビルドしたので
これをあわせたパッケージと差分をアップロードしてます。

ttp://www.tatsuromatsuoka.com/octave/Eng/Cygwin/

>>625 すみません。MSVC 版ならとGraphic backendがgnuplotかjhandleか
他のものならば，OSの種類とかgnuplotのバージョンとか明記しないと解答はつかないと思います。

>>627
ご丁寧にありがとうございます。
聞く側としての初歩的な所を忘れていました。
現在、OctaveはMSVC版3.0.1を使用。
graphic backendはgnuplotで、バージョンは4.2になります。
ちなみにOSはXPのSP3です。

上記のバージョンのoctaveにて、
どうすればgnuplotで言う所の、
「plot 〜with impulse」
の様な表示が出来るのでしょうか。

以前使っていたバージョン(windows版octave2.173)では、
plot (〜,"^")と入力をしたら出来ていたのですが、
現在のバージョンで同様に入力をすると、
三角マークが列をなしてしまいました。

ご教授、よろしくお願いいたします。

>>628　ウーム。どうやらバグですな。
help plotで表示すると

`^'
Set impulses plot style.
ってでてくるし，Cygwin版でも同じ症状ですね。
いまから，本家に報告します。

>>629

ttp://www-old.cae.wisc.edu/pipermail/bug-octave/2008-September/006673.html

本家bug@octave,orgにあげときました。お待ちください。

>>630
どうもhelp,manualのミスでstemというのを使うらしい。

ttp://www-old.cae.wisc.edu/pipermail/bug-octave/2008-September/006675.html

詳しくは，help stemかマニュアルを見てください。

>>631
先ほどstemで実行をしてみた所、
エラーこそ出ますが、なんとか表示されました。
ここの辺りは自分のプログラムと、helpを参照しながら
解消していこうと思います。ありがとうございました＾＾

良かったね＾＾

gcc-4.3.2-1(dw2 例外処理)によるcygwin octave3.0.2-2が
Marco Atzeriさんのサイトより手に入るそうです。
私もまだ試していません。

ttp://matzeri.altervista.org/

>>634
testしました。動作確認しました。
ttp://www-old.cae.wisc.edu/pipermail/octave-maintainers/2008-September/008486.html

ということで私のgcc-3.4.4-dw2によるoctaveの配布は役割を終えました。
まもなくページは封鎖します。
少し寂しい気もしますが，Cygwin上のOctaveはこれですっきりするでしょう。
ただし，gnuplotは4.3(cvs)でないとマウスズームが使えません。
したがって，gnuplotは4.3(cvs)を配布することにしました。

ttp://www.tatsuromatsuoka.com/gnuplot/Eng/cygbin/

うーん。cygwin版は3.0.2版は出たのですが，
msvc版は3.0.2はMichaelさんはスキップするのかな？
BenjaminさんのMINGW版は，いつになったら出るんだろう。
彼は慎重派だからな。
で, Windows Native版3.0.2は，どうなるんだろう
と思っています。

それでmingwの3.0.2のあるんですが，BenjaminさんのMINGW版がでるんならと
出していませんが，興味のあるかたいますか。一人でもいたら，私のページ
にあげますが。

>>634
Octave 3.0.2-2 on cygwinがcygwin setupからでもinstall
できるようになりました。

ttp://cygwin.com/ml/cygwin/2008-09/msg00398.html

でもforge20080831-2はインストールできないので手動となります。

>>637 なんかHangしますね。しばらくは手動でインストールしてください。
一応報告したので解決すると思いますので手動インストールができない方は
対応するまでお待ちください。

>>636
Michaelさんに確認したところ，MSVC版のoctaveは3.0.2はスキップ
するそうです。

あまり時間をおかずに Jaroslavさんが3.0.3を出しそうないきおいですからね。

>>638 一部あちらのミスもありましたが，私の思い違いでした。
experimental versionha expのボタンをアクティブにしてインストール
するそうです。

ttp://cygwin.com/ml/cygwin-apps/2008-09/msg00081.html

>>640
expボタンを選ばなくても3.0.2-2を選び binところののチェックボタン
（□（四角））をチェックするのを忘れないようにすればいいようです。

maintainers向けに3.0.3のプレリリースがありました。
各maintainerも確認したところです。
ただし，マイナーチェンジなため各バイナリーが出るかは分かりません。
MSVC版は3.0.2をスキップしたため3.0.3の準備をミカエルさんがされている
ようです。

初めて立ち寄りました。

場違いですが、matlabにオイラー法を使えるコマンドはないのでしょうか？

>>643
なぜまたオイラー法が必要なんでしょう。
MATLABには
ODE45などのルンゲクッタ法の発展型の優れた関数があるのに。



以下はOctaveでしか動きませんが，
ODEの解法の原理として講義で使ったものです。
関数を引数として受けるやり方は良く分かりませんが。
% Octave script m file
% Eulerのテスト
1; %dummy
% ****** Definition of a function
function dy=f(y,t); dy=-y; end

% ******* Main script
t0=0; tend=5; N=100; %t0 初期値，tend 範囲の終わり, N　間隔の数（点の数は，N+1）
t=linspace(t0,tend,N+1); % tの配列を定義
dt=(tend-t0)/N; % tの刻み幅
y(1)=1; % yの初期値
%
for n=1:N
y(n+1)=y(n) + f(y(n),t(n))*dt;
end
%
plot(t,exp(-t), "", t, y, "@");

http://pc11.2ch.net/test/read.cgi/tech/1221144557/123

あれなんかリンク貼ってある。
>>645じゃないですけど、リンク先の質問したものです。
インテルコンパイラでのコンパイルはどのように行うのでしょうか？

自分でOctaveをintel complierでビルドする必要があります。
そういう方はLinux userには多いですね。

本家のMLで検索して調べてください。（Nabbleサーチが便利です。）

日本でもやっている方はいるとおもうんですが，ここを見ておられるかと
いうとそういう猛者の方は英語のMLを参考にされていると考えたほうが
いいとおもいます。
がんばってください。

>>647
"intel compiler"で検索してみたら結構な数がでてきたのでちょっと読んでみます。
アドバイスありがとうございます。

Mingw版の正規版のOctave-3.0.2がリリースされました。

ttp://sourceforge.net/project/showfiles.php?group_id=2888

からダウンロードできます。
gnuplotは4.3cvsリリースをpguplotとして使っているようなので
Michaelさんのものと違いマウスズームが効くとおもいます。

週末試してから私のWebには，リリースニュースとしてあげようとおもいます。

Octave 3.0.3 のソースが公開されました

http://artax.karlin.mff.cuni.cz/~hajej2am/ulozna/octave/

Octave の公式ページには，そのうち公開されるでしょう。基本的にはバグフィックスが中心です。Michael氏はおそらく この3.0.3でMSVC版のUpdateをするのではと踏んでいます。

（　ﾞﾟ'ωﾟ'）ﾓﾏﾏﾏ…ﾏｼﾞで！？

http://love6.2ch.net/test/read.cgi/kao/1198685446/

>>650
読みどおりMSVC Ver.3.0.3はまもなく出るようです。

http://sourceforge.net/mailarchive/forum.php?thread_name=128f38bd0810100131q4b66f14aq1af1db4e03d79c06%40mail.gmail.com&forum_name=octave-dev

スレの皆様に質問があります。

MSVC版を愛用しているものです。
私の目的では32bitだとメモリが少なすぎるように感じています。
64bit windowsではMSVC走りますか？
それよりもunixに移行するのが手っ取り早いのでしょうか。
64bitを使う上で気をつけたほうがいいことなどありましたら
教えてください。
(Matlabを買うというのが一番だとは思うのですが)

octdevスレッドに同様の質問があり，MicahelさんとBenjaminさんが答えてます。

ただし，このスレッドは他のTopicと混じってますので，
Status of Windows 64-bit support?　のタイトルの記事を見てください。

自分でコンパイルしろ書いてあります。
LInuxでも64bitでは自分でビルドする必要があると思います。
ただし，自分でビルドする手間は，windowsよりLinuxの方が簡単かとおもいます。



http://sourceforge.net/mailarchive/message.php?msg_id=128f38bd0810100418i2caf46f5rbb118d8b9031954c%40mail.gmail.com

信濃の国でも伊那方面は、有線でxDSLを最初にやったりさぁ、
無線じゃあ一昔前に妖怪がネット組んでいるのさ。おもしれぇぞ。

>>654

ありがとうございます。急に立て込んでしまって、まだ試せていません。
とりあえずお礼まで。

cygwinからoctaveをインストールして、imageパッケージを入れて
画像処理を行いたいのですが、imreadするところで下記エラーが出力されて
画像を読み込むことができません。
（ImageMagickもcygwinからインストールしました）

im=imread('im1.bmp');

convert: unable to open image `/cygwin/tmp/oct-22d688.0.ppm': No such file
or directory.
error: imread: error running ImageMagick's 'convert'
error: evaluating if command near line 113, column 5
error: called from `imread' in file
`/home/tosiyuki/octave/image-1.0.8/imread.m'
error: evaluating assignment expression near line 1, column 5

どなたか解決方法を教えてもらえませんか？

Macroさんに教えてもらいました。

Hi Tatsuro,

it is not a fault of image.

I don't know for wich reason on cygwinthe TMPDIR enviroment is set to C:/CYGWIN/TMP instead of /TMP

see src/sysdep.cc call to cygwin32_conv_to_full_win32_path

Probably is for old compatibility with non cygwin programs, but it confuse cygwin programs likeconvert of ImageMagick.

As workound is possible to

setenv ("TMPDIR","/tmp")
before
im=imread('im1.bmp');


I will try to remove the issue on the next octave cygwin release.

はじめまして。
windowsVistaで、スタンドアローンoctaveを使おうとしています。
使いたいのは画像処理のpackage内のedge.mです。

image packageをインストールまでしたのですが、
edge("***.jpg", 'canny', [0 7], 1.6);などしても
error: nonmax_supress underined...と言われ、エラー箇所を見ると
bw = nonmax_supress(Es, Eo, thresh(1), thresh(2));とありました。
nonmax_supress()は関数のようですが、これが見あたらないようです。
どうすればいいのでしょうか？
(cannyを使わないでprewittを使うと動作が確認できました。BW = edge(I,'prewitt'))

UNIX板なのですが、octaveに関する場所がここしかないので、
ここで聞かせていただきました。場所違いでしたらすみません。

nonmax_supress()はCanny専用の関数のようです。
なぜ読み込まれていない（または組み込まれていない）のでしょう？
http://octave.sourceforge.net/doc/f/nonmax_supress.html

http://octave.svn.sourceforge.net/viewvc/octave/trunk/octave-forge/main/image/src/nonmax_supress.cc?revision=HEAD&content-type=text/plain
を見るとnonmax_supress()はC++で書かれていて、
それをoctaveで動的に利用できるようです（http://www.obihiro.ac.jp/~suzukim/masuda/octave/html/octave_71.html）。

C++のソースコードは公開されているので、
それをコピーして手動でoctaveでコンパイルして、
nonmax_supress()を持っておけばいいのだと思うのですが。。。

>>659
image packageは自分で使わないのでよく分かりません
つい最近，octave-3.0.3 (MSVC)がでたのでとりあえずそれを使って
試してください。だめなら[email protected]に投稿します。

http://sourceforge.net/project/showfiles.php?group_id=2888&package_id=40078

>>309
Octave3.01では
wavread.m, wavwrite.m は不要
sound.m の中身を一行
system("wavplay \"%s\"", wavefile);
とダブルクォーテーションに変えてやる必要がある。

検索してたどり着く人用コメントでした

>>659
文法違いませんか。
以下のように使うのでは，edgeの第一引数はMxNの
グレースケールの行列のはずでファイル名ではないですよ。

help imread，help edge　でhelpを読みましょう。

octave:2> IM=imread('test.jpg');
octave:3> image(IM);
octave:4> BM=edge(IM(:,:,1),'Canny');
octave:5> image(BM);

>>665
すみませんm(__)m。書き方が刺々しいですね。
反省m(__)m。

(;'ｰ`)FreeBSDのportsにあるoctaveのコンパイルオプションに-ffast-mathつけるのやめてくれ

理由があるならそれ書いて send-pr しれ

>>666 確かに変ですね。いつか本家で聞いてみます。

さっきのは>>667、>>668の間違い。m(__)m

>>667 >>668
いま，octaveのMLできいていますが，FreeBSDで聞けといわれたらすみませんが
自分でFreeBSDの方で聞いてください。m(__)m

あ，668 の人だけど >668 は >667 宛てね
(FreeBSD 事情なら >668 ということで)

Octave　MLより
3.0.4 call for patches (Jaroslav Hajek-2 Dec 08, 2008; 03:49pm)

http://www.nabble.com/3.0.4-call-for-patches-to20890571.html

ということでOctave-3.0.4のリリースの準備が始まりました。

>>667

最初の書き込みがいけなかったため，該当する返事はまだ得られていません。
（-ffast-mathオプションが使えるかという風にとられたため。
違う趣旨の返事しか来てません。再度趣旨を訂正しましたがこれまでのところ
答えはついていません）

では，
http://www.nabble.com/-ffast-math-option-at-compling-octave-in-FreeBSD-ports---to20878991.html

>>668
--- John W. Eaton wrote:
> I think you should ask the people who packaged Octave for FreeBSD.
>
> Based on the description in the GCC manual, I would generally
> recommend against using -ffast-math.
>
> jwe

ということで，FreeBSDで聞いてくださいとのこと。m(__)m
jwe氏もagainstと書いている通り，おすすめしないといっておられます。
-ffast-mathなんでついてるのかな。

(;'ｰ`)つーか、-ffast-math付けると計算をはしょるから意図しない答えになることがある

>>676
なるほどそうなんですか。私はJaroslav Hajekさんが書かれたことぐらいしか
知らないですが。どちらにしろつけるのは，変ですよね。

portsのMakefileには-ffast-mathついてないんだけど、
/etc/make.confに自分で書いてないかな？

まあ，私が書き込んだことでoctaveのhelpとmaintainersのMLに残るのでFreeBSD
portsを使っている方へのアナウンスにはなっているわけでだし
あとはFreeBSD のportsをメンテされてる方がどう考えているかですよね。
ところでmkoctfile -vで出てくるg++でのコンパイルの表示に
-ffast-mathついてますか?
基本的にはconfigureで指定したoptionがついてますけどね。

(;'ｰ`)うわっ、ごめんよ！
FreeBSDのportsに-ffast-mathついてるのはmath/atlasだったよ

>>680 了解。上にあげときます。

>>680 上に報告しときました。

誰かovtaveで区分求積法を使って積分の計算をやってる例とか貼ってくれない？

>>683
区分求積って台形公式でいいんでしょうか。
I = trapz (Y) ;
I = trapz (X,Y)；
台形公式法（trapezodial method）による，数値積分をする関数。X,Yはともにベクトルであり，引数がYだけ時はX軸の間隔が1であるとしてyの台形の面積の和をとる。通常は，X,Yの両変数を指定し，和をとる。

octave> format long
octave> x=linspace(0,pi,100+1);
octave> I=trapz(x,sin(x))
I = 1.99983550388744
octave> x=linspace(0,pi,200+1);
octave> I=trapz(x,sin(x))
I = 1.99995887647922
octave> x=linspace(0,pi,1000+1);
octave> I=trapz(x,sin(x))
I = 1.99999835506566

精度が必要な時は，適応的自動積分のquadを使ってください。

[v, ier, nfun, err] = quad (f, a, b, tol, sing)

Quadpackという積分パッケージに収録された適応的自動積分の関数。
第1引数は，被積分関数の名前で関数は，
y = f (x)
の形を持ち，xとyは、スカラーである。第２，第３引数は，積分範囲であり，−無限大や+無限大を片側または両方に指定できる。Octaveで無限大は，infまたはInfであらわす。
第４引数のtol は，ベクトルであり，第一要素は，要求したい精度の絶対値である。第二要素は，相対許容度である。
第５引数のsingは，ベクトルで被積分関数が特異点となる独立変数の値からなり，必要に応じて指定する。
第４，第５引数は，オプション引数なので省略は可能である。
左辺にvは，積分結果の値（関数の返り値）。ierはエラーコードで0は，積分の成功を意味する。nfunは関数の評価回数でと積分の困難さの指標になる。Errは，積分結果の評価誤差である。

octave> format long
octave> function y=f(x);y=1./sqrt(1-x.^2);end;
octave> [v, ier, nfun, err] = quad (@f,0,1)
v = 1.57079632679489
ier = 0
nfun = 315
err = 5.44135847491134e-011
octave> pi/2
ans = 1.57079632679490

『octave for windowsメモの筆者』さんすごい！素敵です！

>>683ですが、貼ってくれた方ありがとうございます！

>>689 >>690

大学の講義で使っている資料の一部を貼りました。
著作権の問題があり公開はしていませんが，
著作権問題がない部分については公開したほうがいいですかね。

ちょっと手間がかかるので要望あればぼちぼちやります。

大学の宿題で「観測値と予測値をプロットするoctaveのプログラムを示しなさい。」
っていう問題が出たんですけど授業を聞いてないボクにはさっぱりです・・・。。
誰か教えてくれませんか？

ノート買えばいいじゃん

Octave 3.0.4のテストリリースのソースのリリースがありました。
しばらく，すれば3.0.4が出ることになります。

ttp://www.nabble.com/3.0.4-RC1-to21495212.html

>>692 まあ，>>693さんがいっているように授業受けた人に聞きましょう。

↓に私の講義資料で著作権に関係のない部分をのせといたからこれでも読んで勉強しなさい。

ttp://www.tatsuromatsuoka.com/octave/jpn/08CompAlgo_02.pdf

予測値ってのはもう用意されてるのかね
予測手法も適切に選択するまでが課題だったら随分大雑把で困難な宿題だ

すいません。
勉強します・・・。

質問させてください。

cos(pi*90/180)

を実行すると

6.1232e-17

というネイピア数で表示されてしまうのですが、これを0と表示させるにはどうすれば良いのでしょうか？

つ round

ところでネイピア数は0.271828...e-1だが？

じゃない、0.271828...e1 だ

というか質問の意味がいまわかった。

そのeはネイピア数じゃなくて、浮動小数点数の10のべき乗表記で、
6.1232 × 10 ^ -17
じゃまいか？

>>699

roundで丸めるのではなく、不動小数点の10のベキ乗表記を小数点第何位までにするかという設定（要するに0の極限値の設定方法）があった気がしたので、そのやり方がわからないのです。

>>701

ネイピア数の定義がよくわかってなくてすみませんでした。
おっしゃるとおり6.1232 × 10 ^ -17です。

octave3.0.4のソースのRC2が出ました。

ttp://www-old.cae.wisc.edu/pipermail/octave-maintainers/2009-January/010310.html

octave3.0.4ももうじきです。
一応3.0シリーズのラストと位置づけられているようで，
octave 3.2シリーズの準備がoctave 3.1で行われています。

>>703
しかし，3.0になるときも，次は3.0とでながらかなり細かいバージョンアップを続けた経緯があるのでどうなるかは。分かりませんが。

cygwin版のoctave-3.0.3がでました。ただし，パッケージングミスがあり,
mkoctfileが動きません。/usr/lib/octaveのlibctave.dll.a(oが抜けている)
をliboctave.dll.aにすれば動きました。報告をして，Macroさんが確認された
ので数日中にはfixされるでしょう。

output_precision=16って入力して表示桁数を変えようと思ったら,
output_precision is no longer a built-in variable
と出ました。
その後output_precition=16となりましたが数値の表示桁数が変わりませんでした。
どうすればよいのでしょうか

>>706
output_precision(16)としてください。octaveはver3から殆この手のものは=ではなくて、()で設定します。

>>706
output_precision(16)としてください。octaveはver3から殆この手のものは=ではなくて、()で設定します。

>>706
output_precision(16)としてください。octaveはver3から殆この手のものは=ではなくて、()で設定します。

>>706
output_precision(16)としてください。octaveはver3から殆この手のものは=ではなくて、()で設定します。

携帯から書込みしたら重複投稿してしまいました。ごめんなさい。m(_ _)m

octaveでシミュレーションをやっているんですがどうしても下のようなエラーが出てしまいます
先生や先輩に聞いてもなぜかよくわからないというしまつ・・・orz
どなたか教えていただけないでしょうか？よろしくお願いします

octave:2> test2
error: invalid column index = 2001
error: evaluating binary operator `+' near line 72, column 28
error: evaluating assignment expression near line 72, column 17
error: evaluating for command near line 67, column 1
error: near line 73 of file `/home/soturon/test2.m'


ソース→http://www.uploda.org/uporg1984631.tex.html

712 ソースの場所に行ってみましたが見れませんでした^^;

めでたしめでたし

すいませんでした！
次こそ大丈夫なはずです

ソース→http://ultraup.net/man/src/u_man2312.zip.html
パス：test

>>715 Get Files Thank youとでますが，どこに落ちたかわからん。
サイトがサイトだけに，落ちないのかも，octave for windowsメモで
Googleれば私のサイトは見つかるのでそこのアドレスに添付で
送ってください。

>>715 Get Files Thank youとでますが，どこに落ちたかわからん。
サイトがサイトだけに，落ちないのかも，octave for windowsメモで
Googleれば私のサイトは見つかるのでそこのメアドに添付で
送ってください。

>>716
m(__)m
そこのアドレス==>そこに記載のメールアドレス　

メール欄にsageと入れよう！

>>719 入れましたがどのような意味で？　おじさんはワカラナイ？？

たびたび申し訳ないです＞＜
メールで添付しときました！よろしくです

つーかアップローダーはどこがいいんだろ・・・
うｐろだとか有名なとこはすぐ消えちゃうし

>>721
メールとどきましたこれから見ます。

Uploaderではなく，Yahooのアカウントをお持ちなら，YahooのGeocitiesをフリーで
広告付きだけど50Mbyteのweb spaceを借りれますよ。
私は，自分のドメインと1Gbyteのスペースと持ちたいから
有料の契約してますが，普段使う分には50Mbyteあればいいでしょう。

>>722
Geocities？ブリーフケースのこと？？

>>720
メール欄にsageと入れると書き込みをしてもスレッド一覧の最上位にいちいち浮上しません。
sageと入れないと一部の人からくどく思われる場合がありますので、普通はsageと入れます。
ちなみにsageは日本語の「下げ」をローマ字で書いたものです。他言語のsageの意味はありません。

>>724 m(__)m　ありがとうございます。

>>724 m(__)m　ありがとうございます。

>>712 error: invalid column index = 2001ですね。

66行目　b(:,1)=[0;0];によって
octave:3> size(b)
ans 2 1
で1行2列の行列になってしまうわけです。

66行目　b=zeros(2,2001);
とすれば大丈夫です。
しかし，For Loopの塊ですね。
octaveでは，遅くなりますね。
漸化式likeなので
スパース行列を使えば高速化できるかも知れませんね。
ttp://www.nuce.nagoya-u.ac.jp/e8/Matsuoka/07OctaveNum/07OctaveNum.html
偏微分方程式なのでまともにやったら2重ループなのですが，遅くなるので，
漸化式を係数行列にしてスパースを使って高速化してます。

古い資料なので，octaveのバージョンが古いですが，基本は同じから
参考にしてください。
でもいまのままでも，答えは少し待てば出るので，まあ，
現状でいいかもしれませんがこういうテクニックも
あることもしっておくといいですよ。

>>723 ブリーフケースではありません。Webをおけるスペースです。
まあ，Yahooにいってみてください。
webがおけるのでファイルもおけて，
<html>
<a href='test2.m'>test2.m</a>
</html>
とでもして，index.htmlをつくり同じフォルダーに
test2.mをいれとけばいいです。

日本のyahooの利用規約では公開目的ではなく単にファイルを置くだけなのは禁止されているはず。
国内で利用するんだったらhttp://www.axfc.net/uploader/が良いかもしれないね。
流れてしまう速度はそんなに速くないはずだから。

>>727 おっとこの講義資料の最後の回のものです。

ほんとループが長い！まだ回ってる・・・

>>728,>>729
つまりヤッホーのほうは利用しないで>>729のほうがいいってことですかね？？
なにはともあれありがとうございます

>>730
是非拝見させていただきます

>>729 あっとそうでしたね。

<html>
<body>
<P>数値計算ソフトのソースの例(octave)
<a href='test2.m'>test2.m</a>
</body>
</html>

とでもしとかなきゃだめでしたね。
こうすれば，一応公開目的なのでおけますよね。
ただし，パスワードがかけられませんからね。
まあ，見られてもいい（つまり公開）であればいいのですよね。

グラフを印刷する時、軸のラベルの文字の大きさを変えるのはどうするんですか？

print "hoge.eps" -deps -color -F:20

ってやると、目盛の数字は大きくなるけど、ラベルは小さい文字のままになりま
す。

>>733
xlabel('{\fontsize{25} x}');

使ってるプロットツールによると思うんだけど
文字はそのままでグラフを小さく出力すると印刷時には相対的に文字が大きくなる
例えば俺はgnuplot使ってるから
__gnuplot_set__ size 0.4 0.4;
とかやって全体的に文字サイズを変えてしまうのはどうだろう

__gnuplot_set__ size 0.4 0.4;
この設定は3.0からは使えません。

ttp://www.gnu.org/software/octave/NEWS-3.The following plot commands are now considered obsolete and will be removed from a future version of Octave:

__gnuplot_set__
__gnuplot_show__
__gnuplot_plot__
__gnuplot_splot__
__gnuplot_replot__
Additionally, these functions no longer have any effect on plots created with the Matlab-style plot commands (plot, line, mesh, semilogx, etc.).html　より


The following plot commands are now considered obsolete and will be removed from a future version of Octave:
__gnuplot_set__
__gnuplot_show__
__gnuplot_plot__
__gnuplot_splot__
__gnuplot_replot__
Additionally, these functions no longer have any effect on plots created with the Matlab-style plot commands (plot, line, mesh, semilogx, etc.).

>>736 すいません。m(__)m。また，確認ミス。
とろろで「この設定は3.0からは使えません。」は使えませんはうそで，
正しくは，plot関数などによるplotには影響せず，__gnuplot_set__などは
あくまでも，別のルーチンでgnuplotのコマンドをダイレクトに使用する描画
にのみ影響する。

が正しいです。

といいのは，octaveの次のバージョンから独自のgraphicsエンジンを
使うという構想があり，このようのgnuplotの直接命令は，plotなどの
関数に影響を与えないように変更されています。

今プログラムを作ってるんですが下のようなエラーがどうしても消えません
どこがおかしいのでしょうか？

error: invalid column index = 0
error: evaluating binary operator `*' near line 35, column 14
error: evaluating binary operator `+' near line 35, column 12
error: evaluating assignment expression near line 35, column 4
error: evaluating for command near line 33, column 2
error: evaluating for command near line 29, column 1
error: near line 42 of file `/home/tse550u/soturon/ndynamicsg.m'


ソース：http://www1.axfc.net/uploader/He/so/190986

すいません、>>738のやつパスはｄｗです

新しいOctave入ってるPCで見てみたらマジだ
昔gset sizeでやってた頃に「これは使えなくなるから__gnuplot_set__ でやれ」
ってOctaveに怒られた時からいつかこうなるとは思ってたけど

>>738
error: invalid column index = 0
が怪しい。行列やベクトルは最初に現れたsizeで初期化されますから，
sizeを指定したいときは，
A=zeros(N,M);
などであらかじめ指定してください。

>>740
今は，JhandleというMATLABそっくりな図を描くグラフィックバックエンドも
ありますしね。MATLAB互換路線になっていますから，まあ仕方ないでしょう。
私は，データ吐き出して,
gnuplotで直接図描きますから関係あまりないのですが^^;。

ところで，jhandleは，epsはファイル形式は，postscriptのみだったと
思います。
emfなどはpstoeditをつかって変換するしかないでしょう。

>>741 error: invalid column index = 0
おっと　octaveはindexが1から始まるからでたエラーですねm(__)m。

>>734-737 ありがとうございました。

自前のグラフィックスですかあ。私はgnuplotが結構好きなんで、微妙だなあ。
なんかOctave全体が重くなるのもいやだし。

筆者さん、Octaveの授業もやってるんですね。わたしも筆者さんと同じ分野の人
間です。授業でOctave使おうとしたことがあるんですが、他の教員から「そんな
ブラックボックス使わせても教育にならん」と反対されました。かといって、現
行の授業でやってるようにC言語だと、病的な文法の瑣末なところで引っ掛る学生
がいっぱいいて、シミュレーションそのものを体感させるところまで行かないん
ですよね。

"プログラミング演習"の例題をOctaveで解かせるのは、研究室の学生にさせたこ
とがあります。"プログラミング演習"は過去におそらく3種類出版されているかと
思いますが、わたしは80 年代末に出たBasic版が一番のお気に入りでした。基本
算法の解説が一番詳しかったように記憶しています。現在はFORTRAN版もBasic版
も入手困難なようですね。FORTRAN版は手元にあるんですが、BASIC版はどっか行っ
ちゃいましたorz。今出まわってる改訂版は、うーん、解説は適当だし、Cと
FORTRAN両方のコード載せるのも場所の無駄だし、おまけに今時FORTRAN77だし。

いきなりOctaveやらせると感動がないじゃん
「スゲーこの言語スゲー!!俺こんなに手抜いちゃって良いの?」
っていう感じで

FORTRANでもいいんじゃね

Matlab/Octaveの醍醐味は初級プログラミングのレベルではなくて
普通にループを回すと何時間も何日もかかるループ処理が
行列演算で１発で出るという所なので
Excelで済む程度のデータなら必要無いでしょうね

逆に例えばモンテカルロ計算なんかだと、
一般的なプログラムではループ10万回ぐらい回して
そのループの中で乱数を発生させて、関数に入れて結果を見る事になりますが

Matlab/Octaveでは最初に、10万個の乱数配列を作り
行列演算化した関数に入れて1発で結果が出るみたいな感じなので
発想を根本的に変えないといけません

>>748
そうですね。いかに，行列で問題を記述するかにかかってて，
偏微分の差分法の陰解法では良くガウスザイデル法とかつかいますが，
Octaveなら行列で方程式書いて左除算で済ませる。

係数行列は疎行列だからsparseを使って高速化するとか
発想が違いますよね。

まあ，ブラックボックスっいう方いますが，硬い微分方程式なんぞ
ルンゲクッタで解こうとすると大変で，lsode使うとあっという間に
解けますよね。

適応的自動積分でも同じだとおもいます。

シンプソン公式を自分でCで書いて使うような時代ではないと思ってます。

ブラップボックスでもいいから，数値計算ができることの意味を分からせる
のが大事だとおもいます。

開発版octave 3.1.52がでました。
ttp://www.nabble.com/Octave-3.1.52-available-for-ftp-to21923205.html#a21923205

Octaveの次のMajor Version Upはver 3.2であり，octave-3.1.xxは開発版です。
その3.1.52がでました。

gnuplotは一応defaultのgraphic backendでopenglによるgraphic backend
を機能させるためにはいろいろライブラリィが必要なようです。

詳しくは以下を
ttp://www.nabble.com/Octave-3.1.52-available-for-ftp-to21923205.html#a21923205

そういえば新パソにFreeBSD7.1を入れて、
pkg_add -vr octave-devel
とインストールしたら
GNU Octave, version 3.1.51
でした。leasqrが入って無いのでoctave-forge-optimを入れてないような気がする。
バージョンごとに何が良くなったのか、見えにくいのが何とも…

そういえば、何語なのかわからないけど（チェコ語？）
Octave/Matlabのグラフィックスの例題があります
http://enacit1.epfl.ch/cours_matlab/graphiques.html
gnuplotは置いておいて、jhandleならかなりMatlabに近い事が出来ました

>>751 3.1.xx はdevelopment branchなので，octave-forgeは多分
自分でpkg installで行う必要があると思います。
しかも，octave-forgeは3.0.x用に作られているので結構苦労する可能性は
あると思います。

一般的なものは3.0.xxです。現在3.0.3が，最新版で3.0.4はもうすぐ出るのだと
思います。

(独り言)3.0.4RC3が出てからだいぶたつけど，いつになったら3.0.4の
正式リリース出るんだろう。