Ｃ＋＋マトリックスクラス

C++愛好者の皆さん，マトリックスクラスは何をお使いですか？
４則演算等はもちろんのこと，LU，固有値，特異値，一般逆等
もサポートされているものが望ましいです．

Lapack++やtntは激速だけど，使用が不完全すぎます．　かと
言って，ちまたのwebで素人が配布してるのは遅すぎるし，ド
キュメントが充実していないのがほとんどです．

この際，Lapack並みの速さは求めないとして，ドキュメントの
充実しているC++マトリックスクラスは無いものでしょうか？
教えて下さい．

無い。

...と思う。

無いんですか．．．
では，皆様がお使いのマトリックスクラスは何ですか？
「自作!!」と胸を張って答えて下さる方，是非入手先を
教えて下さい．
ちなみに，僕自身はLapack++を渋々使っています．

こんなページ見つけました．
http://webnz.com/robert/cpp_site.html
でも，どれが良いんだろう？
（VineLinuxのVinePlusにはnewmatというのが入っている．）

>>4
ありがとう。
ちょっと探検してみます。

そんなに行列を使いたいなら、Cでなくてもいいでしょうに。
COBOL使え。Sなんてのもあるぞよ。

　COBOLとかSはマトリックス数値解析に使うような言語なのでしょうか？
Lapack並みの計算速度 or C++やJAVA等によるクラスならではの操作性の
どちらも無い気がするのですが？　無知なもので．．．

　是非COBOLやSでマトリックス数値計算をする利点をご教授下さい．

COBOLで・・・・なんて冗談ですよ。ちなみに＞６じゃないけど。
Ｓならわかる。
でも、その分野ならMATLABそのものじゃないですか。
C++から使いたければ、簡単に関数呼べますよ。

冗談きついです．

> でも、その分野ならMATLABそのものじゃないですか。
> C++から使いたければ、簡単に関数呼べますよ。
MATLABは有料なのが問題なんです．　また，移植性の
ことを考えると，C++のソースのみで完結させたいです．
でも，最近のMATLABは内部でBlas,Lapackを使っている
らしくて，異常に計算速いと評判です．

どなたかC++マトリックスクラスユーザはいらっしゃい
ませんでしょうか？

>ことを考えると，C++のソースのみで完結させたいです
MATLABで作った関数をCに落としちゃうんじゃだめ？

MTL 使ってる俺は厨房か?

俺は行列と言えば、３ＤＣＧで使うぐらいだからな…
というわけでDirectXのD3DMATRIX構造体を使っている。

ヤフーオークションで、幻の人気商品、発見！！！

今は無き「コピーガードキャンセラー」↓
http://page2.auctions.yahoo.co.jp/jp/auction/b18505597

ヤフーオークション内では、現在、このオークション
の話題で、持ちきりです。

>>8
どうやるの？

http://www.adultcross.com/toybox/

>>12
javax.vecmathも使えるかな？
　ttp://objectclub.esm.co.jp/vecmath/index-j.html

この板の住人でJava使ってるやつっているの？

MATLABは，Lapackの上にラッパークラスを被せているだけなので，
本質的にLapackを使うのとは変わらないように見えるが，
Matlab上でメモリアロケーションの工夫は，ほとんど出来ないので，
大規模数値計算プログラミングをするならば，MATLABはほとんど使えない．
せいぜい，1000×1000のフルマトリックス位が精一杯てとこだな．
自作しないでやるなら，newmatが一番完成度高いと思う．
lapackのラッパークラスを自作した方が，何かと都合が良いとはおもうけど．

>>18
嘘言っちゃだめ！！

http://www.townpaage.ne.jp/adult/rankem.cgi?id=sawada

http://sources.redhat.com/gsl/

GSL(GNU Scientific Library)って基本的にＣですよね？ >>21
マトリックス演算を"A*B"の様には書けないと思うのですが．．．
もちろん，ラッパーを作れば出来そうですが，あるんでしょうか？

自作

>>19
どこが嘘なのでしょうか？
MATLABでは，フルマトリックスと，スパースしか作れないので，
有限要素法のようにバンドマトリックスが必要となる場合には，
スパースからバンドへの変換をしなくてはならず，これは，非常に
大規模な有限要素解析の際には，致命的になります．ここは，反論あります？
あと，MATLABは，Lapackの上にラッパーを被せているというのは，間違えないと
思うのですが．確かに，色々拡張機能はありますが，行列計算に関するルーチンに関しては，
LAPACKの拡張と考えて差し支えないと考えていますが，いかが？
1000×1000のフルマトリックスが関の山と書いたのは，それ以上出来ないと
言ったのではなくて，それ以上の規模になったら，プログラム自作が一番良い
ということが言いたかったわけです．
反論は具体的にお願いします．

