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

1 ：デフォルトの名無しさん：2009/01/24(土) 18:32:01

このスレッドは、他のスレッドでは書き込めない超低レベル、
もしくは質問者自身何が何だか分からない質問を勇気を持って書き込むスレッドです。
FORTRAN使いが優しくコメントを返しますが、
お礼はFORTRANの布教と初心者の救済と次期Fortran2008規格でのCOMEFROM文採用をお願いします。

●注意事項
・質問する前にGoogle等の検索サイトで検索しましょう。
・回答者にわかりやすい様に、質問内容はできる限り詳しく書きましょう。
・エラーの場合は起きた状況、環境（OS・コンパイラ・バージョン）、エラーメッセージも詳しく書きましょう。


●前スレ
くだすれFORTRAN（超初心者用）その3
http://pc11.2ch.net/test/read.cgi/tech/1196384126/

くだすれFORTRAN（超初心者用）その2
http://pc11.2ch.net/test/read.cgi/tech/1164121236/

くだすれFORTRAN（超初心者用）
http://pc8.2ch.net/test/read.cgi/tech/1138063703/


●関連スレ
FORTRAN �W
http://pc11.2ch.net/test/read.cgi/tech/1163319215/

672 ：デフォルトの名無しさん：2009/10/28(水) 02:02:02

>>671
EOR　end of record のチェックをするという手もある。自由フォーマットだとうまく行かないが。

あとは、一度一行を文字列として読んで、次にその文字列を内部ファイルとして、それから読み出す
というやり方も考えられる。

データの幅が一定フォーマットなら比較的楽。
空白で切ってるだけで可変長の場合は、ややめんどいかな？

673 ：672：2009/10/28(水) 02:11:12

PROGRAM test
IMPLICIT NONE
! INTEGER, PARAMETER :: iostat_eor = -2, iostat_end = -1
REAL :: a(20)
INTEGER :: io, i, k
CHARACTER (LEN = 136) :: text
! prepare data
DO i = 1, 10
WRITE(9, '( 10(F10.5, X) )') (REAL(k), k = 1, i)
END DO
CLOSE(9)
! (2)
!
DO
READ(9, '(a)', IOSTAT = io) text
IF (is_iostat_end(io)) EXIT ! f2003
a = a / 0.0 ! NaN
READ(text, *, IOSTAT = io) a
IF (is_iostat_eor(io)) EXIT ! f2003
PRINT *, a
END DO
STOP
! (1)
!
DO
READ(9, '(10F10.5)', IOSTAT = io) a
IF (is_iostat_eor(io)) EXIT ! f2003
IF (is_iostat_end(io)) EXIT ! f2003
PRINT *, a(1:10)
END DO
STOP
END PROGRAM test

674 ：デフォルトの名無しさん：2009/10/28(水) 15:56:59

質問です。
シンプソン法を用いてy=exp(x)の計算プログラムを作っています。
以下のプログラムの１５、１６行目のfun()の中には、
何を入れればいいのでしょうか？

675 ：デフォルトの名無しさん：2009/10/28(水) 15:58:29

program simpson1
implicit none
integer::m,i
real(8)::a,b,s,s1,s2,h,x,fun
external fun
1 continue
write(*,*)'INPUT DATA'
write(*,*)'A,B,M ?'
read(*,*)a,b,m
if(m<=0)stop
h=(b-a)/m
s1=0.D0
s2=0.D0
do i=1,m/2-1
s1=s1+fun(i)
s2=s2+fun(i)
end do
s1=s1+fun((m-1)*h+a)
s=h/3*(fun(a)+fun(b)+4*s1+2*s2)
write(6,100) s
100 format(2x,'IR=',D16.6)
end program simpson1
!
function fun(x)
implicit none
real(8)::fun,x
fun=exp(x)
return
end function fun

676 ：671：2009/10/28(水) 16:19:35

>>672,673
ありがとうございます。

「一度一行を文字列として読んでから、その文字列を内部ファイルで読み出す」は応用が利いて便利ですね。
区切りを空白とする以外はフォーマットを限定しない仕様にしたかったので、
この方法がベストだと思いました。NaNのチェックはisnan()を使いました。

最後にもう一点お聞きしたいのですが、EORはformattedで読む場合も関係あるのでしょうか？
また、どういう使い方をするのでしょうか？

677 ：デフォルトの名無しさん：2009/10/28(水) 17:00:51

サブルーチンについての質問したいことがあるのでソース貼らせていただきます。
お忙しいとは思いますが、よろしくお願いします。まず下のソースを見てください。
このプログラム自体は完成しているのですが、
DO 1040 a=1,M
DO 1030 k=1,24
AVE(a)=AVE(a)+(sl(k,a))
1030 continue
WRITE(2,*)AVE(a)/24
1040 continue
の部分をサブルーチンで直して来いといわれたので、
サブルーチンについて調べて色々やりくりしてみたのですが
サブルーチンについての理解が乏しいため上手くいきません。
もしこの箇所をサブルーチンを利用して書き換えるならどのようになるのか教えていただけませんか？
よろしくお願いします。

678 ：デフォルトの名無しさん：2009/10/28(水) 17:02:47

上の続き（ソース）です。よろしくお願いします。
　　　implicit none
integer MM,KK,M,K,a
parameter (MM=400,KK=100)
integer year(MM),month(MM),day(MM),sl(KK,MM),AVE(MM)
character cdummy,CFNAME*110,CFNAME2*80
M=0
CFNAME='/home/maekawa/numeric/kure.txt'
open(1,file=CFNAME,status='old')
1010 M=M+1
read(1,10,END=1020)cdummy,year(M),cdummy,month(M),
@ cdummy,day(M),(cdummy,sl(k,M),k=1,24)
10 format(A5,I2,A1,I2,A1,I2,24(a1,I3))
write(6,*)year(M),month(M),day(M),(sl(k,M),k=1,24)
goto 1010
1020 close(1)
M=M-1
CFNAME2='/home/maekawa/numeric/aki.txt'
open(2,file=CFNAME2,status='old')
DO 1040 a=1,M
DO 1030 k=1,24
AVE(a)=AVE(a)+(sl(k,a))
1030 continue
WRITE(2,*)AVE(a)/24
1040 continue
CLOSE(2)
end

679 ：デフォルトの名無しさん：2009/10/29(木) 02:45:00

>>676
EORはF90以降だったようだ。
昨日は勘違いして７７時代からあった気になっていた。使い方を良くわかっていなかったｗｗ
文法書を見てもらった方が早いが、そもそもFORMATつきでADVANCE='NO'でないと駄目のもよう。

77だとIOSTATでエラー番号を見ればEORを検出できるが、I/Oエラー番号はベンダー依存。
他のI/Oエラーが出ない自信があれば、ERR=nnで飛ぶことも可能。

77だとブランクを数値ゼロに解釈するので、文字列を内部ファイルとして読む場合に
色々前処理をしないとうまくいかないかもしれない。
90なら文字列をTRIMするなり、フォーマットをBNにするなりで対処しやすい。

まとまり無くてスマソ

680 ：デフォルトの名無しさん：2009/10/29(木) 12:12:40

そもそもFortranは「任意のサイズの○○」というのがやりにくい言語だよね
そういうもんだと思って諦めてるけど

681 ：デフォルトの名無しさん：2009/10/29(木) 14:11:01

最初に１回から読みして要素数(行数)確認してから
Allocate...とか、かな。
でもこの辺の取り回しはどの言語も同じだよね。
使っていくうちに増える場合には
一度既存のデータ配列と同じサイズのTempな変数をAllocateしてそこに移して、
元のを開放してからまた同名で要素一個大きいのをAllocateして
TEMPのを移して、TEMP開放、の流れは基本同じだし。

682 ：デフォルトの名無しさん：2009/10/29(木) 14:13:51

終始一貫変化しない任意のサイズ、でなかったら７７だときついな。
データを一度ファイルに落として、そのファイルの行数を数えてから
自動配列で配列生成してそのファイルを読む・・・・・とか。
不可能じゃないけどイヤだなｗ。

683 ：デフォルトの名無しさん：2009/10/29(木) 18:38:33

>>681
Fortranに限らず、データの個数が任意に変わる場合は
ポインタ使ったリスト構造を作るのが普通じゃない？

>>682
結局は、「Fortranでそういうことスンナ」って話になるよねｗ
速度的にも宜しくないし、データ量の多い科学計算には相容れない

684 ：デフォルトの名無しさん：2009/10/30(金) 04:46:30

・・・馬鹿でかい配列を最初に定義してそこで処理するのもなあ。
メモリがもったいない事は今はあまりないけどやはり非効率だし
場合によっては十分大きくとったつもりでも足りなかったり。

685 ：デフォルトの名無しさん：2009/10/30(金) 09:47:58

IBMのFortranのオンラインマニュアルにも、Fortran77と違って
大きな配列を最初に定義しなくて良くなりますよ〜、という事を
Allocateのメリットの例にしてるくらいだからな。

686 ：デフォルトの名無しさん：2009/10/30(金) 13:27:32

FORTRANは言語仕様が堅くて大変だねえ
C/C++ならそういう問題で頭を悩ます必要もないのに

687 ：デフォルトの名無しさん：2009/10/30(金) 14:46:45

>>686
最適化がしやすいようにわざと堅く作ってあるんだよ
C/C++は確かに自由な文法だが最適化がしにくい

688 ：デフォルトの名無しさん：2009/10/30(金) 15:49:03

最適化しにくいし、コンパイラ毎に推奨される書き方の違いは
Fortranより大きいと思う。個人的にはこっちがいやで
Fortranを中心に使っている。

でもやはりというかCのI/Oその他の自由度の高さゆえに
最終的には混ぜこぜコードの出来上がり、になる事も多いけど。
Cはコーディングの自由がありすぎて、結局サイトや教科書の
部分的丸写しが多くなっちゃうから上達しないだけなんだろうけどな。

689 ：671：2009/10/30(金) 18:45:48

>>679
レス遅くなりましたが、ありがとうございます。調べたところ、要するに
「form='formatted' かつ advance='no'　の場合に限り、 eor=文番号 が使える」
ということなんですね。
↓のコードでテストしてみました。

＜ソース＞
real :: f(5)
write(10,*) 1.0 , 2.0 , 3.0
write(10,*) 4.0 , 5.0 , 6.0
write(10,*) 7.0 , 8.0 , 9.0
rewind(10)

do i=1,2
f=0./0
read(10,'(3f)',advance='no',iostat=io,eor=999) f
print *,'io=',io
print *,f
enddo
stop

999 print *,'EOR detected!'
print *,'io=',io
print *,f
end

＜出力＞
io= 0
1.000000 2.000000 3.000000 4.000000 5.000000
EOR detected!
io= -2
6.000000 0.0000000E+00 NaN NaN NaN

690 ：671：2009/10/30(金) 19:00:00

ふと思ったのですが、「レコード」というのはunformattedの場合の概念であって、
単なるASCIIであるformattedでは関係ないはずですよね。
そうすると、EOR（End Of Record）ではなく、EOL（End Of Line）などと
呼ぶほうが適切なのでは、と思ってしまいました・・・。

いずれにしろ、私的には文番号はデバッグ時を除いて使いたくないので、
>>676の方法がベストですね。

691 ：デフォルトの名無しさん：2009/10/31(土) 20:14:11

>>690
IOSTATの番号で見ればいいので文番号は必須ではない。
しかし、いったん１行を文字列に読んで、加工するなりして、
内部ファイルで読むのがいろいろ出来て便利。


EORの呼び名に関しては、御説ごもっともだが、多分ADVANCE＝’NO'によって、
逐次型本来の行単位のアクセスではなく、その内部の要素ごとで、ファイルの現在位置を
止めうるようにしたので、ダイレクトアクセス的なネーミングにしたのではないかと想像？

692 ：671：2009/11/01(日) 21:01:41

>>691
ほんとですね。ありがとうございます。

iostat= と eor= の少なくともどちらか１つを指定していないと、
EORになった時点で↓のようにエラー終了するようですね。
forrtl: severe (268): end of record during read, unit 10, file ・・・

693 ：デフォルトの名無しさん：2009/11/04(水) 05:59:17

IOの挙動がらみで複数の行き先指定すると
その後にgotoスパゲッティになることがあるからムズムズする。

694 ：デフォルトの名無しさん：2009/11/04(水) 16:22:39

メジャーな環境で動く汎用的なコードを書こうという場合、
安心して使えるFortranのバージョンってどれでしょうか？
少し前は「90にしとけ」って話をよく聞きましたが、最近は
95や2003も問題なく通用すると思って良いんでしょうか？

695 ：デフォルトの名無しさん：2009/11/05(木) 01:39:36

>>694
F95が標準の地位を占めたという感じでないか？
更新していないシステムでは、まだF90の環境が結構残っている。
F77はコマンド名として存在していてもF90のコンパイラを固定カラムなどのオプションで
呼び出していることが多い。

F2003はようやくCRAYが完全対応したところで、ベンダーごとに独自の部分的な対応に
とどまっているのが実情。簡単なところから対応しているから大体は重なっているが。

696 ：デフォルトの名無しさん：2009/11/05(木) 05:21:56

もともとがF77で始めたせいで、いまでは固定書式＆左側６文字空白＆72文字制限つきで
でも命令文はF77~95の混在したコードになってもうた。プリプロセスも混ぜてるから拡張子は.fppだ。
まあでも、gfortran と intel fortran はちゃんと処理してくれるから、いいや。
FujitsuのもOKみたいだし。

無駄に行数が増えてしまうけど左側が空いてるとなんか安心するw

697 ：デフォルトの名無しさん：2009/11/05(木) 14:12:15

>>695
ありがとうございます。
F95で書こうと思います。

698 ：デフォルトの名無しさん：2009/11/06(金) 02:03:35

今から始めるんならF95でいいと思うよ。Fortranは古いコードと混ぜてもトラブルは
少ないように規格されてるから過去の資源も使い易いし。

2003でスコープ範囲指定が出来るようになるらしいけど、いつ出るんだろう？
まだ見たことないや。

699 ：デフォルトの名無しさん：2009/11/06(金) 15:29:10

たいていの商用コードはIntelとAMDのx86環境を前提としてるよね。
昔はIBM、日立、富士通、Sun、SGIって色々あったけど、今はそれらは
ほとんど淘汰されてスパコンレベルしか残ってないし。
ま、コードを開発する側としてはやりやすくなったんだろうけど。

700 ：デフォルトの名無しさん：2009/11/07(土) 00:55:13

>>698
>2003でスコープ範囲指定が出来るようになるらしいけど、いつ出るんだろう？
>まだ見たことないや。

変数のスコープか？
PRIVATE、PUBLICは今あるままで、PROTECTED（いわゆるリードオンリー）が加わる程度だと思ったぞ。
他にもあんのかな？

701 ：デフォルトの名無しさん：2009/11/10(火) 09:26:51

COMMON文ってmoduleみたいにスコープを指定ってできないの？

702 ：デフォルトの名無しさん：2009/11/12(木) 13:05:37

名前を共有するんじゃなくて、メモリブロックを共有するための仕組みだからねえ。

703 ：デフォルトの名無しさん：2009/11/12(木) 23:57:09

>>701
おまえふざけてんのか？ｗ

704 ：デフォルトの名無しさん：2009/11/13(金) 03:32:54

お寿司ってハンバーグみたいにソースで食べるってできないの？

705 ：デフォルトの名無しさん：2009/11/13(金) 08:15:17

ああ、普通そうやって食べる。

706 ：デフォルトの名無しさん：2009/11/13(金) 13:00:36

www
両方ともマヨネーズならアリだよな

707 ：デフォルトの名無しさん：2009/11/13(金) 13:04:07

>>706
え？

708 ：デフォルトの名無しさん：2009/11/14(土) 09:03:12

・・・だめ？

709 ：デフォルトの名無しさん：2009/11/14(土) 16:44:02

本物のプログラマはマヨネーズは使わない。
Quiche Eaterはマヨネーズを使う。

710 ：デフォルトの名無しさん：2009/11/18(水) 19:18:00

質問です。
割り付け配列って何ですか？

711 ：710：2009/11/18(水) 19:33:35

710です。たとえば以下のプログラムの変数A,B,Cを割り付け配列すれば
どこを変更すればいいのですか？

program matrix_wa
implicit none
integer::N,I,J
integer,parameter::IDIM=100
real(8),DIMENSION(IDIM,IDIM)::A,B,C
!
write(*,*) 'INPUT MATRIX DIMENSION'
read(*,*) N
write(*,*) 'INPUT MATRIX ELEMENTS OF A'
do I=1,N
read(*,*) (A(I,J),J=1,N)
enddo
write(*,*)'INPUT MATRIX ELEMENTS OF B'
do I=1,N
read(*,*) (B(I,J),J=1,N)
end do
CALL MATADD(IDIM,N,A,B,C)
CALL MATOUT(IDIM,N,C)
stop
end program matrix_wa

712 ：デフォルトの名無しさん：2009/11/18(水) 19:34:17

