- 548 名前:デフォルトの名無しさん [2006/06/29(木) 01:10:31 ]
- >>546
ググルと色々サンプルプログラムが出てくるので参照してみたらいかがかな。 ピボッティングだったら教科書の類にもよく出てきているので、他人のプログラムを 読んでみるといいと思うんこ。 >>547 このループでc(i,j)に次々足してゆくわけだが、一番最初に0.0になっていないと困る。 DO 30 k = 1, m c(i, j) = c(i, j) + a(i, k) * b(k, j) 30 CONTINUE Fortranでは変数の初期値は不定なのでごみが入っていることがある。 最近のコンパイラ(というか正確にはローダーだが)変数をゼロに初期化してくれることが 多いが、これを期待してはいかんざき。変数はみんな初期化する癖を付けるのがよろしいアル カイダ。
|

|