トップページ⇒tech
63コメント12KB

+-*/ 数値解析プログラミング /*-+

■ このスレッドは過去ログ倉庫に格納されています
0001main( ){ int x; x=1+1 printf(x); }NGNG
ただでさえ少ない数値解析プログラマ。
固有(値)の悩みも多いはずだ。
ここで情報交換しよう。

C/C++とFortranユーザ間での罵りあいはやめてね。
(いないかもしれないが)
WinとUnixユーザ間の罵り間もいやだよ。

0002デフォルトの名無しさんNGNG
";"抜けちまった....鬱だ...。
00033NGNG
>>2
このスレもうだめポ。
0004デフォルトの名無しさんNGNG
>>1
まず、おまえのハンドル名のデバグが先だ。
0005main( ){ int x; x=1+1; printf(x); } NGNG
>>4
デバグしますた.....。
0006デフォルトの名無しさんNGNG
printf("%d", x);
0007main( ){ int x; x=1+1; printf("%d", x); } NGNG
>>6
デバグしました...........(えーん)。
0008デフォルトの名無しさんNGNG
>>1のハンドルネームのデバッガか、ここは.....。
0009デフォルトの名無しさんNGNG
#include <stdio.h>
0010 ◆XjKqCmiQ NGNG
>>9
もうダメぽ.....。
0011 ◆XjKqCmiQ NGNG
>>10
ついにハンドルネームの限界に来たか.....。
00121NGNG
あきらめました。
0013デフォルトの名無しさんNGNG
なんだか久しぶりに笑った。
00141NGNG
ハンドルネームが限界長こえると、IDが入るなんて知らなかった。
0015デフォルトの名無しさんNGNG
>>14
> ハンドルネームが限界長こえると、IDが入るなんて知らなかった。
違う
#include と入力したために先頭に「#」が入力された為。

0016デフォルトの名無しさんNGNG
↑違うって。
001716NGNG
>>14 へのレスです
0018デフォルトの名無しさんNGNG
↑違うって。
0019デフォルトの名無しさんNGNG
で、数値解析は?
0020デフォルトの名無しさんNGNG
↑違うって。
0021デフォルトの名無しさんNGNG
シンプレクティック、シンプレクティック、シンプレクティック…
何回練習しても舌噛みそう。
あと、プラストミックとなぜか混ざる。
0022デフォルトの名無しさんNGNG
GRAPEって知ってる?
ttp://www.iam.uni-bonn.de/sfb256/grape/
ひょんな所から、これのデータを吐くプログラムを使うことに成ったのだが、全然分からん。
このデータをjpegやgifなんかのデータに変換したいのだけれども、できるもの?
0023デフォルトの名無しさんNGNG
PrintScreenじゃだめなの?
002422NGNG
>>23
残念だが、駄目。時間発展問題だから少なくとも吐きたいデータは100枚ほどある。
0025デフォルトの名無しさんNGNG
みんなFortranかい?
0026デフォルトの名無しさんNGNG
>>21
>シンプレクティック
トポロジーですか?
0027デフォルトの名無しさんNGNG
モレは流体。
0028 ◆UXqvbfSY NGNG
0029デフォルトの名無しさんNGNG
希少なスレにつきage
0030デフォルトの名無しさんNGNG
マ板向きの>1がいるスレはここですか?
0031デフォルトの名無しさんNGNG
ていうか、MATLABとか使わんの?
0032デフォルトの名無しさんNGNG
>>31
予測計算には使えるが、精度にかんして確信が持てない。
0033デフォルトの名無しさんNGNG
自分よりよっぽど信頼できると思うが。
0034デフォルトの名無しさんNGNG
C++ 用の行列計算ライブラリの良さげなヤツをきぼ〜ん
003534NGNG
↑ちなみに今迄はTNTを使ってたけど、試しにBlitz++を評価中。
0036デフォルトの名無しさんNGNG
>>34
どんな計算に使うんだ?
003734NGNG
>36
いろいろだけど量子計算かな。
ターゲットマシンはPCからSMP,Vectorまでいろいろっす。
0038デフォルトの名無しさんNGNG
>>37
どんくらいの大きさの行列で、精度はどんだけほしいの?
行列の要素の分布具合が悪いと、手作りしないと丸め誤差
が膨らんで手におえなくなる.......。
003934NGNG
# Blitz++ が g++-3.X で使えないという罠?

厨房質問で恐縮っすが
>38
シビアな場合は論外だと思っていたんすが、
ケースに応じて有効なライブラリってあるんすか?

取り敢えずLAPACK程度の使い勝手で複素数が無理なく扱えて、
パフォーマンスがそこそこ良いものが所望っす
#LAPACK程大規模でなくてもいいんですが
0040デフォルトの名無しさんNGNG
># Blitz++ が g++-3.X で使えないという罠?
http://oonumerics.org/MailArchives/blitz-dev/msg00457.php
これのこと?

# oonumerics 死んでるんで, google 経由ではっけん
004140NGNG
ついでに質問. TNT ってダメなの? >> 34
# Blitz は PII 155MHz 64MB RAM だと install すらできなかったから寅馬
004234NGNG
>40
THANX!

>41
悪くはないっす。
ただ、複素数を取り扱うときはそのままではダメだったやうな。
ソース見ればすぐ修正できる(absね)はず。
現バージョンは未見。

Blitz++ 使う気になったんは、ソース見て Expression Template のテク
にクラクラしたせいという、つまらん理由っす。

ところで何処かにC++で複素行列のSDV分解できる(Freeな)ライブラリないすか?
004334NGNG
>41
ちなみに漏れの使ってるDebian GNU/Linux(sid) にはパッケージが
用意されてるっす。
0044デフォルトの名無しさんNGNG
>>39
基本的にうぇる・ぽーずな問題なら使えるものは
転がっていないこともない。(探してみ)
私は病んだ行列ばかり扱うのでシコシコ組むしかない。
004540NGNG
>42,43
なるへそ. そろそろ ver. 1.0 が出るということだったので興味が
あったのですが.
# ちなみに僕も deb 使い.

#量子計算 あんど SVD ということは Schmidt 分解でも求めるので
#しょうか....

>44
質問です. Hermite 行列とか Unitary 行列でも病んでることって
あるんでしょうか?
# 散逸系はやばいかもしれませんけど...











004634NGNG
>44
サンクスコ
004734NGNG
>45
Header Fileだけで構成されているし、ソースも短くコンパクトで
読みやすいので、問題があっても簡単に対応できるという意味では
良いライブラリだと思うっす。だけど、TNTが想定していない型(複素数型も)の
変数を使おうとしたら、ちょっと弄らんといけないかも。
0048デフォルトの名無しさんNGNG
お前ら素直にシミュレーション板逝ったほうがええのんちゃいます?
http://science.2ch.net/sim/
0049デフォルトの名無しさんNGNG

 ま た 大 阪 人 か
0050デフォルトの名無しさんNGNG
>>44
E行列もU行列も、計算機で取り扱うときは要素の分布によって
非常に困るときがイパーイある。引き算による桁崩れもおきるし、
乗除による丸め誤差の拡大とか、割り算でマシーン・ゼロとか。
モレは最近は多倍長で逃げるが、時間かかるなあ。
0051デフォルトの名無しさんNGNG
>>50
Hermit行列 を数値解析の人はE行列と略すの?
005250NGNG
>>51
痛い突っ込みサンクス。そう、H行列だ。
(モレ頭の中ではE(良い)行列なんだよ)
0053デフォルトの名無しさんNGNG
>>48
シュミレーション板っていうのはウォー・シュミレーション専用
と思っていた........。
0054デフォルトの名無しさんNGNG
シュミレーションは趣味かと思っていた.......。
0055デフォルトの名無しさんNGNG
>>53-54
わかっているとは思っていますが、「シミュレーション」です。
005634NGNG
>50
漏れのばやい、対角から離れた要素の値が指数関数的に減少するんで
(関数基底がバレるぅ(^^;)行列サイズが大きくなると、detが常に機械零になるのが
悩みっす。

多倍長は何使われてます?
005734NGNG
>45
今頃気が付いたんすが (^^;、
パッケージが用意されてるのは Blitz++ の方っす。
# TNT は Header だけなんで、適当に置いとけば OK
005834NGNG
試しに LiDIA 入れてみますた…漏れが使うんぢゃないけど
0059デフォルトの名無しさんNGNG
>>56
GMP。いろいろ使ってみたが、これが一番早くて、関数もまあまあ
信用できるレベルだ。三角や指数はつくらないといけないが。
006040NGNG
>三角や指数はつくらないといけないが。
これつらすぎます.
cln はいかがでしょ?
006134NGNG
>59
なるほろ、どーもです。

>60
Taylor展開で収束するまで計算すれば良い(…のか?)けど、
重たそう
0062デフォルトの名無しさんNGNG

おい、おまいら
           ___ 00        /
|__ __  ____         /
|    __|   |   ____  /|
|__ __|   /              | 

を流行らすのじゃ!

http://comic.2ch.net/test/read.cgi/comic/1025843878/551-591


0063デフォルトの名無しさんNGNG
今仕事がなくていやん。
■ このスレッドは過去ログ倉庫に格納されています