確かに嘘書いていました．
MATLABで，1000*1000行列のLU分解は１０秒で済みました．
フルマトリックスの効率が悪くないというのは，大いに納得しましたが，
果たして，数十万*数十万のバンドでは計算できるのかどうか疑問です．
どうなのでしょうか？

プログラム板が妥当かと思った冬の昼下がり

なんでMATLABってあんなに早く計算できるの？

俺もMTL使っている。
厨房だ

>>9
遅レスだけど
有料が嫌ならOctaveやScilabを使えばいいじゃん．

win環境で、2*2行列だけ、昔（学生時代）自作したのを拡張して使ってる。

　　　　　　　　（ ﾟフﾟ） ﾊｲﾊｲﾊｲﾊｲﾊｲ!!
　　　　　　　　（つつ
　-＝≡（　ﾟ ﾟ）　> >≡(　ﾟ ﾟ)
　　　　　　↑　�W�W　 　↑
　　　マーガレット　 フランソワ

　　　　　　　 　  （ﾟフﾟ） ｺﾉｽﾚ････
　　　　　　　　　　⊂Ｕ⊃
　　　　　　  （　ﾟ ﾟ）　∧　　(ﾟ ﾟ　)
　　　　　　　　 　 　 �W�W

　　　　　　　　　　 （ ﾟフﾟ） ｵﾓｼﾛｲ!
　　　　　　　　　　　(つつ
　　　　　 （　ﾟ ﾟ）　　 / |　　(　ﾟ ﾟ)
　　　　　　　　　 　 �W�W

　　　　　　　　　 ⊂（ﾟフﾟ）つ ｻｲｺｰ!!
　　　　　　　　　　　   Ｕ＿
　　　　　　（ ﾟ ﾟ ）　　/　　�W （ ﾟ ﾟ ）
　　　　　　　　　　 �W

　　　　　　　　（ﾟフﾟ ） ﾃﾞﾊﾏﾀｺﾞｼﾞﾂ｡
　　　　　　　　⊂⊂） )))
　　　　（ﾟ ﾟ　）)) < <　　(ﾟ ﾟ　) )))
　　　　　　　　   �W�W　

／￣＼
（´υ｀　）
＼д／

／￣＼
（´υ｀　）
＼ д ／

　／￣＼
│´υ｀　│
　＼ д ／

　／￣＼
│´υ｀　│
　＼ д ／　　　　　　　

マトリックスって何？
誰か分かりやすく教えて！

>>36
赤いカプセルを飲めばわかるよ。

　C++の話題ではないです。参考意見程度に聴いてチョ。

　自分は、マトリクス処理するときは、APLで組んでたけど。。
　もともと言語仕様にマトリクスが含まれてるから、余計な事
（行列演算パッケージのインストール）しなくても使えたし。
　APLで面倒だったのは、コーディングのときかな？特殊な文
字を色々と使うので、そこいらへんを憶えるのが・・・。

俺は自作。
行列の雛形をbasic_matrixテンプレートクラスで作っておいて、
要素の格納形式に応じて、それから派生させたmatrix、skyline_matrixその他。

ネットでいろいろ見たんだけどね。
自分の要求を満たす物はなかったので。

やっぱりvalarrayからつくったほうがいいのかなぁ？

