- 1 名前：デフォルトの名無しさん [04/12/31 00:46:46]
- はっきり言ってマイナーな分野。
やってる人も少ないだろうし市場も小さいと思う。
技術の話でも業界の動向でも何でもいいから
どんどん書き込んでくれ
- 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 ]
- 早稲田が勝ち組だなんて言うつもりか？
|
 |
