くだすれFORTRAN（超初心者用）

1 名前：デフォルトの名無しさん [2006/01/24(火) 09:48:23 ] このスレッドは、他のスレッドでは書き込めない超低レベル、  もしくは質問者自身何が何だが分からない質問を勇気を持って書き込むスレッドです。  FORTRAN使いが優しくコメントを返しますが、  お礼はFORTRANの布教と初心者の救済をお願いします。  467 名前：デフォルトの名無しさん [2006/06/13(火) 12:44:42 ] >>466 ういうい。 CHARACTER DUMMY1*1, DUMMY2*1 READ(*, '(I2, A1, I2, A1, I2)') IH, DUMMY1, IM, DUMMY2, IS IF (DUMMY1 .NE. ':') STOP 'ERROR' ・・・・・・・・ IF (IH .GT. 24 .OR. IH .LT. 0) STOP 'ERROR' 云々 まあ、いきなりストップせず、WRITE 文で入力値を出力して、 どこがおかしいのか書いたほうが親切ではある。 468 名前：デフォルトの名無しさん [2006/06/13(火) 13:25:31 ] alpha 10.0 beta 0.00 .... gamma 4E-3 ↑という感じのファイルを読み込みたいのですけど、 一行に文字と色々なフォーマットの数字が並んでるんですけど、 できますか？ 文字は20文字固定で、余りは半角で埋めようと思っています。 469 名前：デフォルトの名無しさん [2006/06/13(火) 16:37:03 ] >>468 できる。 もし文字と数値の間に必ず空白が空いているなら簡単。 CHARACTER NAME*20, TEXT*100 READ(9, '(a)') TEXT IPOS = INDEX(TEXT,' ') NAME = TEXT(1:IPOS - 1) READ(TEXT(IPOS + 1:IPOS + 20), '(F20.0)') VALUE こんな感じで、一度１行全体を文字列として読み、名前と値の境目を見つけたら そこより前は文字変数に突っ込み、境目から後ろは文字列を内部ファイルとして もちいて、数値として数値変数に読み出せばいい。 （ここでは２０桁の数値を仮定した。F20.0にしたので小数点がないと とんでもない値になるが、この辺もINDEXで尻を見つけるとか、 FORTRAN90でいいならTRIMや何やらを使えば回避できる。） 470 名前：擬似乱数のチェック [2006/06/13(火) 16:55:21 ] [0,1]の一様乱数の標準偏差は0.28...で、0.25ではありませんでした。 お手数をかけました。 471 名前：素人 [2006/06/13(火) 18:03:12 ] -1以上1以下の値ｘを入力してe^xの近似式を出力するプログラム。e^xの近似式はe^x=1+x+x^2/2!+x^3/3!+..+x^n/n!を利用。 (最後の項に絶対値が0.00001より大きい間、項を増やして足し続ける)また、適正なxの値が入力される限り、何度もこの計算ができるようにすること。って課題が出されたんですけどまったく分かりません。誰か助けてください！ って昨日打ったんですけどc++でお願いします！！ 472 名前：デフォルトの名無しさん [2006/06/13(火) 19:15:56 ] >>471 ちょｗｗｗ　おまｗｗｗｗ Ｆｏｒｔｒａｎスレに聞くなよｗ >>450を自分で直せｗｗｗｗ 473 名前：デフォルトの名無しさん [2006/06/13(火) 21:04:58 ] テキストファイルから読み込んだ値を変数に代入したいのですが、 テキストの値の書式が 1,22.489,445,67888,23 のようにカンマ区切りの場合、どのように読み込むのが良いでしょうか？ 474 名前：デフォルトの名無しさん [2006/06/13(火) 21:09:33 ] >>473 書式無しの自由書式。この場合コンマで区切って読んでくれる。 それで問題がある場合は、文字列として読んで、自分でコンマのところで切って 文字列を内部ファイルとして数値に変換する。 475 名前：473 [2006/06/13(火) 21:29:10 ] >>474 integer i,j,k,l,m real r open(10,file='hoge.txt') read(10,*)i,r,j,k,l,m のような感じで良いのでしょうか？？ 476 名前：474 mailto:sage [2006/06/13(火) 21:39:44 ] >>475 いいはず 477 名前：473 [2006/06/13(火) 22:00:59 ] >>476 ありがとう 478 名前：デフォルトの名無しさん [2006/06/14(水) 06:18:06 ] c言語にはfortranでいうcycle機能ってあるの？ 479 名前：デフォルトの名無しさん mailto:sage [2006/06/14(水) 06:43:18 ] continue; 480 名前：デフォルトの名無しさん [2006/06/14(水) 12:09:22 ] すいません。昨日からはじめた超初心者です。 初心者用のサイトを見ながらやってますが、 書いてある通りに行かないので質問させてください。 FTN95でコンパイルする際、 「ftn95.exe プログラムのファイル名.for95」と入れれば.exeファイルができる と説明してあるのですが、 上のように入力しても、「.exe」ではなく「.obj」ファイルになってしまいます。 どうすれば.exeになるのでしょうか？ 481 名前：デフォルトの名無しさん [2006/06/14(水) 12:12:35 ] >>480 コンパイル時に/linkオプション。つまり ftn95.exe プログラムのファイル名.for95 /link 482 名前：デフォルトの名無しさん [2006/06/14(水) 12:20:59 ] >>481 ありがとうございます。 出来ました。 483 名前：452 mailto:sage [2006/06/14(水) 21:43:11 ] >>467 ありがとうございました。 484 名前：デフォルトの名無しさん [2006/06/15(木) 21:53:09 ] real X read(*,*)X として、小数点がつかない 5 などを読み込んだとき、Xにとんでもない値が入ってしまうのを 回避する方法はあるでしょうか？ 485 名前：デフォルトの名無しさん [2006/06/15(木) 23:11:09 ] 行き詰ってしまっているので、失礼します 1行に (整数):(整数)(3文字の文字列)(8ケタの整数)(17文字の文字列)(空白1)(14文字の文字列) というデータが並び、これが何百行と続くデータを読みこみ、 (8ケタの整数)というもののみを、他ファイルに出力したいと思います。 (整数)というのはケタ数が決まっていません。これがエラーでコンパイルできません。ご指導お願いします program read implicit none integer::i,j,k,l open(unit=11,file='input.dat') open(unit=22,file='output.txt') 　do l=1,100 　read(11,'(I8,1X,I15,3X,I8,32X)i,j,k 　write(22,*)k 　end do close(11) close(22) end program read 486 名前：デフォルトの名無しさん [2006/06/15(木) 23:13:19 ] l=1,100　となっているのは、100行読み取りたいとき、ということにしてます 修正です 　read(11,'(I8,1X,I15,3X,I8,32X)')i,j,k でお願いします 487 名前：デフォルトの名無しさん [2006/06/15(木) 23:18:46 ] また失礼します。コンパイルはできるのですが、実行ができないのでした・・ fmt: read unexpected character apparent state: unit 11 named input.dat last format: (I8,1X,I15,3X,I8,32X) lately reading sequential formatted external IO aborted (core dumped) と表示されます。環境はcygwinで、g77でコンパイルをしています。ぐだぐだと失礼しました 488 名前：デフォルトの名無しさん [2006/06/16(金) 00:06:33 ] 事故解決しました。すみません 489 名前：デフォルトの名無しさん [2006/06/16(金) 00:30:48 ] >>484おながいいたします あと、"を文字として書くときってどうすれば良いでしょうか 490 名前：デフォルトの名無しさん mailto:sage [2006/06/16(金) 00:48:22 ] >>484 普通の処理系なら read(*,*) X に対して 5 と入力した場合、小数部を省略したとみなして X は 5. となるはずだが？ >>489 >"を文字として書くとき print *, '"' print *, """" print *, 1H" 491 名前：デフォルトの名無しさん [2006/06/16(金) 00:52:39 ] 4:345:23:346543 346:344:4:455 のような、ケタ数が一致しない複数の数値が:で区切られた値を読み取るにはどうしたら良いのでしょうか 492 名前：デフォルトの名無しさん mailto:sage [2006/06/16(金) 01:59:52 ] >>491 各行に含まれる数値の数が一定なら、以下のように一旦文字列として 読み込んで加工し、内部ファイルを使えば良い。 program hoge implicit none integer i, k,l,m,n, ios character*80 str open(11,file='data.dat',status='OLD') 100 continue read(11, fmt='(A)', iostat=ios) str if (ios .ne. 0) goto 999 do 200 i=1, len(str) if (str(i:i) .eq. ':') str(i:i) = ',' 200 continue read(str,*) k,l,m,n write(*,*) k,l,m,n goto 100 999 continue close(11) stop end 493 名前：デフォルトの名無しさん [2006/06/16(金) 02:03:36 ] >>492 ありがとうございます。内部ファイル・・？ これを参考にして頑張ります 494 名前：デフォルトの名無しさん [2006/06/16(金) 13:46:56 ] >>493 >read(*,*)X >として、小数点がつかない >5 >などを読み込んだとき、Xにとんでもない値が入ってしまうのを >回避する方法はあるでしょうか？ 自由書式の場合は問題ない。 FORMATを与えていた場合に問題が生ずることがある。 小数点があれば与えたFORMATよりも、入力の小数点が優先されるが、 小数点がないと与えたFORMATの桁数にしたがって小数点位置が決められる。 >>491の安易な解決法は：を、に置換して自由書式で読み出すことだｗ >>491の問題は>>492の方法で解決できるが、Fortran90以降を使うと 文字列関数が豊富な分やりやすい。ただし77の範囲でも実現は可能。 だらだらと書くもよし、：までを切り出した後、文字変数に右詰に して返す関数を作るもよし。 内部ファイルと外部ファイルの違う点は、ファイルの読み込み位置が 毎回先頭に戻ってしまう事で、これがちょっとうざい。 外部ファイルとなんでも一緒のつもりで行って昔ちょっとはまったｗ 495 名前：デフォルトの名無しさん [2006/06/16(金) 23:07:38 ] 文字列の配列ってできますか？ character*50 c(10) みたいな感じで。 496 名前：495 [2006/06/17(土) 00:01:25 ] あっ、できました。 また質問いたします。 integer i character*50 c do i=1,4 write(c,'(i4)')i end みたいな感じにして、 cに「1234」のような文字列が入るようにしたいのですが、 上手くいきません。 どうすればよろしいのでしょうか？ 497 名前：デフォルトの名無しさん mailto:sage [2006/06/17(土) 12:03:21 ] >>496 integer i character*50 c write(c,'(4(I1))') (i, i=1,4) write(*,*) c end 498 名前：デフォルトの名無しさん [2006/06/17(土) 12:30:04 ] >>496 別解 DO i = 1, 4 WRITE(c(i:i), '(i1)') i END DO 499 名前：495 [2006/06/17(土) 17:53:25 ] >>497 >>498 どうもサンクス もう一個、 読み込んだファイルの全行数を知ることはできませんか？ 500 名前：デフォルトの名無しさん [2006/06/17(土) 18:02:08 ] >>499 多分1回１行づつ数えるしかないと思うのだが。ほかに方法あるのかな？ iline = 0 DO READ(9, iostat = io, *) IF (io < 0) EXIT iline = iline + 1 END DO こうかな？ 巻き戻し REWIND(9) 501 名前：デフォルトの名無しさん mailto:sage [2006/06/17(土) 18:50:57 ] >>500 あ、できました。どうもありがとう READ(9,*, iostat = io) ですね 502 名前：デフォルトの名無しさん mailto:sage [2006/06/20(火) 23:21:38 ] 今一番新しいfortranの規格ってfortran2003ってやつですか？ これ普及すんのかなぁ 503 名前：デフォルトの名無しさん [2006/06/21(水) 00:32:08 ] >>502 する！ まだ実装したコンパイラないけどなｗ 504 名前：デフォルトの名無しさん mailto:sage [2006/06/21(水) 09:50:12 ] ﾏｼﾞｽｶ まだ77でさえしぶとく残ってるのにね… 505 名前：デフォルトの名無しさん mailto:sage [2006/06/21(水) 11:51:22 ] >>504 ６６だってまだまだ健在だぜ！ ９０で書いてるって言ってる奴も、珍妙な７７や６６の入り混じったハイブリットだったりするし。 506 名前：デフォルトの名無しさん mailto:sage [2006/06/21(水) 12:00:43 ] そこまでいくとｶｺｲｲなｗ 507 名前：デフォルトの名無しさん mailto:sage [2006/06/21(水) 18:54:00 ] fortran90で ｎ個の配列を準備してＲＡＮＤＯＭ＿NUMBER関数を用いて ランダムな値を配列に代入する。 ばらばらの順になっている配列の実数を、小さな順に並び替えるサブルーチンを作成し、 そのサブルーチンをcallする。 結果を出力する。 でばらばらの順になっている配列の実数を、小さな順に並び替えるサブルーチンを作成し、 そのサブルーチンをcallする。 というところがさっぱりわかりません。 どうすればいいでしょうか？？ 508 名前：デフォルトの名無しさん mailto:sage [2006/06/21(水) 21:06:27 ] >>507 とりあえず、「ソート アルゴリズム」でググれば いろいろとソート（並び替え）の方法が出てくるよ。 あとはそのうち簡単なものを選んで実装するのみ。 509 名前：デフォルトの名無しさん mailto:sage [2006/06/22(木) 09:08:07 ] >>507 作ってみた program rndm real,dimension(:),allocatable :: x; integer :: i,n,stat read(*,*) n allocate(x(n),stat=stat) if(stat/=0) stop "could not allocate x" do i=1,n call random_number(x(i)) enddo print *,"before sort"; call printx(x) call srt(x) print *,"after sort"; call printx(x) deallocate(x) contains subroutine printx(x) real,dimension(:) :: x; integer :: i do i=1,size(x); print *,x(i); enddo end subroutine printx recursive subroutine srt(x) integer :: i,iflag; real :: tmp real,dimension(:) :: x iflag=0 do i=1,size(x)-1 if(x(i)>x(i+1)) then iflag=1; tmp=x(i); x(i)=x(i+1); x(i+1)=tmp endif enddo if(iflag/=0) call srt(x) end subroutine srt end program rndm 510 名前：デフォルトの名無しさん [2006/06/22(木) 13:17:26 ] >>509 おまいFortranプログラマーでないだろうｗ 行内複数命令、暗黙の型からの逸脱、再帰関数・・・・・ 他の言語だと推奨されるところだが、Fortranでは嫌われるものばかりｗ 育ちが分かるぜｗ　うっしっしｗ 511 名前：デフォルトの名無しさん mailto:sage [2006/06/22(木) 13:34:23 ] >>510 一応77育ちで最近90に手を出したところ。 行内複数命令はやりたくなかったけど、「改行が多過ぎます」の 書き込みエラーが出たのを回避するため。行数を減らす時に 「implicit none」を消してしまった。 >>507がfortran90と書いていたので、90で取り込まれた仕様を 使ってみた。call random_number(x(i))のところはdoループに しないでcall random_number(x)でも良かったみたい。 512 名前：509,511 mailto:sage [2006/06/23(金) 09:49:52 ] 行内複数命令を使わずに出来るだけ短くしてみた program rndm real,dimension(:),allocatable :: x read(*,*) n allocate(x(n)) call random_number(x) print *,"before sort" print *,x call srt(x) print *,"after sort" print *,x deallocate(x) contains recursive subroutine srt(x) real,dimension(:) :: x iflag=0 do i=1,size(x)-1 if(x(i)>x(i+1)) then iflag=1 x(i:i+1)=(/x(i+1),x(i)/) endif enddo if(iflag/=0) call srt(x) end subroutine srt end program rndm 513 名前：デフォルトの名無しさん mailto:sage [2006/06/23(金) 15:25:04 ] >>512 x(i:i+1)=(/x(i+1),x(i)/) うわ！　こういうswap tips があったのか。 514 名前：デフォルトの名無しさん [2006/06/23(金) 16:04:06 ] >>513 隣接で良ければ x(i:i+1) = x(i+1:i:-1) という手もあるな． 515 名前：デフォルトの名無しさん mailto:sage [2006/06/23(金) 16:59:24 ] >>514 Fortran 90 以降の並列型配列処理の強力さが全開！！ 516 名前：デフォルトの名無しさん [2006/06/23(金) 17:15:11 ] >>515 強力なんだが、暗黙の内にテンポラリ配列に展開するので、 超巨大配列に適用すると、すゲー時間がかかったりメモリーを食ったりする。 サブルーチンに渡すときに７７の時のつもりで先頭要素のアドレスだけ渡している気になると 実はテンポラリにコピー、引き渡し、受け取り、テンポラリから元の配列へコピーとやってたりするｗ 517 名前：デフォルトの名無しさん [2006/06/26(月) 18:15:09 ] すいません、くそ初心者なんですが文の途中に改行するときに 入れる記号を教えてください！ 518 名前：デフォルトの名無しさん mailto:sage [2006/06/26(月) 18:45:14 ] >>517 すいません、くそ初心者なんですが文の途中に改行するときに & 入れる記号を教えてください！ 519 名前：デフォルトの名無しさん [2006/06/26(月) 20:03:11 ] >>517 Fortran90なら>>518にあるように　＆　を尻に付ける． FORTRAN７７なら　６　カラム目に何か文字を入れる．＋とか数字がよく使われるかな． 520 名前：デフォルトの名無しさん mailto:sage [2006/06/26(月) 20:45:54 ] >519 > FORTRAN７７なら　６　カラム目に何か文字を入れる．＋とか数字がよく使われるかな． カーニハン、プローガー、プログラム書法、第2版、木村泉訳、共立出版、p.21では、 継続行を表すのに$が使われています。 その理由が次のとおり書かれています。 （ここで継続行をあらわすのに「$」を使っているのは、「$」がFortranの規格上 公式に認められている文字のうち、構文上特別の意味を持たないただ一つの ものであることによる。こうすれば混乱が生ずる機会を最小にすることができる とともに、欄をまちがったときにそれに気づく可能性が大きい。） 521 名前：デフォルトの名無しさん mailto:sage [2006/06/26(月) 21:14:41 ] なるほど… 90でも$にしたほうがいいの 522 名前：デフォルトの名無しさん mailto:sage [2006/06/26(月) 23:20:08 ] 主プログラムの中で配列の番号指定して それをfunctionに渡してやりたいんだけど・・・うまく渡せないんだけど･･･ どうしたらいいの？（＞＜） program UNKO real chinko(10)　←ファイルから読み込ませる i=5 なんちゃらかんちゃら function sex(ぱらめーた) sex=chinko(i)*manko chinko(i)がゼロになるｗ 523 名前：デフォルトの名無しさん mailto:sage [2006/06/26(月) 23:21:53 ] >>520 そういえばプログラム書法にはそうあったような気がするな。 漏れはＡＣＯＳっ子だったので＆を６カラム目に置かされていたｗｗｗｗ 数字を使う人たちは１，２，３，４・・と行が継続行の数だけ数字を増やして行く流派の人たちだった。 ＋は昔のＦＯＲＴＲＡＮのテキストに多かった気がする。 >>521 ９０では前の行の尻に＆だ。＆以外は許されない。 続いた行頭においてもよかったはずだが。 ９０で固定カラムというのならまぁ＄でもいいがｗｗｗｗ 524 名前：デフォルトの名無しさん [2006/06/26(月) 23:27:57 ] >>522 いくつかやり方がある。 い）必要な要素だけを渡す。 program UNKO real chinko(10)　←ファイルから読み込ませる anal = sex( chinko(5)、ぱらめーた　 ) function sex(chinpo, ぱらめーた) sex=chinpo*manko ろ）配列全部渡す。 program UNKO real chinko(10)　←ファイルから読み込ませる anal = sex( chinko、ぱらめーた　 ) ■７７式　甲 function sex(chinpo, ぱらめーた) real chinpo(10) i = 5 sex=chinpo(i)*manko ■７７式　乙 function sex(chinpo, ぱらめーた) real chinpo(*) i = 5 sex=chinpo(i)*manko ■９０式 function sex(chinpo, ぱらめーた) real, intent(in) :: chinpo(:) i = 5 sex=chinpo(i)*manko 525 名前：デフォルトの名無しさん mailto:sage [2006/06/26(月) 23:29:07 ] 下品すぎます 526 名前：デフォルトの名無しさん mailto:sage [2006/06/26(月) 23:45:31 ] 下品ですいません（＞＜） 一瞬で答えが出ましたねすごいですｗ 必要な要素だけを渡す場合にchinko(5)のところはchinko(i)でもいいんですか？ 具体的には主プログラム中の繰り返し回数iによってfunction中の値も変えたいんです（＞＜） 527 名前：デフォルトの名無しさん [2006/06/26(月) 23:54:08 ] >>526 OK DO impo = 1, 69 anal = sex( chinpo(i) , ....... ) ....... END DO 528 名前：５１８ [2006/06/27(火) 11:40:59 ] やっぱりわかりません泣　どこが違うのでしょうか？ 　　H1 = 0.0639 H2 = 0.0001 PAI = 3.14159 T = 0.0164 OMG = 2.0 * PAI / T D = 0.5 V = 20.0 WRITE(12,10) 10 FORMAT(' T X1(t) X2(T)') DO 100 N=0,300,1 X = FLOAT(N) / 200.0 Y1 = EXP(-H1 * OMG * X) * (D * COS(SQRT(1.0 - H1**2) * OMG * X) $ + (D * H1 + V / OMG) / SQRT(1.0 - H1**2) $ * SIN(SQRT(1.0 - H1**2) * OMG * X) * EXP(-H1 * OMG * X) Y2 = EXP(-H2 * OMG * X) * COS(SQRT(1.0 - H2**2) * OMG * X) WRITE(12,20) X,Y1,Y2 20 FORMAT(3F7.3) 100 CONTINUE 529 名前：デフォルトの名無しさん [2006/06/27(火) 13:59:36 ] >>528 Y1 = EXP(-H1 * OMG * X) * (D * COS(SQRT(1.0 - H1**2) * OMG * X) この行のＤの前の括弧　（　に対応する　閉じ括弧がないから叱られている． 継続行の問題ではない．もしくは，継続行にしているうちに壊してしまったのだろう． >>524　への追加だが，　 ＦＯＲＴＲＡＮ７７での配列の渡し方としては，配列の大きさを渡すのが一番一般的だと思う． 数学ライブラリのたぐいはみんなこうだろう． ■７７式　丙 anal = sex( n, chinko, ..... ) function sex( n, chinko, ..... ) real chinko(n) ....... 530 名前：529 mailto:sage [2006/06/27(火) 14:27:48 ] ついでながら， PAI = 3.14159 これだと，単精度の有効桁にちょっと足りない．６〜７桁はあるのでもう１，２数字つけておくと安心だ． まぁ，本来的には他の数字の有効桁にもよるのだが，常数だからギリギリまで正確に与えた方がいい． PAI = 4.0 * ATAN(1.0) 倍精度なら PAI = 4.0d0 * ATAN(1.0d0) がよく使われる．これだと，機械によらずに有効な桁まで数値が与えられる． 倍精度の時にπを与えるのがかったるいので便利である． 531 名前：デフォルトの名無しさん [2006/06/27(火) 23:30:54 ] すいません、教科書見ながらやってみたのですが全く手に負えません。 この課題お願いします。 ｍ行ｎ列の行列（実数）Ａ，Ｂをキーボードから読み込んで、 行列の積Ｃ＝ＡＢを計算し、結果を画面に表示させなさい。 序盤からなかなかうまいことがいかず・・・お願いします。 532 名前：デフォルトの名無しさん mailto:sage [2006/06/27(火) 23:57:19 ] >>531 m≠n のときはどうするんだｗ 533 名前：デフォルトの名無しさん [2006/06/28(水) 00:39:57 ] >>531 Fortran90ならＭＡＴＭＵＬ　という組み込み関数があるぞｗ 534 名前：531 [2006/06/28(水) 01:06:04 ] >>532 m≠nのときも考えるのだと思います。 一応はいりだけですが・・ DIMENSION A(M,N) WRITE(*,*) 'M,N=' READ(*,*) M,N DO 1 I=1,M DO 2 J=1,N WRITE(*,*) I,J, '要素を入力' READ(*,*) A(I,J) 2CONTINUE 1CONTINUE (つД`)たすけて〜 535 名前：デフォルトの名無しさん [2006/06/28(水) 01:30:42 ] >>534 助けてやってもよいが、行列の積を考えると　ｍ＊ｎ　と　ｍ＊ｎ　の積は取れないと >>532は言っている。 ｍ’　＊　ｎ　と　ｎ　＊　ｍ　のように　積をとる部分では次元がそろってないといけない。 もう少し問題文を正確に教えてくれ。 536 名前：デフォルトの名無しさん [2006/06/28(水) 01:53:09 ] >>534 一応ＦＯＲＴＲＡＮ７７で作ってみた。若干中級レベルの技法を使っている。 まぁ、素直に多重ＤＯ　ＬＯＯＰをまわすのもよし。 ＦＯＲＴＲＡＮ７７では行列の大きさを実行時に決められない。Ｆｏｒｔｒａｎ９０ならできる。 PROGRAM matrix PARAMETER (m = 2, n = 3) REAL a(n, m), b(m, n), c(n, n) C input A WRITE(*, *) ' input Matrix A ' WRITE(*, *) ' size =', n, ' * ', m READ(*, *) ((a(i, j), j = 1, m), i = 1, n) C input B WRITE(*, *) ' input Matrix B ' WRITE(*, *) ' size =', m, ' * ', n READ(*, *) ((b(i, j), j = 1, n), i = 1, m) C C = A * B DO 10 i = 1, n DO 20 j = 1, n c(i, j) = 0.0 DO 30 k = 1, m c(i, j) = c(i, j) + a(i, k) * b(k, j) 30 CONTINUE 20 CONTINUE 10 CONTINUE C print C DO 40 i = 1, n WRITE(*, '(9F8.4)') ( c(i, j), j = 1, n ) 40 END DO STOP END >>536　 537 名前：534 [2006/06/28(水) 02:36:30 ] >>535 ｍ行ｎ列のＡとｎ行ｍ列のＢということだと思われます。 >>536 ありがとうございます。 理解してみようと思います。 538 名前：デフォルトの名無しさん mailto:sage [2006/06/28(水) 03:17:43 ] 宿題ではないのですが、 WINDOWSのスタートアップに入れておいて、 PC起動してから約1時間後に勝手に計算を始めるFORTRANのプログラムをつくりたいのです。 そこで、以下のような概略を考えました。 ・スタートアップでプログラムスタート ・現在時刻を取得して格納・DOループスタート。 このループ中で時刻取得を繰り返し、上で取得した最初の時刻の1時間後過ぎになったらループ終了 ・やりたい計算を実行 ・プログラム終了 と、こんな感じで考えたのですが、PC起動後ずっと時刻取得のループをまわし続けるのが どうも無駄のような気がします。 こういう場合、どういうやり方が上手なんでしょうか 539 名前：デフォルトの名無しさん [2006/06/28(水) 08:52:40 ] >>538 ・案１ Windows のタスク機能を使う。 [スタート] ボタンをクリックし、[プログラム]、[アクセサリ]、[システム ツール] の順に ポイントし、[タスク] をクリックします。 ・案２　ＭＳ／ＤＥＣ／Ｃｏｍｐａｑ／ＨＰ／ＩＮＴＥＬ　ＦＯＲＴＲＡＮ　の場合 ＳＬＥＥＰ　関数を呼んで眠っていてもらう。 540 名前：デフォルトの名無しさん mailto:sage [2006/06/28(水) 17:01:19 ] 行列ﾃﾗｳｻﾞｽ 541 名前：デフォルトの名無しさん mailto:sage [2006/06/28(水) 18:30:52 ] >>540 Ｆｏｒｔｒａｎ９０なら行列気持ちよくいじれる．オナニーできるくらい． 542 名前：デフォルトの名無しさん mailto:sage [2006/06/28(水) 19:14:22 ] いいなー。ところで、 Gauss-Jordanのピボット選択ってどうやるんですか 絶対値が最大のと０になるところを入れ換える。。 バカなんでさっぱりわからん。。 543 名前：デフォルトの名無しさん [2006/06/28(水) 19:30:15 ] >>542 最低限度で良ければ，０が出た時に，そこから未処理分の列の係数をみていって 絶対値最大の行と入れ替えればいいだけだろう。 Ａｘ＝ｂ のｂの対応行も入れ替え。 入れ替えた行を覚えておいて元に戻すのがめんどいかな。 544 名前：デフォルトの名無しさん mailto:sage [2006/06/28(水) 22:53:56 ] 課題が出たのですが、まったくわかりません。 助けてください。 m行n列の２つの実行れい、A、Bをキーボードから読み込んで 行列の席C＝ABを計算し、結果を画面に出力しなさい。 留意点 　汎用性のあるプログラムを設計すること。 　出力設計については結果が見やすいことが肝心である。 　問題の分析、再定義が必要（特に配列の扱いについて） 　行列A、Bの各要素は**.*のように整数２桁、小数１桁の実数である。 と書いてありますが、まったくわかりません。 お願いします。 545 名前：デフォルトの名無しさん [2006/06/28(水) 23:05:25 ] >>544　昨日の>>531とは別人か？ ＦＯＲＴＲＡＮ７７とＦｏｒｔｒａｎ９０でどっちがいいね？ その問題は設定が悪いｗ　　Ｍ＊Ｎのサイズがあいまい。 >>536を見て、もっと詳しい要望を出せ。出力ＦｏｒｍａｔをＦ４．１に変えればおｋ。 ＞行列A、Bの各要素は**.*のように整数２桁、小数１桁の実数である。 これは整数と小数の言葉を間違って使っているｗ 546 名前：デフォルトの名無しさん mailto:sage [2006/06/28(水) 23:20:42 ] あぁー　ホント行列わかんねぇよ 考えてんだけど・・　Doループで頭がクルクルパー・・ >>543の入れ替えた行を覚えておいて元に戻すのがめんどいかな って具体的にどうするの？ 547 名前：デフォルトの名無しさん [2006/06/29(木) 00:48:44 ] すいません>>536についてですが c(i, j) = 0.0 というのはＣの初期値ということでよろしいのでしょうか？ 548 名前：デフォルトの名無しさん [2006/06/29(木) 01:10:31 ] >>546 ググルと色々サンプルプログラムが出てくるので参照してみたらいかがかな。 ピボッティングだったら教科書の類にもよく出てきているので、他人のプログラムを 読んでみるといいと思うんこ。 >>547 このループでｃ（ｉ，ｊ）に次々足してゆくわけだが、一番最初に０．０になっていないと困る。 DO 30 k = 1, m c(i, j) = c(i, j) + a(i, k) * b(k, j) 30 CONTINUE Ｆｏｒｔｒａｎでは変数の初期値は不定なのでごみが入っていることがある。 最近のコンパイラ（というか正確にはローダーだが）変数をゼロに初期化してくれることが 多いが、これを期待してはいかんざき。変数はみんな初期化する癖を付けるのがよろしいアル カイダ。 549 名前：544 mailto:sage [2006/06/29(木) 08:46:21 ] >>545ＦＯＲＴＲＡＮ７７です。 L＊N行列のものと、N＊M行列のものの積で汎用性を持たせろということなので 10*6と6*9や3*3と3*3もできなければならないものと思いました。 >>548ググってみましたがサンプルプログラムが汎用性があるのがわかりません＿|￣|○ 汎用性… 550 名前：デフォルトの名無しさん [2006/06/29(木) 10:48:31 ] >>549 まず『汎用性』というものを具体的に定義しなければいけない。 ＦＯＲＴＲＡＮ７７での汎用性であるから『行列の大きさを自由に変えられる サブルーチン形式』と解釈しておけばいい。 以下例・・・ 551 名前：デフォルトの名無しさん [2006/06/29(木) 11:29:38 ] PROGRAM matrix PARAMETER (nmax = 40) REAL a(nmax, nmax), b(nmax, nmax), c(nmax, nmax) CALL inpmat(nmax, nra, nca, a) CALL inpmat(nmax, nrb, ncb, b) IF (nca .NE. nrb) STOP 'matrix size mismatch' CALL matmlt(nmax, nra, nca, ncb, a, b, c) CALL primat(nmax, nra, ncb, c) STOP END C==== C = A * B ============================================================== SUBROUTINE matmlt(nmax, n, m, l, a, b, c) REAL a(nmax, m), b(nmax, l), c(nmax, l) DO 10 i = 1, n DO 20 j = 1, l c(i, j) = 0.0 DO 30 k = 1, m c(i, j) = c(i, j) + a(i, k) * b(k, j) 30 CONTINUE 20 CONTINUE 10 CONTINUE RETURN END >>551 552 名前：デフォルトの名無しさん [2006/06/29(木) 11:37:24 ] C===================================================================== SUBROUTINE inpmat(nmax, nrow, ncol, a) REAL a(nmax, nmax) WRITE(*, *) 'input Matrix Size ( row & col ).' READ(*, *) nrow, ncol IF (nrow .GT. nmax .OR. ncol .GT. nmax) THEN WRITE(*, *) 'matrix size must be smaller than', nmax, '.' STOP 'Increase nmax!' END IF DO 10 irow = 1, nrow WRITE(*, *) 'input Matrix element of ', irow, ' th row.' READ(*, *) (a(irow, icol), icol = 1, ncol) 10 CONTINUE RETURN END C===================================================================== SUBROUTINE primat(nmax, nrow, ncol, a) REAL a(nmax, nmax) WRITE(*, *) DO 10 irow = 1, nrow WRITE(*, '(1x, 40F5.1)') ( a(irow, icol), icol = 1, ncol ) 10 CONTINUE RETURN END >>552 553 名前：５５０ [2006/06/29(木) 11:40:01 ] ホントは整合配列をきっちりしてやんなきゃないが、面倒なのでいんちきしたｗ すぐ直せるから調べて直してくれｗ 色々エラーチェックとかも適当ｗ じゃあ 554 名前：544 mailto:sage [2006/06/29(木) 14:22:42 ] かなり難しいですね。。 わからない部分がほとんどですが、１つずつわからないところを つぶしていきたいと思います。 わからなくなったら聞きますのでよろしくお願いします。 ありがとうございました。 555 名前：デフォルトの名無しさん mailto:sage [2006/06/29(木) 15:03:41 ] サイズも渡せばいいだけじゃね？ 556 名前：デフォルトの名無しさん mailto:sage [2006/06/29(木) 16:48:13 ] >>555 うみゅ。 渡す配列毎に、使用した縦・横方向と、定義した縦方向の大きさが必要。 557 名前：デフォルトの名無しさん mailto:sage [2006/06/29(木) 21:01:54 ] FORTRAN77流にお行儀悪くやってみたお PROGRAM HOGE DIMENSION A(100), B(100), C(100) WRITE(*,*) 'INPUT L, N, M:' READ(*,*) L, N, M WRITE(*,*) 'INPUT A(',L, '*', N, '):' CALL MTINP(A, L,N) WRITE(*,*) 'INPUT B(', N, '*', M, '):' CALL MTINP(B, N,M) CALL MTMUL(A, B, C, L, N, M) WRITE(*,*) 'C =' CALL MTPRT(C, L,M) STOP END 558 名前：デフォルトの名無しさん mailto:sage [2006/06/29(木) 21:02:47 ] >>557の続きだお SUBROUTINE MTINP(A, M, N) DIMENSION A(M,N) READ(*,*) ((A(I,J), J=1,N), I=1,M) END C SUBROUTINE MTPRT(A, M, N) DIMENSION A(M,N) CHARACTER F*20 WRITE(F,'(''(1X,'',I3,''F7.1)'')') N WRITE(*,F) ((A(I,J), J=1,N), I=1,M) END C SUBROUTINE MTMUL(A, B, C, L, N, M) DIMENSION A(L,N), B(N,M), C(L,M) DO 10 J=1, M DO 20 I=1, L C(I,J) = 0. DO 30 K=1, N C(I,J) = C(I,J) + A(I,K) * B(K,J) 30 CONTINUE 20 CONTINUE 10 CONTINUE END 559 名前：デフォルトの名無しさん [2006/06/29(木) 23:19:03 ] >>554は、たぶん配列の添え字と縦横の関係も分かっていないと思われｗ Ｆｏｒｔｒａｎでは、第一添え字が縦方向、第二添え字が横方向を指定する。 たとえば３＊３の場合 Ａ（１，１）　Ａ（１，２）　Ａ（１，３） Ａ（２，１）　Ａ（２，２）　Ａ（２，３） Ａ（３，１）　Ａ（３，２）　Ａ（３，３） という並びになっている。これはＣ言語とは逆でＦｏｒｔｒａｎ型のほうが数学記号と合う。 また『行　（row　ロウ）』というのは、こういう横方向の並びをいう。 Ａ（１，１）　Ａ（１，２）　Ａ（１，３） こういう縦方向の並びは『列　（column　カラム)』という。 Ａ（１，１） Ａ（２，１） Ａ(３，１） ■やや高度な話題 Fortranでは、２次元の配列は、メモリー上ででは、１列づつ縦書きで日本語とは逆に 左から右へ読む形で１次元に並んでいる。つまり、 Ａ（１，１）　Ａ（２，１）　Ａ(３，１）　Ａ（１，２）　Ａ（２，２）　Ａ(３，２）　Ａ（１，３）　Ａ（２，３）　Ａ(３，３） これはサブルーチンに配列を渡すときに知っていると役に立つ。 ただしFortran90以降では、分散メモリー型の並列機に配慮したせいか、 規格上ではこの種の並びが保障されなくなった。とはいえ、まず信じて大丈夫ｗｗ 560 名前：デフォルトの名無しさん [2006/06/30(金) 02:08:56 ] どなたか、VB.NETでビリヤード程度のゲームのソースを お持ちの方、またサイトを知っておられる方 教えていただけませんか?? 561 名前：1/3 [2006/06/30(金) 04:06:33 ] どなたか↓をC言語にしてもらえませんか？F2C: FORTRAN to C Translator を使いこなせなくて・・・。よろしくお願いします。 C*** 　report2 C*** 　gakusekibangou DIMENSION YA(30),VA(30),YS(30),VS(30),YR(30),VR(30) OPEN(UNIT=16,FILE="dy.d") OPEN(UNIT=17,FILE="dv.d") G=9.8 R=0.3 CANALYTIC SOL. DO 10 K=1,30 T=K EP=EXP(-R*T) VA(K)=G/R*(1.0-EP) YA(K)=G/R*T-G/R/R*(1.0-EP) 10CONTINUE YS0=0.0 YS1=0.0 YS0=0.0 VS0=0.0 YR0=0.0 VR0=0.0 H=0.01 DO 20 K=1,3000 562 名前：2/3 [2006/06/30(金) 04:07:25 ] CSABUN T=0.01*K YS2=2.0*YS1-YS0+G*H*H-R*(YS1-YS0)*H VS2=(YS2-YS1)/H CRUNGE-KUTTA BK1=VR0 BM1=G-R*VR0 BK2=VR0+H*BM1/2.0 BM2=G-R*(VR0+H*BM1/2.0) BK3=VR0+H*BM2/2.0 BM3=G-R*(VR0+H*BM2/2.0) BK4=VR0+H*BM3 BM4=G-R*(VR0+H*BM3) YR1=YR0+H*(BK1+2.0*BK2+2.0*BK3+BK4)/6.0 YR1=YR0+H*(BM1+2.0*BM2+2.0*BM3+BM4)/6.0 IF(MOD(K,100).NE.0) GO TO 30 KK=K/100 YS(KK)=YS2 VS(KK)=VS2 YR(KK)=YR1 VR(KK)=VR1 30 YS0=YS1 YS1=YS2 VS0=VS0 VS1=VS2 YR0=YR1 VR0=VR1 563 名前：3/3 [2006/06/30(金) 04:08:26 ] 20 CONTINUE WRITE(17,100) DO 40 K=1,30 WRITE(17,150) K,VA(K),YS(K),VR(K) 40 CONTINUE WRITE(16,200) DO 50 K=1,30 WRITE(16,150) K,YA(K),YS(K),YR(K) 50 CONTINUE 100 FORMAT(5X,"T",5X,"VA",11X,"VS",11X,"VR") 150 FORMAT(2X,T3,3X,F10.5,3X,F10.5,3X,F10.5) 200 FORMAT(5X,"T",5X,"YA",11X,"YS",11X,"YR") STOP END 564 名前：デフォルトの名無しさん [2006/06/30(金) 10:01:15 ] >>560-563 UZEEEEEEEEEEEEeeeeeeeeeeeeeeeeeewwwwwwwwwww 565 名前：デフォルトの名無しさん [2006/06/30(金) 17:45:48 ] やっとこさコンパイルできたと思ったら実効できない・・・。 run-time error error format/data mismatch 00401000 MAIN# [+03e4] これはどういうエラー？？ 566 名前：デフォルトの名無しさん [2006/06/30(金) 18:28:15 ] >>565 FORMAT文のエラーでないかな。 オプションをつけるとソース行が出ることがある。 トレースバック　というデバッグオプションをオンにしてみそ。 デバッグオプションがわかんない場合は、 １）近くのわかっている人に聞くか、 ２）あきらめて、プログラムにＰＲＩＮＴ　＊、’ｕｎｋｏ　１’　等のＰＲＩＮＴ文をバラまいてエラー行を特定せよｗ 567 名前：デフォルトの名無しさん [2006/06/30(金) 20:46:43 ] すいません、real*8の変数がNaN(Not a Number)になったとき、プログラムを 終了させるにはどうしたらいいですか・

---

---

[ 続きを読む ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef