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

1 名前：デフォルトの名無しさん [2006/11/22(水) 00:00:36 ] このスレッドは、他のスレッドでは書き込めない超低レベル、 もしくは質問者自身何が何だか分からない質問を勇気を持って書き込むスレッドです。 FORTRAN使いが優しくコメントを返しますが、 お礼はFORTRANの布教と初心者の救済と次期Fortran2008規格でのCOMEFROM文採用をお願いします。 ●注意事項 ・質問する前にGoogle等の検索サイトで検索しましょう。 ・回答者にわかりやすい様に、質問内容はできる限り詳しく書きましょう。 ・エラーの場合は起きた状況、環境（OS・コンパイラ）、エラーメッセージも詳しく書きましょう。 ●前スレ くだすれFORTRAN（超初心者用） pc8.2ch.net/test/read.cgi/tech/1138063703/ ●関連スレ FORTRAN �W pc8.2ch.net/test/read.cgi/tech/1163319215/ 910 名前：デフォルトの名無しさん mailto:sage [2007/10/23(火) 20:44:40 ] >>909 通常は拡張子が .f .for .ftn の場合は固定形式、 .f90 の場合は自由形式と解釈される。 # ifort は .f95 をソースとして認識しない。 コンパイル時のオプションに -fixed / -free を指定することによって強制的にそれぞれ 固定形式／自由形式の指定をすることも出来る。 また、拡張子を大文字（.F .FOR .FTN .F90）にするとコンパイル前にプリプロセッサで前処理される。 > 90の機能と何でも混ぜていいわけではないんですね。 １つのソースに固定形式／自由形式を混ぜてはいけない。 ただし、別々のファイルにすれば混在させることも出来る。 911 名前：９０７ mailto:sage [2007/10/23(火) 23:48:50 ] >>910 ありがとうございます。非常に勉強になりました。 912 名前：デフォルトの名無しさん mailto:sage [2007/10/24(水) 03:21:45 ] High Performance Fortranていう仕様があるようですが、どういうところで使われているんですか？ 普通のユーザには縁のないものですか 913 名前：デフォルトの名無しさん mailto:sage [2007/10/24(水) 09:04:39 ] >>912 90年代に並列プログラミングを考慮した規格として提案された。 日本のスーパーコンピュータメーカはコンパイラをちゃんと作ったが、 アメリカがやる気を見せず崩壊した。 地球シミュレータでもHPFは動いている。というかHPFはESに要求された仕様の一つだったような。 914 名前：デフォルトの名無しさん mailto:sage [2007/10/24(水) 11:13:06 ] >>913 崩壊、ですか・・・ 推進評議会のサイトを見ると、「ユーザが最小限の指示文によってデータの分割配置の方法を 指定すれば、残る作業(計算の分割と通信の生成)をコンパイラが自動的に行う」とありますが、 分散メモリシステムで必要となる領域分割と通信を自動でやってくれる、ということですか？ 今はFortran90+MPIという組み合わせが一般的だと思いますが、これよりも良い性能が出るんですか？ 915 名前：デフォルトの名無しさん mailto:sage [2007/10/24(水) 13:22:45 ] >>914 MPIはどちらかというと消極的な理由でデファクトスタンダードになった。 ソフトウェア中心主義の視点からは、MPIはハードウェアの抽象度が低いので、 HPFのほうが上位に立つのだが、現実的にはハードウェアへの依存性が高くて逆転する。 ここに過渡期に実務で使った人の比較意見が載ってる。 ttp://grape.mtk.nao.ac.jp/~makino/articles/future_sc/note031.html#rdocsect36 916 名前：デフォルトの名無しさん mailto:sage [2007/10/24(水) 16:03:35 ] >>915 消極的ですか。確かにそんな気もしますね OOPの可能性にｗｋｔｋするも、性能とのトレードオフに嘆き、手続き型も含めた現実的な方法に帰着する、という話と似てるかも ところでそのサイト見ていつも思うんですが、マ○ノさんてすごい人ですよね ここまでコアな人はなかなかいないと思う 917 名前：912=914=916 mailto:sage [2007/10/24(水) 16:04:56 ] >>915 スイマセン、お礼言うの忘れてました。 どうも有難うございました。 918 名前：デフォルトの名無しさん mailto:age [2007/10/28(日) 19:52:11 ] WRITE文で改行しない標準FORMATを使いたいのですがどうすればよいでしょうか？ 具体的には、 WRITE(*,*) 'x = ', x という文を改行なしで表示したいのです。ここで、xはDOUBLEで、 xの大きさによって浮動小数点表示と固定小数点表示をコンパイラが判断して欲しいです。 よろしくお願いします。 919 名前：デフォルトの名無しさん [2007/10/28(日) 22:10:09 ] >>918 改行しないのはFAQ　このスレでも何度も出ている。 Fortran90　では、ADVANCE='NO'　 数値の大きさでフォーマットを変えるのは　G　Format　これが７７時代からある。 サンプル PROGRAM unko IMPLICIT NONE REAL :: x, y x = 0.1 y = 1.234567e11 WRITE(*, '(A, G15.5)', ADVANCE = 'NO') 'x=', x PAUSE WRITE(*, '(A, G15.5)', ADVANCE = 'NO') 'y=', y STOP END PROGRAM unko 920 名前：デフォルトの名無しさん mailto:sage [2007/10/29(月) 23:32:38 ] COMMONブロックをモジュールに書き換える作業をしております。その際わからないことがあったので質問させてください。 jjoo.hp.infoseek.co.jp/tips/f90/common_to_module.html のサイトを参考に進めてきました。現在のプログラムの概略は次のような感じです（巨大なため全文ではありませんがご容赦くださいませ） module sub integer i,j end modeule sub program main 　use sub 　implicit none call 1 call 2 end program main subroutine 1 implicit none integer i,j (略) end subroutine 2 implicit none integer i,j (略) end この状態でコンパイル及び動作は出来るのですがサブルーチンにおいて変数部を subroutine 1 use sub implicit none (略) end とするとコンパイル時に symbol 'rmax' ai (1) is defined in multiple modules というエラーが出てしまいます。 これはどのように解決したらよいのでしょうか。よろしくお願いします。 921 名前：918 mailto:sage [2007/10/29(月) 23:48:47 ] >>919 欲しい結果が表示されました。ありがとうございます。 922 名前：デフォルトの名無しさん mailto:sage [2007/10/30(火) 00:56:45 ] >>920 それはエラーメッセージに書いてあるがごとく、複数のモジュールの中で同じシンボル名を 重複して使っているという事だろう。 コモンブロックに当たるモジュールは１個だけ作って、それをＵＳＥで呼び出せばいい。 コモンのときは全副プログラムにコモンを置かねばならなかったが、もジュールは違う。 923 名前：デフォルトの名無しさん mailto:sage [2007/10/30(火) 14:05:25 ] CHARACTER*(*) HOGE　のように、配列数が(*)となっているのはどういう意味なんでしょうか 動的に確保するということでしょうか？ 924 名前：デフォルトの名無しさん mailto:sage [2007/10/31(水) 01:18:34 ] >>923 HOGEの長さは HOGEが定数の場合、PARAMETER文で指定された文字列の長さ。 HOGEが副プログラムの仮引数の場合、対応する実引数と同じ長さ。 HOGEが関数副プログラムの場合、その関数を使用する側で予め決められた長さ。 925 名前：920 mailto:sage [2007/10/31(水) 05:54:42 ] >>922 遅くなって申し訳ありません。レス有り難うございました。 なるほど、モジュールは一度読み込めばよいと言うことだったんですか。 そこで、サブルーチン中のuse文を消してみたところ そのエラーは吐かなくなりました。 その代わり x0(n)=2.d0*rn 　1 Unexpected array Reference at (1) と言ったようなエラーが出るようになってしまいました。エラー文では配列の参照がおかしいみたいなことを書いてますがどのようなことなのでしょう。 内容が同じ物をcommon文としてかいた時は問題なく動作したので 配列の宣言がおかしいってこともないと思うのですが・・・ 926 名前：デフォルトの名無しさん mailto:sage [2007/10/31(水) 10:58:31 ] >>923 FORTRAN90以前では、基本的に変数は参照渡しで、実際には先頭番地しか渡していない。 また呼び出し側と呼び出される側での引数の整合性もチェックしていない。 さらに実行時にも、デバッグオプションをつけない限り、はみだしチェックのようなことはやってない。 したがって、副プログラム側では、引数の大きさはどうでもいいことになる。 このためFORTRAN66あたりまでのプログラムでは、配列にも適当な長さを指定するのが伝統だった。 （文字列は、まだ存在しなかったので・・・というか整数型に突っ込んでいた・・・・質問から外れるが） というわけで、FORTRAN77規格で、それを明示的に指定した書式が（＊）だ。（とおもうｗ） 長さ『＊』は、『適当』という意味だ。文字列のみならず数値型配列の最後の要素にも使える。 Fortran90以後は、色々コンパイル時のチェックが入るようになって、PASCALみたいでうぜぇという人達もいる。 最近のスクリプト言語は、型にゆるいのが多いが、昔のFORTRANに先祖がえりしているのである（嘘ｗ 927 名前：デフォルトの名無しさん mailto:sage [2007/10/31(水) 11:01:20 ] >>925 モジュールの宣言は１個で、USEは使うルーチンごと毎回だ。 カスケード的にUSEしているなら、深いほうで１回読んでおけばいい。 928 名前：923 mailto:sage [2007/10/31(水) 18:13:56 ] どうもありがとうございます。 >>924 ＞HOGEが定数の場合、PARAMETER文で指定された文字列の長さ。 HOGEをHOGE="hoge"と変数にしたら駄目でした。定数しか駄目なんですね。 >>926 ということは、サブルーチン仮引数の配列は全部*にすれば、配列サイズの違いに悩まされることもないのですよね？ 例題などで、配列を10などと明示的に指定しているのを多く見かけるのですが・・・ 929 名前：デフォルトの名無しさん mailto:sage [2007/10/31(水) 23:21:08 ] >>928 ＞ということは、サブルーチン仮引数の配列は全部*にすれば、配列サイズの違いに悩まされることもないのですよね？ ＞例題などで、配列を10などと明示的に指定しているのを多く見かけるのですが・・・ いや、別に１０のままでも、本質的には何も変わらない。 どっちにしろ配列の先頭アドレスしか渡さないので、サイズの違いに悩まされることは無い。 あと、多次元配列の場合、最後の次元以外は一応意味があるので、ちゃんと指定しないといけない。 ただ、FORTRAN77では整合寸法を別に引数で送るのが適切とされている。 930 名前：デフォルトの名無しさん [2007/11/01(木) 01:12:57 ] 何が原因かわからないのですいませんがどこが間違ってるのかご指摘お願いします． 縦300横10の配列を読み込んで，ただの列にしたいだけなんです． エクセルを使えば，FORTRAN使うほどのことでもないのでしょうか？ それも含めて教えてください． program hairetu open(30,file='a.dat',status='old') do 10 j=1,10 read(30,*)(b(i,j),i=1,300) 10 continue close(30) open(20,file='b.dat',status='new') do 50 i=1,300 do 40 j=1,10 write(20,*)b(i,j) 40 continue 50 continue close(20) stop end 931 名前：930 [2007/11/01(木) 01:19:49 ] 付け足しです．このようなエラーメッセージが出ています． 0004) read(30,*)(b(i,j),i=1,300) *** B is a function name - a function reference is not allowed in a READ statement list *** Invalid input/output list 2 ERRORS 932 名前：デフォルトの名無しさん mailto:sage [2007/11/01(木) 02:20:00 ] >>930 > どこが間違ってるのかご指摘お願いします． ・配列 b の宣言がない。 ・「縦300横10」のデータなら読み込み時の i, j の終値が逆。 > エクセルを使えば，FORTRAN使うほどのことでもないのでしょうか？ これくらいの処理ならExcelでもなんでも結果が出せればおｋ ということで。 gawk "{ for (i=1; i<=10; i++) print $i }" a.dat > b.dat 933 名前：デフォルトの名無しさん mailto:sage [2007/11/01(木) 02:36:36 ] ちなみにスレ違いだがWindows 2000以降なら以下のようなバッチファイルでも処理できる。 @echo off if exist b.dat del b.dat for /F "delims=" %%r in (a.dat) do ( for %%f in (%%r) do echo %%f >> b.dat ) 934 名前：923 mailto:sage [2007/11/01(木) 13:21:29 ] >>929 ＞どっちにしろ配列の先頭アドレスしか渡さないので、サイズの違いに悩まされることは無い。 そうなんですか。ちょっと混乱してきたのですが ということは、仮引数の配列数は結局は関係ないということですか？ ＞いや、別に１０のままでも、本質的には何も変わらない。 では仮引数の配列数を明示指定するのは、単なる慣習みたいなもの、ってことですか？ 935 名前：デフォルトの名無しさん mailto:sage [2007/11/02(金) 00:14:19 ] >>934 そういえなくも無い。1960年代のプログラムを見ると、任意長の配列がよく大きさ１０で宣言されている。 ただテキストの場合は、本当に１０個の要素だと思うがｗ SUBROUTINE ZERO(n, a) DIMENSION a(10) DO 10 i = 1, n a(i) = 0.0 10 CONTINUE END 見たいなのがよく使われていた。 これとまぁ等価なのが　DIMENSION a(*) で、もっと上品な書き方は整合配列を用いて DIMENSION a(n) とやる方法。後者をお奨めする。だがFortran90で書くことをもっとお奨めするｗ 先頭番地しか渡していなので、呼ぶ側では、配列の途中から渡したりできる。 REAL a(100) CALL zero(100, a) CALL zero(50, a(51)) --->Fortran90での a(51:100) = 0.0 普通のテキストには書いてないが、77までのプログラムでは、こういうのを知らないとまともに プログラムが読み書きできない。昔は誰も教えてくれないので徒弟制度で叱られながら学ぶか、 他人のコードを読んで泣きながら覚えるしかなかった。 漏れの読んだテキストで、そういう事が書いてあったのは、この本くらいか。 ttp://www.star.le.ac.uk/~cgp/prof77.html （他の本にも書いてあったが、単に読んだ時のスキルが足りなくて理解できなかったのかもしれないがｗ） 936 名前：923 mailto:sage [2007/11/02(金) 14:17:22 ] >>935 なるほど・・・よくわかりました。そういうコアな話を聞ける人が周りにいないもので、とても勉強になります。 丁寧な解説ありがとうございました 937 名前：デフォルトの名無しさん mailto:sage [2007/11/04(日) 14:04:31 ] ＞ただテキストの場合は、本当に１０個の要素だと思うがｗ これどういう意味ですか？ 文字列配列も普通の配列に変わりはないと思うんですが 938 名前：デフォルトの名無しさん mailto:sage [2007/11/04(日) 18:15:19 ] >>937 テキストは本当に１０個に固定の要素を扱っているのではないかという事。 深い意味は無い。 939 名前：937 mailto:sage [2007/11/05(月) 01:33:05 ] 気になったので簡単にテストしてみました。 サブルーチン内で普通に文字列代入する（sub1）のと直接要素に代入する（sub2）のとでは違うようですね。 後者は通常の配列と同様な参照渡しなのに対し、前者は少し違った処理をする、ということでしょうか？ implicit integer(i-n) character*3 text1 /'aaa'/ character*3 text2 /'aaa'/ call sub1(text1) write(*,*) text1 call sub2(text2) write(*,*) text2 end subroutine sub1(text) character*2 text text = 'b' end subroutine sub2(text) character*2 text text(1:1) = 'b' end <結果> b a baa 940 名前：デフォルトの名無しさん mailto:sage [2007/11/05(月) 03:17:08 ] >>939 たぶんサブルーチン側ではtext*2で宣言しているので、 sub1 では　text(1:2) = 'b ' の扱いになったんだと思う。 参照渡しなのは変わんないと思う。 text = 'b' は、一文字目は'b'として2文字目以降宣言長まではスペースを補填するのだと思う。 941 名前：デフォルトの名無しさん mailto:sage [2007/11/07(水) 13:28:20 ] fortranの仕様ってどこが決めてるの？ 少なくとも数値計算屋だよね？ 942 名前：デフォルトの名無しさん mailto:sage [2007/11/08(木) 00:39:48 ] >>941 ttp://www.j3-fortran.org/ いちおう、（コンパイラ、ハードウェア）ベンダー、ユーザー、大学人なんかが入っているはず。 なんだかイギリスが生意気にデカイ面している。 Fortran2008で、Crayのハードウェアに最適化されているCoArrayが規格に入ったのが気に食わない。 943 名前：デフォルトの名無しさん mailto:sage [2007/11/08(木) 01:06:31 ] コーエンちゃん早く2003の処理系つくって 944 名前：デフォルトの名無しさん mailto:sage [2007/11/08(木) 14:40:25 ] >>942 やっぱ向こうじゃNECよりCrayが優先されるんだろうね 945 名前：デフォルトの名無しさん mailto:sage [2007/11/08(木) 23:02:52 ] >>944 『向こう』ってなんだよｗ　欧米コンプレックスの４０代の爺さんかよｗｗｗｗ Fortran2008は、本来はFortran95のようにマイナーアップデートのはずなのに、 こういう妙な機能を入れ込ませようとする政治臭の漂うことをやると、信頼性をなくすと思うのだがな。 946 名前：デフォルトの名無しさん mailto:sage [2007/11/10(土) 22:42:08 ] HP Fortran で以下の二つを実現することは可能ですか？ 1) リトルエンディアン形式のバイナリ出力 2) Intel Fortranのopen(10, 中略, form='binary')に相当するやり方 947 名前：デフォルトの名無しさん mailto:sage [2007/11/10(土) 23:14:02 ] >>946 マニュアル嫁 1) リトルエンディアン形式のバイナリ出力 Itanium機では出来る？ ttp://docs.hp.com/ja/B3908-90009/ch08s08.html 2) Intel Fortranのopen(10, 中略, form='binary')に相当するやり方 DIRECT　ACCESSで１バイトづつ読むというテクニックで切り抜けられる。 948 名前：デフォルトの名無しさん mailto:sage [2007/11/10(土) 23:44:30 ] そういやIntel Fortran もようやく ACCESS='STREAM' が使えるようになったな。 これで処理系毎に書き直さなくて済むようになる。 949 名前：デフォルトの名無しさん mailto:sage [2007/11/10(土) 23:58:25 ] >>947 ありがとうございます。 1) そういえば、CPUの種類は確認してなかった・・。 たぶんItaniumだと思うんですが、、。 職場行ったら確認して試してみます。 2) これはちょっと大変そうですね。 ヘッダー部、データ部、それぞれで型の違う変数を扱わなければならず、 変換が大変そうです。 4byteの変数を出力する際、1byteずつ逆から出力する必要があるってことですよね・・ 950 名前：デフォルトの名無しさん mailto:sage [2007/11/11(日) 01:12:25 ] >>949 文字列で読み出して、それをINTEGERやREALやREAL*8とEquivalenceをかけて置けばOK 処理系によって、RECL=1　にしたときに１byteになるものと、１語になるものがある。 たとえばIntelFortranではデフォールトは１語＝4バイトだが、オプションで切り替えできる。 951 名前：デフォルトの名無しさん mailto:sage [2007/11/11(日) 01:14:51 ] Fortran2003でOPEN( ... ACCESS='STREAM') が規格に入ったので、新しいコンパイラなら先取り導入 されている可能性大。 952 名前：949 mailto:sage [2007/11/11(日) 12:31:09 ] >>950 Equivalence使ったことありませんでした。 マニュアルみたら違う変数で同じメモリ領域を使えるようにする、、ってやつなんですね。 つまり整数1byte型x1,x2,x3,x4がy4,y3,y2,y1と同じメモリ領域になるようにし、 逆から順に1byteずつ出力すればいいってことですね。 RECLに関してはマニュアルを読んでみます。 >>951 'DIRECT'、'SEQUENTIAL'、' POSITION=APPEND' マニュアル見たらACCESS指定子に選択できるのは上記3つだけみたいでした。 953 名前：デフォルトの名無しさん mailto:sage [2007/11/11(日) 13:45:48 ] >>949,950 fortran90以降でequivalence文を使うのは非推奨。 subroutine trans_endian(str,inte) character(len=1),dimension(4),intent(in) :: str integer,intent(out) :: inte inte=transfer(str(4:1:-1),0) end subroutine trans_endian 954 名前：1/2 mailto:sage [2007/11/11(日) 15:26:51 ] >>953 とこんな感じで整数をエンディアンの自動判別しながら読める。 正直TRANSFERを使うより、EQUIVALENCEの方がすっきりするよな。 F90はソースコードの見かけに流れされすぎ。 SUBROUTINE check_endian() IMPLICIT NONE CHARACTER(LEN = 4) :: ctmp INTEGER :: itmp ctmp = 'ABCD' itmp = TRANSFER(ctmp, itmp) IF ( CHAR( MOD(itmp, 256) ) == 'A' ) THEN qlittle_endian = .TRUE. ELSE IF ( CHAR( MOD(itmp, 256) ) == 'D' ) THEN qlittle_endian = .FALSE. ELSE CALL abort('endian indeterminable') END IF RETURN END SUBROUTINE check_endian !--------------------------------------------------------------------- 955 名前：2/2 mailto:sage [2007/11/11(日) 15:27:24 ] !--------------------------------------------------------------------- FUNCTION rd_int32() RESULT(ires) IMPLICIT NONE CHARACTER (LEN = 4) :: tmp INTEGER :: i, io, ires IF (qlittle_endian) THEN DO i = 1, 4 ir_record = ir_record + 1 READ(ir, REC = ir_record, IOSTAT = io) tmp(i:i) END DO ELSE DO i = 4, 1, -1 ir_record = ir_record + 1 READ(ir, REC = ir_record, IOSTAT = io) tmp(i:i) END DO END IF ires = TRANSFER(tmp, ires) RETURN END FUNCTION rd_int32 956 名前：デフォルトの名無しさん [2007/11/11(日) 17:22:06 ] >>943 コーエンといえば、これﾜﾛﾀ。 F2003が実用的な意味で出てくるのか甚だ疑問だなｗ ttp://www.nag.co.uk/market/Malcolm_Cohen_Fortran_Jubilee_Presentation.pdf The Irony of FORALL ? “world’s slowest high-performance feature” (HPC vendor 1) ? “months just to get the semantics right...” (HPC vendor 2) ? “even on massively parallel machines, it’s slower than DO” (HPC vendor 3) ? The analysis needed to eliminate the costly array temps... ...would parallelise the obvious DO loop alternative. So at best FORALL is as good as DO, usually it is slower, sometimes much slower. 11 957 名前：デフォルトの名無しさん mailto:sage [2007/11/13(火) 01:50:22 ] readで配列の行数列数を指定せずにファイルを読み込むにはどうしたらいいんでしょうか？ 「end-of-file during read, unit 5」と表示されてしまい、どうやらEND＝100のように文番号を指定する必要があるみたいなのですが、 どのように記述すればいいのか分かりません。 どなたか教えてくれませんか？ 958 名前：デフォルトの名無しさん mailto:sage [2007/11/13(火) 07:28:49 ] >>957 とりあえず切り抜けるには、 READ(nn, *, END =10) A 10 CONTINUE 959 名前：957 mailto:sage [2007/11/13(火) 11:17:16 ] >>958 出来ました！ かなり助かりました。ありがとうございます。 960 名前：デフォルトの名無しさん [2007/11/13(火) 15:23:03 ] はじめまして、初めて書き込みするものです。 Ｃ言語のプログラムをfortran90で書き換えようとしてるんですが、 C言語の do{ 処理文 }while(条件) に対応するfortranの命令文は何を使って書けばいいんでしょうか？ あと、C言語で使われている、間接演算子*aなどはfortranではどう扱えば いいのでしょうか？サブルーチンの仮引数として書いたときに、エラーが 出てしまうので、修正の仕方が分からないで困っています。 アドバイスをいただけないでしょうか？ 961 名前：デフォルトの名無しさん mailto:sage [2007/11/13(火) 21:41:03 ] >>960 C言語の do 〜 while にそのまま当てはまる文はない。 ループの終了条件判定を前へ持ってきても良いのなら DO WHILE (条件) 処理文 END DO 条件判定が後ろにないとダメな（最低１回は実行する）場合は DO 処理文 IF (.NOT. 条件) EXIT END DO 962 名前：デフォルトの名無しさん mailto:sage [2007/11/13(火) 22:27:56 ] >>960 Fortranでは引数は基本的にすべて参照渡しなので 普通に渡して仮引数を書き換えれば実引数の値も変わる。 サブルーチン側で勝手に書き換えられては困る場合は 仮引数に INTENT(IN) 属性を付けておけば書き換えられない ようになる。というかコンパイル時にエラーにしてくれる（はず） 963 名前：デフォルトの名無しさん mailto:sage [2007/11/14(水) 00:07:58 ] >>961 FortranにREPEAT〜UNTIL型の構造が入らなかったのは何故なんだろうね？ FORTRA66時代のDO〜LOOPがREPEAT〜UNTIL型だった呪いかな？ｗ >>962 F2003で引数にVALUE属性が導入されるらしい。 コンパイラによっては今でも独自拡張で出来なくも無いが。 964 名前：デフォルトの名無しさん mailto:sage [2007/11/14(水) 13:38:29 ] >>961 >>962 ありがとうございます！とても参考になりました。 教えていただいたことを使ってもう一度やってみます。 また分からないことが出てきてしまうと思うので そのときはまたよろしくお願いします。 965 名前：デフォルトの名無しさん mailto:sage [2007/11/19(月) 14:10:39 ] 保守代わりに下らない質問 77にはDO WHILE文は無いよね？ 77ではIF(〜) GO TO 〜ってやるしかないよね？ 966 名前：デフォルトの名無しさん mailto:sage [2007/11/20(火) 21:25:14 ] 失礼します、初めて書き込みさせていただきます。 以下のようなプログラムを書いたのですが（相当省略してますが）、 nの値を8以上にするとmの値が正しく出ないのですがなぜなのでしょうか。 （とりあえずaはどうでもいいです） integer :: i, j, m, n real(8) :: a read(* ,*) n do i = 1, 2n-1, 2 m = 1D0 do j = 1, i m = m * j a = a + dble(m) end do end do write (*, *) m mはiまでの階乗を表したくてこのようにしているのですが… 967 名前：デフォルトの名無しさん mailto:sage [2007/11/20(火) 23:02:03 ] >>965 DO WHILEはないが、GOTOの種類ならたくさんあるぞｗ >>966 ｍはintegerだから 基本整数型が32bitの場合、2147483648以上の数は扱えない。 968 名前：デフォルトの名無しさん mailto:sage [2007/11/20(火) 23:55:17 ] >>966 >>967氏を補足するならば、ｉｎｔｅｇｅｒで表される最大値は2**31-1だ。 m=m*jをx=x*DBLE(j)にしろ。ｘはreal(8)な。m=1の代わりに、x=1.0d0。 まぁ整数で何も考えずに階乗を求めると、途中でマイナスになって、 ２の補数表現とかコンピュータの勉強になるので、教育的には非常にいい課題だ。 969 名前：デフォルトの名無しさん mailto:sage [2007/11/21(水) 00:15:09 ] >>967>>968様 全然知りませんでした。おっしゃる通りに直してみます。 本当にありがとうございました！ 970 名前：デフォルトの名無しさん [2007/11/21(水) 02:34:46 ] Silverfrost FTN77 4.03をダウンロードしたのですが、 セットアップのところを押すと、 「C:￥WINDOWS￥SYSTEM32￥AUTOEXEC.NT．システムファイルは 　MS-DOSおよびMicrosoft　Windowsアプリケーションを実行するの 　に適していません。」 と表示されます。どうすればいいでしょうか？教えて下さい。 お願いします。 971 名前：デフォルトの名無しさん mailto:sage [2007/11/21(水) 19:50:57 ] >>970 ttp://support.microsoft.com/default.aspx?scid=kb;ja;324767 972 名前：デフォルトの名無しさん mailto:sage [2007/11/22(木) 00:50:42 ] >>970 特にこだわる理由が無いならFTN95を入れるという選択肢もある。 973 名前：970 [2007/11/22(木) 14:49:08 ] ありがとうございました。 これで家でも使えるようになりました。 974 名前：デフォルトの名無しさん mailto:sage [2007/11/23(金) 21:09:52 ] 学生時代、情報科の先生にプログラミング言語は何教わった？と聞かれて Fortran77ですと答えたら、それ古いぞ今はCだと言われた。 それからだいぶ経つのに、77を教わっている学生さんが居るとは、 教えてる奴はよっぽど77好きなんだろうな。 975 名前：デフォルトの名無しさん mailto:sage [2007/11/23(金) 21:19:46 ] 計算科学あたりでは今でも F77 教えてたりする。 今はもう C すら古いというのに。 976 名前：デフォルトの名無しさん mailto:sage [2007/11/23(金) 22:06:43 ] Visual Fortran10 with IMSL + visual studioって構成でアプリケーション作ったんだけど， IMSLとのリンクが動的にしかできないからIMSLのインストールされていないPCで動かすことができません． INCLUDE 'link_fnl_shared.h' use IMSL_LIBRARIES ならうまくいくのですが， 　　　　　　　INCLUDE 'link_fnl_static.h' use IMSL_LIBRARIES だとリンクエラーになりコンパイルできません． どなたかアドバイスいただけませんでしょうか？ 977 名前：デフォルトの名無しさん mailto:sage [2007/11/23(金) 23:21:36 ] >>974-975 真理は時間や空間には依存しないもの。 流行り廃りがあるような内容は、学問ととして伝授するに値しない。 FOTRANは、普遍的不変なる価値を持つのだ。 ５０年前のコードが動いて、実際使われているんだから、 こんな再利用可能な習得価値のある言語は無い。 978 名前：デフォルトの名無しさん mailto:sage [2007/11/23(金) 23:32:26 ] バグを抑制するための仕様が貧弱だから 保守が大変なんだよ。 979 名前：デフォルトの名無しさん mailto:sage [2007/11/23(金) 23:57:45 ] C言語のことかーーーっ！！！ 980 名前：デフォルトの名無しさん mailto:sage [2007/11/23(金) 23:58:50 ] ちょっと前はとりあえずF77でやっときゃおｋ、だったが今はさすがにF90だろ いつまでも研究室レベルのコードでいいなら別に構わんけど 981 名前：デフォルトの名無しさん mailto:sage [2007/11/24(土) 00:31:30 ] F90からいきなり入った連中は、変数命名の伝統を守らなかったりするので、再教育が面倒。 規格を硬くして厳格に守るべき部分もあるが、自由度を大きくしたまま慣習法によって厳格に守る部分 があってもいい。いざというときだけ破格で逃げる。とにかくバランスを間違えてはいけない。 英米系の慣習法と独仏系の制定法の伝統にもよるのかもしれないが・・・・ 情報学科よりの連中は、PASCALがなぜ廃れたかの教訓を得ていない。 それはそうと本来情報教育の講座が、しっかりした教材を用意すべきだと思うのだが、 連中はFORTRANとかを無闇に馬鹿にして、すぐに忘れ去られる新規言語に飛びつくからﾜﾛｽﾜﾛｽ まぁ、FORTRAN全盛時代の代表的なテキストたる浦昭二のFORTRAN本も暗黙の型を 守っていないのでどうしようもないのだがｗ 982 名前：デフォルトの名無しさん mailto:sage [2007/11/24(土) 00:44:51 ] FORTRAN77 を使ってあまりの酷さに逃げたわ。 何、あの引数の数。構造体が無いから酷い。 それを避けようとしたら COMMON 変数まみれ。 最悪だ。 983 名前：デフォルトの名無しさん mailto:sage [2007/11/24(土) 01:22:57 ] >>982 FORTRANではバッチ処理で物理メモリーを限界まで使うことが多いが、 そうなると自由にメモリーを再利用できるCOMMONは非常に便利。 メモリーが貴重な時代には非常に合理的な仕組みだった。 今風に言うとガベージコレクションを自分でやれる。 その時代時代の最高のハードウェアの制約を前提として限界まで使うにはFORTRANはいい言語。 無限のメモリーやCPU速度を暗黙の前提にしている仮想マシン上の観念言語とは違う。 FORTRAN７７程度で逃げるのは根性無し。FORTRAN66で根性つけて出直せｗ 984 名前：デフォルトの名無しさん mailto:sage [2007/11/24(土) 10:01:32 ] プログラマにはそんな無駄な根性なんて要らない。 985 名前：デフォルトの名無しさん mailto:sage [2007/11/24(土) 10:19:09 ] 必要なのはデスマを生き抜く体力のみ。 986 名前：デフォルトの名無しさん mailto:sage [2007/11/24(土) 19:59:15 ] Fortran77 は、1977年以前に生まれた人間しか使ってはいけないそうな 987 名前：デフォルトの名無しさん mailto:sage [2007/11/25(日) 14:17:12 ] 楽勝で使えるぞな 988 名前：デフォルトの名無しさん mailto:sage [2007/11/26(月) 00:39:37 ] そろそろ次スレ用のテンプレについて議論しよう。 989 名前：デフォルトの名無しさん mailto:sage [2007/11/26(月) 12:29:10 ] テンプレって無視されるのが普通だからどうでもいいような気もするけどね。 コンパイラのリンクとか宿題丸投げ用テンプレでも入れる？ F77やF90の簡単な説明もあると良いかな？自分が習ってるのがどっちか分からない子用に 990 名前：デフォルトの名無しさん mailto:sage [2007/11/26(月) 23:22:19 ] >>989 たしかにF77とF90以降の区別は重要だな。質問のときにはっきりさせろと。 あとは、文字列を数字に直すというFAQかな？ 991 名前：デフォルトの名無しさん mailto:sage [2007/11/27(火) 11:33:52 ] テンプレ議論するには残り少なすぎ・・・ まぁ内部ファイルの質問は多いから入れとこうか 992 名前：デフォルトの名無しさん mailto:sage [2007/11/28(水) 11:46:42 ] >>995が次スレを立てるという事で 993 名前：デフォルトの名無しさん mailto:sage [2007/11/30(金) 07:37:24 ] F 994 名前：デフォルトの名無しさん [2007/11/30(金) 08:34:03 ] バイナリデータを書き込むと先頭に余分な6バイトがついてしまいます。 純粋なデータのみを書き込みたいのですが、どのコマンドを使えばよいのでしょうか？ 995 名前：デフォルトの名無しさん mailto:sage [2007/11/30(金) 09:40:05 ] >>994 エスパー力が試されているな。 実行環境とどんな型のデータをどんな風に 書き込んでいるのか分からない。 昔Compaq/AlphaでVisual Fortranを使った時の話 open(io,file="hogehoge",access="direct",form="unformatted",recl=2,statu="new") としてバイナリファイルを作った時、2バイトずつ書き込もうとしても、 コンパイル時にオプション(-assume byterecl)を与えないと8バイトずつになってしまった。 余分な6バイトがつくというのは、こういう状況か？ 996 名前：995 mailto:sage [2007/11/30(金) 09:40:56 ] ×statu="new" ○status="new" 997 名前：デフォルトの名無しさん mailto:sage [2007/11/30(金) 11:29:16 ] >>994 コンパイラ・処理系などは何だ？ 998 名前：デフォルトの名無しさん mailto:sage [2007/11/30(金) 14:00:28 ] くだすれFORTRAN（超初心者用）その3 pc11.2ch.net/test/read.cgi/tech/1196384126/ 999 名前：デフォルトの名無しさん mailto:sage [2007/12/01(土) 12:56:23 ] 1000 名前：デフォルトの名無しさん mailto:sage [2007/12/01(土) 12:57:05 ] ( ´ー｀) 1000 ! 1001 名前：１００１ [Over 1000 Thread] このスレッドは１０００を超えました。 もう書けないので、新しいスレッドを立ててくださいです。。。

---

---

[ 新着レスの取得/表示 (agate) ] / [ 携帯版 ]

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