06/05/12 14:43:10
FFTで高速なサブルーチンはありませんか？
117:デフォルトの名無しさん
06/05/12 16:23:38
fftwとか
118:デフォルトの名無しさん
06/05/22 16:30:47
幾つか座標を渡すと、
曲率半径を計算してくれる
ライブラリありますか？
119:デフォルトの名無しさん
06/05/22 19:50:13
数値計算プログラムって一見動いてエラーはでないけど、
微妙に答えが間違ってたりしてデバッグがむずいなーなどと思うことがある。
120:デフォルトの名無しさん
06/05/22 20:27:45
>>119
つ【ﾃｽﾄ問題】
121:デフォルトの名無しさん
06/05/23 14:51:17
>118
三点なら簡単に求まるけど、それ以上の点から導きたいわけ？
122:118
06/05/23 15:11:18
三点の場合教えて下さいでつ。
（後で使うときには拡張して複数点にするかもしないけど）
123:デフォルトの名無しさん
06/05/23 15:33:32
３点から同じ距離の点を方程式たてて解くか
２点づつの垂線の交点座標から求めるか
そして、その方法は複数の点には簡単には拡張出来ない。
「円の方程式と最小二乗法 」
「最小２乗法による円弧推定」をキーワードで検索すれば出る
124:デフォルトの名無しさん
06/05/24 22:41:23
Athlon64用のオープンソースな多倍長精度ライブラリってないか？
125:デフォルトの名無しさん
06/05/25 00:40:59
Athlon64専用ってわけじゃないけどgmpじゃだめ？
126:デフォルトの名無しさん
06/05/25 14:58:12
>>125
サンクス。試してみる。
127:デフォルトの名無しさん
06/05/30 22:32:32
>>116
GPUFFTW
URLﾘﾝｸ(gamma.cs.unc.edu)
これは単精度なのかな
しかし速そうだ
128:デフォルトの名無しさん
06/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:デフォルトの名無しさん
06/06/13 18:58:20
>>128
そのまま復帰すると、同じ演算を行って同じ例外が発生します。
C ならばエラー処理を行いたいあたりで適当に setjmp しておいて、
ハンドラ内から longjmp で脱出するのが定石だと思います。
130:128
06/06/13 19:07:13
>>129
>そのまま復帰すると、同じ演算を行って同じ例外が発生します。
そうか．なるほど
ありがとうございます．
131:デフォルトの名無しさん
06/07/20 02:16:23
上げますよー
132:デフォルトの名無しさん
06/10/02 19:25:51
age
133:デフォルトの名無しさん
06/10/20 01:29:32
age
134:デフォルトの名無しさん
06/12/04 10:57:53
4次ルンゲ＝クッタ法といのがありますが、1次や2次、3次のルンゲ＝クッタ法というものも存在するのでしょうか？
もしあるのであれば、どのように計算するのでしょうか？
135:デフォルトの名無しさん
06/12/07 00:49:40
URLﾘﾝｸ(grape.c.u-tokyo.ac.jp)
ここのルンゲクッタの章を読んで
136:デフォルトの名無しさん
07/02/12 13:13:43
ベクトル解析に向いている言語ってなに？
137:デフォルトの名無しさん
07/02/12 21:35:25
>>136 ベクトル解析に向いている言語ってなに？
試行錯誤しながら数値実験していくのならば下が使える。
URLﾘﾝｸ(www.nasuinfo.or.jp)
評価版でも 100 次位まで十分に動く。気に入ったら買ってやってくれ。5 千
円だと言っているから電卓代だと思える
138:デフォルトの名無しさん
07/02/12 21:49:17
MATLAB買え。
139:デフォルトの名無しさん
07/02/18 00:31:33
>>112
URLﾘﾝｸ(tharkad.pa.msu.edu)
ってのもある。
あと、
blitz++, free pooma, Matrix Template Library, なんてのもあるよ。
Boostにも行列系ライブラリはある
って11ヶ月前のものにレスしてる^^;
140:デフォルトの名無しさん
07/03/07 18:05:55
拡散方程式で、拡散係数がテンソルなのですが、対角項、
\frac{\partial^2 f}{\partial x \partial y}はどのように計算すればよいのでしょうか？
141:デフォルトの名無しさん
07/03/07 21:59:27
Mathematica買え
マジで感動
値段分の価値はある
142:デフォルトの名無しさん
07/04/17 18:39:14
速度ではC++はFortranに勝てないのか？
143:デフォルトの名無しさん
07/04/17 21:00:55
処理系や内容による。
144:デフォルトの名無しさん
07/06/04 08:09:42
科学計算向けだと、どのCPU（Opteron or Core 2 Duo）とどのコンパイラの組み合わせがお薦めですか？
145:デフォルトの名無しさん
07/06/04 21:23:14
AlphaとCompaq Fortran
146:デフォルトの名無しさん
07/06/05 03:55:46
OpteronかCore 2 Duoで検討しています
147:デフォルトの名無しさん
07/06/06 01:20:29
OpteronならBorland C++、Core2 DuoならLSI C-86
PathScaleとかIntelは不許可
148:デフォルトの名無しさん
07/06/07 01:54:49
C++Builderで検討してみます。
149:デフォルトの名無しさん
07/06/29 16:19:54
40000*40000の行列の逆行列を求めたいのですが、どのライブラリがお薦めですか？
150:デフォルトの名無しさん
07/07/08 13:55:29
intelのMKLのScaLAPACKはCでも使えますか？
151:デフォルトの名無しさん
07/07/09 14:14:51
何か面白いシミュレーションを見たい！
152:デフォルトの名無しさん
07/07/09 14:19:01
見せてくれよ！
153:デフォルトの名無しさん
07/07/09 14:53:06
>>150
使える
154:デフォルトの名無しさん
07/07/16 10:12:19
rubyで科学計算したら研究室追い出されますか？
155:デフォルトの名無しさん
07/07/16 13:31:25
>>154
特に計算速度が要求される場合、スクリプトはいけない。
156:デフォルトの名無しさん
07/07/16 15:38:07
>>154
つ「とんでる力学」牧野 淳一郎
雑誌「パリティ」での連載をまとめたもの。rubyを使った物理シミュレーションを行っている。
「地球流体電脳倶楽部」で行列のライブラリとか結構開発されているよ。あとは画像処理するための前処理とか。
さすがに、rubyで直接シミュするのは例外だとは思うけど。
157:デフォルトの名無しさん
07/07/16 20:21:19
URLﾘﾝｸ(c-au.2ch.net)
URLﾘﾝｸ(c-au.2ch.net)
158:デフォルトの名無しさん
07/09/09 11:34:26
MATLABで科学技術計算やってますｽｲﾏｾﾝ
159:デフォルトの名無しさん
07/09/09 12:01:12
この金持ちめ
160:デフォルトの名無しさん
07/09/10 00:38:43
python + f2py + Fortran でやってる
MATLABは後輩が買ってたけど使ってないな…
161:デフォルトの名無しさん
07/09/16 19:49:46
金持ちならマテマテカじゃん！
162:デフォルトの名無しさん
07/11/07 12:01:41
定数での割り算は最適化で定数のかけ算に変換されるものなのでしょうか？
163:デフォルトの名無しさん
07/11/08 00:18:38
-S オプションつけてコンパイルするんだ
164:デフォルトの名無しさん
08/02/24 17:31:30
＞＞１６３
どのコンパイラの話？
165:デフォルトの名無しさん
08/05/13 19:52:53
MATLAB+Mathematicaを学部レベルで導入出来ている大学は勝ち組。
166:デフォルトの名無しさん
08/05/14 23:31:41
早稲田が勝ち組だなんて言うつもりか？