- 1 名前：デフォルトの名無しさん [04/12/31 00:46:46]
- はっきり言ってマイナーな分野。
やってる人も少ないだろうし市場も小さいと思う。
技術の話でも業界の動向でも何でもいいから
どんどん書き込んでくれ
- 102 名前：デフォルトの名無しさん mailto:sage [2006/02/16(木) 00:17:19 ]
- って、複素数と読み間違えてた・・・
ってか、整数？
FFT は整数だと近似演算になるぞ。
- 103 名前：デフォルトの名無しさん mailto:sage [2006/02/16(木) 00:53:11 ]
- >>100
ttp://www.jjj.de/fft/fftpage.html
- 104 名前：100 mailto:sage [2006/02/16(木) 20:34:05 ]
- >>103
どうありがとう
いくつかありますね。見てみます
>>102
整数が近似演算なのは了解済み
計算量が多く、演算時間が十分にとれないため実数演算と許容できる精度で近似がとれていればOK
- 105 名前：デフォルトの名無しさん mailto:sage [2006/02/17(金) 06:23:15 ]
- 整数でやるとはやくなるもん？
個人的には、SSEのついてるCPUならfftwで
./configure --enable-float --enable-sse
でつくった短精度FFTが相当速いと思ってるんだけど、
暇があったら比較して。
フラグはFFTW_PATIENTで。
- 106 名前：100 mailto:sage [2006/02/17(金) 18:22:23 ]
- >>105
>>個人的には、SSEのついてるCPUならfftwで
基本的には合意
ただ、今回はマルチプラットフォームでPPCなどでも走らせるので
SSEはおろか．．．コプロもない
そんなもんで数値計算するのが悪いんだが、猫の手も借りたい用途
ただ、個人的には猫の手は？
- 107 名前：100 mailto:sage [2006/02/17(金) 18:23:54 ]
- >>105
あと、元データが32bit整数なので、実数であればdoubleが必要
- 108 名前：デフォルトの名無しさん mailto:sage [2006/02/18(土) 03:35:15 ]
- doubleも３２かもしれない
- 109 名前：デフォルトの名無しさん mailto:sage [2006/02/18(土) 03:40:16 ]
- スレ違いかもしれないけど
２値マスク画像の輪郭パスを取り出すアルゴリズムってなんかある？
単一オブジェクトの輪郭追跡だけじゃなくて、複数オブジェクトあり、中抜きありな
マスク画像をオブジェクト単位の輪郭パスに分離するようなもの
- 110 名前：デフォルトの名無しさん mailto:sage [2006/02/18(土) 09:17:03 ]
- つ[画像処理スレ]
- 111 名前：ハーピィ mailto:sage [2006/02/24(金) 11:49:52 ]
- E・∇・ﾖﾉｼ ＜111ｹﾞｯﾄ♫
- 112 名前：デフォルトの名無しさん [2006/03/14(火) 02:36:16 ]
- Linux板から流れてきますた。
これまでFortranメインで数値計算(差分法、積分方程式法)をやってきたのですが、
C++で新しくコードを書くことになって、大規模な計算の核部分は既存のFortranコードを
呼ぶことにしているのですが、ちょっとした計算や行列、ベクタの準備や修正などでは、
例えばMATLABライクな使い方ができれば良いなと思っています。
とりあえずgmm++やliboctaveを試してみたのですが、他にもなんか使いやすい
C++の基本行列計算や数学関数の簡易ライブラリがあったら教えてくらさい。
- 113 名前：デフォルトの名無しさん mailto:sage [2006/03/15(水) 21:44:30 ]
- >>112
使ったことないのでアレなんだけどgsl
www.gnu.org/software/gsl/
のラップで
GSLwrap
gslwrap.sourceforge.net/
GSL--
cholm.home.cern.ch/cholm/misc/gslmm/
とかどうでしょ？
- 114 名前：112 mailto:sage [2006/03/18(土) 01:51:40 ]
- うわ、これいい感じですね。使いやすそうです。
試しみます!良い情報ありがとうございますたー。
- 115 名前：http://www.vector.co.jp/soft/win95/util/se072729.html mailto:http://msdn2.microsoft.com/ja-jp/library/h2k70f3s.aspx [2006/03/18(土) 20:46:14 ]
- TextSS のWindowsXP(Professional)64bit化おながいします
もしくは64bitにネイティブ対応したテキスト置換ソフトありますか？
- 116 名前：デフォルトの名無しさん [2006/05/12(金) 14:43:10 ]
- FFTで高速なサブルーチンはありませんか？
- 117 名前：デフォルトの名無しさん mailto:sage [2006/05/12(金) 16:23:38 ]
- fftwとか
- 118 名前：デフォルトの名無しさん [2006/05/22(月) 16:30:47 ]
- 幾つか座標を渡すと、
曲率半径を計算してくれる
ライブラリありますか？
- 119 名前：デフォルトの名無しさん [2006/05/22(月) 19:50:13 ]
- 数値計算プログラムって一見動いてエラーはでないけど、
微妙に答えが間違ってたりしてデバッグがむずいなーなどと思うことがある。
- 120 名前：デフォルトの名無しさん [2006/05/22(月) 20:27:45 ]
- >>119
つ【ﾃｽﾄ問題】
- 121 名前：デフォルトの名無しさん mailto:sage [2006/05/23(火) 14:51:17 ]
- >118
三点なら簡単に求まるけど、それ以上の点から導きたいわけ？
- 122 名前：118 mailto:sage [2006/05/23(火) 15:11:18 ]
- 三点の場合教えて下さいでつ。
（後で使うときには拡張して複数点にするかもしないけど）
- 123 名前：デフォルトの名無しさん mailto:sage [2006/05/23(火) 15:33:32 ]
- ３点から同じ距離の点を方程式たてて解くか
２点づつの垂線の交点座標から求めるか
そして、その方法は複数の点には簡単には拡張出来ない。
「円の方程式と最小二乗法 」
「最小２乗法による円弧推定」をキーワードで検索すれば出る
- 124 名前：デフォルトの名無しさん mailto:sage [2006/05/24(水) 22:41:23 ]
- Athlon64用のオープンソースな多倍長精度ライブラリってないか？
- 125 名前：デフォルトの名無しさん mailto:sage [2006/05/25(木) 00:40:59 ]
- Athlon64専用ってわけじゃないけどgmpじゃだめ？
- 126 名前：デフォルトの名無しさん mailto:sage [2006/05/25(木) 14:58:12 ]
- >>125
サンクス。試してみる。
- 127 名前：デフォルトの名無しさん mailto:sage [2006/05/30(火) 22:32:32 ]
- >>116
GPUFFTW
ttp://gamma.cs.unc.edu/GPUFFTW/
これは単精度なのかな
しかし速そうだ
- 128 名前：デフォルトの名無しさん [2006/06/13(火) 17:10:39 ]
- fenv.hまわりの関数の使いかたが分からないので質問させてください．
環境はLinux 2.6.11，gcc 3.3.5，glibc 2.3.2です．
NaNの発生を検出しようと思いまして以下のようなコードを書きました．
例外が発生したところでfpe_handlerを1回呼んで復帰してほしいのですが，
以下のコードではfpe_handlerが何度も呼ばれるようでhogeが無数に表示されます．
例外フラグをクリアする必要があるのかなと思いまして(A)を追加しましたが動作は変わりません．
恐らく根本的に使い方が間違っているのだと思いますが，この辺の資料が少なく分からないでいます．
参考になるコードなどどなたか御存知ないでしょうか？
#include <stdio.h>
#include <signal.h>
#define __USE_GNU/* required for feenableexcept */
#include <fenv.h>
/* void print_flags () {printf ("fetestexcept: %x\n", fetestexcept(FE_ALL_EXCEPT));} */
void fpe_handler (int signum) {
/* feclearexcept (FE_ALL_EXCEPT); /\* (A) *\/ */
printf ("hoge\n");
}
int main () {
double d0 = 0; double d1 = 0; double d2 = 0;
signal (SIGFPE, &fpe_handler);
feenableexcept (FE_INVALID | FE_DIVBYZERO | FE_OVERFLOW); /* 浮動小数点例外トラップするようにする */
d2 = d0 / d1;
return 0;
}
- 129 名前：デフォルトの名無しさん mailto:sage [2006/06/13(火) 18:58:20 ]
- >>128
そのまま復帰すると、同じ演算を行って同じ例外が発生します。
C ならばエラー処理を行いたいあたりで適当に setjmp しておいて、
ハンドラ内から longjmp で脱出するのが定石だと思います。
- 130 名前：128 mailto:sage [2006/06/13(火) 19:07:13 ]
- >>129
>そのまま復帰すると、同じ演算を行って同じ例外が発生します。
そうか．なるほど
ありがとうございます．
- 131 名前：デフォルトの名無しさん [2006/07/20(木) 02:16:23 ]
- 上げますよー
- 132 名前：デフォルトの名無しさん [2006/10/02(月) 19:25:51 ]
- age
- 133 名前：デフォルトの名無しさん [2006/10/20(金) 01:29:32 ]
- age
- 134 名前：デフォルトの名無しさん mailto:sage [2006/12/04(月) 10:57:53 ]
- 4次ルンゲ＝クッタ法といのがありますが、1次や2次、3次のルンゲ＝クッタ法というものも存在するのでしょうか？
もしあるのであれば、どのように計算するのでしょうか？
- 135 名前：デフォルトの名無しさん [2006/12/07(木) 00:49:40 ]
- ttp://grape.c.u-tokyo.ac.jp/~makino/kougi/system_suuri4_1998/all/all.html
ここのルンゲクッタの章を読んで
- 136 名前：デフォルトの名無しさん [2007/02/12(月) 13:13:43 ]
- ベクトル解析に向いている言語ってなに？
- 137 名前：デフォルトの名無しさん [2007/02/12(月) 21:35:25 ]
- >>136 ベクトル解析に向いている言語ってなに？
試行錯誤しながら数値実験していくのならば下が使える。
www.nasuinfo.or.jp/FreeSpace/kenji/sf/fastTour/sfFastTour.htm
評価版でも 100 次位まで十分に動く。気に入ったら買ってやってくれ。5 千
円だと言っているから電卓代だと思える
- 138 名前：デフォルトの名無しさん mailto:sage [2007/02/12(月) 21:49:17 ]
- MATLAB買え。
- 139 名前：デフォルトの名無しさん mailto:sage [2007/02/18(日) 00:31:33 ]
- >>112
tharkad.pa.msu.edu/~asteiner/hal/html/index.html
ってのもある。
あと、
blitz++, free pooma, Matrix Template Library, なんてのもあるよ。
Boostにも行列系ライブラリはある
って11ヶ月前のものにレスしてる^^;
- 140 名前：デフォルトの名無しさん [2007/03/07(水) 18:05:55 ]
- 拡散方程式で、拡散係数がテンソルなのですが、対角項、
\frac{\partial^2 f}{\partial x \partial y}はどのように計算すればよいのでしょうか？
- 141 名前：デフォルトの名無しさん mailto:sage [2007/03/07(水) 21:59:27 ]
- Mathematica買え
マジで感動
値段分の価値はある
- 142 名前：デフォルトの名無しさん [2007/04/17(火) 18:39:14 ]
- 速度ではC++はFortranに勝てないのか？
- 143 名前：デフォルトの名無しさん mailto:sage [2007/04/17(火) 21:00:55 ]
- 処理系や内容による。
- 144 名前：デフォルトの名無しさん [2007/06/04(月) 08:09:42 ]
- 科学計算向けだと、どのCPU（Opteron or Core 2 Duo）とどのコンパイラの組み合わせがお薦めですか？
- 145 名前：デフォルトの名無しさん [2007/06/04(月) 21:23:14 ]
- AlphaとCompaq Fortran
- 146 名前：デフォルトの名無しさん mailto:sage [2007/06/05(火) 03:55:46 ]
- OpteronかCore 2 Duoで検討しています
- 147 名前：デフォルトの名無しさん mailto:sage [2007/06/06(水) 01:20:29 ]
- OpteronならBorland C++、Core2 DuoならLSI C-86
PathScaleとかIntelは不許可
- 148 名前：デフォルトの名無しさん mailto:sage [2007/06/07(木) 01:54:49 ]
- C++Builderで検討してみます。
- 149 名前：デフォルトの名無しさん [2007/06/29(金) 16:19:54 ]
- 40000*40000の行列の逆行列を求めたいのですが、どのライブラリがお薦めですか？
- 150 名前：デフォルトの名無しさん [2007/07/08(日) 13:55:29 ]
- intelのMKLのScaLAPACKはCでも使えますか？
- 151 名前：デフォルトの名無しさん [2007/07/09(月) 14:14:51 ]
- 何か面白いシミュレーションを見たい！
- 152 名前：デフォルトの名無しさん [2007/07/09(月) 14:19:01 ]
- 見せてくれよ！
- 153 名前：デフォルトの名無しさん [2007/07/09(月) 14:53:06 ]
- >>150
使える
- 154 名前：デフォルトの名無しさん [2007/07/16(月) 10:12:19 ]
- rubyで科学計算したら研究室追い出されますか？
- 155 名前：デフォルトの名無しさん mailto:sage [2007/07/16(月) 13:31:25 ]
- >>154
特に計算速度が要求される場合、スクリプトはいけない。
- 156 名前：デフォルトの名無しさん [2007/07/16(月) 15:38:07 ]
- >>154
つ「とんでる力学」牧野 淳一郎
雑誌「パリティ」での連載をまとめたもの。rubyを使った物理シミュレーションを行っている。
「地球流体電脳倶楽部」で行列のライブラリとか結構開発されているよ。あとは画像処理するための前処理とか。
さすがに、rubyで直接シミュするのは例外だとは思うけど。
- 157 名前：デフォルトの名無しさん [2007/07/16(月) 20:21:19 ]
- c-au.2ch.net/test/-/mass/1183191945/i#b
c-au.2ch.net/test/-/soc/1152718893/i
- 158 名前：デフォルトの名無しさん mailto:sage [2007/09/09(日) 11:34:26 ]
- MATLABで科学技術計算やってますｽｲﾏｾﾝ
- 159 名前：デフォルトの名無しさん mailto:sage [2007/09/09(日) 12:01:12 ]
- この金持ちめ
- 160 名前：デフォルトの名無しさん mailto:sage [2007/09/10(月) 00:38:43 ]
- python + f2py + Fortran でやってる
MATLABは後輩が買ってたけど使ってないな…
- 161 名前：デフォルトの名無しさん mailto:sage [2007/09/16(日) 19:49:46 ]
- 金持ちならマテマテカじゃん！
- 162 名前：デフォルトの名無しさん [2007/11/07(水) 12:01:41 ]
- 定数での割り算は最適化で定数のかけ算に変換されるものなのでしょうか？
- 163 名前：デフォルトの名無しさん mailto:sage [2007/11/08(木) 00:18:38 ]
- -S オプションつけてコンパイルするんだ
- 164 名前：デフォルトの名無しさん [2008/02/24(日) 17:31:30 ]
- ＞＞１６３
どのコンパイラの話？
- 165 名前：デフォルトの名無しさん [2008/05/13(火) 19:52:53 ]
- MATLAB+Mathematicaを学部レベルで導入出来ている大学は勝ち組。
- 166 名前：デフォルトの名無しさん mailto:sage [2008/05/14(水) 23:31:41 ]
- 早稲田が勝ち組だなんて言うつもりか？
|
 |
