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

1 名前：デフォルトの名無しさん mailto:sage [2007/11/30(金) 09:55:26 ] このスレッドは、他のスレッドでは書き込めない超低レベル、 もしくは質問者自身何が何だか分からない質問を勇気を持って書き込むスレッドです。 FORTRAN使いが優しくコメントを返しますが、 お礼はFORTRANの布教と初心者の救済と次期Fortran2008規格でのCOMEFROM文採用をお願いします。 ●注意事項 ・質問する前にGoogle等の検索サイトで検索しましょう。 ・回答者にわかりやすい様に、質問内容はできる限り詳しく書きましょう。 ・エラーの場合は起きた状況、環境（OS・コンパイラ・バージョン）、エラーメッセージも詳しく書きましょう。 ●前スレ くだすれFORTRAN（超初心者用） pc8.2ch.net/test/read.cgi/tech/1138063703/ くだすれFORTRAN（超初心者用）その2 pc11.2ch.net/test/read.cgi/tech/1164121236/ ●関連スレ FORTRAN �W pc8.2ch.net/test/read.cgi/tech/1163319215/ 357 名前：デフォルトの名無しさん mailto:sage [2008/06/08(日) 16:42:46 ] 丸め誤差ってコンパイラオプションによっては違った結果になるよね あんま気をつけたことはないけど、精度命な人はその辺も考えてるのかな 358 名前：デフォルトの名無しさん mailto:sage [2008/06/08(日) 18:45:04 ] 文を改行後も続ける場合、&を行末に付加するって仕様、早く直らないかなぁ。 面倒で仕方ない。 359 名前：デフォルトの名無しさん mailto:sage [2008/06/08(日) 20:57:27 ] >>358 ？ 昔のように頭の方に継続記号を入れろってことか？ 360 名前：デフォルトの名無しさん mailto:sage [2008/06/08(日) 21:34:06 ] 考えようによってはCみたいに文末ごとに ; を付けるよりよっぽどお手軽だぞｗ 361 名前：デフォルトの名無しさん mailto:sage [2008/06/08(日) 22:06:47 ] >>358 Fortran文法は行番号があって，しかも予約語がない end = 1 + end print *, end, end; end とか do 10 do= 10, 10 10 end do end いや聞くだけでわくわくする文法ですね 362 名前：デフォルトの名無しさん mailto:sage [2008/06/09(月) 00:42:13 ] PASCALとかCとかって終端文字；を書くのがめんどくさくて死にそう FORTRANはやっぱイイナと思うよ。 363 名前：デフォルトの名無しさん mailto:sage [2008/06/10(火) 12:04:01 ] Fortran77の質問です。初心者なので次の課題がでてよく分かりません。お願いします。 　 「フィボナッチ数列を配列ｆ（）とDo文を使って計算させ、30項までの結果を8i8の書式で出力せよ。 「日本で現在の流通通貨は1万円札から1円玉まで9種類ある（2千円札を除く）。 5桁の金額 kin　を最も少ない枚数でそろえるには、それぞれ何枚になるか出力する。 貨幣9種の額はファイル　kahei.txt　から読み込む。出力の書式は○○yen　○○mai.」 364 名前：デフォルトの名無しさん mailto:sage [2008/06/10(火) 17:28:54 ] >>360 文が複数行に渡るなんてよくあるし、Cみたいな;のほうが良いと思う俺はFortraner失格ですか？ 最近遊びでVSのC#を触ったが、色んな意味で感動したw 他の言語触るのはなかなかイイもんですな 365 名前：デフォルトの名無しさん mailto:sage [2008/06/10(火) 18:11:18 ] >>363 フィボナッチらしき（間違えている）ものを作った。 宿題に合わせて、間違えを修正しろ。 program main integer f(10) f(1)=1 f(2)=1 do 10 i=3,10 f(i)=f(i)+f(i) 10 continue write(6,'( 8i8 )') (f(j),j=1,10) end 366 名前：デフォルトの名無しさん mailto:sage [2008/06/10(火) 18:48:57 ] >>363 とりあえず入出力はつくった。 計算式をよく考えろ。 program main integer kahei(9) integer kin open (1,File='kahei.txt') do 10 i=1,9 read(1,*) kahei(i) 10 continue print *, "input kingaku" read *, kin do 20 i=1,9 write(6,'( i5, "yen ", i1, "mai" )') kahei(i), kin/kahei(i) 20 continue end kahei.txt 10000 5000 1000 500 100 50 10 5 1 367 名前：デフォルトの名無しさん mailto:sage [2008/06/11(水) 00:45:49 ] >>363 しばらくぶりに７７で書いた。すでに>>365-366があってがっくしｗｗｗ PROGRAM FIBONA INTEGER F(0:30) F(0) = 0 F(1) = 1 DO 10 I = 2, 30 F(I) = F(I - 1) + F(I - 2) 10 CONTINUE WRITE(*, '(8I8)') (F(I), I = 0, 30) STOP END PROGRAM KANE INTEGER KAHEI(9), NUM(9) OPEN(9, FILE = 'KAHEI.TXT') DO 10 I = 1, 9 READ(9, '(I8.8)') KAHEI(I) 10 CONTINUE WRITE(*, *) 'INPUT KINGAKU' READ(*, *) KIN DO 20 I = 1, 9 NUM(I) = KIN / KAHEI(I) KIN = MOD(KIN, KAHEI(I)) 20 END DO DO 30 I = 1, 9 WRITE(*, '(I6, A, I5, A)') KAHEI(I), ' YEN ', NUM(I), ' MAI' 30 CONTINUE STOP END 368 名前：デフォルトの名無しさん mailto:sage [2008/06/11(水) 11:32:44 ] >>365-367 どうもありがとうございました。助かりました。 369 名前：デフォルトの名無しさん mailto:sage [2008/06/12(木) 00:12:47 ] >>367 「キルヒアイス、お前は優しいな。」 20ループと30ループ合体できない？ DO 20 I = 1, 9 NUM = KIN / KAHEI(I) KIN = MOD(KIN, KAHEI(I)) WRITE(*, '(I6, A, I5, A)') KAHEI(I), ' YEN ', NUM, ' MAI' 20 END DO 370 名前：デフォルトの名無しさん mailto:sage [2008/06/12(木) 02:28:49 ] >>369 それも思ったが、I/Oと主演算は分離すべきかなと思って分けた。 天皇陛下万歳！ 371 名前：デフォルトの名無しさん [2008/06/13(金) 01:34:43 ] すいません、初心者なんですが windows上で動くFORTRANのエミュってありますか？ 372 名前：デフォルトの名無しさん mailto:sage [2008/06/13(金) 01:44:23 ] >>371 VMwareにLinux入れてg77使えばおｋｗ 373 名前：デフォルトの名無しさん mailto:sage [2008/06/13(金) 01:46:34 ] >>371 ある。７７も９５も。 ぐぐれｗ 374 名前：デフォルトの名無しさん mailto:sage [2008/06/13(金) 04:53:55 ] >>373 お勧めの物はないのでしょうか？ 375 名前：デフォルトの名無しさん mailto:sage [2008/06/13(金) 08:01:12 ] fortran のエミュってどういう意味なんだろう？ 単に IDE があってコンパイルができればいいなら silverfrost ftn95 でググれば見つかるよ。 376 名前：デフォルトの名無しさん mailto:sage [2008/06/16(月) 02:07:35 ] >>374 G95とかどう？ ttp://ftp.g95.org/ Self-extracting Windows x86かそのへん 377 名前：デフォルトの名無しさん mailto:sage [2008/06/17(火) 01:57:48 ] >>376 cygwin入れるのがだるい。 378 名前：デフォルトの名無しさん mailto:sage [2008/06/17(火) 02:42:03 ] >>377 まあマジレスすると Self-extracting Windows x86 の方はMinGWだから Cygwin環境が無くても単独で動くわけだが。 379 名前：デフォルトの名無しさん mailto:sage [2008/06/17(火) 09:13:45 ] しかし Cygwin のインストールがだるいって相当だな 380 名前：デフォルトの名無しさん mailto:sage [2008/06/17(火) 19:07:53 ] 倍精度計算で、D○.○のようにして出力すると 例えば　0．12345D-3　のように表示されてしまうのですが これを1.2345D-4 と表示する事は可能でしょうか？ 試行錯誤しているのですが全然分からず教えて頂きたいです 381 名前：デフォルトの名無しさん mailto:sage [2008/06/17(火) 19:41:00 ] 昔はcygwin入れるのは当たり前って感じだったが今はな・・・ 最近影薄くて心配になるけど、細々と続いてるようね 382 名前：デフォルトの名無しさん mailto:sage [2008/06/17(火) 19:51:44 ] >>380 ESでできなかったかな。 383 名前：デフォルトの名無しさん mailto:sage [2008/06/17(火) 23:44:40 ] >>380 フォーマット指定するだけじゃん・・・ 384 名前：デフォルトの名無しさん mailto:sage [2008/06/17(火) 23:55:47 ] >>382-383 ESとはなんでしょうか･･ ごめんなさい、どのようにしたら良いのかさっぱりで・・ 385 名前：デフォルトの名無しさん mailto:sage [2008/06/18(水) 00:43:11 ] >>384 フォーマットの指定で"D○.○"の所を"ES○.○"にする "D○.○"→0.100D+00 "ES○.○"→1.000E+00 "ES○.○E3"→1.00E+000 試してないから間違ってるかも 386 名前：デフォルトの名無しさん mailto:sage [2008/06/18(水) 01:20:56 ] >>380 Fortran90でよければ、上にみんなが書いているようにESでおｋ。 FORTRAN77ならば、1P D○．○ 0P　でいける。 P　指定子は破壊力がでかいので、使ったらすぐ　0P　で元に戻すが吉。 ま、マニュアル嫁 PROGRAM ppp IMPLICIT NONE REAL(8) :: d d = 1.234567d-2 PRINT '(D15.6, ES15.6, 1P, D15.6, 0P)', d, d, d STOP END PROGRAM ppp 実行出力 0.123457D-01 1.234567E-02 1.234567D-02 Press any key to continue . . . 387 名前：デフォルトの名無しさん mailto:sage [2008/06/19(木) 11:57:12 ] P指定子ってすごいよな・・・ 今考えるとアリエンけど、昔はああいう仕様は全然平気だったんだろうな 388 名前：デフォルトの名無しさん mailto:sage [2008/06/19(木) 17:08:36 ] 世の中の処理系のほとんどが対応してると言っていいFortranのバージョンって、95でしょうか？ ちょっと前は、95はまだ早いから無難に行くなら90にしとけ、って声を聞きましたが。 389 名前：デフォルトの名無しさん mailto:sage [2008/06/19(木) 20:44:00 ] >>387 P指定子ってなんでありえないの？？ よく使うのだが… 390 名前：デフォルトの名無しさん mailto:sage [2008/06/19(木) 23:07:51 ] 昔、>>380の様に考えて出力用にP指定子を使ったが、クリアしていなかったので READのところでFフォーマットの値がスケールされて読まれてしまい、とんでもない結果が出て悩んだｗ たしかEとかベキ記号がついているときはスケールされないので、すげー悩んだ。 Fortran90の勉強を始めたとき、ESフォーマットがあったので、あぁアリガテーと思った。 391 名前：デフォルトの名無しさん mailto:sage [2008/06/20(金) 12:47:18 ] >>389 マジで言ってんの・・・？ まぁ気をつけてれば問題ないけどさ 392 名前：デフォルトの名無しさん mailto:sage [2008/06/23(月) 16:22:50 ] 学校の課題で次のものが出たんですけど、初心者なのでわかりません。お願いします。ちなみにFORTRAN77です。 今日の月日と曜日を入力し，これを用いて今年の任意の月日の曜日を出力させよ． 　各月の日数はﾌｧｲﾙから読み込む．ﾌｧｲﾙ名や変数名は自分で付ける． 曜日は数値で代用してよい．日曜=1，月曜=2，･･･土曜=7． 目標の月日が何日後か（前か）を出す．それから曜日を計算する． today 6 23 2 when 1 1 -174days 3you today 6 23 2 when 12 31 191days 4you 393 名前：デフォルトの名無しさん mailto:sage [2008/06/23(月) 17:25:50 ] >>392 「今年」と銘打っているのに、曜日の入力を要求するんだな。 今日の月日と曜日に不整合があった場合、どうするんだろう。 まずは次の課題を考えなさい。 今日の月日を入力し、一年の経過日数（正月＝1）を出力するプログラムを作成せよ。 394 名前：デフォルトの名無しさん mailto:sage [2008/06/23(月) 20:34:09 ] わかってるところとわからんことを書かないと、アドバイスのしようがないのでは？ 全部かいてもらう気なら別だが。 395 名前：392 mailto:sage [2008/06/23(月) 20:41:37 ] >>394 上にも書きましたが初心者です。 全部わからないので全部書いてください。 396 名前：392 mailto:sage [2008/06/23(月) 20:54:35 ] 上の人は別人です まだFORTRAN77を触り始めて間が無いので表示するプログラムが限界です 勉強のためにコメントをつけてくだされば嬉しいです って言いたいと思う 397 名前：デフォルトの名無しさん mailto:sage [2008/06/23(月) 21:17:09 ] >>396 本質的には >>395 と言っている事が一緒だね。 さらに解説ヨロってのも凄いが。 398 名前：デフォルトの名無しさん mailto:sage [2008/06/23(月) 21:36:47 ] まぁとりあえず、作ってみようぜ 399 名前：デフォルトの名無しさん [2008/06/23(月) 23:13:44 ] 学校の課題なんだから必要な構文の説明あるだろ 初心者じゃなくて怠け者だろ 400 名前：デフォルトの名無しさん [2008/06/24(火) 02:33:14 ] この問題がわかりません↓ 0.0 から 1.0 までの 10 個の値、 0.12, 0.32, 0.09, 0.83, 0.58, 0.39, 0.72, 0.44, 0.71, 0.24 を使って、起点 0.0、幅 0.2 としたときの頻度分布 を求めるプログラムを作成せよ。 結果はこのように表示されるようです。 s < 0.2 2 0.2 <= s < 0.4 3 0.4 <= s < 0.6 2 0.6 <= s < 0.8 2 0.8 <= s 1 if文を使うのはわかるのですが、頻度分布をどう求めていいのかわかりません。 fortran77です。 どうかよろしくお願いします。 401 名前：デフォルトの名無しさん mailto:sage [2008/06/24(火) 09:21:12 ] >>396 初心者には、入出力まわりが一番難しいのだが。 402 名前：デフォルトの名無しさん mailto:sage [2008/06/24(火) 10:36:12 ] >>400 PROGRAM HIST INTEGER NUM(5) REAL S(10) DATA S / 0.12, 0.32, 0.09, 0.83, 0.58, $ 0.39, 0.72, 0.44, 0.71, 0.24 / DO 10 I = 1, 5 NUM(I) = 0 10 CONTINUE C DO 20 I = 1, 10 IF (S(I) .LT. 0.0) THEN WRITE(*, *) 'INPUT ERROR!' STOP ELSE IF (S(I) .LT. 0.2) THEN NUM(1) = NUM(1) + 1 ELSE IF (S(I) .LT. 0.4) THEN NUM(2) = NUM(2) + 1 ELSE IF (S(I) .LT. 0.6) THEN NUM(3) = NUM(3) + 1 ELSE IF (S(I) .LT. 0.8) THEN NUM(4) = NUM(4) + 1 ELSE IF (S(I) .LE. 1.0) THEN NUM(5) = NUM(5) + 1 ELSE WRITE(*, *) 'INPUT ERROR!' STOP END IF 20 CONTINUE C 403 名前：デフォルトの名無しさん mailto:sage [2008/06/24(火) 10:36:47 ] WRITE(*, *) 's < 0.2 ', NUM(1) WRITE(*, *) '0.2 <= s < 0.4 ', NUM(2) WRITE(*, *) '0.2 <= s < 0.4 ', NUM(3) WRITE(*, *) '0.6 <= s < 0.8 ', NUM(4) WRITE(*, *) '0.8 <= s ', NUM(5) STOP END s < 0.2 2 0.2 <= s < 0.4 3 0.2 <= s < 0.4 2 0.6 <= s < 0.8 2 0.8 <= s 1 Press any key to continue . . . 404 名前：デフォルトの名無しさん mailto:sage [2008/06/24(火) 11:20:33 ] >>392 あまりチェックしてない。 曜日求め部分をよく練っていない。 うるう年対応してしまったｗ　今年限定ならデータいじればよし。 FUNCTION ILEAP(IY) ILEAP = 0 IF (MOD(IY, 4) .EQ. 0) ILEAP = 1 IF (MOD(IY, 100) .EQ. 0) ILEAP = 0 IF (MOD(IY, 400) .EQ. 0) ILEAP = 1 RETURN END C FUNCTION IDAYS(IY, IM, ID) INTEGER MDAYS(12) DATA NFEB /28/ DATA MDAYS /31, 0, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31/ MDAYS(2) = NFEB + ILEAP(IY) IDAYS = ID DO 10 I = 1, IM - 1 IDAYS = IDAYS + MDAYS(I) 10 CONTINUE RETURN END C 405 名前：デフォルトの名無しさん mailto:sage [2008/06/24(火) 11:23:50 ] PROGRAM CALEND 1 WRITE(*, *) 'TODAY' READ(*, *) IM0, ID0, IDW0 IF (IM0 .LE. 0) STOP 'NORMAL END' WRITE(*, *) 'WHEN' READ(*, *) IM1, ID1 IDAYS0 = IDAYS(2008, IM0, ID0) IDAYS1 = IDAYS(2008, IM1, ID1) IDIFF = IDAYS1 - IDAYS0 IDW = MOD(IDIFF, 7) + 7 + IDW0 - 1 IDW = MOD(IDW, 7) + 1 WRITE(*, *) IDIFF, 'DAYS', IDW, 'YOU' GOTO 1 STOP END TODAY 6 23 2 WHEN 1 1 -174 DAYS 3 YOU TODAY 6 23 2 WHEN 12 31 191 DAYS 4 YOU TODAY 0 0 0 NORMAL END Press any key to continue . . . 406 名前：デフォルトの名無しさん mailto:sage [2008/06/24(火) 12:03:44 ] >>404　乙 「各月の日数はファイルから読み込む」となっているね。 最近みかける宿題の出題者は、 無駄にファイル処理をさせようとしていない？ 407 名前：デフォルトの名無しさん mailto:sage [2008/06/24(火) 12:09:14 ] >>406 問題よく読んでなかったｗ 後で直す。 408 名前：デフォルトの名無しさん mailto:sage [2008/06/24(火) 13:48:19 ] >>402 ありがとうございました！！ NUMを使うんですね。 解決できました。感謝です。 409 名前：デフォルトの名無しさん mailto:sage [2008/06/24(火) 15:03:32 ] >>408 > NUMを使うんですね。 ピントがずれてないか？ 優しい人たち乙。 私なら「宿題は自分でやれ」で放置だな。 410 名前：デフォルトの名無しさん mailto:sage [2008/06/24(火) 15:47:09 ] >>409 初心者はピントがずれているもんだよ　orz 411 名前：デフォルトの名無しさん [2008/06/24(火) 20:20:12 ] 入力データ: test.datの中身 0,10,100,a,b,500, プログラムソース character dmm(10)*10 open(10,file='test.dat',status='old') read(10,*) (dmm(i),i=1,6) do i=1,6 write(*,*) dmm(i) end do stop end 出力結果 0 10 100 a b 500 compaq visual fortran 6.0では たとえば入力データの"b"->"*"に変更した場合のように 文字列にアスタリスクが入っている場合、エラーで落ちてしまいます。 cygwinのg77では問題なく表示されました。 コンパイラによって読めない文字列ってあるのでしょうか？ 412 名前：デフォルトの名無しさん mailto:sage [2008/06/24(火) 23:55:38 ] >>411 漏れのCVF6.6Cはｂを＊にしてもエラーでない。 単なるVer.6のバグでないか。 おぼろげな記憶ではVer.6→Ver.6.5は有料うｐデートで、結構大きな変化があったはず。 Ver.6.0→Ver.6.1　と　Ver.6.5→Ver.6.6は無料。 413 名前：407 >>392 mailto:sage [2008/06/25(水) 01:25:17 ] FUNCTION IDAYS(IM, ID, MDAYS) INTEGER MDAYS(12) IDAYS = ID DO 10 I = 1, IM - 1 IDAYS = IDAYS + MDAYS(I) 10 CONTINUE RETURN END C********************************************* FUNCTION IDW(IDIFF, IDW0) K = MOD(IDIFF + IDW0, 7) + 7 - 1 IDW = MOD(K, 7) + 1 RETURN END C********************************************* PROGRAM CALEND INTEGER MDAYS(12) OPEN(9, FILE = 'month.dat') READ(9, *) (MDAYS(I), I = 1, 12) C 1 WRITE(*, *) 'TODAY' READ(*, *) IM0, ID0, IDW0 IF (IM0 .LE. 0) STOP 'NORMAL END' WRITE(*, *) 'WHEN' READ(*, *) IM1, ID1 IDAYS0 = IDAYS(IM0, ID0, MDAYS) IDAYS1 = IDAYS(IM1, ID1, MDAYS) IDIFF = IDAYS1 - IDAYS0 WRITE(*, *) IDIFF, 'DAYS', IDW(IDIFF, IDW0), 'YOU' GOTO 1 STOP END 414 名前：407 >>392 mailto:sage [2008/06/25(水) 01:36:07 ] month.dat 31, 29, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31 >>406 今回は、うるう年を外部データとして与えようという親心なのかもしれない。 >>396 プログラムの解説 FUNCTION IDAYS(IM, ID, MDAYS) 一月一日から特定の月日までの日数を計算する。 FUNCTION IDW(IDIFF, IDW0) 今日からの日数差と今日の曜日から、その日数差の日の曜日を求める。 この計算式は、考えるのが面倒なのでちょっと複雑。 もう少しシンプルな書き方があるかもしれない。 FORTAN７７のMOD関数は負数を引数とすると、負の数を返すので＋７とかが必要。 Fortran90ならMODULO関数というのが新設されていて便利。 だがいずれにせよ７で割り切れるときは０を返すので、１〜７にするために-1/+1の 人工的な操作をした。自分で決める曜日を０〜６にすれば話はもっと楽になる。 PROGRAM CALEND まず各月の日数をファイルから読み込む。 今年はうるう年なので、うるう年の日数データを用意した。 あとは、I/Oの後、上で用意した関数を使ってやればいい。 ここでループから抜けるために、今日の月が０以下の数の場合終了するようにした。 415 名前：デフォルトの名無しさん mailto:sage [2008/06/25(水) 10:25:23 ] >>414 詳しい説明ありがとうございました。 416 名前：411 mailto:sage [2008/06/25(水) 13:24:19 ] >>412 Ver.6.0から6.1へアップデートしてみました。 今度は*のデータの後にブランクが入るようになりました。 つまり、>>411の入力データでaを*にすると、 * b と表示されて、500のデータは読めていない。 おっしゃる通り、Ver.6.0とか6.1固有のバグなのかもしれません。 こういうデータ処理には使うなってことでしょうか。 ありがとうございました 417 名前：デフォルトの名無しさん mailto:sage [2008/06/26(木) 02:19:25 ] >>416 俺の予感では、 read(10,*) (dmm(i),i=1,6) をフォーマットつけて read(10,’（A)') にすればおｋな気がする。根拠は宇宙からのメッセージ！ｗ 418 名前：デフォルトの名無しさん mailto:sage [2008/06/26(木) 14:08:07 ] FortranってCの/*〜*/みたいな範囲コメントってできないんでしょうか？ 419 名前：411 mailto:sage [2008/06/26(木) 14:16:04 ] >>417 そんなこと・・・あるかもなって思って試してみたけど、 だめだわ。エラーで落ちる。 cygwinのg77でもダメ。 その宇宙からのメッセージ、あんまあてに出来んぞ ま気持ちだけは嬉しかったわ。ありがとな（´・ω・`） 420 名前：デフォルトの名無しさん [2008/06/27(金) 00:29:58 ] >>418 fortranは、第1行目に「c」或いは「*」を入れるか、 「!」の後がコメントになる その他のコメントは聞いたことないし見たことない f90以降にそんなフォーマットがあるのかもしれんけど、おれは見たことない 421 名前：デフォルトの名無しさん mailto:sage [2008/06/27(金) 01:31:33 ] >>418 Fortranは行を基本単位とする言語なので、行毎に指示をするのが自然な形だ。 CやPascalのように、行が意味を持たない代わりに終端符号などをいちいち付ける言語とは 根源的な思想の違いがある。 しかしプリプロセッサを使えば実は範囲コメントも実現出来る。 ただ、プリプロセッサは本来の規格の外にあるので処理系依存の機能になる。 自分の処理系のFortranマニュアルでｆｐｐとかを調べてみるべし。 良いニュースとしては将来Fortran2008で、プリプロセッサが規格に入ることになっているので、 それまで、待つという手もある。 １０年くらい待つことになると思うがｗｗｗｗｗｗｗ 422 名前：418 mailto:sage [2008/06/27(金) 10:29:44 ] >>420,421 ありがとうございます。 行をまとめてズバッとコメントアウトしたいときとか不便だなぁと感じるのですが、仕様上諦めるしかないんですね・・・ 世の中のコンパイラがFortran2008に対応するのいつごろだろう 2003だってまだ新しい感じですよね 423 名前：デフォルトの名無しさん mailto:sage [2008/06/27(金) 12:39:40 ] >>422 え、それこそ1文字目にCや*を入れるだけでしょ。1行の前半分をコメントにするのは難しいけど。 まさか、エディタでしこしこ100行くらいの1文字目にC入れるのがしんどいとかだったら言語の問題じゃなくて あんたの未熟さの問題だけど。 424 名前：418 mailto:sage [2008/06/27(金) 17:39:39 ] >>423 確かに私は未熟者ですが・・・。 手動でしこしこやらずに、そういうことをサクッと一発でやる方法ってないですかね？ viでしこしこ書いてるんですが、sedなどの外部コマンド呼ぶのもアレだし emacsや他のIDEならそういう便利機能ってあるんでしょうか？ 425 名前：デフォルトの名無しさん mailto:sage [2008/06/27(金) 18:00:20 ] 例えばVisual Studio + Intel Fortranの場合、 Shift押しながら複数行選択し、コメントアウトボタンを押せば行頭に!つくよ。 426 名前：418 mailto:sage [2008/06/27(金) 18:14:37 ] >>425 ありがとうございます。VSはさすがに便利ですね。 今の開発環境は、単に端末からサーバに入ってvi、なんですが Windowsに移行しちゃおうかしら 427 名前：デフォルトの名無しさん mailto:sage [2008/06/27(金) 19:47:11 ] do while文は最適化を妨げるから使うなって聞いたのですが、本当ですか？便利なので使いまくってますが。 428 名前：デフォルトの名無しさん mailto:sage [2008/06/27(金) 23:03:57 ] >>426 viなら置換コマンドを使え。 10行目から20行目までの先頭に ! を挿入 :10,20s/^/!/ 30行目から40行目までの先頭の ! を削除 :30,40s/^!// 行頭の ! をすべて削除 :%s/^!// 429 名前：デフォルトの名無しさん mailto:sage [2008/06/27(金) 23:05:25 ] >>426 emacs では C-c ; でリージョンをコメントアウトできる 430 名前：418 mailto:sage [2008/06/27(金) 23:37:45 ] >>428,429 ありがとうございます。勉強になりました。 431 名前：デフォルトの名無しさん [2008/06/28(土) 08:28:25 ] >>418 コメントアウトしたいブロックを「#if 0」〜「#endif」で挟む。 UNIX系OS上の処理系なら、ソースファイルの拡張子を「.F」に、 Fortran90なら「.F90」にしておけば、コンパイル時に 自動でプリプロセッサを通してくれるはず。 432 名前：デフォルトの名無しさん mailto:sage [2008/07/02(水) 16:58:39 ] Fortran90のサフィックスは普通.f90ですが、2003や2008だとどうなるんですか？.f03や.f08？ 433 名前：デフォルトの名無しさん mailto:sage [2008/07/02(水) 20:15:58 ] >>432 サフィックスは固定形式と自由形式の違いの認識じゃない？ f90 以降は、.f90だと思う。 434 名前：デフォルトの名無しさん mailto:sage [2008/07/03(木) 09:47:29 ] >>433 そうでしたね・・・お恥ずかしい質問失礼しました。 435 名前：デフォルトの名無しさん mailto:sage [2008/07/05(土) 17:34:50 ] Fortran90以降の（自由形式の）ソースファイルとして認識されるファイルの拡張子は処理系によって .f90, .f95, .f2k, .f03, .f08 などがある。 固定形式の方は .f, .for, .ftn などがある。 436 名前：デフォルトの名無しさん mailto:sage [2008/07/10(木) 11:41:15 ] www.uploda.org/uporg1533847.f.html このファイルなんだけど、 stv=c1*cos(t)-(0.5d0)*(exp(-c*t/2.0d0)*(c*cos(t*e1)-(e)*sin(t*e1) の部分で絶対にエラーを吐かれます 何がいけないのかまったくわからないのですが文字数制限とかでしょうか？ 437 名前：デフォルトの名無しさん mailto:sage [2008/07/10(木) 11:46:53 ] 事故解決しました 438 名前：デフォルトの名無しさん mailto:sage [2008/07/11(金) 20:28:27 ] nはサブルーチンの引数で0<n<100の整数とします。 GOTO n 1 CALL hogehoge01(...) GOTO 999 2 CALL hogehoge02(...) GOTO 999 ... 99 CALL hogehoge99(...) 999 PRINT *, '終わり' みたいな事がやりたのですがコンパイルエラーになります。 現在は IF THEN FLSE で場合分けしているのですが分岐の数が多くて あまりスマートではありません。何か上手い方法は無いでしょうか？ 439 名前：デフォルトの名無しさん mailto:sage [2008/07/11(金) 21:23:35 ] >>438 使っている言語を詳しく 440 名前：デフォルトの名無しさん mailto:sage [2008/07/12(土) 00:31:01 ] >>439 Intel Fortran Compiler 10.1 です 441 名前：デフォルトの名無しさん mailto:sage [2008/07/12(土) 00:36:14 ] すいません、言語でしたね。Fortran 90 を使っています。 Fortran 95 で Intel Compiler がサポートしている機能なら使っても問題ないです。 442 名前：デフォルトの名無しさん mailto:sage [2008/07/12(土) 14:38:07 ] >>438 割り当てGOTO文の行き先はASSIGN文で与えないといけない。ただし、常数じゃないと駄目だったはず。 GOTO n, (1, 2, 3, 4, ....) と書く方法もあるが。 まぁFortran90だったら、素直に　SELECT　CASE　にしとけ。 今時Fortran66時代の割り当てGOTOを持ち出し、しかもIF-THEN-ELSEより割り当てGOTOを スマートだと感じ、エラーメッセージを貼らないで質問し、 おまけにhogehogeなどというUnix界隈の嫌な語を使うｗｗｗｗ はっきり言えば、チミはセンスがズレてるから、ちゃんとFortran90以降を基礎から習え。 443 名前：デフォルトの名無しさん mailto:sage [2008/07/24(木) 10:03:56 ] 昔のFORTRANって引数に渡した定数をサブルーチン側で変更すると 呼び出し側にも影響が出たって本当ですか? function(1, 2) とかで呼んだ場合定数の1と2が入れ替わるみたいな話を耳にしたのですが。 444 名前：デフォルトの名無しさん mailto:sage [2008/07/24(木) 17:17:05 ] フォートラン90なんですが、途中結果をターミナルに出力させた後に pause 文を入れて、出力結果を確認してリターンキーを押さないと プログラムが再開しないようにしてあります。 これを、リターンキーを押さなくても、ある程度の時間（１0秒とか）が 経つと勝手に再開するようには出来ないのでしょうか。。 445 名前：デフォルトの名無しさん mailto:sage [2008/07/24(木) 19:30:07 ] >>444 OS とコンパイラは何？ 446 名前：444 mailto:sage [2008/07/24(木) 20:38:42 ] WinXpのcygwinで g95でコンパイルしてます 447 名前：デフォルトの名無しさん [2008/07/24(木) 21:40:55 ] fortranを使って双子素数（連続した奇数が素数になっているものの組）を出力する課題に取り組んでいるんですが 完全にお手上げ状態です。 一応素数の出力はできてるみたいなんですが・・・ *素数の出力 INTEGER n,i,k * READ(5,*) n i=2 WRITE(6,*) i * DO 10 i=3,n,2 k=3 30 IF (k**2 .GT. i) GOTO 20 IF (mod(i,k) .EQ. 0) GOTO 10 k=k+2 GOTO 30 20 WRITE(6,*) i 10 CONTINUE * STOP END こんな感じです。 違う方法でやったほうがいいんでしょうか… 448 名前：デフォルトの名無しさん mailto:sage [2008/07/24(木) 22:46:10 ] 解読できんw 449 名前：デフォルトの名無しさん mailto:sage [2008/07/24(木) 23:05:23 ] >>446 CALL SLEEP(10) 　! 10秒間待つ >>447 そのプログラムを修正するのなら、 　・調べた奇数が素数だったのか素数で無かったのかを記録する変数を用意して 　　「DO 10 i=…」 のループで毎回結果を記録しておく。 　・素数が見つかったら上記の変数を調べて、１つ前の奇数も素数だった場合は 　　両方の値を出力する。 450 名前：デフォルトの名無しさん mailto:sage [2008/07/24(木) 23:34:50 ] >>443 その通り。ただ定数を直接書くのではなくパラメータ指定した定数を書かないと はっきり分からんだろう。これはサブルーチンの引数が参照渡しなので起こりうる。 SUBROUTINE sub_test(a, b) IMPLICIT NONE REAL, INTENT(IN OUT) :: a, b REAL :: tmp tmp = a a = b b = tmp RETURN END SUBROUTINE sub_test !=========================================== PROGRAM test IMPLICIT NONE REAL, PARAMETER :: a = 1.0, b = 2.0 EXTERNAL :: sub_test PRINT *, 'a=', a, ' b=', b CALL sub_test(a, b) PRINT *, 'a=', a, ' b=', b STOP END PROGRAM test Fujitsu-Ｌａｈｅｙ Fortran95 実行結果　(Ｉｎｔｅｌ　Fortran は実行時エラー) a= 1.00000000 b= 2.00000000 a= 2.00000000 b= 1.00000000 Program Completed Press Enter to Continue. >>444 過去レス>>326と重複している。その回答を読め。 451 名前：デフォルトの名無しさん [2008/07/25(金) 00:13:07 ] >>449 本当に申し訳ないのですが素数だったか素数じゃなかったかを判別するプログラムが作れません… まるなげになってしまって申し訳ないですが具体的にどうなるか教えていただけないでしょうか？ （面倒だったり自分でやれとおっしゃるならいいですが） 452 名前：デフォルトの名無しさん mailto:sage [2008/07/25(金) 09:09:52 ] 449 ではないが。 素数を画面に書かせることはできてるわけでしょ？ それなら （１）素数を画面に書くかわりに、適当な変数たとえば prime_prev に代入 （２）あらたに素数を見つける。その素数と prime_prev との差を調べる （３）差が2なら双子素数。2でないなら、（２）で見つけた素数を prime_prev に代入して 　　　ふたたびステップ2へ。 でいいんじゃね？ 453 名前：デフォルトの名無しさん mailto:sage [2008/07/25(金) 11:54:47 ] 普通素数はエラトステネスのふるいをつかって求めるから、 その表で連続している奇数を捜せば良いのではないかい？ 454 名前：名無し [2008/07/28(月) 17:48:07 ] FORTRAN77を用いてモンテカルロシミュレーションのプログラムを作成しようと思っています。 しかしFORTRAN77を用いたプログラミング自体始めたばかりであり、プログラムの作成方法についてピンときません。 もしよかったならどのようにしてプログラムを作ればよいのか教えてもらえませんか？ 大まか過ぎる質問ではありますが・・・。 455 名前：デフォルトの名無しさん mailto:sage [2008/07/28(月) 23:27:16 ] その質問だと回答してくれる人少ないんじゃない？ まず 1. 計算機にやらす内容を手作業でできるかどうか（時間は別として） ここでつまづくならFortranの問題じゃなく、 モンテカルロって何？ 何を入力して、何を出力するの？ などがわかってないのでちょっと別問題。 1.をクリアしてるなら、次に 2. その処理をコーディングしていくと、どこでつまづくのか を考えよう。 その部分が明確になっているなら、このスレッドの方は回答くださるかと思います。 大まかな質問なので大まかな回答しかできないけど RANDOM_NUMBERやRANDとかいった乱数発生関数があるのでこれを利用することになるんだと思うよ。 ちなみにどんな関数があるかはコンパイラに依存するので・・。 456 名前：デフォルトの名無しさん mailto:sage [2008/07/28(月) 23:57:35 ] >>454 分子シミュレーション講座1　モンテカルロ・シミュレーション ISBN978-4-254-12691-4、神山新一・佐藤明、朝倉書店 付録　7.7　基本的なFORTRANの計算プログラム集 図書館で借りろ。 457 名前：デフォルトの名無しさん mailto:sage [2008/07/29(火) 00:49:18 ] >>454 FORTRANでプログラム出来るとはありがたいこと。 まず天皇陛下に感謝をささげよ。 次に先輩から下敷きになるプログラムを貰うこと。 これで将来の運命が決まる。 筋の悪いウンコプログラムを貰ったらもうアウトｗ 筋が良くても長大な完成版を貰ったらやはりアウトｗｗ 筋が良くて簡素で、自分で拡張してゆける物を貰えたら、おのずから上達して成功する。 つまり、すべては運なので神に祈れｗ

---

---

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

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