!==============================================
SUBROUTINE MATOUT(IDIM,N,X)
!=============================================
implicit none
integer IDIM,N,I,J
real(8),DIMENSION(IDIM,IDIM)::X
write(*,*)'C='
do I=1,N
write(*,*) (X(I,J),J=1,N)
end do
return
end subroutine MATOUT
!===============================================
SUBROUTINE MATADD(IDIM,N,A,B,C)
!=============================================
implicit none
integer::IDIM,N,I,J
real(8),DIMENSION(IDIM,IDIM)::A,B,C
DO I=1,N
do J=1,N
C(I,J)=A(I,J)+B(I,J)
end do
end do
write(*,*) 'C='
do I=1,N
write(*,*) (C(I,J),J=1,N)
end do
return
end subroutine MATADD

713 ：デフォルトの名無しさん：2009/11/21(土) 10:18:39

ちょっとお聞きしたいのですが、
３つの物理量の実測（とびとびのデータ）からある答えの物理量
がわかるデータがあるとき、
３つの物理量がある入力値の場合の答えの物理量を出すときに、
データの線形補完をしたいのですが
そういうやり方とかないでしょうか

714 ：デフォルトの名無しさん：2009/11/21(土) 12:07:51

>>713
言いたいことがよく分からないのだけど、補完みたいなことをしたいの？

715 ：デフォルトの名無しさん：2009/11/21(土) 21:21:33

回帰分析とかいう分野の事?

716 ：デフォルトの名無しさん：2009/11/21(土) 21:25:40

>>713
最小二乗法とか？

717 ：デフォルトの名無しさん：2009/11/22(日) 00:11:48

>>714-716
ありがとうございます。
あまりうまく伝えられなくて申し訳ないです。。

詳しくないのですが多分回帰の範囲と思います
3変数から一つの解を出す方法なのですが、
データから多項式で式が一発で出せるようなのがあればいいのですが、
そういうのはないでしょうか。。

718 ：デフォルトの名無しさん：2009/11/22(日) 02:38:52

それを求めたいのか、求めるプログラムを作りたいのか。
前者ならR言語が使えるらしい。

719 ：デフォルトの名無しさん：2009/11/22(日) 12:13:04

>>717
y=a*x1 + b*x2 + c*x3

のような式において、x1、x2、x3、ｙのリストから係数a、b、cを
導きたいなら、Excelの関数で１発でできる。重回帰分析でググれ。

720 ：デフォルトの名無しさん：2009/11/22(日) 12:48:30

>>713
> データの線形補完をしたいのですが

最小二乗法使っとけ。

線型でいいということは測定できる変数x、y、z とそれらから分かる量wの間に

w = a*x + b*y + c*z

という関係があると仮定できるってことでしょ？で、係数a、b、cは分かんない。
この分かってない係数を測定データから決定することで、測定データの中には無
いx、y、zの組のときにwがどんな値になるかを知りたいのだと邪推してみる。
>>713 の頭の中を覗けるわけではないので間違ったこと言ってるかもしれんけど。

それをどうしてもFortranで書かなければならないのなら仕方が無いが、そうでも
なければ >>718 が言うように R でも何でも使ってればいい。

でも最小二乗法ならwのデータも要るのか？元質問だとwは測定して無いっぽいな。
こういうときはどうすんだ？

元質問がグダグダなので無責任に色々と言い放ってみた。グダグダじゃなくても
責任ある発言なんてできないけど。

721 ：デフォルトの名無しさん：2009/11/23(月) 15:09:24

Windows 用の Fortran 95 コンパイラを買おうと思っているのですが，
NAG Fortran Builder 5.2 と Intel Visual Fortran 11.1 日本語版の
どちらにしようか迷っています。すみませんが，以下の疑問それぞれに
ついて判るかた，教えてください。

【1】 以下の項目について，どちらが優れているのでしょうか？
　　(A) コンパイルしてできた実行ファイルの処理速度の速さ
　　(B) 開発環境の使い勝手の良さ（入力補完など）

【2】 Intel Visual Fortran でも NAG Fortran Builder のように
　　GUI アプリケーションが作れるのでしょうか？

【3】 Intel Visual Fortran と NAG Fortran Builder の２つ以外に
　　日本語版の Fortran 95 コンパイラはあるのでしょうか？

なお，私のパソコンは，
Windows が XP Home Edition SP3，
CPU が Intel Pentium 4 (2.40GHz)，
RAM が 512MB です。