.　　　　 l ::::|l　.i 　 ::::::::::;;;',;::;ﾞ､:;;',､:::::;;;;;l　,,､-､'""､;;:;::::;;;',.　ヾ;;ﾞ､　,r<_,.,､､zﾞ'､､ﾞ;､ヾ;ﾞ､;;;;;;l;;;!;;::::::::::::
　　　 　 ',::::l ! :.l. . .:::::::::::;;;;',;:;;,ﾞ､;;',ﾞ､::;､rl'"　 ,,.ﾞ､',..,ﾞ'､;;;;;;;;ﾞ､ 　 ヾ;､　,,,;;ｼ''!"~~ﾞﾞ';,ヾx､.ﾞ､;;;!;;l;;;:;::::::::::
　　 　 　 ',:::l ',:;;',.:ﾞ､;;;::::::::;;;;;',;;;;, ヾ;r'､::;;;;!.r'',,､zｪｭ､.　ヾ､;;;;ﾞ､ 　　 ヽ''"　 l,　 ,,､:':'､　ヾ:;､;!;l;;;';;::::::::::
.　　　　　 ﾞ､:!. ',:;',､::ﾞ､;;;;:::::;;;;',ﾞ､;;／ヽ.ヾ;;;'ill!'''"ﾞi　~`　　 ﾞヾ､ﾞ､.　 　 　 　 lr''ﾞ!::li!::;i,　　ヾj//;;;;:::::: 　
　　　　　　 ヾ　ﾞ､:',ﾞ､:ﾞ'､;;;;::::;;;ﾞ､ヾ;､　ヽ;lrヾ､　 .,i,　　　　　　 ｀`　　　　 　 l;;::ゞﾝ;::;j　　 ﾂ;/;;;;;;::::;;
　　　　　　　 ﾞ､　ﾞ､.!:ヽ､ヾ;;;;;;;;;ﾞ;､ ヽ､;ll!.　l　 ,ｒ';::l,　　　　　　　　　　　　　 l;;;. .::: ,;;j　　',;ｲ:;;;;;;;;::;;;;
　　　　　　　　 　 ,'ヽ､;:::l';,､ﾞ､r､;;ﾞ;､. ill!.　 ';;:'ﾞ'i,,j;::l,　　　　　　　　　　　､,.　ゞ--一ﾞ-　,//:;;;;;;;;;;;;;;;
　　 　 　 　 　 　 i　::::;:::!;;;;;;;lﾞ'',､.ﾞ､ヽll!.　　'i;:: ﾞ''ﾞ:. l､　　　　　　　　　　　　　　 　 　 ,/'/;;;;;;;;;;;;;;/;;
　　 　 　 　 　 　 ! :::: :::l;;;;;;;;;! l　 .ﾞ､. ヾ､.　 ﾞ､;;;;:､‐'''''"　　　　　　　　　 　 　 　 　／ /;;;;;;;;;;;;;ｼﾂ'
.　　　 　 　 　 　 |.:::;:: :::l;;;;;;;;;l　 ､〈ﾞ､　　　､　'´ 　 　 　 　 　 ,　　　　　　　　　　´　,.';;;;;;､-'" ´
.　　　　　　 　 　 l::::; ::i:::l;;;;;;;;;ヽ, ヽ､ﾞ､　　　　　　 　 　 　 　 ;l　　　　　　　　　　 ／''" / 　　藍より青し 公式ホームページ　　　　
　　　　　 　 　 　 l::::;;:::',::',;;;;;;;;;ﾞ､ﾞヽ､. ﾞ､.　　　 　 　 　　　　　ﾞ'　　　　　　 　 　 　 　 ,ｲ 　　　http://www.aiyoriaoshi.com/index2.html　　
　　　　　　　　　　',:::;;;:::',::',;;;;;;;;;;r'"~ﾞ''ｰ'､　　　　　　　　　　　-‐ '''''"　　　 　 　 　 ,.';;;! 　　 　２話　- 夕餉 - 画像
　　　　　　　　　　 ﾞ､:;;;;::;',:ﾞ､'''"　　　　　 ヽ､　　　　　　　 　 　 ､:::'　　　　 　 　 ／!;;;i!l 　　　http://wqio6xt.hp.infoseek.co.jp/aiao/idx022.html　　

人間は、地球をほろぼすウイルスだ！！
今いる世界は現実ではないことはないのだ！

ublas

http://www.boost.org/libs/numeric/ublas/doc/overview.htm
http://www.kmonos.net/alang/boost/classes/ublas.html

（＾＾）

（＾＾）

何のためのC++だ
自分の要求にあった派生クラス作れ

（＾＾）

　　 ∧＿∧
　　（　　＾＾ ）＜ ぬるぽ（＾＾）

━―━―━―━―━―━―━―━―━[JR山崎駅（＾＾）]━―━―━―━―━―━―━―━―━―

