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

1 名前：デフォルトの名無しさん [2006/01/24(火) 09:48:23 ] このスレッドは、他のスレッドでは書き込めない超低レベル、  もしくは質問者自身何が何だが分からない質問を勇気を持って書き込むスレッドです。  FORTRAN使いが優しくコメントを返しますが、  お礼はFORTRANの布教と初心者の救済をお願いします。  181 名前：デフォルトの名無しさん [2006/02/19(日) 10:26:10 ] COMMON /M/!,N,LZ COMMON /P/P,GAMMA,GPH とゆうCOMMON文があるのですがCOMMON文自体を理解していなく 何の処理をしているのかいまいちわかりません ｆ９０ではMODULE文にできるらしいのですがどうのように 書き直せばいいのかどなたかご教授お願いします。 182 名前：デフォルトの名無しさん mailto:sage [2006/02/19(日) 12:47:24 ] >>181 COMMONブロック P の先頭から順に実数型変数 P, GAMMA, GPH が割り当てられる。 （同じプログラム単位中でP, GAMMA, GPHの定義がない場合） COMMONブロック M は変数名の並びがコメントアウトされているので無効？ （エラーになる？） モジュールにするとこんな感じ module hoge implicit none private !integer, public, save :: N, LZ real, public, save :: P, GAMMA, GPH end module hoge 但し、他のプログラム単位中でCOMMONブロックM, Pが他の用途に使われていたり P, GAMMA, GPHの並びが違っていたりすると正常に動作しないことがある。 183 名前：デフォルトの名無しさん [2006/02/19(日) 16:57:35 ] COMMON文はグローバル変数としての機能もあるが、その概念を外延すると 理解しにくいこともある。むしろ同じメモリー領域の手動割付と考えると分かりやすい。 昔メモリーが少なくて貴重だった時代には、手動で半動的にメモリー割付を してメモリーを有効活用していた。こういう場合同じCOMMONブロックに全然 対応の付かない変数が書かれていたりする。 COMMONブロックでは、書いた変数の名前の順に先頭から順番にメモリーが 割り振られてゆく。 COMMON /P/P,GAMMA,GPH この場合Pという名前のブロックに、（P,GAMMA,GPHを単精度４バイトとして） P,GAMMA,GPHに対応する４バイトの領域が３つとられる。 このコモンブロックを別のサブルーチンでは、COMMON/P/A,B,C　と受けることも 出来る。この場合P,GAMMA,GPHの内容はA,B,Cと一致する。（変数の番地が 一緒になる） またこのコモンブロックを別の形で受けてもいい。 COMMON/P/A(3) で受けた場合、P,GAMMA,GPHの内容はA(1),A(2),A(3)に対応する。 重要なのは順番で変数名ではない。これはMODULEの場合と異なる。 COMMON /P/GAMMA,GPH,P で受けたとすると、元のP,GAMMA,GPHはそれぞれGAMMA,GPH,Pに対応する。 このように違う名前でCOMMMONを受けていなければ、単純にMODULEに直せる。 そうでない場合は慎重に作業する必要がある。 184 名前：デフォルトの名無しさん [2006/02/19(日) 23:37:30 ] この式のプログラムできませんか？さっぱりわからんのです。。。 F(P)=(9/(128*PAI**2))*(1/P)|EXP(-P*(1/TAN(X/2)))*(1+COSX**4)dX [|]はインテグラルで、範囲は０からπ/２です。シンプソンのプログラムはあるんですが、もうさっぱりです！ どなたかお願いします。。。 185 名前：デフォルトの名無しさん [2006/02/20(月) 02:13:21 ] >>184　 tan(X/2)=1/Yとおくと積分範囲は∞〜１ dX=-2/(1+Y^2)dY cos(X)=(Y^2-1)/(Y^2+1) 1/(1+cos(X)^4) = (Y^2+1)^4/(Y^8+6Y^4+1) 代入して (1/p)∫exp(-p/tan(X/2)) / (1+cos(X)^4) dX = -(2/p)∫(Y^2+1)^3 exp(-pY) / (Y^8+6Y^4+1) dY 積分範囲を１〜∞に反転させれば頭のマイナスが消える。 とりあえず原点付近の１／ｔａｎの発散を避けるようにしてみた積分範囲が無限大にまでになるが、 被積分項は急速に小さくなるので適当なところで打ち切れるだろう。 適当なので間違ってたらスマソ。 186 名前：デフォルトの名無しさん [2006/02/20(月) 12:18:04 ] fortranは、コマンドライン引数使って実行できないのですか？ 187 名前：デフォルトの名無しさん mailto:sage [2006/02/20(月) 12:23:48 ] >>186 iargc() getarg 188 名前：デフォルトの名無しさん mailto:sage [2006/02/20(月) 12:35:30 ] >>186 規格に入ったのはFortran2003から、Fortran95までの規格には存在しない。 しかしたいていはベンダー拡張ルーチンでできることが多い。 マニュアルを調べてみるといい。POSIX互換の形式をとることが多いので、>>187が 書いたようなコマンド名になることが多い。 189 名前：デフォルトの名無しさん mailto:sage [2006/02/20(月) 15:16:01 ] >>187, 188 なるほど、勉強になります。ありがとうございます。 190 名前：デフォルトの名無しさん [2006/02/20(月) 18:40:38 ] >>182 183 ありがとございます m(__)m 191 名前：デフォルトの名無しさん [2006/02/20(月) 23:15:53 ] VAXVMS上で動く，Fortran IVで書かれたプログラムを Fortran77以降のものに自動変換してくれるプログラムとか 無いですか？ 192 名前：デフォルトの名無しさん mailto:sage [2006/02/21(火) 00:07:10 ] >>191 見たこと無い。 193 名前：デフォルトの名無しさん mailto:sage [2006/02/22(水) 13:44:48 ] 宿題を頼んだ人は事後報告をするように。 194 名前：デフォルトの名無しさん mailto:sage [2006/02/22(水) 19:59:20 ] >>193 宿題に困るレベルの人が総括（事後報告）できると思う？ 「ありがとうございます」で納得してあげてちょ。 195 名前：デフォルトの名無しさん mailto:sage [2006/02/23(木) 00:17:12 ] 解答もらったらそれっきり音沙汰ない奴もおるからのぉ。 196 名前：デフォルトの名無しさん [2006/02/24(金) 01:07:22 ] さて、遠慮せず次の宿題を持ってきたまへ。 197 名前：デフォルトの名無しさん mailto:sage [2006/02/24(金) 01:31:52 ] 宿題たまってるけど眠いのでもう寝ます。 お休みなさい。 198 名前：デフォルトの名無しさん [2006/02/24(金) 07:10:10 ] Fortranで書かれた数値計算のプログラムをC#に移植しようとしています。 至る所で引数に配列の要素数をもつ関数が定義されているのですが、 C#のGetLengthやLengthをつかってその都度判定するような仕様にしては問題あるのでしょうか。 つまり、GetLengthなどの関数を実行する時間を節約するために要素数を引数で与えているのでしょうか。 それともFortranでは配列数を取得するような関数が無いためなのでしょうか。 199 名前：デフォルトの名無しさん [2006/02/24(金) 12:39:09 ] >>198 FORTRNでは、FORTRAN９０になるまで配列の大きさを知る命令は存在しない。 FORTRAN９０で関数SIZEが導入された。 FORTRAN７７以前では、配列の大きさは副プログラム（サブルーチン、関数）に、 引数として明示的に渡すか、COMMON等で渡すか、あるいは”知っている”もの としてサブルーチン側で勝手に使うかした。 FORTRANでは変数はアドレス渡しであり、FORTRAN７７以前の書式を使うと 配列であっても渡しているのは先頭要素のアドレスだけで、配列の大きさもの情報も 次元も渡していない。（FORTRAN９０で導入された新しい書式を使うと配列の サイズや次元の情報も渡される。） ここで気をつけなければならないのは、FORTRAN６６時代のプログラムでは、 サブルーチン側の引数の仮配列の大きさを実際の長さとは無関係に１０などにして 内部では好き勝手な長さとして使用する書き方が標準的に採用されていたことだ。 （FORTRAN77では、サイズを＊で定義するやり方に相当するものかな？） こういう書法の伝統を引いている場合、>>198の移植方法では大混乱をきたすだろう。 200 名前：199 [2006/02/24(金) 12:53:08 ] 参考　（ｔabが消えるので７カラム目までの空白は脳内で補完してください） ＊FORTRAN66時代にありがちなプログラム REAL A(20, 20), B(200) CALL ZERO(200, B) CALL ZERO(400, A) CALL ZERO(20, A(1, 2)) CALL ZERO(20, A(1, 3)) CALL ZERO(20, A(1, 4)) END SUBROUTINE ZERO(N, X) REAL X(10) DO 10 I = 1, N X(I) = 0.0 10 CONTINUE END これは昔よく使われたやり方。サブルーチンに渡されるのは配列の先頭要素の アドレスだけなので、次元が整合していなくても問題はない。また配列の途中 だけを渡して部分的に利用しても問題はない。（FORTRANの配列はCとは 逆順で列方向に対して連続的） CALL ZERO(200, B)　これは素直な呼び方。 CALL ZERO(400, A)　これは２次元配列を１次元的に還元して渡している。 CALL ZERO(20, A(1, 2))　これは２次元配列の２列目だけを渡している。 CALL ZERO(20, A(1, 3)) この呼び方は列ベクトルの使用法として良く使われる。 CALL ZERO(20, A(1, 4))　同様に３，４列目を渡している。 ＊FORTRAN77時代 SUBROUTINE ZERO(N, X) REAL X(N) または REAL X(*) 201 名前：199 mailto:sage [2006/02/24(金) 13:05:48 ] １９９，２００で書いたこととは無関係に一般論的な結論を言わせてもらうが、 fortranでは宣言した配列の大きさとは無関係な大きさで配列を使うことが ままあるので、引数として渡している大きさを配列の宣言長に直すのは、 とても危険である。 FORTRANのライブラリなどでやたらと配列の寸法を引数として求められるのも、 配列は先頭要素のアドレスしか渡しておらず、そのほかの情報を一切渡していない からである。 またFORTRANでは、この性質を大いにプラスに利用しているので （というかこれが無いとまともな実用的なプログラムが書けない）これを非難するのは お門違いである。 FORTRAN９０では、配列の次元や寸法も渡しているがその性で、>>200で書いた ような列ベクトルを渡す時に、余分な処理が必要になって速度が低下する。 こういう時にFORTRAN77式の先頭アドレス渡しの簡潔さが心地よく感じられる。 202 名前：198 [2006/02/25(土) 00:14:20 ] >>199 丁寧なご解説、痛み入ります。 なるほど、配列に関する色々な情報ではなく、配列のアドレスを渡しているだけなのですね。 そういえば他の箇所で多次元配列が1次元の配列として処理されていたり、 首をかしげたおぼえがあります。 このあたり、移植の際には気をつけないと駄目ですね。 ありがとうございました。ご指摘の点を踏まえて頑張ってみます。 203 名前：デフォルトの名無しさん mailto:sage [2006/02/25(土) 01:49:11 ] 10年以上前の話だけど、 なんとなく入学した大学で、入学直後からFORTRAN77の宿題を大量に課せられた。 もちろん僕にはFORTRANが一体何なのかもよくわからなかった。 とりあえず、テキストをマネしながらキーボードを叩いてはみたものの、 なんでこんな勉強しなきゃいけないのかわからなかった。 正直言って、楽しくなかった。 同学部の友達が居なくて、先生も冷たくって、僕は１人だった。 結局、大学には通わなくなり、僕はその場所を去った。 その後、僕はプログラムとは無縁の人生を送ることになった。 このスレにFORTRANの宿題を持ってくる人・・・ それに対して丁寧に解説をされている方々・・・ ログを上から読んでみたけど、なんだかスゲー楽しそうじゃん？ なんだか羨ましいよ。マジで。 あの頃、この2chのFORTRAN（超初心者用）スレがあったら・・・・。 なんて、思ってしまったけど、仮定法過去完了で人生を語るのはやめようと思う。 皆さんお元気で。 ひっそりと心の奥底で応援いたします。 204 名前：デフォルトの名無しさん [2006/02/26(日) 01:20:51 ] >>203 昨日の涙を今日の笑顔に変えて、でっかい夕日に向かって走るんだ！ それはそうとFORTRANの教科書や教育課程はろくなものが無い。 面白くなくて訳分からなくなるのが標準だ。 やさしく分かりやすく書くということを、内容を整理して書くことではなく、 本文と無関係な下手糞な挿し絵をちりばめる事だと思っている世界だｗ 205 名前：デフォルトの名無しさん [2006/02/27(月) 00:55:37 ] さぁ宿題をもってきたまえ 206 名前：デフォルトの名無しさん mailto:sage [2006/02/27(月) 13:35:32 ] >>203 僕はＦＯＲＴＲＡＮスレで助けてもらったおかげで 大学院を無事修了できました。 ここの住人の方々には本当に感謝しています。 僕もがんばって、知恵で人を助けられるような人間になりたい 207 名前：デフォルトの名無しさん mailto:sage [2006/02/27(月) 14:16:38 ] >>204 うむうむ 「作法」が執筆された時に 77 があったら良かったのに。 ratfor で見た目がちょっと違うというだけで、 初学者には敬遠されてしまうから。 これも仮定法過去完了だぬ。 208 名前：デフォルトの名無しさん [2006/02/27(月) 16:04:00 ] >>206 おめ！ これからも山崎渉を・・じゃなかったFORTRANスレを応援してくださいｗ 209 名前：デフォルトの名無しさん mailto:sage [2006/02/28(火) 13:20:06 ] >>204 『FORTRAN 入門』浦昭二　が名著だと思う。 うちの大学ではどこでも転がっている標準本。 Fortran90 への対応を希望する。 210 名前：デフォルトの名無しさん [2006/02/28(火) 19:56:27 ] >>209 オレもそれのFORTRAN77版で習ったんだけどさ、 あれが名著と言われるようでは日本もお終いだと考えるべきだと思う罠ｗｗｗ 前半例題で後半が文法という形式以外は例題もいまいちだしいいとこなしだと思うよ。 211 名前：デフォルトの名無しさん mailto:sage [2006/02/28(火) 20:08:01 ] 日本語で書かれたFortran90以降の機能が詳しく解説されてる本ってない？ 近所で探してみたけど、どれも全くの入門書止まりだったり FORTRAN77に毛が生えたような本しか見つからなくて。 212 名前：デフォルトの名無しさん mailto:sage [2006/02/28(火) 20:38:07 ] >>211 これが結構いい。 コードも比較的Fortran90の機能だけで書こうとしている。 入門 Fortran90実践プログラミング SOFTBANK BOOKS 東田 幸樹 (著), 山本 芳人 (著), 熊沢 友信 (著) ttp://www.amazon.co.jp/exec/obidos/ASIN/4890524932/249-3009705-2437946 しかし絶版だったｗｗｗｗ 213 名前：デフォルトの名無しさん mailto:sage [2006/02/28(火) 21:19:16 ] >>210 そうかな？ 文法事項を段階的にかつ網羅的に適用した、 適切な入門例題だと思うが。 科学技術計算を意識した例題が多いというのも、 FORTRAN という言語の性質上、好ましい。 >>211 そうだよな。 Fortran90 以降の名著が少ない（あっても絶版）というのは、 Fortran 市場を反映しているようで、残念でならない。 『FORTRAN 入門』浦昭二を90向けに改訂するなら、喜んで手伝うが。 214 名前：デフォルトの名無しさん mailto:sage [2006/02/28(火) 21:21:08 ] >>209-210 あの本はFORTRAN77版だけじゃないんですか？ 詳細ｷﾎﾞﾝﾇ．（"・ω・`）ﾉ 215 名前：デフォルトの名無しさん mailto:sage [2006/03/01(水) 03:02:46 ] >>214 あの本はFORTRAN６６版がある。 アマゾンとかで検索してFORTRANというだけの７７のついていないのがそれ。 それを元にFORTRAN77版を書いているので、 わりとすっきりしていたFORTRAN６６版に比べるとFORTRN７７版は いまいちな気がしないでもない。 FORTRAN90の教科書だったら、FとかLahey ELFのようにFORT７７以前との 互換性を捨て去ったバージョンで書くのがすっきりしていいと思う。Fなんかは 色々入門書が出ていたと思う。もっとも肝心のFがNAGのフロントエンドを入れた せいで、当初の過去互換性を切り捨てた簡潔な言語になるという理想から 逸脱して迷走しているのがずっこける。 www.fortran.com/F/ www.fortran.com/ FORTRAN９０で導入された新しい機能を用いれば、昔PASCALなんかで やっていた正則な構造化（モジュール化）されたプログラム書法が身につくと思うのだが、 いかんせん現実問題としてはFORTRAN77やFORTRAN６６のプログラムをポンと渡されて、 動かせと言われたりするから、まあどうしたものか分からん。 216 名前：デフォルトの名無しさん [2006/03/02(木) 15:37:49 ] PROGRAM GGGG IMPLICIT NONE INTEGER*8 :: N,K,NH,M,I REAL*8 :: GANMA2(0:10000),GANMA3(0:10000),GANMA4(0:10000) REAL*8 :: BETA2,BETA3,BETA4,SECHYP,EAMP0,ER(0:2047),EI(0:2047) REAL*8 :: COS,SIN,ERK,EIK,T0,T,PI,FRQ,FQ,H,C,CD,CDD,WL,ROU,TAU0 COMPLEX*16 :: E(0:4094),EK M=11 N=2**M NH=N/2 PI=4.0D0*DATAN(DBLE(1.0D0)) T=5.12D-12 TAU0=0.22D-12 C=3.0D+8 CD=-42.0D-1 CDD=0.0D0 WL=1.55D-6 H=5.0D-5 BETA2=-CD*WL*WL/(2.0D0*PI*C) BETA3=WL*WL*WL*(WL*CDD+2.0D0*CD)/(4.0D0*PI*PI*C*C) BETA4=0.0D0 T0=TAU0/(2.0D0*DLOG(DSQRT(DBLE(2.0D0))+DBLE(1.0D0))) EAMP0=10.0D0 217 名前：デフォルトの名無しさん [2006/03/02(木) 15:41:26 ] DO K=0,N-1 SECHYP=DBLE(1.0D0)/DCOSH((T/T0)*(DBLE(K)/DBLE(N)-DBLE(0.5D0))) E(K)=CMPLX(EAMP0*SECHYP,0.0D0) END DO DO K=0,NH-1 ROU=DBLE(K)*T0/T FRQ=PI*ROU/T0 GANMA2(K)=-H*BETA2*PI*PI*ROU*ROU/(T0*T0) GANMA3(K)=-H*BETA3*PI*PI*PI*ROU*ROU*ROU/(T0*T0*T0)*2.0D0/3.0D0 GANMA4(K)=-H*BETA4*FRQ*FRQ*FRQ*FRQ*1.0D0/3.0D0 FQ=ROU/T0 END DO DO K=0,20,1 EK=E(K)*CEXP(CMPLX(0.0d0,(GANMA2(K)))) EK=EK*EXP(CMPLX(0.0d0,(GANMA3(K)))) EK=EK*CEXP(CMPLX(0.0d0,(GANMA4(K)))) E(K)=EK END DO END PROGRAM GGGG 218 名前：デフォルトの名無しさん [2006/03/02(木) 15:47:14 ] DO K=0,20,1 ERK=ER(K)*DCOS(GANMA2(K))-EI(K)*DSIN(GANMA2(K)) EIK=EI(K)*DCOS(GANMA2(K))+ER(K)*DSIN(GANMA2(K)) ERK=ERK*DCOS(GANMA3(K))-EIK*DSIN(GANMA3(K)) EIK=EIK*DCOS(GANMA3(K))+ERK*DSIN(GANMA3(K)) ERK=ERK*DCOS(GANMA4(K))-EIK*DSIN(GANMA4(K)) EIK=EIK*DCOS(GANMA4(K))+ERK*DSIN(GANMA4(K)) ER(K)=ERK EI(K)=EIK END DO END PROGRAM GGGG 連投すみません、216,217と計算が実数か複素数かだけの違いなのですが 217ではなぜか範囲を広げるほど虚数の値がかけ離れていきます 桁落ちでもないのですが、どなたか考えられる原因を教えてくださいm(__)m 219 名前：デフォルトの名無しさん [2006/03/03(金) 12:54:26 ] >>218 範囲を広げるという意味が良くわからんので、適当に書いてみるｗ EK=E(K)*CEXP(CMPLX(0.0d0,(GANMA2(K)))) EK=EK*EXP(CMPLX(0.0d0,(GANMA3(K)))) EK=EK*CEXP(CMPLX(0.0d0,(GANMA4(K)))) この式で、EK(K)は最後に掛けたほうがいい。なぜなら、それ以外の部分は 全部位相因子なので幾ら掛けても、絶対値は１になるはずの量だから。 まずはそこでチェックしてみたらどうか？ つまり EK = EXP(CMPLX(0.0d0,(GANMA2(K)))) EK = EK * EXP(CMPLX(0.0d0,(GANMA3(K)))) EK = EK * EXP(CMPLX(0.0d0,(GANMA4(K)))) ここでABS(EK)＝１．０ｄ０のはず（誤差を抜かして）　 もう一つの式でも同様にして、その段階で虚部を比較すれば、 より限定的な情報が得られて問題解決に近づくだろう。 EK = EK(N) * EK　 （ここで後ろのEKはpahase_factorとかの変数名にしたほうがいいと思うが まぁ好みの問題だ） あとCEXPは総称名EXPに統一したほうがいいだろう。 精度の規定値が気になるなるが、KIND=8をつければ強制できるはず。 220 名前：デフォルトの名無しさん mailto:sage [2006/03/03(金) 13:52:48 ] > ERK=ERK*DCOS(GANMA3(K))-EIK*DSIN(GANMA3(K)) > EIK=EIK*DCOS(GANMA3(K))+ERK*DSIN(GANMA3(K)) とあるけど、一行目で、ERKを書き換えてしまってはダメじゃないの？ 221 名前：デフォルトの名無しさん [2006/03/03(金) 15:40:33 ] >>220 実数、虚数にわけてる計算はあっているので、ただしいはずです・・・ 222 名前：デフォルトの名無しさん [2006/03/03(金) 16:13:32 ] >>219 ありがとうございます、試してみます EK=EXP(CMPLX(0.0d0,(GANMA2(K)),8)) といった感じで精度をあわせてしまえばいいんですよね？ 223 名前：デフォルトの名無しさん [2006/03/03(金) 18:07:37 ] >>221 いや、１行目でERKが書き換わっているので、正しいとは思えない。 その2行で１つの計算のはずだから、虚部を求める時に新しくなった実部を 使ってはいけない。 これで正しいとすると、二つが一致しないのは当然。 ただこれが正しいことはまずありえないと思う。 224 名前：デフォルトの名無しさん mailto:sage [2006/03/03(金) 18:11:08 ] >>221 いや、式はあってるんだけど、1行目でreal partを計算する際に、 元の変数のreal partを書き換えてしまっているのでは？ 2行目の右辺に出てくるERKはもはや元の変数のreal partではないないと思うけれど。 225 名前：デフォルトの名無しさん mailto:sage [2006/03/03(金) 18:11:51 ] かぶったw 226 名前：デフォルトの名無しさん mailto:sage [2006/03/03(金) 18:15:59 ] 式があっているってのは、実部虚部への分け方はあっている、ってことね。 そのプログラムは変だと思うよ。 227 名前：デフォルトの名無しさん [2006/03/03(金) 22:14:06 ] >>223~226 気づかなかった；；；サンクスです。 228 名前：デフォルトの名無しさん [2006/03/07(火) 22:19:56 ] もりあげようぜ！ｗ 229 名前：デフォルトの名無しさん mailto:sage [2006/03/07(火) 23:26:03 ] 盛り上げようにも大学が春休みで宿題がない予感。 230 名前：デフォルトの名無しさん mailto:sage [2006/03/10(金) 22:52:46 ] 欠測のある平均値の求め方を教えてください。 10行ごとの各列の（欠測を除外した）平均値と欠測数を求める問題です。 時間,A,B,C 1,23.5,23.2,22.7 2,24.5,25.2,22.1 3,25.4,24,27.5 4,26.4,,22.4 5,23.5,22.4, と2000行、続きます。 231 名前：デフォルトの名無しさん [2006/03/11(土) 00:02:09 ] >>230 適当なのでもっとうまい方法もあるやも知れぬ。参考までに。 INTEL　FORTRAN　Ver.9 PROGRAM vipp IMPLICIT NONE INTEGER :: i, j, k REAL :: x(10) i = 0 DO i = i + 1 x = 0.0 READ(9, '(i5, 10F15.5:)', ADVANCE = 'NO', EOR = 98,END = 99) k, (x(j), j = 1, 10) 98 PRINT *, j, ":", SUM(x(1:j)) / REAL(j) END DO 99 PRINT *, "EOF" STOP END PROGRAM vipp 実行結果 C:\Documents and Settings\All Users\Documents\FORTRAN\vip\Debug>vip 3 : 23.13333 3 : 23.93333 3 : 17.63341 3 : 16.26667 3 : 15.30000 EOF 232 名前：２３１ mailto:sage [2006/03/11(土) 00:09:25 ] >>230 ごめｗ　列ごとの平均かｗ　 風呂から上がったらまた考えるｗ 233 名前：230 mailto:sage [2006/03/11(土) 00:25:20 ] >>232 トライ、ありがとうございます。 温度なので、23度あたりになるはずです。 5列分の平均ですが、 A 24.66 , B 23.7 , C 23.675 内部ファイルとLEN 関数をつかって、 うまくできないかとチャレンジ中です。 234 名前：デフォルトの名無しさん [2006/03/11(土) 00:28:05 ] くたばれFORTRAN 235 名前：デフォルトの名無しさん [2006/03/11(土) 00:51:46 ] がんばれFORTRAN 236 名前：２３１ mailto:sage [2006/03/11(土) 00:52:36 ] >>233 単純に思いつく方法は １．手抜き方式、０．０が読み取られたときを欠損値とするというやり方 ２．内部ファイルに読み込んで、INDEX関数で”，”を区切りとして切り取るやり方 汎用性からみると２かな？ 第３の方法を求めて FORMAT文がF95で進化したので、空データの読み取りが検知できないか今調べているｗ 237 名前：230 mailto:sage [2006/03/11(土) 00:56:23 ] >>233 使用コンパイラ　Compaq VFでは、 空白を欠測ではなく0と認識してしまったので、ダメでした。 文字列として読み込んで、^@ の処理を考えます。 238 名前：デフォルトの名無しさん mailto:sage [2006/03/11(土) 01:28:02 ] こんな感じッスかね？ program hoge implicit none integer, parameter :: N = 10, NT = 2000 real, parameter :: NG = -1.0E10 real :: A(N), B(N), C(N), t, aa, ab, ac integer :: ka, kb, kc, i, j, ios character(len=80) :: tmp open(unit=11,file="huge.dat",status="OLD",blank="NULL") loop_j: do j = 1, NT/N A = NG; B = NG; C = NG loop_i: do i = 1, N read(unit=11,fmt="(A)",iostat=ios) tmp if (ios < 0) then ! エラー処理は省略... exit loop_j end if read(unit=tmp,fmt=*,iostat=ios) t, A(i), B(i), C(i) end do loop_i ka = count(A <= NG); aa = sum(A, mask=(A > NG)) / (N-ka) kb = count(B <= NG); ab = sum(B, mask=(B > NG)) / (N-kb) kc = count(C <= NG); ac = sum(C, mask=(C > NG)) / (N-kc) print *, "欠測 = ", ka, kb, kc print *, "平均 = ", aa, ab, ac end do loop_j close(11) end program hoge 239 名前：デフォルトの名無しさん mailto:ｹﾞｪｪ [2006/03/11(土) 01:33:55 ] test 240 名前：230 mailto:sage [2006/03/11(土) 01:56:46 ] >>238 ありがとうございます！ 無事に解決できました。 10個とも欠測の時間帯の平均値が NaN と表記されて、 いたれりつくせりです。 count(A <= NG)　と　mask=(A > NG) がミソですね。 0割になりそうでならないのが少し不思議です。 241 名前：２３１ [2006/03/11(土) 02:05:07 ] 出遅れたかｗｗｗｗｗｗｗｗｗｗｗ PROGRAM vipp IMPLICIT NONE INTEGER :: i, j, k CHARACTER(LEN = 136) :: buff REAL :: a, x(10) = 0.0 INTEGER :: n(10) = 0 DO READ(9, '(a)', END = 99) buff k = INDEX(buff, ',') READ(buff(:k), '(i10)') j buff = buff(k + 1:) DO i = 1, 10 k = INDEX(buff, ',') IF (k == 0 .AND. LEN_TRIM(buff) == 0) EXIT ! end by "," IF (k /= 1) THEN ! k == 1 is missing data : ",," READ(buff, '(F15.0)') a x(i) = x(i) + a n(i) = n(i) + 1 END IF buff = buff(k + 1:) IF (k == 0) EXIT ! end of record END DO END DO 99 PRINT *, "EOF" PRINT *, n PRINT *, x / REAL(n) STOP END PROGRAM vipp >>241 242 名前：２３１ mailto:sage [2006/03/11(土) 02:06:24 ] C:\Documents and Settings\All Users\Documents\FORTRAN\vip\Debug>vip EOF 5 4 4 0 0 0 0 0 0 0 24.66000 23.70000 23.67500 NaN NaN NaN NaN NaN NaN NaN 243 名前：230 mailto:sage [2006/03/11(土) 02:26:44 ] >>241 わかりやすい注釈までいれていただいて、ありがとうございます。 k = INDEX(buff, ',') buff = buff(k + 1:) LEN_TRIM(buff) == 0 など、文字列処理の強力さはさすがですね。 配列手続といい、文字列処理といい、 f90 の強力さに、あらためて惚れなおしました。 244 名前：230 mailto:sage [2006/03/11(土) 02:38:50 ] >>241 今、調べなおしてみると、 LEN_TRIM 以外は、f77 の機能のようですね。 f77 の文字列処理の先見の明に驚いてしまいます。 245 名前：２３１ mailto:sage [2006/03/11(土) 02:39:14 ] >>243 場合わけはあまりエレガントでなかった。もう少し冗長でも分かりやすくすべきかも。 buff = buff(k + 1:) これで長さが足りない分に空白が入るかどうか確かめておいた。 ちゃんと空白が入っているようなので、ごみが尻のほうに残ることは無い。 文法書にもそう書いてあった気がするのだが、記憶がアイマイミー。 246 名前：230 mailto:sage [2006/03/11(土) 02:54:36 ] >>245 INDEX を用いた切り分けや、テクニカルな文字列処理の手本 を示していただいただけでも、感謝感激です。 本当にありがとうございました。 247 名前：２３１ mailto:sage [2006/03/11(土) 03:38:39 ] >>246 喜んでもらえてうれしいお（＾ω＾）ﾉｼ >>238も勉強になった。漏れは固定FORMATの愛好者なのでCVS形式はほとんど使ったこと無い。 自由書式の時にCVS形式の空データは読み飛ばすとは知らなかった。０が入ると思っていた。 F77はINDEX関数がとても便利だったが、TRIM関数が無かったせいで初心者時代に苦労した。 ある程度習熟してからは、不定長の文字引数で尻から先頭に向かって空白以外が出るまで 数える関数などを自分で用意できるようになったが、そこに到達するまではかなり時間がかかった。 内部ファイルの使い方も共同研究の外国人の書いたプログラムを読んでいて初めて知った。 それまでこんな便利なものがあるとはまったく知らなかったｗ 後になって改めて文法書を読むと、文字変数を入力装置にできるとはさりげなく書いてあるのだが、 それがとても便利なものだということはまず分からない。 この辺のノウハウを徒弟制度や口伝に頼っているFORTRAN界の未来は暗いなと思ったりしてｗ 248 名前：デフォルトの名無しさん [2006/03/11(土) 08:27:30 ] ところでFortranに分割コンパイルっていつ入るんですか？ 249 名前：デフォルトの名無しさん mailto:sage [2006/03/11(土) 10:30:55 ] >>248 単純に複数のプログラム単位をコンパイル＆リンク出来るという意味でなら FORTRAN IIで副プログラム呼び出しとBinary Symbolic Subroutine Loaderが 使えるようになった1958年から。 複数のプログラム単位のコンパイルを自動化出来るという意味でなら FORTRAN IIモニタの出来た1959年から。 近代的なモジュールの機能が使えるようになったという意味でなら Fortran 90の規格が出来た1991年から。 ということでよろしいかな？ 250 名前：デフォルトの名無しさん [2006/03/12(日) 21:48:13 ] 非線形連立代数方程式をとくFORTRANのコードをC++に移植しています。 解の一部がどうしても両プログラムで一致せず、 丸一日かけて色々と試行錯誤したところ、Fortranでreal型を使っていたところ、 c++でdouble型にしてしまっていたことが原因でした。 このような場合、c++への移植の際にもfloat型を使うなどしてもとのコードと精度を合わせた方がよいのでしょうか。 なんとなくより大きな範囲を扱える変数型の方がよさそうな気がするのですが、 移植元が倍精度を使っていないので迷っています。 移植元のコードは1980年頃に書かれたもので、計算時間との兼ね合いから単精度を選択したと捉えたほうがよいでしょうか。 また、計算誤差の発生しないdecimal型があると知りました。 一般に数値計算をする場合はこのような型を使うのでしょうか。 精度が高い代わりに計算速度が遅いと聞きましたので、 一般にはどのような基準で使い分けて折り合いをつけていくのか興味があります。 すみませんが、どなたかご存知の方がいらっしゃいましたらお教えいただけますか。 よろしくお願いいたします。 251 名前：デフォルトの名無しさん [2006/03/12(日) 23:43:23 ] >>250 精度が出てるのを確かめられたならdouble型にすべし。 floatにする必要はないと思う。 C++の組み込み型にdecimalがあるとはきいたことないが、 だた10進数の内部表現ができるってだけで(0.1とかが正確に表せるのかな？) 数値計算自体に誤差が出ないわけではない。はず・・・たぶん 252 名前：デフォルトの名無しさん [2006/03/13(月) 11:13:40 ] ところでFORTRANにリアルタイムイベントハンドラはいつ入るんですか？ 253 名前：デフォルトの名無しさん [2006/03/13(月) 13:47:57 ] >>250 1980年代のコードである場合、浮動小数点には今とは違う事情がある。 現在の浮動小数点はIEEEの形式にほぼ統一されているが、かつてはメーカーごとに 実数の内部形式は異なっていた。そのために、同じワード長でも仮数部により多くの ビットを割り当てているものもあれば、指数部により多くのビットを割り当てているもの もあって、機種によって有効桁の範囲が違っていた。 （IBMは仮数部に、CRAYは指数部により多くのビットを割り当てていた。 さらに言えば、ＮＥＣのＡＣＯＳのあるシリーズのように１ワード＝３６bitとか、 CDCのいくつかのシリーズのように１ワード５６？bitとかでワード長自体が 違うものもあった。） というわけで、書かれた対象機種が分かれば必要精度に関するヒントになる。 ただ数値演算プロセッサが倍精度で作られている以上 基本的に倍精度にしておくのがよろしいと思う。 まぁSSEで高速化を優先するならまた話は別だが・・・ 254 名前：250 [2006/03/13(月) 20:34:42 ] >>251, 253 ありがとうございます。 頂いたアドバイスに従って とりあえず、倍精度で移植を進めてみようと思います。 >現在の浮動小数点はIEEEの形式にほぼ統一されているが、かつてはメーカーごとに >実数の内部形式は異なっていた。 そういえば、コンピュータの精度に応じてコメントアウトすべしというような注釈のついた関数がありました。 そういうことだったんですね。勉強になります。 255 名前：デフォルトの名無しさん [2006/03/14(火) 01:24:29 ] >>250 DECIMAL型について言えば、あれは商業用の機能で、数値計算的にはあまり 利益がないと思う。特に計算スピードはがくっと落ちるはずである。 専門ではないから正確ではないが、商業の世界では利息で生じた小数点以下の 端数をあるところで切り捨てる（または四捨五入？）する慣例がある。浮動小数 点では２進法を使っている関係で丸め誤差から小数点以下の端数の切り捨てで １０進法の時と微妙なズレが生じ、大量の金額に対して計算すると慣例とのズレ が無視できなくなる。これに対する対策がdecimal型すなわち１０進型である。 一般的には、１桁ごとに１０進法の数字を１個割り当てて計算するので極めて 計算効率が悪い。１桁につき最低４ビット必要。同じビット数で表される数の 範囲が狭くなってしまう。 （確か８０８７コプロにもBCD演算機能があったはずだが、固定小数点だったような 気がする。） 数値計算的には、倍精度、倍々精度という方向に行くことはあっても、 １０進演算の方向に行くことは特殊な場合を除いてないと思う。 実際FORTRANで１０進演算はできないが、利用者からの不満は聞かない。 一方COBOLやPL/Iなどでは１０進演算が出来るようになっている。 256 名前：デフォルトの名無しさん [2006/03/14(火) 01:53:35 ] >>254 ＞そういえば、コンピュータの精度に応じてコメントアウトすべしというような ＞注釈のついた関数がありました。 昔は精度に関する機種依存性が結構あって移植の問題になっていた。 Fortran90で数値の精度に関する関数が大量に導入されたのは、 この辺でみんな苦労していたからだと思う。ただ９０年代に入ってから IEEE形式に集約されて行ったので、ややありがたみが減った。 最近は０割りが発生しても強行突破して計算を続けてくれるので便利だと思う。 投機的実行の為にわざわざIF文で排除したはずの０割りを実行して糸冬了とか 脳の血管切れそうな事態が起きなくなった。 Fortran2003（一部はFortran95の拡張規格）でIEEE754に関する様々な関数が 導入されているので、将来はよりマニアックな要請に対応できるはずである。 257 名前：デフォルトの名無しさん [2006/03/17(金) 01:29:56 ] PROGRAM TEST IMPLICIT REAL*8 (A-H,O-Z) IMPLICIT INTEGER (I-N) * DIMENSION G(50) * WRITE(*,*)'個数(入力数+1)' READ(*,*)N * WRITE(*,*)'整数' DO 10 I = N, 0, -1 READ(*,*)G(I) 10 CONTINUE * DO 500 I = N, 0, -1 WRITE(*,*)I,G(I) 500 CONTINUE * DO 510 J = N, 0, -1 WRITE(*,*)J,G(J) 510 CONTINUE * END これ実行するとＧ(０)の値がちゃんと出力されないんですが、何故ですか？ ＤＯ５００とＤＯ５１０で、出力が違う値で出てきてしまいます。 あまりにも簡単すぎる質問で申し訳ないですが、誰か教えてください。 258 名前：257 [2006/03/17(金) 01:31:55 ] frtで実行するとこんなのが出ます。 ↓ jwe0019i-u The program was terminated abnormally with Exception Code EXCEPTION_A CCESS_VIOLATION. error occurs at _MAIN__ line 19 loc 0040108c offset 0000008c _MAIN__ at loc 00401000 called from o.s. error summary (Fortran) error number error level error count jwe0019i u 1 total error count = 1 259 名前：デフォルトの名無しさん mailto:sage [2006/03/17(金) 01:55:44 ] >>257 配列の使い方を復習しましょう。 ttp://www.g.dendai.ac.jp/lecture/fortran/fort07.html 260 名前：257 [2006/03/17(金) 02:28:35 ] >>259 なるほど・・・初歩の初歩ですね。 とてもよくわかりました。 助かりました。ありがとう。 261 名前：デフォルトの名無しさん [2006/03/17(金) 02:29:05 ] >>257 数学科や計算機学科での自然な数は０から始まるのだが、 FORTRANの自然数は１から始まる古典的な自然数体系を取っている。 したがってG(50)と宣言したときはG(1)〜G(50)までの領域が取られており、 G(0)はout of rangeになる。 エラーメッセージにACCESS_VIOLATIONと出ているが、これはは大抵 配列のはみ出しと思っていい。 G(0)がどうしても使いたければG(0:50)という形で宣言すればいいが、 まあ素直に１から始まる配列になれたほうがいいだろう。 エラーメッセージから見る限り富士通のコンパイラを使っていると思うが、 デバッグ用のオプションに配列はみだしをチェックするサブスクリプトチェックの項目が あるだろうから、デバッグ段階ではそれをオンにしておけばはみだしが発生したとき そのソース行を教えてくれる。エラーが取れたら、チェックを外せばよい。 262 名前：259 mailto:sage [2006/03/17(金) 11:56:35 ] >>261 めちゃめちゃ親切に説明して頂いて助かります。 富士通のコンパイラまでわかるなんて、凄いです。 エラーメッセージの意味が理解できなかったので、まだまだでした…精進します。 263 名前：262 mailto:sage [2006/03/17(金) 11:58:40 ] 自分、257でした。 259さん、すみません。 264 名前：http://www.vector.co.jp/soft/win95/util/se072729.html mailto:http://www.microsoft.com/japan/windowsxp/64bit/default.mspx [2006/03/18(土) 19:33:30 ] 64bitに対応したトリップ検索プログラムありますか？ TextSS　のWindowsXP(Professional)64bit化おながいします もしくは64bitにネイティブ対応したテキスト置換ソフトありますか？ 265 名前：デフォルトの名無しさん [2006/03/18(土) 23:10:21 ] この板でも爆撃喰らうんだな。ＶＩＰでは日常茶飯事だがｗ 266 名前：デフォルトの名無しさん [2006/03/26(日) 07:13:26 ] xxxx01.txt xxxx02.txt .... みたいな連番になっているファイルを順番に読み込むようにしたいのですが、 読み込むファイル名の部分をどのような書式にすれば良いのでしょうか？ （（文字列+整数）をまとめて文字列として認識する方法が分かりません） 267 名前：デフォルトの名無しさん mailto:sage [2006/03/26(日) 11:30:18 ] >>266 内部ファイルを使うのが一番手っ取り早いだろう。 CHARACTER*10 FNAME INTEGER I DO I = 1, 20 　WRITE(FNAME,'(A, I2.2, ''.txt'')') 'xxxx', I 　OPEN(10, FILE=FNAME, STATUS='OLD') 　： 　ファイルの読み込み 　： 　CLOSE(10) END DO 268 名前：266 [2006/03/26(日) 20:41:47 ] >>267 どうもありがとうございます。助かります。 申し訳ないですが、もうひとつ質問、 連番になっていなくても、ディレクトリにあるファイル全部を 読み込むような手段はあるでしょうか？ 269 名前：デフォルトの名無しさん mailto:sage [2006/03/26(日) 23:22:58 ] >>268 Fortran95 までの標準仕様ではディレクトリの操作や検索は出来ない。 しかし、大抵の処理系にはファイル操作関係の拡張手続きが付いているので それらを使えば可能。 以下は Compaq Visual Fortran での例： 　use DFLIB 　implicit none 　character(len=100) :: dir, files, fname 　integer(kind=4) :: handle, fnlen 　type(FILE$INFO) :: info 　： 　dir = "D:\DataDir\" 　files = trim(dir) // "*.txt" 　handle = FILE$FIRST 　do 　　fnlen = GETFILEINFOQQ(files, info, handle) 　　if ((handle == FILE$LAST) .or. (handle == FILE$ERROR)) exit 　　fname = trim(dir) // info%NAME 　　open(unit=10, file=fname, status="OLD") 　　： 　　ファイルの読み込み 　　： 　　close(10) 　end do 　： 詳しくは処理系のマニュアル等を参照のこと。 270 名前：266 mailto:sage [2006/03/27(月) 20:27:30 ] >>269 どうもありがとうございました。 連番に直してから読み込むことにします。 271 名前：デフォルトの名無しさん mailto:sage [2006/04/05(水) 17:37:56 ] format 文について、教えてください。 配列Aのサイズ変化にあわせて、 format を可変にしたい場合どうすればよいのでしょうか？ program main integer, parameter :: NUM=4 real*8 A(NUM) integer n do n=1,NUM A(n)=0.1*n end do ! write (6,'(1x, NUMf5.1)') (A(n),n=1,NUM) ! これはエラー write (6,'(1x, 4f5.1)') (A(n),n=1,NUM) ! これなら通る end 272 名前：デフォルトの名無しさん mailto:sage [2006/04/05(水) 18:08:54 ] write文で DO型並びを使わず愚直に DO文をかく。 write(6, '(1x, $)') do n=1, NUM write(6, '(f5.1, $)') A(n) end do write(6, '(/)') FMT記述子を文字列変数にするという手もあるが、 こちらの方がシンプル。 273 名前：271 mailto:sage [2006/04/05(水) 20:36:03 ] >>272 ありがとうございます。無事にできました。 改行制御編集記述子と斜線編集記述子を組み合わせれば、 できるのですね！ 274 名前：デフォルトの名無しさん mailto:sage [2006/04/05(水) 20:58:19 ] 仕様上はFORMAT書式に対応する項目がなくなった時点で書式制御が終了するので 出力する項目に対して多めに指定してやっても良い（ちょっと強引だが） write(6, '(1X, 100F5.1)') (A(n),n=1,NUM) 処理系依存だが、<> でくくった中に整数式を書ける処理系もある。 write(6, '(1X, <NUM>F5.1)') (A(n),n=1,NUM) ついでに、コロン編集記述子を使えば入出力項目がなくなった時点で 以降の書式を無視して終了させられる。 write(6, '(1X, 100(F5.1, :, ", "))') (A(n),n=1,NUM) 275 名前：デフォルトの名無しさん [2006/04/05(水) 21:01:15 ] >>272 ＄　指定子ってFortranの規格外でないかい？ 276 名前：271 mailto:sage [2006/04/05(水) 21:13:08 ] >>274-275 情報、ありがとうございます。 過大定数の仮代入やコロン編集記述子は、アリですね！ <> でくくった中に整数式を書ける処理系 Compaq Visual Fortran ではできました！すごい！ 処理系依存はできるだけ避けたいと思っているので、 考えどころですね。 改行制御編集記述子は規格外でしたか！ なかなか使い勝手のよさそうな機能なのに。 277 名前：デフォルトの名無しさん mailto:sage [2006/04/06(木) 10:53:52 ] $の代わりにこっちでどうかな？ write(*,'(f5.1)',advance='no') A(n) 278 名前：271 mailto:sage [2006/04/06(木) 12:35:24 ] >>277 ありがとうございます。 こちらでもできました！ 停留入出力という技もありましたね！ F90 の書式仕様の理解は、手を抜いていましたが、 まじめに取り組むべきだと思いました。 279 名前：デフォルトの名無しさん [2006/04/07(金) 19:33:01 ] 教えてください！！ 初めて書き込みさせていただきます。 Ｆｏｒｔｒａｎのプログラムの課題が出ているのですがわからなくて・・・ 簡単な内容なんだと思うのですが、 みなさんのお力貸していただけないでしょうか？？ 内容は 配列データａ(i) (i=0,Num-1)　が与えられたとき、 このデータの平均値および標準偏差を計算するプログラムを書け。 というものです。 お願いします！！ 280 名前：デフォルトの名無しさん mailto:sage [2006/04/07(金) 20:31:05 ] >>279 平均値と標準偏差を求める式は知っているな？ それをそのまま書けばよい。 avg = sum(a(0:Num-1)) / Num v = sum((avg - a(0:Num-1))**2) / Num s = sqrt(v) 281 名前：デフォルトの名無しさん [2006/04/08(土) 00:17:30 ] >>280 ありがとうございます！！ 書いてみます！

---

---

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

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