　　　 　∧＿∧
ﾋﾟｭ.ｰ　(　　＾＾ ） ＜これからも僕を応援して下さいね（＾＾）。
　　＝〔~∪￣￣〕
　　＝ ◎――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉

　__∧＿∧_
　|（　　＾＾ ）|　＜寝るぽ（＾＾）
　|＼⌒⌒⌒＼
　＼ |⌒⌒⌒~|　　　　　　　　　山崎渉
　　 ~￣￣￣￣

　__∧＿∧_
　|（　　＾＾ ）|　＜寝るぽ（＾＾）
　|＼⌒⌒⌒＼
　＼ |⌒⌒⌒~|　　　　　　　　　山崎渉
　　 ~￣￣￣￣

　　 ∧＿∧
　　（　　＾＾ ）＜ ぬるぽ（＾＾）

　　　 　∧＿∧　 ∧＿∧
ﾋﾟｭ.ｰ　（　 ・３・） (　　＾＾ ） ＜これからも僕たちを応援して下さいね（＾＾）。
　　＝〔~∪￣￣￣∪￣￣〕
　　＝ ◎――――――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉&ぼるじょあ

今すぐ会いたいの。ひろみはあなたなしじゃもうだめみたい。
私、ここにいるから・・・・・探しに来て、くれる？
７日間会費フリー、１０分間無料になってるの、だから来て。
あなたに会えなくちゃ、寂しくて死んじゃいそうだから
待ってます。来てくださいね！
→→→http://www.gals-cafe.com

cpplapackhttp://cpplapack.sourceforge.net/

エロサイトのサンプル
http://homepage3.nifty.com/coco-nut/sample/sample.html

堤さやかちゃんの引退記念作です。
これは絶対見るしかないでしょう。
甘えたしゃべりかた、小さな身体に大きなオッパイ、そしてこの顔。
どれをとっても特Ａ級！こんな子がＡＶ女優だったなんて信じられませんね。
無料ムービー観てちょ
http://www.exciteroom.com/

スレンダーなロリ体形ながらもオマンコ方はしっかり発達しちゃってます。
大きなオチンチンを喉の奥までくわえ込んでのフェラは大人顔負けですぞ。
挿入シーンは少ないながらも小さなオマンコにしっかり挿入されているシーンが勃起モノ！
若い娘はいいですな。

無料ムービーはこちら！
http://www.cappuchinko.com/

http://life.fam.cx/a011/

　　　 (⌒V⌒)
　　　│ ＾ ＾ │＜これからも僕を応援して下さいね（＾＾）。
　　⊂|　　　　|つ
　　　（＿）（＿）　　　　　　　　　　　　　　　　　　　　　　山崎パン

　　　　.,. -──-､
　　 ／. : : : : : : : : ＼
　 /┛┗ : : : : : : : : : ヽ　 l⌒)OOo　　　　　∩＿＿＿∩
.　{:┓┏ '⌒'　 '⌒'i: : :}　 {　　）　 )　　　　　| ノ --‐'　 ､_＼
　{:: : : : |　ｪｪ　 ｪｪ |: : :}　 ｜　::::::/　　　　/　,_;:;:;ノ､　　● |
.　{ : : : :|　　　,.　　 |:: :;!　 /　　::::|，,･_ 　 |　　　　( _●_)　 ミ
.　ヾ: :: :i　r‐-ﾆ-┐| ::ﾉ /　　　　）’,∴ ･¨彡､　　　|∪|　　ミ
　　 ゞイ! ヽ 二ﾞノ ｲゞ　 /　　 /￣∵ ’　　／　　　 　ヽノ￣ヽ
　　/　＿` ｰ一'´￣　 /　　　/　　　　 ／　　　　　　　／＼ 〉
　　(＿＿＿）　　　 /　　　 /　　　　/　　　　　　　 /

test

boostのublas。これ最強。

そうかぁ？
あれにもいろいろ不満が・・・

http://etc.2ch.net/test/read.cgi/charaneta/1075974266/745

:::::::::::/　　　　　　　　　　 ヽ::::::::::::
:::::::::::|　　ば　　じ　　き　　i::::::::::::
:::::::::::.ゝ　か 　 つ 　 み　 ﾉ:::::::::::
:::::::::::/　 だ　　に　　は　ｲ:::::::::::::
:::::　 |　　な。　　　　　　　ﾞi　　::::::
　　　＼_ 　　　　　　　 ,,-'
――--､..,ヽ__　　＿,,-''
:::::::,-‐､,‐､ヽ.　)ノ　　　　　 ＿,,...-
:::::_|/ ｡|｡ヽ|-i､　　　　　 ∠＿:::::::::
／. ` ' ● '　ﾆ ､　　　　 ,-､ヽ|:::::::::
ニ　＿_ｌ＿_＿ノ　　　　 |・ |　|, -､::
/￣　＿　　| i　　　　　ﾟr ｰ'　>>1 |::
|(￣`'　　）/ /　,..　　　 i　　　　　'-
`ー---―' / '（__ ）　　　ヽ 、
====（ i）==::::／　　　　　　,/ﾆニニ
:/ 　　 　ヽ:::i　　　　　　　/;;;;;;;;;;;;;;;;

荒れてるな

blitz++

>>1
Java製だけど

Jakarta Commons Math
http://jakarta.apache.org/commons/math/

というのがある。

まだまだ不完全で改良の余地はあるかな。
早く目的にあったものを作ってほしいなら
自分でもこのオープンソースソフトウェアの開発にコントリビュータ
として参加してみるのもいいと思う。

C++ のmatrixクラスか。使い勝手が良くて効率が良いライブラリはない。

なんせ、演算子のオーバーロードを使うと、関数値としてマトリックス
を返すときに、要素が全部コピーされるという絶望的な状況に陥るし。

それを避けるために、閉包や数式テンプレートのテクニック
使うと、呪文みたいに難しくなる。

あと、演算子＊＋−の両側のオペランドが作業用の一時マトリックス
かどうかをチェックしながら、グローバルオブジェクトの参照を返す
というテクニックもあるけど、マルチスレッドでは無力。

結局、今のところ決定的なテクニックはない。FORTRANやCみたいに
シコシコ引数を渡すのが一番効率的。

あれてるなぁ．使ってるのはcvmat．

>>9
遅レスだけど
有料が嫌ならOctaveやScilabを使えばいいじゃん．

ublas + OpenTissue

TNT

演算子のオーバーロード遅いって言ってるやつへ
遅いがまじめに返すと、演算子のオーバーロードをしたら遅くなるって言ってるやつ多いが、
演算子は普通の関数と処理速度は変わらんから、組み方が悪いだけと言える。
読み取らせたいだけなら、参照渡しで渡せば良いだけだしな

演算子のオーバーロードがどのぐらい利いてくるかちょっと試すためにdouble型の要素を3つ持つ、
ベクトルの足し算が出来るの作ってみるとさ、確かに、普通にコンパイルすると、
多少こっちの方が遅かった。
でも、最適化を有効にしてコンパイルすると全く速度が変わらなかった。

結局アルゴリズムが同じなら、最適化されればどっちも同じ速度らしい。

ttp://www.esnips.com/doc/c30b4188-1e9b-48ea-9c7a-bd902b44fcea/inst_all_new
例
#include <instMatrix.h>
using namespace inst;
#define MAT(rows,cols) MATRIX<Float,rows,cols>
#define VEC(rows) MATRIX<Float,rows,1>
const VEC(n) X0 = {-1.0f, -1.0f, -1.0f};
const MAT(n,n) W =
{0.2f, 0.0f, 0.0f,
0.0f, 1.2f, 0.0f,
0.0f, 0.0f, 0.2f};
const VEC(n) T = {-0.9f, 0.0f, 0.9f};
VEC(n) NET(Int t) { return W * X(t-1) - T; } // 便利
VEC(n) X(Int t) { if(0==t) return X0; else return Each(satlins, NET(t)); }

ＰＳ．俺の名前空間は、INST_DEBUGを予め定義してください
ＰＳ２．その名前空間内のリストは、（test用のみだから）iteratorはまだ実装していないから、
本格的なリスト操作なら、stlを使ってください。

>>76

ベクトル演算くらいで偉そうにグダグダ言うな。この馬鹿！
やるならマトリックスの乗算で議論しろや。

あっ、言っておくが、RVOのテクニック使っても、(A＋B)*C*(D+E)の
ような(演算子オーバーロードによる)演算では要素コピーはほんの
一部しか避けられん。

じゃ、数式テンプレートと言い出すだろうが、行列のサイズが大きく
なるにつれてコンパイル時間が異常にかかってしまう。サイズ10×10
以上では事実上不可能。

76じゃないが、
＞(A＋B)*C*(D+E)
速度を気にする場所でこんな書き方するほうがおかしい。
temp1 = D;
temp1 += E;
temp2 = A;
temp2 += B;
temp2 *= C;
temp2 *= temp1;
こうすりゃ一緒じゃねーの？
あと数式テンプレートってなんのこと？
まさか値がコンパイル時にわからないもの同士の計算をコンパイル時に出来ると思ってないよね？

めっさ勘違いしてた。
たまたま自前のベクトルクラス作ってたんだが、
＞(A＋B)*C*(D+E)
こう書いた方が速かった。
（環境はVC2008EE、operator +,*は当然インラインメンバ関数）

下手に一時変数を自分で書くよりも、
メモリからのロード順や計算順序を最適化しやすいのかもしれない。
逆アセ見たら、無駄な一時オブジェクトを発生させることもなく
計算用レジスタに直接ぶちこまれて直接計算してた。
というわけで
＞RVOのテクニック使っても
＞要素コピーはほんの一部しか避けられん
大嘘こいてんじゃねえ。

↑馬鹿

計算結果が正しいかチェックしてみな

馬鹿はお前だ

>>81
式テンプレートすら知らんお前がしゃしゃりでてくんな糞野郎