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

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

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

●前スレ
くだすれFORTRAN（超初心者用）その4
http://pc12.2ch.net/test/read.cgi/tech/1232789521/
くだすれ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/

>>1
乙！ (≧∀≦)b

お前ハイテンションだな
この手の顔文字って結構イタい

>>1
おつ

質問です。
integer,parameter::dp=selected_real_kind(20)
a=1._dp
みたいに変数(この例ではdp)で精度を決めてるんだけど、指数部の表現ってどうやるか教えてください。
例えば、0.001を表現する時に倍精度なら1.d-3ってな感じでやってますが、変数を使った場合どうなるでしょう？
今は、1._dp/1e3 という風に表していますが、ちょっと美しくないです。
ご存知の方いらっしゃいましたら、教えていただけませんでしょうか？

>>4
1.0e-3_kd

>>3
やっぱ今はやりの顔文字はこれだよな (^p^)b

●関連スレ【 FORTRAN �W 】の URL は変更されていますよ〜。
http://pc12.2ch.net/test/read.cgi/tech/1163319215/

>>7
専ブラなんで気づかなかった、すまん。

>>8
いえいえ，俺もお気に入りメニューに入れてたんだけど，
なんで新しいレスが来ないのかなと長時間思ってた。
んで，昨日初めて気が付いたｗ

はじめまして、超初心者なので質問させて下さい。

今扱っているデータは1秒ごとに取得した時系列データなのですが
データ取得時に装置の調子が悪くて時々データ抜けが起きてしまっています。
１分ごとに平均値を取りたいのですがデータ抜けのせいで
nを１から１０個目までとして平均値を取ってしまうとずれが出てきてしまいます。
ずれを出さないために以下のようにデータを補完したいのですがやりかたが全く判りません。

00:00:04 x4 y4 z4
00:00:05 x5 y5 z5
00:00:06 x6 y6 z6
00:00:10 x7 y7 z7
00:00:11 x8 y8 z8
00:00:12 x9 y9 z9
↓
00:00:04 x4 y4 z4
00:00:05 x5 y5 z5
00:00:06 x6 y6 z6
00:00:07 x7' y7' z7'
00:00:08 x8' y8' z8'
00:00:09 x9' y9' z9'
00:00:10 x7 y7 z7
00:00:11 x8 y8 z8
00:00:12 x9 y9 z9

データの抜けはランダムに起こっています。
フォートランは先生から貰った既存のプログラムを
データにあわせて改造するくらいしかした事がありません、
色々調べたり試したりしてみていますが
もしやりかたが判る方がいらっしゃいましたら教えて頂けると助かりますm(_ _)m

>>10
FORTRAN以前の話だけど、補間はそれでいいの？
00:00:04 x4 y4 z4
00:00:05 x5 y5 z5
00:00:06 x6 y6 z6
00:00:07 x6 y6 z6
00:00:08 x6 y6 z6
00:00:09 x6 y6 z6
00:00:10 x7 y7 z7
00:00:11 x8 y8 z8
00:00:12 x9 y9 z9
の方が自然だと思うけど。

>>11
10です、前の値で埋める方法も判ればそれはとても嬉しいです、
それもやってみようとしてうまくいきませんでした...
また、前後の値のみだとそれほど信用高いというわけでも無いので
フーリエか最小二乗法で補完出来ればなお良い、という感じです。

３秒くらいの抜けは良い方でデータ取得時の状態が酷いときには
20秒くらい抜けてしまっていたりして、どうしたものかと思い悩んでいます。

訂正12>それもやってみようとしてうまくいきませんでした...

やってみたけれど、結果で20秒とか30秒とかの間
同じ値がずらりと続いていておかしな事になった時間帯があり
あまりにも無惨だったのでやめた事を思い出しました。

>>10
補間なのか補完なのかもよくわからん。
（無いデータを補完しておくのか、抜けたのを内挿するのか）
のちのちのデータ解析にも関わるので、一度指導教官とよく話しあった方がいいぞよ。

詳しいことが分からんからなんとも言えないが、データ読み込み時に、データ補間／補完もできるが、
生データは無加工のままとっておいて、後処理で補完した方がいい気もする。

まぁコメントがつけられるデータ形式なら、補完行にコメントを付けておくこともできるが・・・

自宅PCが起動しなくなったトラブル＆携帯が規制かかってて今まで書き込めませんでした、
すみませんm(_ _)m

>>14
無いデータを補完です、データ抜けでその時間帯のデータが無くなっているので。
生データはちゃんととっておいています、フォートランでデータ加工時に
別名をつけていっているので途中経過も判るようになっています。

先生からは「とりあえず調べてやってみて」と言われまして
色々調べてはいるもののちょっと混乱しているのでこちらに助けを求めに来ました。
もしヒントだけでも教えて頂けるようでしたらお願い致しますm(_ _)m

>>15
PROGRAM test
IMPLICIT NONE
INTEGER, PARAMETER :: nmax = 100
REAL :: x(nmax, 3), x1, x2, x3
LOGICAL :: qdata(nmax)
INTEGER :: i, k, k0, kh, km, ks, io, ndata

x = 0.0
qdata = .FALSE.
DO i = 1, nmax
READ(10, '(i2, 1x, i2, 1x, i2, 3F5.0)', IOSTAT = io) kh, km, ks, x1, x2, x3
IF (io == -1) EXIT
IF (i == 1) k0 = 3600 * kh + 60 * km + ks - 1
k = 3600 * kh + 60 * km + ks - k0
x(k, 1) = x1
x(k, 2) = x2
x(k, 3) = x3
qdata(k) = .TRUE.
END DO
ndata = k
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
DO i = 1, ndata
IF ( .NOT. qdata(i) ) x(i, :) = x(i - 1, :)
END DO
DO i = 1, ndata
PRINT '(i3, ":", L5, 3(1x, F8.2))', i, qdata(i), x(i, :)
END DO

STOP
END PROGRAM test

データファイル
00:00:04 1.0 2.0 3.0
00:00:05 4.0 5.0 6.0
00:00:06 7.0 8.0 9.0
00:00:10 10.0 11.0 12.0
00:00:11 13.0 14.0 15.0
00:00:12 16.0 17.0 18.0

実行結果
1: T 1.00 2.00 3.00
2: T 4.00 5.00 6.00
3: T 7.00 8.00 9.00
4: F 7.00 8.00 9.00
5: F 7.00 8.00 9.00
6: F 7.00 8.00 9.00
7: T 10.00 11.00 12.00
8: T 13.00 14.00 15.00
9: T 16.00 17.00 18.00

とりあえず、１秒単位であることを仮定。
読み込むときデータが抜けたら、論理型配列 qdata を　.FALSE.　にしておいた。
抜けたデータは１個前をコピー。故に連続して抜けた場合は同じのが続く。

こんなんでおｋ？

ああごめん>>17は>>16が書いている。名前欄の>>15は間違い。

>>17のデータは2chだと半角スペースがつめられてしまうのであれだが、
実際は空白が入って小数点が揃うように並べてある。
そうでないと読み込みのところのformatがくるってしまう。

データ読み込みとデータ補間は、それぞれサブルーチンにすればすっきりするだろう。
コメント行で切れ目を入れておいた。

スレ汚しすまん。

昨晩また規制に引っかかってしまって書き込み出来ませんでしたm(_ _)m

>>16,17,18
ありがとうございます、
抜けたところに前のデータを挿入して補完する方法は以前試したのですが
データ落ちが酷いところでは数十秒くらい同じ値が続く事になってしまって
あまり宜しくないという事で、今回別の方法で補完するやり方を探しているところです。

現在は3次スプライン補完という方法があるらしく色々調べているのですが
結構難しくてまだ理解出来ていません(>_<)もう少し頑張ってみます。

データ補間してフーリェかけたり平均とったり・・・とか自爆なのを
今でもときどき見るけど、気を付けてね。
走査平均対象の10データ中の有効なものの数でエラーバー決めて解析すればそれでいいと思うけどねぇ。

平均をとる操作と線形あるいは高次補間って同じ事
((非負の係数で)合計が１になる係数をかけて足し算する事に帰着する)だから
とある補間(補完？)した後に平均値をとっても・・・という気はする。
一度、教官に確認した方がいいよ。単なるテストというか簡単な試練を与えられただけなら
それはそれで教育的ともいえるけど・・・危険な香りがする。

>>17 抜けたデータは１個前をコピー。故に連続して抜けた場合は同じのが続く。
王道だね。できればかけたデータの後半はその後にくる有効なデータで、なら
たぶんもっといいけど面倒だしね。
if (.NOT. qdata(i)) then のあとに、

  j = i
  do while ((.NOT. qdata(j) .OR. (j .LE. ndata))
    j = j + 1
  enddo
  ntrue = j

で次の有効なデータの番地がわかるから・・・あれ、これだとntrueは一個ずれるかも。

しかも .OR. じゃなかったね・・・

スレチだけどexcelのグラフでスムージングをオンにするのはマジでやめて欲しい
妙に滑らかなラインが気持ち悪くて仕方ない

最小二乗法かと思ったら移動平均らしいな＞Excelのスムージング
ちゃんと理解して使ってる人は少ないと思うがｗ

>>20
ご指摘ありがとうございます、うまく説明出来なくて申し訳ないのですが、
全体の解析の結果から数秒〜数分単位でグラフにして見たい部分とかが判った時に
予め近似曲線が求まっていれば楽でしょう、という感じで現在に至ります。

なんか、何から何まで良くわからない状態で本当に申し訳ありません・・・

先生の蔵書の中から参考資料を見つけて頂けました、これを使ってまた少し頑張ってみます。
自分でも至らない点や確認出来た点が多かったのでとても助かりました、色々と助言をありがとうございましたm(_ _)m
また判らない事があったらこちらに駆け込むかもしれませんがその時は宜しくお願いします。

>>26
新Ｍ1か、四年生かな頑張れー。

超初心者なのですが、Intel Visual Fortranを使っていらっしゃる方はいますか？
私にはC言語を白黒画面のコンソール画面でコンパイル＆リンクして実行する程度の簡単な経験しかありません。
殆ど統合開発環境を使った事がないのですが、関連ホームページを調べると、魅力的なCPU毎の最適化機能とか、IMSL数値計算ライブラリを使えるようなので興味を持っています。
私が使い始めることになったバージョンは10.1です。

すいません。つまらない問題なのですが、
複数ファイルに同一の結果を出力する場合どうすればよいのでしょうか。

open(10,file='A.dat',form='formatted')
open(11,file='B.dat',form='formatted')

write(10,101)int((t2-t1)/60),dble((t2-t1)-int((t2-t1)/60)*60)
write(11,101)int((t2-t1)/60),dble((t2-t1)-int((t2-t1)/60)*60)

みたいな表現をコンパクトにしたいのです。

>>29
DO iunit = 10, 11
WRITE(iunt, 101) ....
END DO

ありがとうございます。
write(10,101)・・・
write(31,101)・・・
write(*,101)・・・

みたいな場合はその方法だと無理ですよね
やはり一行ごとに処理するのが基本なのでしょうか

>>31
そう多くないなら１行づつが基本だと思う。

でも一応、
INTEGER.PARAMETER :: iunits(3) = (/10, 31, 6/)

DO i = 1, SIZE(iunits)
WRITE(iunits(i), 101) ....
END DO

とか、あるいは内部サブルーチンにして
CALL wr_info( (/31, 10, 6/) )
......
CONTAINS
SUBROUTINE wr_info(iunits)
INTEGER, INTENT(IN) :: iunits(:)
DO i = 1, SIZE(iunits)
WRITE(iunits(i), 101) ....
END DO
RETURN
END SUBROUTINE wr_info

内部サブルーチンなら、変数は共通。 (/ /) は配列生成子で、これを使えば引数でファイル番号を
並べて指定すればよく、わざわざ配列をメインルーチンで用意する必要がなくなる。

６番は歴史的理由でコンソールになっている（はず）。

６番＝標準、７番エラー出力・・・だっけ忘れた。
なんかStdoutみたいな明示するキーワードあった気がするけど忘れてもうた。

>>24 最小二乗法
でもxもyも観測量なのに普通のxは誤差無し確定値とみなすやつをへーきで使ってたりするけどねぇ。

>>24
エクセルなんか使うなよ・・・人文科学かな？

環境系の同僚に「このあたりのコーディングはだれか詳しい人に助っ人頼んだ方がいいよ」って言ったら
「いや自分でやるのが大事ですから。シミュレーションはエクセルで私がやります」って言われた
本人は、今いいこといったオレ、見たいな風情なのでなにもいえなかった

放っとけ
構うとお前さんに仕事が飛んでくるぞ

>>36 「いや自分でやるのが大事ですから。シミュレーションはエクセルで私がやります」

キリッ！
やる夫みたいだw

program sample444
real :: x(3),y(3),z(3)
integer :: i
open(1, file="aiueo.txt")
do i=1,3
read(1, *) x(i), y(i), z(i)
print *, x(i), y(i), z(i)
end do
close(1)
end program


上記のコマンドでコンパイルしてプログラムを作成しましたのですが
次のようなエラーが出てしまいます！ファイルの入出力だけがうまく進まず
もうここから全然進みません；
どなたかわかるかたよろしくお願いします；
ちなみにtxtファイルは3×3配列の数値がちゃんと入っています；
traceback: not available, compile with -ftrace=frame or -ftrace=full
Fortran runtime error: End of file

>>39
OSとコンパイラの種類
読み込むデータの例も書いた方が良いよ

>>40
わかりました。すいません；
んと、windows vistaで、コンパイラはg95です。
読み込むデータはたとえば
25　50　59
34　23　11
45　43　65
とかです。
どうかよろしくお願いします。

とっさに思いつく限りでは・・・

元データの最後の改行位置は？
open で status="unknown" を追加してみては？

くらいかなぁ。Good Luck!

読み込むデータがInt型なのにrealで読み込もうとしてるのは良くないと思う
確認してないから他にも間違いがあるかもしれないけど

>>39 Fortran runtime error: End of file
だから、>>42 のとおりデータの最後の行の右端に改行コードがないんじゃない？
>>43
大事な事だよね。
整数を実数でよんでそのあと整数化すると１小さくなることがあるし。

>>42
ふむふむ
改行位置で変わったりするのですか？
テキストにも書いていなかったので少し勉強してみます。
ありがとうございました。

>>43
なるほどー
確かにそのとおりですね；
修正しときます！
ありがとうございました。

>>44
改行コードが何か全くしりませんでした；
元データをただ改行するだけじゃだめなんですか？
もう少し勉強してみます；
ありがとうございました！

連投失礼します！
皆さんの助言でなんとかファイル読み込みできるようになりました！
ですが、1列目しかデータが反映されません；
これも改行コードとやらによるものなんでしょうか？
何回もすいません；

あまり関係ないかも・・・
Win機上では改行コードが LF/CR かな２つ無いと読めない・・・
ぎゃくにGCCやG95はそれぞれを改行コードとみなしちゃう、とか
・・・今でもあるのかなぁ・・・・？

qkc だっけ、漢字コード変換のソフトがおまけ機能で
改行コードもWin機向けに変えてくれるから、
もし元データがLinux機で作られてたのなら・・・試してちょ。
unix2dos dos2unix っていうのでLinux上でも改行コード変換はできるけど。

>>33

キーワードの方は知らんけど、標準エラーは 0 だ。

# 理由は知らんが規制が掛かって暫く
# 書き込めなかった。

thanks

>>48 理由は知らんが規制が掛かって

そういえば今回は ac.jp ドメインもいくつか規制食らってたなw

>>50
うち ocn なんだけどね。学校とは縁が切れて、博士号持ちのワープワやってるw

そうなんだ。

大学院進んだ人とかポスドクっておかねが多少もらえても雇用統計では「自発的失業者」らしいね。
それを聞いて、博士号取得をくじけた人をしっているけどそれは正解でもあるとおもった。
・・・博士号もち仲間よりw

>>51
うお、うちの会社でアルバイト（年300万研究補助）しね？

これはひでえなｗ＞自発的失業者

！？
大学院進んだ人は普通に「学生」みたいな区分だと思ってた。
ポスドクも任期付／非常勤だからパートと同じような扱いかと・・・(´；ω；｀)ﾌﾞﾜｯ

sage忘れたorz

泣くんじゃないッ！(´；ω；｀)

fortran90の勉強はじめましたっ！
fortran77は全然知りません :DD

>>58
遺物と向かい合う必要がない限り、知らんで良い。
でも Fortran90 って良書はあるのかな？
色々出てるみたいだけど、倍精度の説明さえロクに
書いていないのは今も謎だ（個人的に）。

これから開催される CUDA Fortran のセミナーや講習会，知りませんか？
ググっても3月に東工大のGPUコンピューティング研究会が開催したのしか
検索されてこなくて。
もっと早くチェックしてればよかったなぁ‥‥

>>52
生活保護対象にならないようにそういう扱いにするというのは聞いたことがある。
にしてもひどい話ではある・・・(´；ω；｀)
>>60
CUDAの試作機とかどっかタダで使えるところないかなぁ・・

nvidiaとか本気で流行らせたいんなら、試作機をリモートでいいから無償提供とかやりゃいいのにね
言語なんていかに裾のが広がってくかが大事なのに

>>62
? ただで（Cだけど）コンパイラ配ってるし、普通のグラボでも使えるし、Fortranコンパイラにも協力しているしで
十分やってると思うけどね。

ぼくのマシンはインテルのへっぽこアクセラレータしかないお・・・でも触ってみたいんだよう、
という人はいると思うんだ。すでに大規模＆本格的に使ってる人(おおくないとおもうけど)が、
かるいデバッグ作業するとかなら普通サイズのぐらぼマシンに自分で設定して、でもいいとおもうけど、
将来使うかもという程度の興味出始めな人には
横着なだけかもしれないけどコンパイルと実行環境が整った場所があるとやはりうれしい。

最近になって、がっつり数値計算をする必要がでてfortranを
勉強しだしたんですが、要素数が不明の２次元データをファイルから読んで、
その要素分だけallocateしたいのですが、
とりあえず最後まで読んで要素の数を数えてから配列をつくろうかな、と思ったんですけど、
fortranにおいてend of fileの扱いをどうしたらいいのかわからずに困っています
使用ＯＳはvista,コンパイラはg95です。初歩的とは思いますが、よろしくお願いします。

open(1,file='data.dat',status='old')
icount = 0
998 continue
read(lun,*,end=999) fjlasdjfkds
icount = icount + 1
goto 998
999 continue

とか、かな。
でもこれ、ファイルの最後の最後が改行コードでおわるか
EOFで終わるか両方無しでいきなり終わるか・・・でいろいろ挙動が変わるかも。

fortran 90を使用しています．
二次元配列に初期値をparameterとして代入するにはどのようにしたらよいのでしょうか．

implicit none
integer,parameter :: i(2,2) = (/1,2,3,4/)

だとエラーを吐きます．コンパイラーはPGIです．
宜しくお願いします．

>66
マジでありがとうございます！！
とりあえずreadの一個目にファイル変数渡して
できました。改行からのＥＯＦも即ＥＯＦも問題なく動きました。
とりあえず、ＥＯＦなしのデータは使わない？と思うので。
あと、個人的な好みなのですが、goto文はあまり使いたくないので、
とりあえずdoループで書いたのですけど、endはやはり、
ちょっとぐぐってみたかんじ、ジャンプさせるしかないようにみえるんですが、
exitで抜けさせることはできないのでしょうか？
もしご存知のかたがいらっしゃいましたら、ご教授おねがいします。
なんにせよ、とりあえず進めそうです！66さんありがとうございます。

>>67
ttp://www.str.ce.akita-u.ac.jp/~gotou/zyouhou/f90/sec6.html#6.4.4.

>>69
ありがとうございます．
地道に代入することにします

>>68

いまコンパイラないから試してないけど、これで走らないかな。

open(1,file='data.dat',status='old')
icont = 0

do
read(1,*,iostat=io) aaa
if(io /= 0 ) exit
icont = icount + 1
end do

>>67
INTEGER, PARAMETER :: n(2, 2) = RESHAPE([1, 2, 3, 4], [2, 2])

[ = (/
] = )]
ただし Fortran2003

>>71
iostat つかっての方がきれいで良いね。

>>73
他のエラーを捉えられなくなるのが難点か。

エラー番号って規格で決まってるの？

レス遅くなって、すみません。
iostatの読み込みで、−１がたったらexitするように書いたら
望むとおりの挙動をしてくれました。
なんとなく勝手がわかってきたきがします。
重ね重ね、ありがとうございます！

>>74
大半は決まってない。
EOF=-1 は、まず大丈夫だが、これが規格なのか、慣例なのか分からん。
たしか EOL=-2 だったか。

Fortran2008の次期規格で、この辺の一部の定数が、INTRINSIC定数として定義されるらしい。

>>77
EOF=-1は慣例で規格としては決まってない。
だから、EOF=-1と決めつけてプログラムするのは危険。
まぁ、そんな危険はほとんどないが。

>>78
EOF=-1 は慣例なのかー。
知らんかったぜ。

すみません、質問させて下さい。

数万行ある（x,y）データを600行ごとに値を読んでサブルーチンで計算させて
そのサブルーチンで出て来た値をまた別のサブルーチンで計算させて最後に出て来た結果を
出力ファイルに書き込むようなメインプログラムの作り方がどうにもわかりません。
判るかたがいましたらご教授頂けると助かります。
f77を覚えたてです、まるっきりの初心者で恐縮ですが宜しくお願いします。

追記。
数万行のデータのうち始めの600行に対してサブルーチン１、サブルーチン２を実行、その結果を出力ファイルに書き込む
次の600行を同じようにサブルーチン１、２実行、結果を同じ出力ファイルの最後の行から付け足して書き込む
これを入力データの行数分だけ行いたい、という感じです。

age

入力データは (f7.0,f9.2) X(I),Y(I) です。
X(I)は0から始まって1,2,3,...600,601,...20000,20001...
という風に１ずつ増えています。

おまえらソースコードの管理どうしてる？

open(10,file=infile)
open(20,file=outfile)
N=600
30 continue
DO 20 I=N,N+6001
DO 20 I=0,N+1
20 READ(10,*) X(I),Y(I)
CALL SUB1(X,Y,G,H,N,MP)
DO 40 I=N,N+600
40 WRITE(20,*) X(I),Y(I),G(I),H(I)
WRITE(6,*) X(I),Y(I),G(I),H(I)
N=N+600
goto 30
close (20)
close (10)
STOP

サブルーチンを一つだけにして結果をファイルに出力させてみました。
上記のようにしてみたのですが、どうしてもうまく600行毎に計算できません。
もし判るかたがいらっしゃいましたら助言を頂けると助かります。

DO 20 I=N,N+6001
DO 20 I=0,N+1
訂正
DO 20 I=N,N+601

>>84
脳内

ぬるぽ

ガッ

>>85
すくなくとも　

N=1

DO 20 I = N, N + 600
もう一個のDO 20 は要らない。

>>90
ぬるぽいんたーえくせぷしょん

>>91
FORTRAN スレなんだから、ABEND だろ！

>>90
すみません、色々試していてコメントアウトしたところとかを載せてしまっていたので
86で訂正しています。ややこしくて申し訳有りません。
今のところサブルーチン２つまわすやりかたのほうで下のような感じにしています。

N=0
30 continue
DO 20 I=N,N+601
20 READ(10,*) X(I),Y(I)
CALL SUB1(X,Y,G,H,N,MP)
DO 40 I=N,N+600
40 WRITE(20,*) X(I),Y(I),G(I),H(I)
HH=1
DO 30 I=N,N+600
XX=X(N)+HH*I
YY=SUB2(G,H,X,Y,XX,N,MP)
KH=0
DO 50 K=N,N+600
IF(ABS(X(K)-XX).LE,1.d0) THEN
GOTO 50
ENDIF
50 CONTINUE
WRITE(20,fmt)XX,YY
30 CONTINUE
N=N+600
goto 30

>>80続き（改行が多すぎるのでエラーになってしまったので）

PARAMETER (MP=36000, NW=600)
DIMENSION X(0:MP),Y(0:MP),G(0:MP),H(0:MP)としています、
MPはデータの行数で、ファイル毎に少し違うのですが最大行数を入れてあります。

600行を読み込んで、SUB1とSUB2を計算させて、ファイルに書き出して次の600行をまた・・・
という感じで操作させたいのですが、何故かおそらく始めの600行のSUB1で出て来るHやGを使って
次の600行以降が計算されているような変な結果になってしまいまして。

始めの600行に関しては結果として正しい値が出て来ているのでサブルーチン部分には問題は無いと思います。
本当にすみませんがよろしくお願いします。

あわわわまた間違えていました!!!30 CONTINUEが2回ありますが
>>93の頭の30 continueは15 continueの間違いで、最後のgoto 30がgoto 15になります!!!
ほんとうに申し訳有りません。。。

どうもGOTO文やDO．．CONTINUE文で混乱しているようだけど，
まだFORTRANを覚え立てなら，わざわざ古い書き方で覚えずに，
Fortran90形式のDO．．END DO構文で慣れた方がいいんではないかい？

特にGOTO文を使うのは間違いを誘発しやすいから，使わないようにするのが世界的な流れだし．

ところでend doとendoはどっちが望ましいの？
個人的にenddoの方が一語にまとまってて好きだけど

>>90
そうなんですか、理解していませんでした(>_<)

f77とf90だとかなり違うのでしょうか？
例えば参考にしているサイト陣がf90だったりf77だったりそれぞれ違っていたと思うのですが、
取り混ぜてしまっていたらまともなコードにはならない・・・という考えで合っていますか？

教科書は今までf77のものを使っていたので、なんにしてもf90のものを探してみます。

(>_<) …

良いこと思いついた
焚書坑儒して77を滅ぼそう

>>94
色々、妙な点が多いプログラムではあるので、可能性は色々あるが、
サブルーチン側でX(１）〜X(６００）までで計算してるんで内科医？

そもそも、古い数値がいらないなら、６００個の配列に毎回書き換えて読み込めばいいし、
サブルーチン側でも６００個を受け取って、１〜６００までの要素で計算すれば良い。

全要素を引渡しているようだが、必要なのかな？

>>100
７７もいいとこあるんだよ。

大学のPCのOSがLinuxでfortranを学び始めました。
家でも練習みたいのをしたのですが、windows vistaでもできるのでしょうか。
例えばメモ帳に書いて、フリーのコンパイラでコンパイルするというのを想定しているのですがこれは無理なんですか。

>>102
知らないから僕ではなく他の人に聞いてみてよ

>>102
http://pc12.2ch.net/test/read.cgi/tech/1163319215/2-3

>>104
そのスレの３つ目は・・・・w

>>102

これが所謂ゆとりというものなのか

>>101
助言ありがとうございます、毎回書き換えてみる方法を色々試してみました。
下のようにメインプログラムを書き換えてみたのですが、
READ(10,*) X(0),Y(0)
do N=1,NW,MP
DO I=1,MP
READ(10,*,err=99) X(I),Y(I)
end do
CALL SUB1(X,Y,G,H,N,MP)
DO I=0,INT(X(MP-1)-X(0))
XX=X(0)+real(I)
K=1
do while ((X(K-1)-XX)*(X(K)-XX) > 0.0d0)
K=K+1
end do
YY=SUB2(G,H,X,Y,XX,K,MP)
WRITE(20,2100) XX,YY
end do
X(0) = X(MP)
Y(0) = Y(MP)
end do
close (20)
close (10)
STOP

まぁまぁw
Win機だと gfortran が一番かな？
タダだし、Linux版その他あってWin機でテストしたコードをよそに持っていっても
規格やら互換の問題がとても少ない・・・とおもう。OpenMPも使えるし。
コマンドプロンプト内の処理しかなく、Linuxコンソールと基本同じ使い勝手。
インスト(といってもファイルを展開するだけなのだが)楽チン。環境変数設定はめんどいw

107続き
最後のほうのデータが抜けてしまって、現在はどうしたものかと悩んでいます。
取り敢えずサブルーチンのDIMENSIONの設定とかをもう少し見直してみます。

>>108
すいません。。PCは大学入るまで使うことが少なくてよくわからなかったのです。
とりあえずいわれたものをダウンロードしてみます。ありがとごいました

>>110
g95だがこのへんとか
ttp://d.hatena.ne.jp/arakik10/20090213/p1

あと最近の本だとこれがわかりやすい。
ttp://www.amazon.co.jp/dp/4627847211

タダのFortranコンパイラを、と訊かれると
自分が使っている＆インストール簡単だった、というそれだけの理由でgfortranを
人にはいつも薦めてるけど・・・じつはg95との違いをよく知らないのよねw

ぐぬーライセンスとか内部のオトナの事情はさておき、１ユーザーとしてみると何が違うのかな？

gfortranでは装置番号を環境変数で指定できない

>>85,>>92 → >>107
同じ人が書いたとは思えない！初心者でも1日でこんなに上達するんだな。

>最後のほうのデータが抜けてしまって
表現があいまいすぎて意味分からん。もしかして女の子？
助言を求めているなら、具体的に書かないとなんとも言いようがない。

×　>>85,>>92 → >>107
○　>>85,>>93 → >>107

あえて７７で書いてる人っているの？

ifortマンセー

>>114
わかりにくくて申し訳有りません(>_<)データは30000行入っていて欲しかったのですが
どうやら最後のほうや途中がデータ計測の時に抜けてしまったようでした。
抜け自体はそれほど多く無いので線形補間というのをすれば問題は無いそうです。
計算式は渡されて、どうやるかはまだ良く判っていませんが、調べてみます。

>もしかして女の子？
はい、最近分野変えで物理系の研究室に転入してきたばかりなのですが
今まではプログラミングには縁もゆかりもありませんでした。

>107は昨晩学校を徘徊して他の研究室に残っていた先輩からもアドバイスを貰いました。
>93については先輩曰く「動くようなプログラムには見えないぞ」との事で
つまるところ>>96と同じような指摘を受けたりなんだりして、手直ししてみました。

自分ではまだまだ勉強が浅いので次に違うプログラムを作る時には
また色々戸惑うだろうなぁとは思っています。
これからまたしっかり勉強します、ありがとうございました。

すいません。学校の課題でy=(x+1)(x-3)=0でNEWTON法を使い、収束判定が0.0000001で終わるプログラムを
「READ」「WRITE」「if」「GO TO」「SUBROUTIN」「DIMENSION]だけ使って作るんですけど・・・
SUBROUTINEとDIMENTIONが全く分かりません。誰か教えてください。

>>118
匿名掲示板だから思い切って聞きますが、処女ですか？

>>120
横レスですが違います

>>120
他人ですが・・・うるさーい(>_<)

＼　　　/　.::::::::::::::::::::::::;;:;;::,ｯ､::::::　　 )　　く 　 ホ 　す
　 ＼　l　 ,ｯｨrj,rｆ'"'"'"　　　 lﾐ:::::::　く　　 れ　 モ　 ま
　　　　 Y　　　　　　　　　　 ,!ﾐ:::::::　ヽ　　な　 以　 な
｀ヽ､　　|　　　　　　　　　　 くミ::::::::　ﾉ 　 い　 外　 い
　　　　 |､__　 ｬー--＿ニゞ　｀i::::,rく　　　か　 は
｀`''ｰ-　ゝ､'ｌ　 ￣ﾞ´彑,ヾ　　　}::;! ,ﾍ.）　　！　 帰
　　　　　 ﾞｿ　　　""""´｀　　　　　〉 L_　 　 　 っ
　　　　　 /　　　　　　　　　　i　　,　　/|　　　　て　　　　r
≡＝-　〈´　,,.._　　　　　　　 i　　't-'ﾞ　|　,へ　　　　　,r┘
,､yx=''"　｀ｰ{ﾞ　＿, -､　　　　;　　l　　　ﾚ'　　ヽr､⌒ヽ'
　　　　　　　 ﾞ､`--─ﾞ　　　　　 /!　　　　　　　　 `､
　　_,,､-　　 　 ﾞ､ ー''　　　　／　;　　　　　　　　　　 `､
-''"_,,､-''"　　　 ﾞ､　　　　／;;'　,'　　/　　　　　　　　 、＼
-''"　　　 ／　　　`ｰ─''ぐ;;;;' ,'　　ノ　　　　　　
　　 ／／　　　 /　　　　　ヾ_､=ﾆﾞ

本当に誰か助けてください(TAT)
流れは
xを読み込む
↓
y1=(x+1)(x-3)とy2=2x-2のy1、y2をそれぞれ算出
↓
H=y1/y2を算出
↓
Hが0.0000001より小さいか判別
↓
小さかったらHを出力して終了、大きかったらx-H=Aを算出
↓
Aをxに代入して0.0000001より小さくなるまで繰り返し計算する

こんな感じのプログラムを作りたいんですが・・・全然分からなくてできません・・・
誰か本当に助けたください

>>124
ｶﾞﾝｶﾞﾚ！

>>90あたりからの流れ好きだよ俺は

>>124
エラーチェック無し、暴走考慮せず、記述通りで無駄あり、Fortran90だ。
PROGRAM unko
IMPLICIT NONE
REAL :: x, y1, y2, a, h
x = 0.0
DO
y1 = (x + 1.0) * (x - 3.0)
y2 = 2.0 * x - 2.0
h = y1 / y2
IF (ABS(h) < 1.0e-7) EXIT
a = x - h
x = a
END DO
WRITE(*, *) 'result=', h, ': x=', x
STOP
END PROGRAM unko

実行結果
result= 0.0000000E+00 : x= -1.000000

>>124
プログラム出しておいてなんだが、たぶん問題が間違ってるか、問題を写し間違っているか、
質問者が問題を理解しそこねている。

この問題だと、単にｙ１のゼロ点近傍に近づいたときに終了し、y2の条件があまり意味を持たない。
あまりにイミフｗ

>>124
ごめんｗ　>>119でNewton法で質問していたのか。Newton法なら意味わかる。
前レスは忘れろｗｗｗ

>>119
＞「READ」「WRITE」「if」「GO TO」「SUBROUTIN」「DIMENSION]だけ使って作る

これは出題者がこうしろと言ってるの？

>>129
ありがとうございます。そして悲しいことにFORTRANは77なんです(TAT)


>>130
そうなんです(汗)
講義で教えられてないのを使うと大目玉をくらうんです。

>>118
よく解らんが、不等間隔のデータを読み込んで、等間隔のデータに補間して出力したいのかな。
＞ do while ((X(K-1)-XX)*(X(K)-XX) > 0.0d0)
この条件式の書き方は昔風で今は流行らない。素直にIF文にした方がいい。

こんな感じの事をしたいのかもしれない。
PROGRAM unko
PARAMETER (MP = 600)
DIMENSION x(0:MP), y(0:MP)
READ(10,*) X(0),Y(0)
DO N = 1, 999999
DO I = 1,MP
READ(10,*,err=99) X(I),Y(I)
END DO
CALL SUB1(X,Y,G,H,MP)
K = 1
DO I = 0, INT( X(MP - 1) - X(0) )
XX = X(0) + REAL(I)
IF ( ( X(K-1) < XX ) .AND. ( XX <= X(K) ) ) THEN
YY = SUB2(G,H,X,Y,XX,K,MP)
WRITE(20,'(2F15.7)') XX,YY
ELSE
K = K + 1
END IF
END DO
X(0) = X(MP)
Y(0) = Y(MP)
END DO
99 CONTINUE
STOP
END

>>131
７７なら７７と最初から言え・・・といいつつ本来の質問を見ていなかった朕。
一応、昔風に芋臭く書いてみた。DIMENSIONをどう使ったらいいのか分からん。
微分もプログラム的にやることにしてみた。ちょっと無理っぽいｗ　　あと、ABSを使った。

PROGRAM main
WRITE(6, *) 'INPUT START VALUE ='
READ(5, *) a
1 x = a
CALL newton(x, h)
a = x - h
IF (ABS(h) .GT. 0.0000001) GOTO 1
WRITE(6, *) 'h=', h, ' x=', x
END

SUBROUTINE newton(x, h)
DIMENSION c(2)
c(1) = -1.0
c(2) = 3.0
y1 = (x - c(1)) * (x - c(2)) ! f = (x + 1)(x - 3)
y2 = (x - c(1)) + (x - c(2)) ! df/dx = 1(x - 3) + (x + 1)1
h = y1 / y2
RETURN
END
実行例１
INPUT START VALUE =
5.0
h= 0.0000000E+00 x= 3.000000
実行例２
INPUT START VALUE =
-3.0
h= 0.0000000E+00 x= -1.000000

>>133
ありがとうございます。(T∀T)
やっぱりDIMENSIONは無理やり使わないと使えないですよね(汗)
ちょっと疑問になったんですが、ABSのところをそのまんまhにしてもエラーが
おきたりしないですか？

>>134
エラーは発生しない。
それとお礼を言われる程のことはしていない。

>>135
嘘教えんなよｗｗｗ誰だよｗ

>>134
hが負になる時の事を考えると、ABSが必要。
どうしてもABSを使いたくないなら、
IF (h < -0.0000001） GOTO 1
IF (h > 0.0000001) GOTO 1
とし給へ。

>>136
＞嘘教えんなよｗｗｗ誰だよｗ

＞エラーは発生しない。
これは嘘ではないと思うが。

＞それとお礼を言われる程のことはしていない。
それともこっちが嘘と言っているの？

プログラマーって短気が多いんだなぁw

>>137
キモイです．．．．．(>_<)

>>136
いやエラーは出ないだろ。

>>132
これは書き方が判りやすくて助かります！おっしゃる通りのようなプログラムです、
どうやら抜け落ちたデータの補間もこれで出来るらしい？です。

107に書いたやりかたで取り敢えず動くのでやってみたのですが
50000行くらいの大きなデータだと、結果が出るまでもの凄く時間が掛かってしまいました。
変なところにバッファしてるんじゃないの？と言われたのですが良く判らず
10000ずつ計算させて切り貼りしています。どうしたものかと悩んでいます。

＞ do while ((X(K-1)-XX)*(X(K)-XX) > 0.0d0)
先輩が教えてくれたのですが、イマイチ意味が判っていませんでした、
普通にIFでやっても良いんですね、ありがとうございます。

>>141
計算機でやる分には、時間がかかってもそのままやらせればいいと思う。
手で切り貼りする作業時間を含めれば、結局トータルでの時間は、計算機でやった方が早いのでは？

時間がかかる原因は、色々ある。
しかしプログラムのどの部分に時間を取られているのか分からないとなんと見えない。
細切れのI/Oが多いと時間はかかる。だが、５００００行くらいなら、６００行づつ読んでもそれほどでもないような。

元のプログラムだと、毎回最初のデータから条件を満たす位置を探し始めているが、
Xが単調に増加しているのなら、前回の位置から探し始めれば早く適切な場所が見つかると思う。
この辺に、時間がかかっている気がしなくもない。

>>142
一桁間違えていました！(>_<)500000行くらいのデータです。
数時間の測定だと100000行以内くらいで済むので５〜６分で計算出来るんですが
200000行くらいだと20〜30分くらい掛かってしまいます。
500000行くらいのデータをやろうとしたら数時間掛かってしまっていました。

Xは単調増加なので前回の位置から探すというやりかたをちょっと考えてみます、
ありがとうございます。

>>143
N=
1 * 10**5, 5. .6min 　　5 / 1^2 = 5
2 * 10**5, 20.. 30min 　　20 / 2^2 = 5
5 * 10**5, 60..180min 〜125/ 5^2 = 5

計算時間が要素数の自乗で増えているが、処理内容から言ってそうなるはずのものではない。
計算回数が自乗になりそうな所を探せば解決する気がする。

データを毎回頭から探すような部分はｎC_2〜O（n^2)だから、可能性は高い。

あと、それとは別に、今デバッグモードで実行しているなら、リリースモードで実行することにより
２〜１０倍くらいはやくなる可能性はある。その辺はコンパイラに依るので先輩などに聞いてみるとよいかな。

>Xは単調増加なので前回の位置から探すというやりかたをちょっと考えてみます、
>>142が言っている「元のプログラム」って>>85や>>93のことではないの？

>>107も>>132も，すでに前回の位置からREADするようになっているでしょ．

>>143ちょっとまて！
>>107でやってみたって言ってるけど，SUB1の引数にNが入るのっておかしくないか？
>>107バージョン CALL SUB1(X,Y,G,H,N,MP)
>>132バージョン CALL SUB1(X,Y,G,H,MP)

データの補間がしたいと言うことから推察するに，スプライン補間かな？（SUB1は行列を解いてる？）
5番目の引数が配列のサイズを渡しているんだとすれば，>>107ではNに相当大きな配列サイズが指定されて，
無意味なメモリ領域を読み取ってそのまま解いているんじゃないか？
SUB2ではその無意味な計算結果は使われないから，それで正常な結果が得られたとしても納得できそう．
ちょっと行き過ぎた推測かも知れないけど．．．

いちど，サブルーチンの引数を確認してみては？

>>146
毎回最初からと言っているのは、READではなくて、補間位置を決める、２つのデータに挟まれている位置を
求める部分の話。

>>107での２重DO−LOOP構造では、DO内でK=1にしているので毎回X(0)から探し始めている
DO I=0,INT(X(MP-1)-X(0))
XX=X(0)+real(I)
K=1
do while ((X(K-1)-XX)*(X(K)-XX) > 0.0d0)
K=K+1
end do
..YY=SUB2(G,H,X,Y,XX,K,MP)
end do
..
end do

一応>>132では（Xの単調増加を仮定して）外に出してある。

>>146
読み込んだデータを延々と配列にため込んでいるので、毎回現在のポインタ位置のような意味でNを渡す
必要があるのではないかと思う。

必要があるといいですね

前半は納得した．
後半「読み込んだデータを延々と配列にため込んでいるので、」というのは違うと思う．
どちらも0:MPの配列を使い回してる．
>>107の抜粋
　READ(10,*) X(0),Y(0)
　do N=1,NW,MP
　 　 　DO I=1,MP
　 　 　 　 　READ(10,*,err=99) X(I),Y(I)
　 　 　end do
　 　 　CALL SUB1(X,Y,G,H,N,MP)
　 　 　．．．（略）．．．
　 　 　X(0) = X(MP)
　 　 　Y(0) = Y(MP)
　end do

>>132の抜粋
　READ(10,*) X(0),Y(0)
　DO N = 1, 999999
　 　 　DO I = 1,MP
　 　 　 　 　READ(10,*,err=99) X(I),Y(I)
　 　 　END DO
　 　 　CALL SUB1(X,Y,G,H,MP)
　 　 　．．．（略）．．．
　 　 　X(0) = X(MP)
　 　 　Y(0) = Y(MP)
　END DO

どちらにしてもサブルーチンの中身が分からないと，X,Yの配列を使い回すのが正しいのかすら推測の域を出ないな．
質問者のレベルから考えると，引数の意味を把握しているか自体が心配．
（このスレ，勉強になるなぁ）

>>149
おｋ把握した。
古い方のプログラムと混乱していた。

古い方は、データをため込んでいるので、現在位置をサブルーチンに渡す必要があるから、
Nを渡していると理解していた。

とすると>>107では、Nの値が狂っているなｗ　N=1で固定しなきゃないかな？ｗ

>>144->>150
今理解するのにちょっと時間が掛かっています(>_<)皆様有り難うございますm(_ _)m

>>146さんの言う通りサブルーチンの中身はスプライン補間になっていました。
SUB1で行列を解いていて、SUB2で間隔1毎の値を計算している？感じです。
引数の意味はまだ良く理解出来ていません、SUB1のNがおかしいのでしょうか？
もう少しいぢくってみます。。。

>>151
メール欄に sage と書き込んで、スレが上のほうに上がらないようにした方がいいですよ。
スレが上に上がると、通りすがりの冷かしが面白半分なレスを書き込んでくる確率が上がるので（＾＾

i==2 ? x = 0 : x = 2
みたいなのは無理なんですか？

>>153
x=2; IF (i==2) x=0

>>154
なるほどサンクスです
でもやっぱ?使えた方が見やすいですね

>>153
これって何の言語ですか？

>>156
Ruby java

>>156
おしっこC言語。３項演算子。

しかし、昔風の書き方の気が・・・
最近はFortranでもCでも、もっとわかりやすく書くように指導されているはずだが・・・

配列の仮引数に、スカラーの実変数をsize=1の配列として渡すのって、
どうやればいいの？
下の例だと、rank mismatchでcompile通んないし、
reshape(s,1)もダメって言われた。
一回、別の配列に値をコピーするしかないのかな？

Program test
integer s

s = 1
call inc(1,s)
print *, s

Contains

Subroutine inc(n, a)
implicit none
integer, intent(in) :: n
integer, intent(inout) :: a(n)

integer i

do i = 1,n
a(i) = a(i) + 1
end do
End subroutine inc

End program test

>>159
INTENT(IN)　属性なら、
call inc(1,[s]) ないし　call inc(1, (/s/))
でいけるが、INTENT(OUT）で値が返る場合はダメだなあ。

integer s(1)
とサイズが１を明示すれば通るけど(intel のやつ)・・・それじゃダメなんだよねぇ・・・・
思いつきで（邪道とされるw）equivalence 使って

Program test
integer sdummy(1)
integer s
equivalence(sdummy, s)
s = 1
call inc(1,sdummy)
・・・・・
でも通るけど・・・意味ねぇなw

あ、始めて気づいたけど、
integer s(1)
s=1
って通るんだね。コンパイラにも寄るんだろうけど。

>>162
たぶんそれ、F90の配列初期化の一括代入になってるだとおもう。

ありがとう。

そっか普段77で書いてるから不慣れというかまったく気付かなかった。
たまに90で書くし初期化で使うのに・・・
Intel (たぶん余所のも) のはF77固定書式で書いても、
その事をコンパイルオプションで明示しないと、
F90の機能をちょっとなら混ぜても
文法が干渉しない限りは有効とみなすんだよね。
・・・危ないような気もする。

>>160
その場合でも結局、一時的な配列に値をコピーしてから渡すって意味になってるんだよ。

>>164
FORTRAN77なら型チェックなんかないし、
配列にスカラー変数を渡すのは正しい操作なんだよね。
内部的にはメモリアドレスを渡してるだけなんだから。

>>161のequivalenceは、今回のprogram文では意味ないけど、
変数sがサブルーチンの仮引数のときに、
無用なコピーを作らず渡せる唯一の方法かもしれない。

>>165
素直にテンポラリを間にかますのが一番いいと思うけどさｗ

EQUIVALENCE　ありなら、なんでもアリだなｗ
POINTERという手はどうだ？
Cとの互換POINTERを使えば、RANKとかのチェックは逃れられる。

Program test
USE, INTRINSIC :: iso_c_binding
IMPLICIT NONE
integer :: s

s = 1
call inc(1, c_loc(s))

Contains

Subroutine inc(n, p)
implicit none
integer, intent(in) :: n
TYPE(C_PTR), intent(in) :: p
integer :: a(n)
integer i

CALL C_F_POINTER(p, a)

do i = 1,n
a(i) = a(i) + 1
end do
End subroutine inc

End program test

>>159みたいなことは御法度なんじゃないの？
そのために90以降いろいろと拡張されたわけでしょ？

質問です。
LAPACKを使うプログラムなんですが、
Intelのコンパイラだとコンパイル出来るのですが、
gfortranだとSyntax errorが出てしまいます。
こういうことってプログラムに間違いが無くても起こりますか？

>>168
起こる。

>>167
そうなんだけど、その辺は７０年代的Pascalマンセー的な発想で、型チェックギチギチにしたら
確かにエラーは減ったが、窮屈で困ってしまったというのが、８０年代以降というか。

実際、Fortran的にはF95でELEMENTAL型の副プログラムが定義出来るようになって、
スカラーで副プログラムを定義しておけばRANKによらずに配列引数が取れるようになっている。

アルゴリズム部分とデータ構造を、もっと独立したユルユルの関係にしたいというか。

>>169
そうですかー…
ありがとうございます

>>168

その LAPACK はだれがどのコンパイラでコンパイルしたの？

>>170
気持ちは分かるなぁ・・・
この辺の善し悪しは開発の規模にもよるんだろうな

>>172
LAPACKはLinuxのパッケージ管理ソフトで入れたんですが、
だからたぶんgccでコンパイルされたものだと思います。

Fortranerの中には数値計算な業界の人も多いと思うけど、ぶっちゃけ聞きますが飯食えてますか？
ニッチでそこそこやっていけてるのか、それとも潰しのきかない悲しい境遇なのか

受験者数と受験者の受験番号と点数が分かっている。
10
17 73
5 64
22 67
6 66
7 90
42 88
9 74
11 79
2 63
13 72
（1行目は受験者数、2行目以降の左は受験番号、右は点数）
以上のデータdata.txtを読み込み、以下に示す結果をresult.txtに出力するようなプログラムを作成せよなお、この試験では59点以下はD、60点以上69点以下はC、70点以上79点以下はB、80点以上はAと判定することにする。
----------------------------------------------
受験者　　　x人
A　　　x人
B　　　x人
C　　　x人
D　　　x人
平均点　　　　　x
標準偏差 x
中央値 x
順位　　受験番号　　　点数　成績
1 25 85 A
2 41 82 A
（以下省略）
----------------------------------------------
・xには数字が入る。
・中央値はMOD関数を用いること。

以上の問題が分かりません。誰かよろしくお願いします。

e^x,sinx,logxのテーラー展開を求めるプログラムを教えてもらいたいのですが…
e^xはできたので、これを参考に教えてください！

subroutine expotn (x)
c
real Tn, ans, err
integer n
c
EPS = 1.0e-6
n = 0
Tn = 1.0
ans = Tn
write(6,10) x
10 format('EXP(', f10.5, ') ')
50 if( abs(Tn) .gt. EPS) then
n = n + 1
Tn = Tn * x / n
ans = ans + Tn
err = ans - exp(x)
write(6,11) n, ans, err
11 format(I3,' exp(X) = ',f12.7,' err = ',f15.10)
go to 50
end if
write(6,*) ' Exp( ',x,' )= ',ans
return
end
c

MinGWにくっついていたG95コンパイラを用いて、
Fortran77のソースコードをコンパイルしようとすると
以下のようなエラーが出ます。
コンパイルするときのコマンドは以下のとおりです。
>g95 F:\Direct_vs_Ewald.f -freal-loops
問題のソースコードは以下の部分です。
volume=a(1)*(b(2)*c(3)-b(3)*c(2))-a(2)*(b(1)*c(3)-b(3)*c(1))
& +a(3)*(b(1)*c(2)-b(2)*c(1))
この処理の部分で&の下に1が出現します。
Error: Unclassifiable statement at (1)
なぜ、このようなエラーが出るのかわかりませんでした。
()の閉じ忘れは確認したところありませんでした。
&記号が二つの文をひとつにする働きがあることは調べたら出てきましたが、
G95コンパイラにはこのような機能は無いということでしょうか？

すみません。いったん寝ます

>>177
SIN(x) = x - x^3 / 3! + x^5 / 5! - x^7 / 7! +.....

subroutine sinetn(x)
c
real Tn, ans, err
integer n
c
EPS = 1.0e-6
n = 0
Tn = x
ans = Tn
write(6,10) x
10 format('SIN(', f10.5, ') ')
50 if( abs(Tn) .gt. EPS) then
n = n + 1
Tn = - Tn * x**2 / ( (2 * n + 1) * 2 * n )
ans = ans + Tn
err = ans - sin(x)
write(6,11) n, ans, err
11 format(I3,' Sin(X) = ',f12.7,' err = ',f15.10)
go to 50
end if
write(6,*) ' Sin( ',x,' )= ',ans
return
end
c

>>177
LOG(1 + x) = x - x^2 / 2 + x^3 / 3 - x^4 / 4 +......
収束半径が小さいので注意。一般の範囲でやりたいなら、LOGの性質を使って範囲を調整汁。

subroutine logtn(x)
c
real Tn, ans, err
integer n
c
EPS = 1.0e-6
n = 1
Tn = x - 1.0
ans = Tn
write(6,10) x
10 format('Log(', f10.5, ') ')
50 if( abs(Tn) .gt. EPS) then
n = n + 1
Tn = - Tn * ( x - 1.0 )
ans = ans + Tn / n
err = ans - log(x)
write(6,11) n, ans, err
11 format(I3,' Log(X) = ',f12.7,' err = ',f15.10)
go to 50
end if
write(6,*) ' Log(',x,' )= ',ans
return
end
c

>>178
g95で試してみたが、６カラム目に＆がくれば問題なくコンパイルできた。
FORTRAN77の継続行指定は６カラム目に無いと駄目。

ずれると、
Error: Unclassifiable statement at (1)
が出る。

エラー行と数字の１が出ているはずで、それがエラー位置を指しているので、よく見よ。


ちなみに、Fortran90では＆は継続行を意味するが、１行目の尻に書く。
詳しくはマニュアルを読んでくれ。

純粋なFORTRAN77では＆は規格外の文字なので本来は使っていけない。
（今となっては余り気にしないが。）

>182
ありがとうございます。
Tabでインデントしていたので気づきませんでした。
半角スペース5個をつけた後、&をつけるようにしたら、そのエラーは消えました。
Warning (121): COMMON block 'coor' is 1600004 bytes at (1) and 160004 bytes at (
2)
というエラーは出ましたが、-wオプションをつけることでなんとかコンパイルはできるようになりました

ファイルをdirect accessで読む場合、
open (10,file='hoge.dat',access='direct',recl=4)
read(10,rec=10000) j
みたいに欲しい場所をすぐ読めますが、sequentialで同じようなことは可能でしょうか？

今まではこんな風に不要な部分は破棄していたんですが、効率が悪いと思いまして。
open (10,file='hoge.dat',access='sequential')
do i=1,9999
read(10) jtmp
enddo
read(10) j

>>184
末尾だったらOPEN時に指定出来るが、任意位置には行けない。
というか行きようが無い。
直接型はレコード長が全要素で共通だから特定の位置までのバイト数がわかるが、
レコード長が不定長の逐次型では原理的に知り得ない。

ところで読み飛ばすときjtmpに読み込まなくても空READだけでおｋ。
あと後ろに近いなら、尻から始めてBACKSPACEで逆行する技もある。
それが負荷を減らすかどうか知らんｗｗｗｗｗ

>>185
やっぱり不可能ですか。
将来ウンTBのデータを扱う時代になったら、データハンドリングはどうなるんでしょうね。
ファイルはダイレクトアクセス、データ構造はレコード長まで完全に把握するというのが
一般的になるんでしょうかね・・・

>>176
受験者 10 人
A 2 人
B 3 人
C 4 人
D 1 人
平均点 68.00000
標準偏差 19.24578
中央値 69.50000
順位　　受験番号　　　点数　成績
1 7 90 A
2 42 88 A
3 11 79 B
4 9 74 B
5 13 72 B
6 22 67 C
7 6 66 C
8 5 64 C
9 2 63 C
10 10 17 D

PROGRAM test
PARAMETER(nmax = 999)
INTEGER id(nmax, 2), idistr(4)
CHARACTER grade

OPEN(10, FILE = 'data.txt' , STATUS = 'old' )
OPEN( 9, FILE = 'result.txt', STATUS = 'unknown')

DO 10 i = 1, nmax
READ(10, *, END = 99) id(i, 1), id(i, 2)
10 CONTINUE
STOP 'error: too many data! increase nmax'

99 ndata = i - 1
sum = 0.0
DO 20 i = 1, ndata
sum = sum + id(i, 2)
20 CONTINUE
ave = sum / ndata

sum = 0.0
DO 30 i = 1, ndata
sum = sum + (id(i, 2) - ave)**2
30 CONTINUE
rms = SQRT(sum / ndata)

CALL sort(nmax, ndata, id)

DO 40 i = 1, 4
idistr(i) = 0
40 CONTINUE

DO 50 i = 1, ndata
IF ( grade( id(i, 2) ) .EQ. 'A') idistr(1) = idistr(1) + 1
IF ( grade( id(i, 2) ) .EQ. 'B') idistr(2) = idistr(2) + 1
IF ( grade( id(i, 2) ) .EQ. 'C') idistr(3) = idistr(3) + 1
IF ( grade( id(i, 2) ) .EQ. 'D') idistr(4) = idistr(4) + 1
50 CONTINUE

IF (MOD(ndata, 2) .EQ. 1) THEN
xmean = id(ndata / 2 + 1, 2)
ELSE
xmean = ( id(ndata / 2, 2) + id(ndata / 2 + 1, 2) ) / 2.0
ENDIF

WRITE(9, *) '受験者 ', ndata, '人'
WRITE(9, *) 'A ', idistr(1), '人'
WRITE(9, *) 'B ', idistr(2), '人'
WRITE(9, *) 'C ', idistr(3), '人'
WRITE(9, *) 'D ', idistr(4), '人'
WRITE(9, *) '平均点 ', ave
WRITE(9, *) '標準偏差', rms
WRITE(9, *) '中央値', xmean

WRITE(9, *) '順位　　受験番号　　　点数　成績'
DO i = 1, ndata
WRITE(9, '(i5, 2i11, a5)') i, id(i, 1), id(i, 2), grade(id(i, 2))
END DO

STOP
END
C

SUBROUTINE sort(nmax, n, in)
INTEGER in(nmax, 2)
DO 10 i = 1, n
DO 20 j = 1, n
IF ( in(j, 2) .LT. in(i, 2) ) THEN
itmp1 = in(i, 1)
itmp2 = in(i, 2)
in(i, 1) = in(j, 1)
in(i, 2) = in(j, 2)
in(j, 1) = itmp1
in(j, 2) = itmp2
END IF
20 CONTINUE
10 CONTINUE

RETURN
END

FUNCTION grade(k)
CHARACTER grade
IF (k .GE. 80) THEN
grade = 'A'
ELSE IF (k .GE. 70) THEN
grade = 'B'
ELSE IF (k .GE. 60) THEN
grade = 'C'
ELSE
grade = 'D'
END IF
RETURN
END
考えなしに書いた。なんなん(・∀・)！？

残差の二乗和を計算するときの精度と実行速度について悩んでます。
速度よりは精度の方を優先したいのですが、
以下の3つのソースで、最もバランスの良いものはどれでしょう？

xとaはそれぞれ倍精度実数として、

1.
(x-a)**2

2.
(x-a)*(x-a)

3.
(x-a)*(x+a)

よろしくお願いします。

指数部が整数の時は同じ掛け算をする(はず)。というわけで精度＆実行速度の点で、ぼほ１＝２と思ってよろし。
でもたとえば、a=1.0 でデータxを３要素で値が(1e-6,1,1e6)とすると、
第3要素のおかげで他の２項の寄与が消えちゃう。
ので、x^2 + 2ax + a^2
の３項にわけで個別に和をとって最後に足す、という事をすると
2ax の部分が残る分マシ。参照値 a とデータ x の分布次第だね。
x が a と同じオーダーなら気にすることは無い。

３は、 x*x - a*a の方が一般には精度という点で推奨される
同じく参照値 a と比べたときの x の値の分布次第、でもあるが
速度については積の演算はCPU毎にコストが違うので、どちらが良いかは一概にはいえない。

>>181
ちょっとまた質問です。
例えば、x-1.0=tと自分で置いたら

Tn = x - 1.0 → Tn = t

write(6,10) x →　write(6,10) t

Tn = - Tn * ( x - 1.0 ) ← Tn = - Tn * t

とおいて、より簡単なプログラムにしてもいいですか？
なるべくなら、式が簡単なほうが分かりやすいですし。

>>193
それでもよい。
なお、そのテイラー展開の収束半径は、|x|＜１のはずなので、その範囲から外れるような値も
計算できるようにしたいなら、LOG（X＊１０＾N)＝LOG（X)＋N＊LOG(１０）　などの公式を使って
スケールしてやる必要がある。

仮にt=1.5の計算をするときには範囲のプログラムとして>>194のようなプログラムを施すということでよろしいのでしょうか？

>>195
もう一個スケール処理するサブルーチンを書いて、それがテイラー展開のサブルーチンを呼べばいいんでないかい？
サブルーチンからサブルーチンを呼ぶ。

>>196
つまりは例として出ているLOG（X＊１０＾N)＝LOG（X)＋N＊LOG(１０）という関数のサブルーチンをlog(x)の範囲のところに呼ぶというわけですか？

こんにちは。sinxのテーラーのプログラムを以下のように書きましたが、うまくいきません…
おそらく、何か足りないのでしょうが…教えてください。
プログラム(1/2)

program sinTelor
c
real x
integer which
c
c begin
c EPS = 1.0e-5
write(6,*) 'Taylor Expansion of mathematical function'
write(6,*) ' .Sine function'
read (5,*) which
end

プログラム(2/2)
c

subroutine sine (x)
c begin
real Tn, ans, err
integer n, nstop
c
nstop=30
EPS = 1.0e-6
n =1
Tn = x
ans = Tn
write(6,10) x
write(6,*) ' sin(x)= ',sin(x)
10 format('SIN(', f10.5, ') ')
50 if( abs(Tn) .gt. EPS .and. n .lt. nstop) then
n = n + 1
Tn = - Tn * x **2 / ( (2 * n -1) *( 2 * n-2 ))
ans = ans + Tn
err = ans - sin(x)
write(6,11) n, ans, err
11 format(I3,' sin(X) = ',f12.7,' err = ',f15.10)
go to 50
end if
write(6,*) ' sin( ',x,' )= ',ans
return
end
return
end

>>198>>199
それは君が馬鹿だからです

>>198
まずこの　綴りの間違いを何とかしろ！ｗ
program sinTelor

CALL sine( which )

を書け。サブルーチンを呼び出さずに何をする気だｗ

兄貴がサブを呼ぶと覚エロ！

>兄貴がサブを呼ぶ
解説プリーズ

すみません（汗）
FORTRAN77で、Gauss-Jordan法のプログラムを
「READ」「WRITE」「GO TO]「IF]「SUBROUTINE」
「DIMENSION」を使って作りたいんですが・・・
全く分りません。誰か教えてください！

>>203
図書館行ってアルゴリズム集から完コピ。マジオススメ

>>203
誰か、F77で掃き出し法のサブルーチン作ってテンプレに書いておけよ。

>>203
goto は使うな！

>>206
それ７０年代の話ｗｗｗ

>>207
え！？

>>203
dimensionも使うな！

微笑ましいよね。
きっと課題に「xxx」「yyy」を使え、と明示してあるんだろうね。

77でgotoを使わないのは、結構苦しいと思うが。

全く無しはしんどいよね。最小限にはすべきだけど。
do while
があれば goto を使う機会が激減する・・・と書こうとしたけど、
お題の203さんの課題では do は入ってないね・・。
if goto の組み合わせで何とかしろ、ということなのか。
これはこれでしんどいｗ

>>203
まじめにアドバイス
Fortranがわからないの？Gauss-Jordan法がわからないの？

「知らない数値解法」で「知らない言語」で書け、
という課題が悩ましいのは自分もそうだったのでわかるけど
丸投げだとあまり回答もらえないですよ。

前者なら、どの程度の基礎的なプログラムなら書けるのかを提示して
どの部分を調べたらよいかをアドバイスもらい、

後者なら、ネットで調べたり、その解法のどの部分のプログラム化が？なのかを提示
したほうがいいと思います。
正直、その解法で、何を解きたいのかよくわからないし。

>>213
放っておいてあげれば良いと思うよ。

A〜Zまでの文字のデータを1つ1つファイルにしておき
小文字で入力された文字列をファイルから読み込んで
大文字で表示するにはどのようにすれば良いのでしょうか？

入力された文字をそれぞれのファイルから
読み込む方法がよくわかりません。お願いします。　Ｆ９０

質問があいまい＆意味不明．．．

ファイル内に記述された特定の文字だけを大文字に変換したいのかな？
もしそうなら、Fortranではなくawkやperlを使った方が断然楽だよ。

いちおう、ファイルを読み込むのも、入力文字を読み込むのもREADを使う。

>>215
完全イミフ。
日本語でおｋ

表示ができるようになったのですが、ファイルの行数を明示しないと
End Of Line になってしまうのですが、そこで
integer io を宣言して行末を発見して読み込みを止めさせようとしたが
うまくできないのです

なお　data614.txt は株価のデータで
価格（整数）,年,月,日
にしてあります。

program hyoujisurudake2
implicit none
integer :: i
integer :: io=0
integer :: xprice,xyear,xmonth,xday
!*4を*8にしたりして試していたけどできなくて、こうしたらできた
open(10,file='data614.txt',status='old')
i = 0
do
!だめだ read(10,iostat=io) ,xprice, xyear, xmonth, xday
read(10,*) ,xprice, xyear, xmonth, xday
! if(io <= -1 ) exit
if(i>=3406) exit
i=i+1
!2 print* , i,xprice
print* , i, xprice, xyear, xmonth, xday
end do
close(10)
stop
end

多分続けて投稿です
上記の株価データですが本当はこうなっているんです
これ
年/月/日,株価

具体的には
1999/12/31,12345
2000/1/1,23456
2000/1/2,23456
てな具合ですがここで大きな問題がある　それは
読み込むときに「,」切りはいいとして
「/」での切り方がわからない・・・
fortran だと何文字をどこそこに格納して　とかいうのはとくいですが
パーサーが無い？　ので「,」で切ってから「/」で切るという処理は
どうやるのでしょうか

>>218
IOSTATがうまくいかないのは、Fortmat書式を抜いているから。


>>219
INDEX関数使えば出来る。

>>219
program hyoujisurudake2
implicit none
integer :: i, io, ip1, ip2, ip3
CHARACTER(LEN = 80) :: text
integer :: iprice, iyear, imonth, iday
open(10,file='data614.txt',status='old')
i = 0
do
read(10, '(a)', iostat = io) text
if(io <= -1 ) exit
i=i+1
ip1 = INDEX(text, '/')
ip2 = INDEX(text, '/', BACK = .TRUE.)
ip3 = INDEX(text, ',')
READ(text( 1:ip1 - 1), *) iyear
READ(text(ip1 + 1:ip2 - 1), *) imonth
READ(text(ip2 + 1:ip3 - 1), *) iday
READ(text(ip3 + 1: ), *) iprice
print * , i, iprice, iyear, imonth, iday
end do
close(10)
stop
end

1 12345 1999 12 31
2 23456 2000 1 1
3 23456 2000 1 2
続行するには何かキーを押してください . . .

この入力データで問題となるのは、FORTRANの書式指定無しの場合、スラッシュ『/』は改行を意味する制御文字に
なっている点だ。これを知らないと、何故か入力がなされなくて、地獄へ落ちる！

>>203
gotoは１関数に１ラベルなら使ってOKだよ！
ソースが見やすくなるからね

>>203
PIVOT処理するのが面倒なのでしてないｗ　

SUBROUTINE GJ(a, b, n) ! Gauss-Jordan method
DIMENSION a(n, n), b(n)
DO 10 irow = 1, n
pivot = a(irow, irow)
IF (pivot .EQ. 0.0) STOP 'pivot is 0.0!'
DO 5 icol = 1, n
a(irow, icol) = a(irow, icol) / pivot
5 CONTINUE
b(irow) = b(irow) / pivot
DO 20 jrow = 1, n
IF (irow .NE. jrow) THEN
q = a(jrow, irow)
DO 30 icol = 1, n
a(jrow, icol) = a(jrow, icol) - q * a(irow, icol)
30 CONTINUE
b(jrow) = b(jrow) - q * b(irow)
END IF
20 CONTINUE
10 CONTINUE
RETURN
END

>>219
2000年元旦に市場で何があったんだ・・・w
>>221
/ , を空白に置き換えてから内部入力で、というやりかたもできそうだけど、どうかな？
text(ip1:ip1) = ' '
(略)
READ(text,*) iyear, imonth, iday, iprice
手元にコンパイラがないから・・これで勘弁

>>222
俺もそうだよw

224だけど、どうかな？

>>224
>/ , を空白に置き換えてから内部入力で、というやりかたもできそうだけど、どうかな？
その方がいいな。　↓コンマは区切り文字になっているから残してもいいのだが、INDEXの代わりにSCAN関数使ってみたかったので消した。
program hyoujisurudake2
implicit none
integer :: i, io, ip1, ip2, ip3
CHARACTER(LEN = 80) :: text
integer :: iprice, iyear, imonth, iday
open(10,file='data614.txt',status='old')
i = 0
do
read(10, '(a)', iostat = io) text
if (io <= -1 ) exit
i = i + 1
CALL cutslash(text)
READ(text, *) iyear, imonth, iday, iprice
print * , i, iprice, iyear, imonth, iday
end do
close(10)
stop
CONTAINS
SUBROUTINE cutslash(text)
CHARACTER(LEN = *), INTENT(IN OUT) :: text
INTEGER :: k
DO
k = SCAN(text, '/,')
IF (k == 0) EXIT
text(k:k) = ' '
END DO
RETURN
END SUBROUTINE cutslash
end program hyoujisurudake2

>>226 コンマは区切り文字になっているから残してもいい

本当だ！
いまの今まで知らんかった・・・。 CSV みたいに','が区切り記号で使われてる
ファイルもらったときも御丁寧に置き換えてた・・。

write文での制御編集記述子についての質問です。

100e18.10
と指定してファイルに書いていたのですが、
これを
100es18.10e2
と指定すると、一部がアスタリスクになってしまいます。
マニュアルを読むかぎり
100が書式の繰り返し数
18が小数点上の桁、小数点以下の桁、指数を含めたすべての桁数
10が小数点以下の桁数
2が指数の桁数
と理解したのですが、どこか誤りがあるでしょうか。
教えていただけないでしょうか。

218 219 です
>220
さん　ありがとう

>221
さんコードまで書いてもらってありがとうございます。
うちのPCだとコンパイルがうまくいかないようですが、
それは多分調整のレベルなのでしょう。

ほんとにこんなに早くレスがくるとは思ってもみず　返事が遅くなりました

>>224
> 224 名前：デフォルトの名無しさん [sage]： 2010/07/04(日) 06:34:23
> >>219
> 2000年元旦に市場で何があったんだ・・・w

すみません、これ分かりやすくするためのものでして・・・。
「文字の数が変わっちゃうのでこまっていた」ということでした。


> >>221
> / , を空白に置き換えてから内部入力で、というやりかたもできそうだけど、どうかな？
> text(ip1:ip1) = ' '
> (略)
> READ(text,*) iyear, imonth, iday, iprice
> 手元にコンパイラがないから・・これで勘弁

いろいろ工夫が必要なのが分かって着ました。
高級言語の割には今時のRubyとかよりずっとマニュアル車っぽくていいようなかんじも・・

ありがとうございました。

これで本題の数値解析に移れます

>>226
> 226 名前：221 [sage]： 2010/07/04(日) 14:20:20
> >>224
> >/ , を空白に置き換えてから内部入力で、というやりかたもできそうだけど、どうかな？
> その方がいいな。　↓コンマは区切り文字になっているから残してもいいのだが、INDEXの代わりにSCAN関数使ってみたかったので消した。
> program hyoujisurudake2

218　です　いろんなやり方があるんですね。
それだけ皆さんは苦労してきたと想像します

>>227
> 227 名前：デフォルトの名無しさん [sage]： 2010/07/05(月) 04:43:04
> >>226 コンマは区切り文字になっているから残してもいい
>
> 本当だ！
> いまの今まで知らんかった・・・。 CSV みたいに','が区切り記号で使われてる
> ファイルもらったときも御丁寧に置き換えてた・・。

219　です
すると "," は標準で区切り文字になっているが、それを解除？することは
難しそうですね。きっとできるでしょうが、とっても深い世界

はい

>>228
Fortranの場合、Formatで指定した桁数からあふれると、アスタリスクが出力される。
今の場合、多分指数部分が溢れているんじゃないかと思う。
E-10とかが出ると、ヤバイ。E2の2には符号の分も含まれる。

>>232
＞すると "," は標準で区切り文字になっているが、それを解除？することは
＞難しそうですね。きっとできるでしょうが、とっても深い世界
それはFORMATをあらわに指定すれば出来る。
上の例でも、ファイルの１行を丸々文字列として読み込んでいる。
read(10, '(a)', iostat = io) text
この場合は、スラッシュもコンマもただの文字として読み込まれる。

読み込まれる。

>>235
ひがむなって。百姓・町人みたいだぞ。

>>234
ありがとうございます。
>>228
にある18を25などの大きな値にしてみたのですが、それでも一部でアスタリスクになってしまうことがありました。

一行にかける文字数などなにか別の制約があるのでしょうか。

>>237
全体の幅を増やせば仮数部は大丈夫になるが、その場合でも指数部の幅は変わらないので、指数部が問題なのでは？

>>237
いいよいいよ。

部分ピボッティング付きのガウス−ジョルダン法。これで対角要素が０でもダイジョブのはず？

SUBROUTINE GJ(a, b, n) ! Gauss-Jordan method with partial pivoting
DIMENSION a(n, n), b(n), temp(n)
DO 10 idiag = 1, n
CALL pivoting(a, b, n, idiag)
pivot = a(idiag, idiag)
DO 5 icol = 1, n
a(idiag, icol) = a(idiag, icol) / pivot
5 CONTINUE
b(idiag) = b(idiag) / pivot
DO 20 jrow = 1, n
IF (idiag .NE. jrow) THEN
q = a(jrow, idiag)
DO 30 icol = 1, n
a(jrow, icol) = a(jrow, icol) - q * a(idiag, icol)
30 CONTINUE
b(jrow) = b(jrow) - q * b(idiag)
END IF
20 CONTINUE
10 CONTINUE
RETURN
END

SUBROUTINE pivoting(a, b, n, k)
DIMENSION a(n, n), b(n)
temp = 0.0
ipiv = k
DO 10 i = k, n
IF (ABS(a(i, k)) .GT. temp) THEN
temp = ABS(a(i, k))
ipiv = i
END IF
10 CONTINUE
DO 20 i = k, n
temp = a(ipiv, i)
a(ipiv, i) = a(k, i)
a(k, i) = temp
20 CONTINUE
temp = b(ipiv)
b(ipiv) = b(k)
b(k) = temp
RETURN
END
メインルーチン：ただしｆ９０機能を使用ｗ
　　　PROGRAM test
PARAMETER( n = 5 )
REAL a0(n, n), a(n, n), b0(n), b(n)
CALL RANDOM_SEED()
CALL RANDOM_NUMBER(a0)
CALL RANDOM_NUMBER(b0)
a = a0
b = b0
CALL GJ(a, b, n)
PRINT *, SQRT( SUM( (MATMUL(a0, b) - b0)**2 ) / N )
END

>>238
ありがとうございました。ようやくおっしゃる意味がわかりました。
指数の部分を2桁にしてことが原因ということですね。
納得しました。明日試してみます。

>>221
> 221 名前：デフォルトの名無しさん [sage]： 2010/07/03(土) 01:44:30
> >>219
> program hyoujisurudake2
> implicit none
> integer :: i, io, ip1, ip2, ip3
> CHARACTER(LEN = 80) :: text
> integer :: iprice, iyear, imonth, iday
> open(10,file='data614.txt',status='old')
> i = 0
> do
> read(10, '(a)', iostat = io) text
> if(io <= -1 ) exit !★ここです
> i=i+1
> ip1 = INDEX(text, '/')
> 中略
> end do
> close(10)
> stop
> end
-----
> 3 23456 2000 1 2
１６行目の if(io <= -1 ) exit !★ここですが「-1」を「0」とか他の値にしても出ないかEnd of File　のエラーになるのですが、これはUbuntu910 + gfortran だからでしょうか？
いくつを入れてもうまくいかないです

>>243
自己レスです
すみません　データーの最後の行に「 / / , 」という空白でータ行が入っていた！
だから最後まで読むけど　終わりを過ぎてしまう　というエラーが出ていた　おろかしー

3元以上、9元までの連立一次方程式を解くプログラムをg77で作成しているんですが、なかなかエラーなしにすんなりとソースプログラムを作ることができません。
ガウス・ジョルダン法を利用して作れと指示があり、一応作ってみたのですが、エラーがいたるところに出て、エラーの内容もわからないため、訂正をお願いしたいです。
なにせFORTRANを使い始めて間もないため、超初心者でもここまでひどいソースは書かないとおもわれるほどのひどいソースかもしれませんが、どうかお願いします。
Cygwin 1.7.5-1 の　fortranコンパイラを使用しています。以下、作った未完成ソースプログラムです。
real i,j,N,k,aij,bij
read(5,*) N
DO 100 i=1,N
DO 110 j=1,N
read(5,*) aij
110 continue
100 continue
DO 120 i=1,N
read(5,*) bi
120 continue
call GJ(aij,bi)
130 write(6,*) bi
end
subroutine GJ(aij,bi)
DO 130 k=1,N
akj=akj/akk
if(i=k)
aij=aij
else
aij=aij-akj*ajk
if(i=k)
bk=bk/akk
else
120 bi=bi-bk*aik
DO 130 i=1,N
return
end

GJのアルゴリズムは確認していないけど、とりあえず配列の作り方を確認してみ。
配列変数には()が必要。

あと、エラー内容をコピペしてみて。

エラーをがんばって見て作りなおしてみましたがエラーがでます。
以下、ＧＪのサブルーチンを抜かしたソースです
　　　　implicit none
integer s,t,n,k
real a,b
read(5,*) n
dimension a(n,n),b(n)
DO 100 s=1,n
DO 110 t=1,n
read(5,*) a(s,t)
110 continue
100 continue
DO 120 s=1,n
read(5,*) b(s)
120 continue
連投になって申し訳ないですが、次にエラーのコピペします。

>>246
先ほど言い忘れましたが、レスありがとうございます。
全部コピペは入りきらないので最初のエラーかつ一番多い文の一つをコピペしました。
内容は、辞書で調べながらやったのですが、宣言が間違っていると言ってるように思います。
114.for: In program `MAIN__':
114.for:2:
integer s,t,n,k
1
114.for:5: (continued):
dimension a(n,n),b(n)
2
Invalid declaration of or reference to symbol `n' at (2) [initially seen at (1)]
お願いします。

> implicit none
> integer s,t,n,k
> real a,b
> read(5,*) n
> dimension a(n,n),b(n)

ここの部分が間違っている。

implicit none
real a(n,n),b(n)
integer s,t,n,k
read(5,*) n

と書くのが正しい配列の定義の仕方。このばあい、属性を意味する"dimension"は書かなくてok

あと注意点だけど、変数定義の途中に実行可能な文を入れてはダメ。
>>247の場合dimensionの前にreadがはいってるけど、readは変数をすべて定義した後にやる。

あ、間違えた。

配列は定義する時点でサイズが決定していないとだめなんだ。
配列サイズをreadする場合には、動的なメモリの割付けが必要で、allocatable属性を指定しなければならない。
readでnを呼んだ後、allocate文で実際にメモリを割り当てる。

implicit none
integer s,t,n,k
real,allocatable a(:,:),b(:)
read(5,*) n
allocate( a(n,n), b(n) )

>>250
多分質問者は、F7７でかつDIMENSIONとGOTOを使えと言われている人なので、その辺を察してあげないと・・・

>>251さんの言う通りで…質問してるのに申し訳ないです。
allocatableというのは聞いたことがないので、使わずにやれということなんだと思います。
>>248,>>249を参考に最初の5行を訂正してみました。

implicit none
integer s,t,n,k
real a(n,n),b(n)
dimension a(n,n),b(n)
read(5,*) n

エラーは>>248と全く同様に出力されました。
ちなみにa,b,n,iに関して全く同様のエラーが何回も出力されています。
日をまたいでレス遅れてしまって申し訳ないです。
よろしくお願いします。

>>252
訂正:a,b,n,iに関して→a,b,iに関して
連投失礼します。

宣言部は全て最初にやらないといけないんだよなぁ
コード真ん中のスコープ内でローカルに宣言できたら便利なのにといつも思う

あ、でも次のFortranで実装するって話を聞いたことがあるような
次の次だっけかも

allocatableがつかえないと、余分にメモリを食っちゃうね。
まず、配列変数の定義の時点で
real, dimension(100,100) :: a
real, dimension(100) :: b
と、ある程度大きいメモリ容量を割り付けておいて、すべての変数に0を代入して初期化。
a(:,:) = 0.0
b(:) = 0.0

そのあとnやa,bの要素をreadすればokかな？
ここで、nは単なる繰り返し制御のための変数になるけど、まぁ仕方ないか。

>>255
real, dimension(100,100) :: a
real, dimension(100) :: b
をそのまま利用すると、Fortran 90 feature at (^) unsupported　(対応してない?)といわれ、
realとdimensionを分けてa,bを定義すると
Invalid form for DIMENSION statement at (^)　(dimensionの宣言形式が違う?)
とdimensionを定義する箇所にエラー
a(:,:) = 0.0
1 2
Unrecognized statement name at (1) and invalid form for assignment or statement-
function definition at (2)
a,bに関して同様のエラーがでます。。。
もはやcygwinのfortranコンパイラがダメなきがしてきました。。
何度も申し訳ないです。

いやいや、fortran90形式の方法を教えてしまった。
90形式と77形式を混ぜてコンパイルしても問題にならないことが多いから、
いいかと思ったけどだめだったみたいだね。ごめん

implicit none
integer i, j
integer s,t,n,k
real, dimension(100,100) a
real, dimension(100) b
do i = 1, 100
do j = 1, 100
a(i,j) = 0.0
end do
b(i) = 0.0
end do
read(5,*) n

Fortran77形式で書いてみた。古い書き方は全然使わないから、間違っているかも。

変数宣言の仕方って何が推奨なんですか？
ttp://www.nag-j.co.jp/fortran/FI_4.html　を見ると、数に関しては

整数型　 integer
実数型　 real
倍精度実数型　 double precision
複素数型　 complex
倍精度複素数型　 complex(kind(0d0))

となっていますが、実数型はreal*4とすべしと言ってるサイトもありますね。
double precisionなんて長ったらしいし、*で済むのならそれがスマートだと
思うんですが、どうですか？

>>258
そのページのすぐ下に書いてあるが、REAL*4とかREAK*8とかは、FORTRANの規格の規定にはない書き方。
ただし、事実上通用する書き方に成っている。
だから規格に厳格に乗っ取りたいという流派なら、使わないという選択になる。
またFortran90以降では、REAL(KIND=4)や省略形REAL(4)などがあるので、そちらを使うのがよりスマート。
そういう文脈を踏まえて、７７時代の慣用的な書きかたをしようというなら、REAL*8とかで書くのが主流派だと思う。

>>257
７７にはIMPLICIT NONEは無い。
またreal, dimension(100,100) a　この書き方もない。
REAL　A
DIMENSION　A(100,100)
と分けて書かねばならんはず。

*の数字は変数が占めるトータルのバイト長らしいから、複素数の場合は注意だな
やはりkind指定が無難かな

以下、intel fortran 11.0

complex :: c1
complex*8 :: c2
complex(4) :: c3
complex(kind=4) :: c4
complex(kind(0.0)) :: c5

print *, c1
print *, c2
print *, c3
print *, c4
print *, c5

end

＜出力＞
(0.0000000E+00,0.0000000E+00)
(0.0000000E+00,0.0000000E+00)
(0.0000000E+00,0.0000000E+00)
(0.0000000E+00,0.0000000E+00)
(0.0000000E+00,0.0000000E+00)

Fortranの命名規約って何か習慣的なものはあるんでしょうか？
例えばJavaなどではキャメルケースが一般的ですが

ない

>>261
キャメルケースってボーランド記法のことか。
Fortranは昔は大文字しか使えなかったし、今も大文字小文字の区別がないので、そういう工夫はあまりしない。
また１行の長さに制限があるので、余り冗長な名前も使わない。

７７時代までは極めて明快だった。
暗黙の型は厳守。型を守ってないものは、文法的に許されてもエラーの元だし、社会的に誰からも相手にされないので。
名前は６文字までなので（IBMがこだわったので）単語の母音を抜いた子音並べが基本。

この手の命名法は、LAPACKなどライブラリに厳然と残っているので、知っておくと便利。
頭文字がAとかRは単精度実数、Dは倍精度実数、Cは単精度複素数、Zは倍精度複素数（本来は規格外だが）。

F90になってから混迷を深めているが、いちいち変数の型を調べるのに、宣言部分を見る気もしないので、
捨て変数に関してはIMPLICIT NONEを掛けても暗黙型の頭文字をつかっているのでは。
LOOP変数のｉ，ｊとか。

変数の命名規則は、開発の仕方や規模によるところが多いんじゃないかな。
他人のコードを頻繁に見るような場合は、パッと見で分かるようなのにしてくれないと死ぬ・・・。

構造型(体?)を使ってみたいのですが、良い参考書等あるでしょうか？
FORTRANは本自体少ない上、FORTRAN90以降の新機能について
詳しく言及している本は少ない気がしますが。

>>265
ない！手持ちの本で十分だ。がんばれ。

>>266
その手持ちの本に構造体書いてりゃな。
書名は何？

>>265
（構造体は必要）ない！手持ちの本で十分だ。がんばれ。

>>267
F90の本で構造体が書いてないものなんかあるのか？
戸川隼人のしょぼいザ・Fortran90/95にも１章費やして書いてあるぞ。
と思ったら、パラパラと見た感じ牛島省の本には無いな。京大土建屋の英単語つづりデタラメ本だから仕方ないかｗ
同じ京大の富田の本には２ページくらいおざなりに書いてある。
しかし、富田本は昔の大型計算センターのテキストみたいな分かりにくさの本だな。
イラストとしてキャラクター・グラフィックスで書いた３Dの山とリサージュ図形があれば完璧だったろう。
薄くていいが定価で買う気はしない。京都の古本屋で百円で売ってそう。

構造体で注意すべきなのは、Fortran90ではALLOCATABLE属性は使えなかったので、POINTERを使うしかなかったが
Fortran95ではALLOCATABLE属性のメンバーを持てることくらいか。

>>268
マジレスするとｗｗｗ
Fortran2003でのオブジェクト指向は、構造体の中に副プログラムをCONTAINする形になるので、避けられない。

使ったことないからよくわからないけど、Fortran の構造型の命名って制限ある？
ピリオドで区切るから・・ .EQ. とかの比較演算子とかぶりそうな気がふとしたので。

>>270
区切りは％だぜ。

. % のどちらでも実質、両方正しいよ。
. 区切りはDECのだった気がするけど、まあ実質業界標準。

>>272
ねーよｗ
Fortran90以降の構造体とDECの構造体は別物だぞ。
.　が使えるわけがない。

>>269
＞構造体で注意すべきなのは、Fortran90ではALLOCATABLE属性は使えなかったので、

Fortran90では構造体の中のメンバはALLOCATABLE属性に出来ないって意味だよね？
それともALLOCATABLE自体がFortran95以降なんだっけ？

F90でも普通にallocatable出来る

>>274
F90では構造体のメンバーがALLOCATABLE属性を取れなかった。
これが大変不評でF95で取れるようになった。

なったね

たね

ね

すみません
初歩的な質問だと思うのですが、
以下のように、行を一つのベクトルと考え、自分の入力したベクトル間の距離を出そうとしたのですが、上手くいきません
どこがいけないのかアドバイスお願いします
最初に定義している文字が多いのは気にしないでください

program aaa
implicit none
integer :: i,j,m,k,o,p,q,zenntai
real :: c1x,c1y,c2x,c2y
integer,parameter :: n=5
real,dimension (n,2) :: a
real,dimension (n) :: b,c

do m=1,n
write(*,*)'a='
read(*,*) (a(m,k),k=1,2)
end do

write(*,*) 'a='
do i=1,n
write(*,'(1x,2F7.3)') (a(i,j),j=1,2)
end do

c1x=a(1,1)
c1y=a(1,2)
c2x=a(2,1)
c2y=a(2,2)

つづく

つづき
do o=1,n

call kyori(b(o),c1x,c1y,a(o,1),a(o,2))
write(*,*) b(o)

end do

do p=1,n

call kyori(c(p),c2x,c2y,a(p,1),a(p,2))
write(*,*) c(p)

end do

stop
contains

subroutine kyori(d,s,t,u,v)
implicit none
real,intent(in) :: s,t
real,intent(in),dimension(n) :: u,v
real,intent(out),dimension(n) :: d

d=((abs(s-u))**2.0)+((abs(t-v))**2.0)

return

end subroutine kyori
end program aaa

>>280
どう上手くいかんのかもう少し分かりやすく言えよと
subroutineにしないでやっても同じなわけ？

すみません

サブルーチンにしないでやっても結果は同じでした。
a=1　 2
1　 3
1　 4
1　 5
1　 6
とすると、とりあえずb(o)だけ書きます
b(o)=0.
1.
1.
394825
2.4050422E+22
c(p)もb(o)と同様にpが上がると凄い勢いでその値が増えていきました

>>283
２点ばかり考えられることを。
１。
入力がフォーマット付きなので、空白を区切りとして数値を二つ与えることはできない。
コンマでしきれば大丈夫だったか？
a=1　 2　＝＞１，２　もしくは　　　　　　　１　　　　　　２


２。
サブルーチン中のここがおかしい。
real,intent(in),dimension(n) :: u,v
real,intent(out),dimension(n) :: d

スカラーを渡しているのだろうから
real,intent(in) :: u,v
real,intent(out) :: d
のはず。

sqrtはなくてもいいの？

>>285
そういう内容に踏み込む場合、REALなんだから二乗する前にABSはとらなくてもいいはずとか、
突っ込みどころいっぱいありすぎ。
そもそもループ変数にoとかpとか使ってる時点で、こいつラリッてるのか、酔っ払ってんのか的気持ちｗ

質問です。
数値計算をしようとfortranでプログラムを書いて、実行してみたんですが
古いPen4のPCだとCPU使用率が100%までいくのに、新しいi7のPCだとCPU使用率が13%以上になりません。
とりあえず、元のプログラムを6個のプログラムに分けて実行してますがいい方法とは思えません・・
どうしたら１つのプログラムでCPUを有効に使えるようになりますか？

コンパイラはg95です

っ OpenMP & gfortran

>>288
ありがとうございます。やってみます

はじめまして。
少々悩んでいることがありまして質問させて頂きます。
行列積のプログラミングを組もうと思っているのですが、
入力部分に数字のみならず記号も含まれている場合どうすればよいかわかりません。
色々試してはみたのですが、結局無理でした。
例えば、
1 2
3 4
という係数行列と
a
b
というベクトルをかけて
a+2*b
3*a+4*b
というふうに出力したいと思っています。
どなたか心ある方、ほんの少しのアドバイスでもよろしいのでお願いいたします。

関数とサブルーチンってどっちが有利？
たとえば

function hoge(a,b) result(c)
c=a+b
return
end function

subroutine hoge(a,b,c)
c=a+b
return
end subroutine

と書いた場合、どっちがメモリ的に
あるいはスピード的に
あるいは可読性的に有利？

>>290
fortranに文字計算させるの？
数値計算ってそういう物じゃないと思うんだけどなあ

subroutine matvec(inmat,invec,outvec)
implicit none
intent(in):: inmat, invec
intent(out):: outvec
real mat(2,2), invec(2), outvec(2)
inteter i,j
outvec(1) = mat(1,1)*invec(1)+mat(2,1)*invec(2)
outvec(2) = mat(1,2)*invec(1)+mat(2,2)*invec(2)
return
end subroutine

・・・・なにか質問の答えとはずれてるな・・・w

function と subroutine の違いって戻り値があるか無いかじゃないのか？

>>294
そうなんだけど結果的には同じように扱えるよね？
上の例で言うと変数cがそれになると思うんだけど

可読性という点から見ると関数は、Fortranだと
配列と区別しにくいから少し解りづらいかも

質問です。

現在、ファイル名が
11016.csv
11151.csv
：
24141.csv
というものが、97個ほどあります。そのファイルを順にOPENして読み込み、作業したいのですが、下のようにopen文を書いたのですが、初めの11016の地点が、何度も出力されてしまいます。連番でないファイルの場合、
この書き方は適当ではないのでしょうか？それとも全体のプログラムを示していないですが、enddoを書く場所が間違っているなどのミスでしょうか？

また、read文なのですが、色々詮索した結果、このように行変更してしまいますと、2行以降の変数は読めていませんでした。ファイル内のデータを削除し、プログラムの変数を減らして一行に収めると、正しく出力されました。
なぜなのでしょうか？どうかアドバイスください。もし、より全体のプログラム内容だったり、読み込んでいるファイル内の形式が必要でしたら、
再度載せます。
----------------------------------------------------------------------
INTEGER :: year, mon, day, data,sum, no,id,pre,wind,winddirection,hightemplowtemp,maxrain,maxwind,code,tmpsunlight
CHARACTER(10) :: name
　　 REAL :: lon, lat
CHARACTER*5 sssss

do ispot = 11016,24141
write(sssss,"(i5)") ispot
open(50, file=''//sssss//'.csv', status='old',iostat=io)
if (io < 0) cycle
do i = 1,12000

read(50,*,iostat=io) year,mon,day,id,lon,lat,code,&
name,pre,wind,winddirection,temp,hightemp,lowtemp,&
tmpsunlight,maxrain,maxwind

----------------------------------------------------------------------

マルチポスト？...ではないか。向こうの質問内容は意味不明だ。
http://oshiete1.watch.impress.co.jp/qa6168042.html

cycle文を使う場所が間違っている気がする。
あとは、装置番号を50ではなくispotにして、readしたあとはclose(ispot)するとかかな。

この手のFortranのIOは面倒で、真面目に考える気になれない...

うん、変だよね・・・Cycle は使わんでもよさそうだし・・。
いろいろ副作用やら曖昧さをさけた流儀にすると・・

character*9 sssss
logical lexist
do ispot=11016,24141
write(sssss,'(i5.5,''.csv'')') ispot
inquire(file=sssss,exist=lexist)
if (lexist) then
write(*,*) 'now handling a file : ', sssss
open(unit=50,file=sssss,status='old')
do ..
read(50,*) ....
enddo
close(50)
endif
enddo

かな・・・。コンパイラ無しな状況で書いてるから自信はないね。

教えていただいたOPENファイルの書き方でうまくいきました！
すごいです。こんな方法あるんですね。 勉強になります。
助かりました。ありがとうございました。

えがったえがった

2GBを超えるunformattedなファイルを read(num) (itmp(i),i=1,n) という感じで（nは既知）
読むとエラーになったんですが、 readにオプション iostat=io を加えると大丈夫でした。
2GB以上はiostatが必要ということでしょうか？

なお、エラーメッセージは以下の具合です。
Intel Fortran・・・「input statement requires too much data」
PGI Fortran・・・「attempt to read/write past end of record」

iostat無しだとエラーを受けて終了してしまい、
iostatありだとエラーをエラーコードとしてioに保存してそのまま読み込みつづけただけじゃないの？

>>302
よく考えたら2GB以上だと仕様が変わるなんてあるわけないですよね。
ファイルにはn個データが入っているはずなんですが、もしかしたらデータがおかしいのかもしれません。

もう一度調べてみます。ありがとうございました。

>>301
書式無し、でサイズ既知、なら・・・
integer itmp(n)
で宣言した後、
read(num) itmp
と、添字なし表現で読めない？

>>304
read(num) itmp
read(num) itmp(1:n)
などと試しましたが同じでした。やはりデータが疑わしいようです。

・・あと思いつくのは・・
無書式だとメモリ上の０・１の値をそのままファイルにぶち込むだけ(な事がほとんど)だから
CPU や OS や コンパイラがファイルを生成したものと違う組み合わせで読み込むと
失敗することがあるよ。

人からもらったデータなら、元にお願いしてプレーンテキストにしてらうのが吉。
情報量は変わらないのだから適当な圧縮をかけるとファイルサイズはそう変わらないはずだ。

Endian とか Padding とかをキーワードにして、いろいろ試してね

文字列として読み込んでみるとか、

何かの値で初期化して
IOSTATつきで読みこんでどこでエラーが出ているか確認し、
そのデータ部をバイナリエディタで見てみるとか

プレーンってformattedで出すってことだよね
バイナリより数倍大きくなるよね

そういえばunfortmattedとバイナリって同じ意味なの？

>>309
今時２ギガが１０ギガになったくらいで騒ぐこともあるまい・・。
とりまわしが面倒なら分割すりゃいいんだし。

binary つかったことない・・
Unformatted は306みたいにメモリ上のならびを
そのままぶち込む・・と思ってたけど違うのかな？

>>310
最近の数値計算じゃテラのオーダーっすよ。
つーかこのスレにいる奴らもそっち系の人ばっかでしょ？

>>311
unformattedってデータとデータの間に区切りとしてなぜか8byte(?)ぐらい
埋め込まなかったっけ？？asciiで書き出すよりファイル容量でかくなって
ビビったような

Padding はありそう・・・とうわけでチェック！
単精度・倍精度・文字列・整数・理論、を100要素のものを無書式で書き落とすと
サイズは(バイトで)、
800,400,100,400,400＋８(intel fortran)または＋１６(gfortran)
だった。OSはLinux、こまかいDistriは知らんw
おまけがつくのはそれぞれのコンパイラの仕様かな？
ともあれ、整数・実数・文字列はなんとなく想像のうちだけど、論理(logical)は
本来１ビットだけど、８ビット単位にPadding しているようだ(まあ、その方がメモリ内での
とりまわしが効率的なんだろうね)・・・即興の実験終わり。

313さんのはどの種類の変数だったの？

>本来１ビットだけど、８ビット単位にPadding しているようだ
ああ、おお間違い。４バイトね。
中身をみたら当然だけど０パディングばかりでスカスカだった。

ああ、もしかすると、301 さんのは
データを書き出したコードのコンパイラと読み込みのが別、とか・・かな。
CPU OS コンパイラ の３つ次第でどうにでも変わるのが無書式だし。

単純に同じシステム内で処理をする分には(コーディングも実行速度も)
一番速いけど、他人とのデータのやりとりではやっぱり危険だね。
自分はなるべくプレーンなテキストでやりとりするけど・・・みんなはどうしているんだろう？
HDF とかなの？

>>313
例えば、１行の read 文で書き出すデータの総バイト数を表す数字を
データの両側に付けるのは仕様だった筈。ここが C と違う。

integer a
real*8 b

read(unit=9)a,b

とすると、12バイトのデータ（aとb）の両側に余計にデータがつく。

open 文でエンディアンの指定も出来たように記憶しているが、最後
に触ったのが３年くらい前なので忘却の彼方。

＞データの両側に付ける
direct accessのためにそういう仕様になってるんだっけか
ただ、その数字自体のバイト数がコンパイラ依存なんだよな
なんで規格でちゃんと決めなかったんだろ

＞open 文でエンディアンの指定
まじ？聞いたことないけど

>>318
> ＞open 文でエンディアンの指定
> まじ？聞いたことないけど

http://gcc.gnu.org/onlinedocs/gfortran/CONVERT-specifier.html

規格の事は知らんけど、こんなのがあったぞ。
コンパイラ依存だったりして。

>317 データの両側に付ける
ためになった。
経験上たいていの無書式書き出ししたファイルサイズが
ちょびっとだけ大きかったから
何か付帯情報がくっついているんだろうな、くらいには思ってたけど
「両側」に「総バイト数」というのは知らなかったでした。
自分のつくったファイルをみるとたしかにそうなってます！
・・・問題は >>318 さんのとおりコンパイラごとに違うという事か。めんど〜

なんか扱いにくいね＞unformat
結局binaryが一番ってことか

ACCESS="STREAM"

fortran90 を使って、n 行の文字列データが格納されている
test.dat というファイルの1行1行を読み込んで、
標準出力に出力しようと思っています。
そこで、(抜粋ですが)

------------------------------------------
character(len = 50) :: chara

open(1, file='test.dat', status='old')

do i = 1, n
read(1, *) chara
write(*, *) chara
end do
------------------------------------------

と書いてみたのですが、文字列データに/や*などが含まれると
文字化けしてしまいます。また、空白が出力されません。
一体、何が原因と考えられますでしょうか。
なお、OSはubuntuで、コンパイラは、gfortranです。
お分かりの方がいらっしゃいましたらお教えいただけると嬉しいです。
よろしくお願いします。

手元にコンパイラ無いから検証なしw

read(1,'(a)') chara
write(*,'(a)') chara

で試してちょ。たぶんうまくいくとおもうけど。
長さ指定無しのA は便利。

>> 324 さん
早速お答え頂きありがとうございました！
試してみたところ、文字化けせず上手く出力されました！
お教え頂きありがとうございます！！

みんなラクダ記法とアンダーバー記法どっちで書いてる？

現在、30年分の日別のデータがずらっとあります。
その生のデータを初めに、5日ごとの半旬平均したあとに30年平均したい場合、
どう書いていいのか想像つかないのですが。５日ごとの半旬平均する場合、
欠損値があれば、数えずに平均したいんです。
1半旬（1〜5日間の平均）、2半旬（6〜10日間の平均）、3半旬（11〜15日間の平均）、4半旬（16〜20日間の平均）、
5半旬（21〜25日間の平均）、6半旬（26〜30日間の平均）が12カ月分で、まず1年を72半旬にしたいのです。
6半旬は、月によっては、（26〜28日間の平均）や（26〜31日間の平均）の場合もあります。

例えば、
1/1　10
1/2 　8
1/4 　6
1/5 　2
1/6 　1
1/9 　2　
1/10 　9
： 　：


とあった場合、

1月1半旬　　（10+8+6+2)/4 = 6.5
1月2半旬 (1+2+9)/3 = 4
：　　　　　　　　　　：

といったような感じなのですが。わかりにくいかもしれませんよね。わからにくいところが
ありましたら、ご連絡ください。どうかよろしくお願いします。

>>326
正直悩み所だわそれ
ごちゃ混ぜになってるｗ

>>328
そっかー．最近はjavaっぽく変数名はラクダ，サブルーチン名はパスカルにしてるんだよねー

たびたびすみません・・・。
323 ですが、Fortran90 で、文字列や数字がごちゃまぜになっている
ファイルから数字だけを抜き出したいのですが、どうすればよろしい
のでしょうか？
例えば、(1234567890)という行から、1234567890だけを抜き出したいと
思っています。

試しに、(1234567890 が書いてある行を対象に

character(len = 1) :: chara
real(kind = 8) :: x
read(1, '(a)') chara, x

とやってもエラーになり駄目でした。
(chara で左括弧だけを取り除きたかったのですが・・・。)
ご存知の方がいらっしゃいましたらお教えいただけるとありがたいです。
よろしくお願いします。

324

数字だけ見つけるには・・・
ichar(x) (長さ１の文字列をアスキー表の整数に変換する奴)で
ichar(x) 48 から 57 の範囲をとるからそれで見分けはつくので、
一度、
 read(1, '(a)') chara
 do j = 1, 50
  x = chara(j:j)
  idummy = ichar(x)
  if ((idummy .LT. 48) .OR. (idummy .GT. 59)) chara(j:j) = ' '
 enddo

な感じで数字以外を空白に置き換えてから・・・かな？
上の例だと指数部のE (やD)とか小数点(.)や正負(+-)も空白に置き換えるので
必要に応じて変えてちょ。

あと、
write(x,'(i1)',ERR=199) idummy
みたいにI/O エラーで判別するやりかたもあるけど、
これは分岐が汚くからichar使うのが無難かなぁ

>>327
エクセルでやった方が幸せになりそう。

>>327
yy/mm/dd みたいにせめて桁はそろえてほしい・・・
331 的やりかたで年月日の３つの数字を取り出して
・・・先は長いｗ

Intel Fortranでコンパイル時にOpenMPを指定するだけでオーバーフローになってしまうのですが
何が原因でしょうか？

F77のプログラムで、普通にコンパイルすれば問題なく作動します。

コンパイルにしくじるんだよね？
unlimit か ulimit して制限取っぱらっても変わらないだろうなぁ。
というかOS なに？

あとは・・・Private 属性でアホみたいに大きいサイズの配列があると・・・とか、かなぁ。

構造体のメンバを一覧でとってくる方法ってありますか？
type(member)
integer id
integer age
real(8) height
end type member

type(member) person
person%id = 1
person%age = 26
person%height = 171.d0

などと宣言して

do ???????????????
print *, @@@@@@@
enddo

のようにすると
1
26
171.00000
のなるようにメンバ変数でループをまわしたいのです

>>335
/AUTO_ALLOCATEでうまくいきますた。

はじめまして。

T検定を書いているんですが、
disperは、不偏分散で、valueは、分散の推定値です。

value(doy) = ((ndata(doy)-1)*disper(doy)+(ndata2(doy)-1)*disper2(doy))/&
(ndata(doy)+ndata2(doy)-2)
t(doy) = abs(sunlight(doy)-sunlight2(doy))/sqrt(value(doy)*((1/ndata(doy))&
+(1/ndata2(doy))))

これのt(doy)の出力が、+inf,時々NaNになります。
直後に

write(6,*) doy,lon,lat,ndata(doy),ndata2(doy),value(doy),t(doy),&
sunlight(doy),sunlight2(doy)

で出力した際の一部は
1 141.678 45.4147 7 10 3.97499 +Inf 1.025 0.82
2 141.678 45.4147 7 10 3.57839 +Inf 0.42 0.46000004
3 141.678 45.4147 7 10 5.2759204 +Inf 2.48 0.24000001
です。エクセルで計算した場合、Tの値はでると思うんですが、何が原因の可能性として
考えられるのでしょうか？
宣言文は、全てREAL,dimension(72)です。

338

解決しました。すみません。

みんな構造体使ってないのかな

使ってるよ

収束しねーよ・・・。(´；ω；｀)

発散しないのなら、変なattracter に
はまった•••となると
それが系の本質かもしれないから
頑張れ！
ただのタイプミスの可能性も残るけどね

Fortran90/95の構文色分けが一番見やすいエディタって何ですか？

vim

真面目な答え

emacs

>>345-346
ありがとうございます。

ははー。90以上だと、やっぱりそこらへんですか。
色々試してるんですけど、Fortran対応ってなってるものでも
実質77対応のが多いんですよね…

秀丸＋定義ファイルとかもやってみようかな

誰かFortran開発という観点でのエディタ・IDEを評してくれ
メジャーな言語だとそういった話題が豊富なんだが、Fortranは寂しい限り

The first argument (X) to the intrinsic SQRT must be of REAL or COMPLEX type, not INTEGER(KIND=3)

平方根計算しようとしたらこんなのでるの

>>349
SQRT関数の引数は実数か複素数でないといけません。整数はダメです。

はじめまして。

前からなぜうまくいかないのか、わからず、違う強引な方法で行っていたのですが、
今回はそうはいかなく、質問させていただきます。
if文なのですが、今、データのなかの日にちを1〜365としています。
それを1〜59、335〜365を冬。60〜151を春、152〜243を夏、244〜334を秋と、if文で
場合分けしたいのですが、うまくいきません。どう書けばよいのでしょうか？

do doy = 1,365

if( 60 <= doy <= 151) then
ave1 = ave1 + sunlight(iyear,doy)
ndatax_sp = ndatax_sp + ndata(iyear,doy)
else if(152 <= doy <= 243) then
ave5 = ave1 + sunlight(iyear,doy)
ndatax_su = ndatax_su + ndata(iyear,doy)
else if(244 <= doy <= 334) then
ave9 = ave1 + sunlight(iyear,doy)
ndatax_au = ndatax_au + ndata(iyear,doy)
else
ave13 = ave1 + sunlight(iyear,doy)
ndatax_wi = ndatax_wi + ndata(iyear,doy)
endif
end do

>>351
とりあえず教科書か↓のページを開いて IF文と論理式について復習するんだ。
http://www.str.ce.akita-u.ac.jp/~gotou/zyouhou/f90/sec3.html#sec3

Fortran90以降が使えるのなら IF文ではなく、SELECT CASE文で書いた方が
きれいに書けるかも。

if…else if…else ifって何個も分岐するならselect case文のほうがいいに決まってるけど、
分岐が2、3個の場合でもselect case文でやるべきか結構迷う。
2、3個程度ならパッと見はifのほうがすっきり見えるからねぇ。

>>352,353

ありがとうございました。case文で書きました。
もう一度復習してみます。

エディタ色々試し中です。
Win版geditを使ってみました。
Fortran95に対応してますが、
拡張子見て勝手にFortranモードになってくれるわけではなくて
一々手動でモード切り替えしないといけないこと以外は
結構良いですね

知恵を貸してください。

3次メッシュコード、年齢層、要素の値の3列
で構成されるファイルが要素別に10ファイルあって、
こいつらを1つのファイルにまとめようとしています。
(メッシュコード、年齢層、要素1、…、要素10)

10ファイルに出てくるメッシュコードと年齢層の組み合わせは
ファイルごとに異なっていて、重複する場合もあれば1ファイルのみで出現することもあります。
また、あるファイルでは数百行しか格納していないのに対して、
一方では数十万行分のデータを格納しているものもあります。

素人考えで12次元の配列にぶちこんで出力すれば良いかなとも思ったのだが、
3次メッシュコードだけで60万行になってしまうし、
要素1〜10の値がまったく含まれないコードも多々ある訳で。
効率的に要素1〜10の値が1つでも含まれるメッシュコードと年齢層の組み合わせだけを
出力できる良い方法はないでしょうか。

教えてください!
関数 idate に関してですが、

gfortran -c -fdefault-integer-8 -O2 -std=legacy Libtad.f

In file Libtad.f:1339
call idate (month,day,year)
Error: Too many arguments in call to 'idate' at (1)

となり、エラーになってしまいました。
ググったら idate はinteger 3 個でいいと書いてあって
なぜエラーになるのか解かりません・・・。
もともとのソースコード(tinker/Libtad.f)は

integer year,month,day
(略)
call idate (month,day,year)

と書いてあります。
gamess というソフトをインストールしようとして、コンパイルしたら、
何故かそのライブラリだけがError を起こして上手くインストールできないのです・・・。
(コンパイルは入っていたスクリプトで行っています。)

OS:Scientific Linux 5.5
gfortran: gcc4.1.2

>>357
http://gcc.gnu.org/onlinedocs/gfortran/IDATE.html
IDATE - Get current local time subroutine (day/month/year)
IDATE(VALUES) Fills VALUES with the numerical values at the current local time.

>>358
ありがとうございます。
idate の部分上手く解決しました。

教えてください
あるtxtファイルから必要なとこだけを抽出したいんです。例として、

120 30
121 30
121 31
120 31
120 30

123 32
124 32
124 33
123 33
123 32

↓から

120 30
121 30
121 31
120 31
120 30

123 32
123 32

のように、ある範囲（ここではX=123以下かつY=32以下）以外のものをはじきたいのですが、
空行をそのまま残す方法がわかりません。
どなたかご教授願います。

5個ずつのブロックで入ってるのなら，
5個読みこんだら空行を書き出すようにしたら?
たとえば

do i = 1, n
read(10,*) a(i, 1), a(i, 2)
if ( a(i, 1) > 132 .or. a(i, 2) > 32) write(20,*) a(i,1 ), a(i, 2)
if (mod(i, 5) == 0) write(20,*)
enddo

とか．

モジュールを外部ファイル化して
メインファイルとモジュール用のファイルをコンパイルするのって
f90 main.f90 mod.f90
でいいの？

>>361
お答え頂きありがとうございます。

私の書き方が悪かったのですが、１つのブロックの中の数はバラバラです。
条件によっては書き出し後でそのブロックの中の数が変わることもあります。

>>361さんのおっしゃる通り、空行まで読むといったん条件を通して
新たなブロックと空行を書くようにすればできそうな感じはするのですが、
「空行まで読むと」という部分はどのようにすればいいのでしょうか？

ichar() で空行を読み込むと "32" が返されるから，

character(1) :: a(100)

open(10, file='hoge.txt')
do i = 1, foo
read(10,'(a)') a(i)
enddo
close(10)

でもう一回同じファイルを開いて do loop を回して範囲外のものをはじきつつ
if( ichar(a(i) == 32) で空行を書き込むとかどうかな?

違ってたらごめん． 正直Fortran で完結させないでシェルスクリプトなり
使った方がラク．

>>360
こんな感じで良いの？

program extraction
implicit none
integer::i,j,k,m,is,large_number=10000
character(len=32) c
open(10,file='hoge.dat')
do i=1,large_number
read(10,'(a)',iostat=is) c
if(is>0) stop "error"
if(is<0) exit
m=index(c,' ')
if(m/=1) then
read(c(1:m-1),*) j
read(c(m+1:),*) k
if(j<=123.and.k<=32) then
write(*,*) j,k
end if
else
write(*,*)
end if
end do
close(10)
end program extraction

>>362
モジュール先に作るべきなんじゃなかったっけ？
f90 -c mod.f90
f90 -c main.f90
f90 mod.o main.o
となるのが無難な気がする．

program hoge
implicit none
integer :: x, y, ios
character(len=100) :: str
open(unit=11, file="input.txt", status="OLD")
open(unit=12, file="output.txt", status="UNKNOWN")
do
read(11, "(A)", iostat=ios) str
if (ios /= 0) exit
read(str, *, iostat=ios) x, y
if ((ios /= 0) .or. ((x <= 123) .and. (y <= 32))) then
write(12, "(A)") trim(str)
end if
end do
close(11)
close(12)
end program hoge

index使わなくてよいのか。
勉強になった。
ありがとん。

>>364>>365>>367
皆様ありがとうございます。

残念ながら、本日はFORTRANを使える環境におらず、
月曜日にもう一度皆様からご教授いただきましたように
プログラムを見直し、計算してみようと思います。

また月曜日にご報告させていただきます。
本当にありがとうございます。それではよい週末を〜

結果をご報告いたします

>>367様のプログラムで計算し、出力ファイルを見ると
全データの中から、すべての空行が出力されておりました。
数値は入っておりませんでした。

>>365様のプログラムでは、
同じように空行のみが出力されておりましたが、数がかなり少なかったので、
範囲内(j<=123.and.k<=32)にある空行が出力されたのではないかと思います。

もう一度じっくり考えてみたいと思いますが、無理そうなら諦めます；；
お二方ともありがとうございました。

>>366
それでできなかった
モジュールのほうはコンパイルできるんだけどメインがモジュール読み込んでなくてコンパイルできない
モジュールのほうをコンパイルしたときにoファイルだけ生成されてmodファイルが生成されないことと何か関係あるのかな？

＞モジュールのほうをコンパイルしたときにoファイルだけ生成されてmodファイルが生成されない
これはおかしいね
modファイルを作らないコンパイルオプションになっているとか？（そんなオプションあるのか知らないけどｗ

Fortran が複雑になって行っている気しかしない。
77 の頃に戻ってくれんかな。まぁ幾らなんでも
77 ではヒド過ぎるから、90 程度の拡張で。

一体何処に向って行っているのだろう。学習のコ
ストが上がるのは良い事ではないと思うのだが。

>>372
コンパイルソフトが何か分からないんだけど
コンパイルソフト何使ってるか調べるコマンドってある？

学習コスト＜学習による恩恵　であれば肯定されるべきだと思うが
言語策定者は「こう使えばシアワセになれるんよ！」ってところを広くアピールしてほしいね
まぁそういうことするのは物書きの仕事かもしれないけど
でもFortranは少ないからなぁ

test

>>373
そうは言っても、もはや世界がFortranだけで閉じていないから、インターフェース周りを強化しないと
ライブラリのラッパーすらかけない。

伝統的数値計算ライブラリならともかく、MPIとかOSのAPとかIは同じ引数が色々な型を取ったりするから
そういうのに対応しようと思うと、多態型とかOOPっぽい機能が必要になるし、Cの呼び出し規格との
共通化も必要になる。

数値FormatがIEEE745に統一された以上、数値計算主役のFortranでは、その代数系を操作する
命令も欲しくなる。組み込み関数ももっと色々欲しくなる。

CPUも並列が当たり前になりつつある今、それに向いた命令も欲しい。


とはいえ他の言語のように、F９０を削ったコア部分と標準Moduleみたいな形に分割したほうがいい気はするけど。
過去互換性をかなり保ったまま肥大化したから、命令の直交性なんかが悪くて、新規学習者に負担が
重すぎる気がする。

>>377
学習負担そんなに重いかな？
Fortran使う分野は特殊なんだし、ユーザもそれなりのレベルを前提としていいと思うけど

あと、
＞MPIとかOSのAPとかI
は
＞MPIとかOSのAPIとか
だよね？

>>378
>＞MPIとかOSのAPIとか
>だよね？
んだ。すづれい。

>学習負担そんなに重いかな？
90/95を77を横目に見つつ66以前を無視する形で学ぶ分にはいいかもしれんが、
これに2003が加わるときついような。やや95と挙動が変わっている。
CLASSの作り方とかも不細工だしｗｗ

計算時間長いプログラムをデーモンとして走らせたいんですけどどうやったらいいんでしょうか．
コンパイラはpgiです．

>>380
まあこの場合、コンパイラはPGIでもIntelでも何でもいいんだが
http://www.codereading.com/notebook/moin.cgi/IgnoreTheHangupSignal

Fortranで用意されている関数やサブルーチンを詳しく網羅しているサイトってないですか？
やはり本家を見るべきですか？（その本家がどこかも知りませんが・・・）

サブルーチンは自分で組み込むもんじゃないのか

内蔵してるサブルーチンでしょ
system_clockとか

intel fortranは拡張機能で、内部手続きを別の手続き
に渡すことを許している。他のコンパイラでは、もし
2003に対応していたら、手続きポインタでさして、わ
たせば、一応機能しているように見える。が、偶然の
可能性もある。2003の手続きポインタがさせる先は仕様
でどうなっているかを教えていただけないだろうか。　　

Fortranである変数や配列のアドレスを知る演算子はありますでしょうか？

非標準だが、LOCまたは%LOC関数が使えることがある。

>>387
ありがとうございます
非常に助かりました！！
ちゅっちゅっ

>>372
富士通だとデフォルトでは生成されない
frt -c -Am -M./ (なんとか.f90)
でいけるのではなかろうか。いずれにせよ、
frt -help | less
とかやって、"module"で検索を掛けるのが吉。

>>389
そのレスを見る前に昔の先輩が書いてたプログラムを見てたらmakeファイルに使用コンパイラがfrtってなってたから
ググってもしやと思ったらなんというタイミング・・・
多分そのとおりだと思う
月曜に試してみてまた報告しにくるわ

1990 4 2 401 45.2448 141.41 稚内 65 27 0 0 0
1990 4 2 402 44.5618 142.3524 北見枝幸 67 18 0 0 0
1990 4 2 404 44.2142 141.4218 羽幌 55 41 0 0 0
1990 4 2 405 44.3442 142.5806 雄武 95 36 0 0 0
1990 4 2 407 43.4612 142.2224 旭川 55 24 0 0 0
1990 4 2 409 44.01 144.17 網走 25 66 0 0 0
1990 4 2 412 43.0324 141.1954 札幌 53 36 0 0 0
1990 4 2 413 43.1236 141.4718 岩見沢 83 29 0 0 0
1990 4 2 418 42.583 144.233 釧路 83 31 0 0 0
1990 4 2 426 42.093 142.4648 浦河 55 60 0 0 0
1990 4 2 428 41.5154 140.0742 江差 53 65 0 0 0
1990 4 2 433 42.5354 140.4542 倶知安 67 34 0 0 0
1990 4 2 435 44.2036 143.2136 紋別 57 49 0 0 0

という形式のデータを読んだ場合(1961~2008年あります）、

do i = 1,137000
read(50,*,iostat=io) year,mon,day,ID,lat,lon,name,aveavecloud,sunsunlight,&
prepre,maxmaxhourpre,maxmax10minpre
if ( io < 0 ) exit
pre(ID,year,mon,day) = prepre / 10.0
sunlight(ID,year,mon,day) = sunsunlight/10.0
avecloud(ID,year,mon,day) = aveavecloud
maxhourpre(ID,year,mon,day) = maxmaxhourpre/10.0
max10minpre(ID,year,mon,day) = maxmax10minpre/10.0
ndata(ID,year,mon,day) = ndata(ID,year,mon,day) + 1
write(6,*) year,mon,day,ID,ndata(ID,year,mon,day)
end do
とし出力させると、ndata(ID,year,mon,day)がIDが402と409だけが２となります。
この原因はなんでしょうか？

すみません。
391での説明で追加です。

ndata(ID,year,mon,day)がIDが402と409だけが２となります。と書きましたが
毎回ではないと思います。402と409ではないときもあるかもしれません。この下に２となると止める
プログラムを書いているので、1990年の4月にとまり確認できたはじめのケースです。
以下どうなっているのかはまだ確認中です。

１つのソースファイル内で固定形式と自由形式を混ぜるとコンパイルがエラーになりますが、
「この行からこの行までは自由形式で」みたいな指定をしてコンパイルさせる方法が
あったりしませんかね？

大学の授業で「エラトステネスのふるい」で３桁の素数を求める問題が出たのですがその解答例で

index=1
do n=2,int(max**0.5)
if(index(n)==0)cycle
index(n**2:max:n)=0
end do

↑のようなものだったのですが（回答の一部抜粋です）indexとcycleの役割がいまいち
わからないので教えてください。maxには書いてませんが上で=999と定義しています。


あと別の質問ですが大学では初めに「g77 -o 〜　〜」でコンパイルするように教わり
別の教授からは「intel -o 〜　〜」と教わったのですが具体的に何が違うのですか？


質問長くなってすいません。回答お願いします

>>390だけど
今日試したら無事コンパイルできて実行できた
プログラム書くのがだいぶ楽になったわ　レスくれた人たちありがとう

>>394
> indexとcycleの役割

index は（おそらく整数の）配列だろう。
「ふるい」にかけた後 index(n) == 1となる n が素数。

cycleはDOループの中で次の行からend doまでをスキップする。

>393
ないんじゃない？
ファイル分けてコンパイルするのが一般的じゃないかしら．

>394
エラトステネスのふるい自体はwikipediaにあるので説明なしで．
index=1 ! 配列index(999)に整数値１を代入，最初は全部素数だとしておく．
do n=2,int(max**0.5)
if(index(n)==0)cycle ! index(n)=0のとき（すなわち素数でない場合）enddoまでスキップ．
!この場合はこの１行↓をスキップするだけ
index(n**2:max:n)=0 !n**2からmaxまでnおきにindexの値を0とする．
! do i=n**2,max,n
! index(i)=0
! enddo !と同義
end do
このプログラムは，既に組み込み関数として用意されているindexとmaxを整数型として利用しているので，あんまり良くない回答例だと思います．

g77とintelの違いはコンパイラの違い．

indexやmaxどころかINTEGERやIF、THENまで変数等に使えるのがFORTRANの良いところじゃないかwww

>>396,>>397
ありがとうございます
indexってネットで検索したら関数として乗っていたのでてっきりそっちかと思ったら
ただの配列だったんですかｗ
おかげで謎が解けました。すっきり

interfaceブロックの中には、型定義だけでなく、implicit noneも
及ばないのな。なんで、こんな仕様にしたんだろう。
Fotran2003からimportで取りこめるようになったが、2003完全サポート
なコンパイラは少ない。

完全サポートって完全準拠って意味？
もうすぐ2011になろうってのに

コンパイルすると「do文にパックド命令が出てきます」って表示されるんだけどどういう意味？
コンパイラは富士通で一応ちゃんと実行できるんだけどなんか怖い

infoseekのweb終了で「フォートランナーのためのFortranあれこれ」っていうサイトが見れなくなってしまいました。
かなり好きなサイトだったのですが・・・ここ見てらしたら是非復活お願いします＞作者さん

.not. と .and. と .or. の優先順位ってどうでしたっけ？

優先順位： （高） .NOT., .AND., .OR. （低）

>>405
ｔｈｘ！

自分は怖いから
くどい位に括弧付けてる

>>407
オレも。

Cの i++ みたいなインクリメント演算って i=i+1 とするより早い（専用の処理をするから）と聞いたのですが
Fortranではそういうのないんですか？

>>409
そりゃ相当おバカなコンパイラだけ。
普通は最適化で同じコードになる。

>>410
そうですか(^^;)
今のご時世気にする必要ないんですね

>>409
何もしなくてもコンパイラが勝手に宜しくしてくれるとかないの？
構文解析だけの問題のような気がするけど。

意味的には、
(i++) の式の値はインクリメントする前の i の値で、
(i=i+1) の式の値は代入された加算後の i の値なので、
意味が異なるから構文解析だけの問題ではない。

でも ++ や += が使えたら便利だとは思う

>>414
ようこそC++の世界へ

まあ、
i = i + 1 くらいだったら、表記法を増やさないほうをとるが、
hoge(3,4)%piyo(6)%i = hoge(3,4)%piyo(6)%i + 1
となってくると、C/C++がうらやましくなったりはする。

C++って配列演算ってあるっけ？

std::valarrayとかか？

配列演算自体はできるみたいね。
速度的にはどうなんでしょ。ググってたらこんなん見つけましたが

http://ir9.jp/prog/ayu/datlog/tech_cpp/1038031395/1038031395_03.html#R560

Fortranを使い続けることの意義を見出したいと思う今日この頃・・・
徒労かなｗ

>>419
たとえ速くとも、こんなん書きたくないよ、
というのが、正直なところ。ライブラリで、
プログラマビリティを上げるのは限界が、
あって、豊かな配列操作を組み込みでもっ
ているfortranにはかなわないでしょう。

C/C++のポインタや配列はエイリアスの問題を持っていて最適化の妨げになっているので
restrictポインタを使わないとFortranにはかなわないよ

Windows版のifort (Intel Visual Fortran)でプロファイルを取る方法を教えてください。

検索すると、
ifort -pg hoge.f90
./a.out
gprof ./a.out gmon.out
というのがよく掲載されているのですが、Windows版には、下記のように
-pgの対応物がないようです。

map_opts -tw -lf -opts -pg

Intel(R) Compiler option mapping tool

mapping Linux options to Windows for Fortran

'-pg' Linux option maps to
--> no equivalent Windows option found

あるデータをバイナリ化にするプログラムについての質問です。
あるデータは、2列の33行のデータです。1列目は年、2列目は実数が入っています。
----------------------------------------------------------------------
implicit none

INTEGER :: i,iyear
real*4 :: val

open( 10, file='a.dat', status='old')
open( 90, file='b.bin', status='unknown', access='direct', &
recl=1)

do i=1,33
read( 10, * ) iyear, val
write( 90, rec= i ) val
enddo

stop
end
------------------------------------------------------------------------
これはコンパイルはできるのですが、実行できませ。
At line 12 of file bainari.f90 (Unit 90 "b.bin")
Traceback: not available, compile with -ftrace=frame or -ftrace=full
Fortran runtime error: Writing more data than the record size (RECL)
というエラーがでます。
なにか考えられる原因はありますでしょうか？？
どうかよろしくお願いします。

仕様でどうなっているかは忘れた。
正格でないかもしれないが、
いっかいのwrite文でかきだせる長さ
はreclで指定した長さ以下でなければ
ならない。

ここで問題となるのはreclの単位。
intel機の場合、1バイト、いっぽう書こう
としてるvalは4バイトなので、ランタイム
エラーが出る。Recl=4としないといけない。

とはいえ、元のコードが間違つているわけでもない。
SUNなどの環境ではreclの単位は4バイトである。
（コンパイルオプションで、1バイトにできるが＾-＾）

ありがとうございます。
ちょっとわからないんですが、僕の使っているパソコンはレッツノートでＣＰＵは
Intel centrino2と書いています。cygwinを使っています。
直し方ありますか？

>424,425

すみません。ＲＥＣＬ＝4となおしたらできました。
ありがとうございます。描画して合ってるか確かめてみます。

初めて訪問させていただきます。
intel fortran にて
error #7836
実引数がスカラーの場合、実引数が形状引継ぎ配列やポインター配列でない配列の要素、またはそれらの要素の部分文字列でなければ、対応する仮引数はスカラーになります。
これはどのようなエラーでしょうか？
debugモードだと発生し、releaseモードだと発生しません。

宜しく、宜しくお願いします。

構文チェックを厳しくするオプションがDebugのときは付くのだと思うのだが。
該当箇所(呼び出し元と呼ばれる手続き)を見せれる?

渡ししている変数の型と、読んでいる手続きの引数の型が
わかればいいです。

>>429
解決しました。
問題の箇所は、変数Aが
　親ルーチンで配列
　子ルーチンは実働が無いので変数
　孫ルーチンで配列
にて宣言されていたためです。
リリースモードではスルーされていました。
ありがとうございました。

>>414
それはアセンブラのINC命令とかADD命令に一対一対応していた頃は意味があったかもしれないが、
いまさらただの表記の慣れの問題にすぎんだろ。

>>416
F2003でASSOCIATE命令が出来たから、がまんできるべ？

>>420
F95以降ELEMENTAL　FUNCTIONをユーザー定義できるようになったから、スカラー用の関数を書けば
多次元配列を引数に出来るので便利。
これはLISPとかのMAPと同じようなもんだが、マクロと実行命令ではスピードが違うから
配列操作が多いならFortran以外ないだろう。

話にあんまついて行けてない・・・
多分レベル差がありすぎるんだと思う
誰かF2003あたりの本書いてください、オナシャス！

２階常微分方程式の問いを考えていたのですがどう手をつけたらいいかわかりません
問題分は
ttp://iup.2ch-library.com/i/i0207521-1292859670.jpg
です。
できれば解説付きでお願いします

せめて画像回転ぐらいしようや
読ませる気が無いとしか思えん

>>433です
貼る画像間違ってました
http://iup.2ch-library.com/i/i0207930-1292946928.jpg
お願いします

ルンゲ・クッタ法を教えて下さい、なら Fortran の話題じゃないんじゃね？
それこそ、ググ(ry

MODULE Integrate
IMPLICIT NONE
CONTAINS
! ルンゲ・クッタ４次積分子
! y : (x1,...,xn=t) at t=t0+dt
! x : (x1,...,xn=t) at t=t0
! f : (x1,...,xn) -> (f1,...,fn)
! (形式)
! f(f_vec, x_vec)
! real(8), dimension(n) :: f_vec, x_vec
SUBROUTINE Runge_Kutta4(y,x,f,dt)
real(8), dimension(:) :: y,x
real(8), dimension(size(x)) :: xx, k1,k2,k3,k4
real(8) :: dt
interface
subroutine f(deriva,x)
real(8), dimension(:) :: deriva, x
end subroutine
end interface
call f(k1,x)
xx=x+0.5d0*dt*k1
call f(k2,xx)
xx=x+0.5d0*dt*k2
call f(k3,xx)
xx=x+dt*k3
call f(k4,xx)
y=x+dt/6.0d0*(k1+2.0d0*(k2+k3)+k4)
END SUBROUTINE
END MODULE

そういやルンゲクッタの証明ってどうやんだっけ・・・
たいていのテキストには式しか書いてないよな

dx(t)/dt = f(x)
の関係を使って、xの高次の時間微分をfのxによる微分に置き換える
(いわゆるチェーンルールを使って)
係数を知ったうえで、4次までの展開になっていることを確認したことが
あるが、それでもかなり面倒な計算になる。

>>437さんありがとうございます
さっそく明日やってみたいと思います

>>437さんのがうまくできなかったんですが
どういう手順で考えていくかだけでも教えてもらえませんか？

しばらく見てなかったので申し訳ない。使用例を書きます。
例として、調和振動子を考える。これは、位置x, 速度vについての微分方程式
　dx/dt = v (1)
　dv/dt = -x (2)
である。形式的に自明な式
　dt/dt = 1 (3)
を追加して、ベクトル z = (x, v, t)を定義すると、
(1)〜(3)はあるベクトル関数fによって
　dz/dt = f(z)
とかける。このfをFortranで具体的に書くと以下のようになる。
module ho
subroutine f(dz_dt,z)
real(8), dimension(:) :: v, x
dz_dt(1) = z(2)
dz_dt(2) = -z(1)
dz_dt(3) = 1D0
end subroutine
end interface
program test
use ho
real(8) :: z(3), dt
z(1) = 1d0 ! x(0)
z(2) = 0d0 ! v(0)
z(3) = 0d0 ! t = 0
dt = 0.01D0
do while (x(3) .lt. 6.28)
call Runge_Kutta4(x,x,f,dt)
write(*,*) z
end do
end subroutine

三重積分の数値計算プログラムについて相談
台形公式を用いたプログラムで分割数nを大きくしていけば精度は上がると思ってるんだけど
n=100とn=1000とでは値が何故か全然違うんだ
被積分関数f(x,y,z)=x^2+y^2+z^2で積分区間をx,y,z共通で[0,1]としたとき値Sは1.0になるわけだけど
n=100のときS=1.000005054で、n=1000のときS=6.25e-2になってしまう
数学板とで迷ったけどアルゴリズム的なことだったらこっちかなと思って聞いてみる
誰か分かる人いたら教えてください

メインプログラムはこんな感じ

n=1000（あるいは=100）

xs=0.0
xe=1.0
ys=0.0
ye=1.0
zs=0.0
ze=1.0
dx=(xe-xs)/float(n)
dy=(ye-ys)/float(n)
dz=(ze-zs)/float(n)
s=0.0

do iz=1,n
z1=zs+dz*float(iz-1)
z2=zs+dz*float(iz)
do iy=1,n
y1=ys+dy*float(iy-1)
y2=ys+dy*float(iy)
do ix=1,n
x1=xs+dx*float(ix-1)
x2=xs+dx*float(ix)

h1=f(x1,y1,z1)
h2=f(x2,y1,z1)
h3=f(x1,y2,z1)
h4=f(x2,y2,z1)

h5=f(x1,y1,z2)
h6=f(x2,y1,z2)
h7=f(x1,y2,z2)
h8=f(x2,y2,z2)

h=(h1+h2+h3+h4+h5+h6+h7+h8)/8.0
s=s+dx*dy*dz*h

end do
end do
end do

とまあこんな感じ
今はf(x,y,z)=x^2+y^2+z^2で計算してるけど最終的にはもうちょっと複雑な関数になる

「情報落ち」だな

>>446
？

ググれ
このごろじゃ計算誤差について何も教えてないのか

最後に合計を見るだけじゃなく
実際に何を足してるのか出力してみろ
何が起きてるのか分かる

>>448
把握した
今このプログラムは実数が単精度になってるから倍精度にして
h1〜h8まで出力してみて確認してみる

あぁああ〜
乱流モデルの差し替え、うまくいかねーーーーーーOTL
なぜあぼーとするんだ・・・
だれか助けて・・・

じゃあDNSで

>>443だけど単精度から倍精度に変えて情報落ちしてるであろうSの計算方法を変えたら問題解決できますた

>>448
昔のほうがしっかりしたカリキュラムだったって話を聞くが本当なのか
ゆとり教育ってやつか

>>453
昔は電子計算機という前提のカリキュラムだったが、最近は学生全員にノートパソコン買わせてるのに
プログラミングも教えないからな・・・

なにそれひどい

>>442さんありがとうございます
おかげである程度理解できました！

FORTRANの課題が出て〆切が月末なんだ。
三問ほど有償で解いてやんよって人いないか。

とりあえず問題晒してみそ
誰か解いてくれるかも

>>435なんですが
１：０〜２０の区間を適当に分割
２：初期値入力
３：２階常微分方程式を一階常微分に直す
　　dy/dx=z dz/dx=f(x,y,z)=-sin(y)としました。
４：ルンゲクッタ法
　　y=y+(k1+2*k2+2*k3+k4)/6
k1=h*f(x,y)
k2=h*f(x+h/2,y+k1/2)
k3=h*f(x+h/2,y+k2/2)
k4=h*f(x+h,y+k3)
　　これの()内が３つあるバージョン。分割した回数繰り返す。
５：ファイルにデータ書き込み出力

こういう流れで作りたいんですがどこが間違ってるのか出力すると初期値を
繰り返したデータになってしまいます。２〜５あたりだけでもいいので
わかる方いたら教えてもらえませんか？
イメージとしては↓のサイトのようにやったんですがうまくいきませんでした
http://www.akita-nct.ac.jp/~yamamoto/lecture/2005/5E/differential_equations/text/html/node3.html

お言葉に甘えて失礼します。
とりあえず一問目。

「ある動物園の一日の総入園者数（1行目）と入園者の年齢（２行目以降）が
収納されているファイルからデータを配列に貯えて、平均年齢、及び最高年齢と最低年齢を
求めて出力するとともに、各年代（10歳未満、10代、20代、・・・70代、80歳以上）ごとの
入園者数を出力するプログラムを作れ。
また、平均年齢は小数点以下第１位まで示すこと。」

データはここ。
http://www5.plala.or.jp/selene-kw/

宿題代行スレでないのは百も承知だが、ここ以外にFORTRAN師探せなかったんだ。
なにかしら反応もらえると助かる。
よろしくおねがいします。

>>459-460
やはり FORTRAN77 か？　９０でもいいのか？

９０でも７７でもいいです
といってもあまり違いが分かってるわけじゃないのですが

>>460
ここは宿題スレだからスレの趣旨には合ってる

>>460
予想してたより簡単で拍子抜けしたｗ
暇ができた時にまだ誰も解いてなかったら答えるよ

でけた
年代別のところは途中で端折ったから自分で補え

program main

implicit none

integer :: i, n, decade(9)
integer,allocatable :: age(:)

open (99,file="ex5-3.d", status='old', action='read', position='rewind' )
read (99,*) n
allocate ( age(n) )
read (99,*) ( age(i) , i=1, n )
close (99)

decade = 0
do i=1,n
if (age(i)<10) decade(1) = decade(1) + 1
if (age(i)>11 .AND. age(i)<20) decade(2) = decade(2) + 1
end do

write (*,'(a8, f5.1)') " ave = ", REAL(sum(age)) / REAL(n)
write (*,'(a8, i5)') " min = ", minval(age)
write (*,'(a8, i5)') " max = ", maxval(age)

write (*,'(a8,i5)') " 0-10 : ", decade(1)
write (*,'(a8,i5)') "11-20 : ", decade(2)

end program main

ありゃ、タブは2chじゃ無視されるんだな
適宜字下げしてくれ

あと、連続したスペースも無視されるようだから a8 の中身も適宜修正してくれ
intelコンパイラ、f90形式でうまくいくのは確認した

スペースを＆＃１６０；←これを半角で入れるとちゃんと字下げされるんだけど
コピペした時に削らないとコンパイルできないので置換する必要がある

早っΣ（；゜д゜）　親切な回答ありがとうございます！！
Ｃ言語の代行が企業だと50万かかると検索して出たから
有償とかそれなり覚悟してたんだ。

3問中の1問は自力で解くことにします。
もう1問も類似したプログラム見つけるなりしてみる。
どうしても無理げだったら聞きにこさせてもらうかも。

本当にありがとうございます！

>>435を
>>459みたく解ける方いませんか？
もうここに頼るほか解く手段が見つかりません

>>470
435の画像が見られないので何ともいえないけど、442では無理なのか？

>>470
今どんな状態か報告すれば親切な人がアドバイスしてくれるかもしれないよ

>>435が消えてたので再度貼ってみました

http://iup.2ch-library.com/i/i0224579-1295190437.jpg

>>471
437,442のやり方でやってみたんですが自分は７７か９０しかわからないせい
か、単純に勉強不足が原因でわからなくてコンパイルとかできませんでした

私事だけどこの課題が解けないと講義の単位が取れないかもしれないので
かなり焦っています。
親切な方ぜひよろしくお願いします。

構造体で配列を定義するとき、配列サイズを以下のように定数で指定すると
「構造体ではattribute指定できない」という旨のコンパイルエラーになりました。

type test
integer, parameter :: n = 10
integer :: a(n)
end type
type(test) :: test1
end

parameterを使わず、integer :: a(10) とすればもちろんコンパイルできますが・・・。
やり方・考え方がおかしいんでしょうか？ご助言お願いします。

>>473です
今日TAの人に助言をもらい完成させることができました。
考えてくれた方がもしいたらありがとうございます

(x)=exp(a*x) をm次まで　x=0の近傍でTaylor展開した値　g(m)、
f(x)との間の誤差を評価するプログラムをfortranで作成せよ。


またプログラムを実行し、ある１種類のaの値に対し、（xの値も2種類程度でよい）

a x exp(a*x) m g(m)


の順でデータを並べてファイル　ftaylor.dat に出力せよ。
また　f(x)-g(m)をmの関数として　gnuplotで表示し、誤差がどのように小さくなってい くかを
確認せよ。(特定の関数で |f(x)-g(m)|がフィットできるかどうか見てみてください
これをどのようにしていけばいいかおしえてください.

f(x)=exp(a*x)をx=0でテイラー展開（マクローリン展開）すると

g(m) = sum_(n=0)^m [ { f'(n)(0) / n! } x^n ]
　　　 = 1 + a*x + a^2*x^2 / 2! + a^3*x^3 / 3! + … + a^m*x^m / m!

になる．実際はg(x,m)という関数なりサブルーチンを作って，
f(x)-g(x,m)が，mを増やすごとに0に近付くことを確認すればいい．

>>474
定数を構造体の前に置く

integer, parameter :: n = 10

・・・

type test
integer :: a(n)
end type

・・・

>>478
ありがとうございました。

>>476
とりあえず作った
ファイル出力は自力でよろしく
使い方もソース読んで自分の頭をひねりましょう

http://madexpress.ddo.jp/kudafort/5-476.f90

構造体って数値計算だと遅くなるから使わないって聞きますが本当ですか？
クラスも構造体も使えないとなると相当生産性悪いような気がしますが

>>481
構造体使うとメモリアクセスが非連続になるからね
この辺の事情はCでも同じじゃないの？

反復計算は配列を何度も読み出すから
これが構造体になるとキャッシュヒット率はガタ落ちになる

なるほど
でも物理系の数値計算って最終的に Ax=b を解く問題に帰着するケースが多いし
計算時間の大部分はこれが占めますよね
そうするとA, bさえ連続な配列にしておけば、残りのデータは構造体でもいい気がします

＞物理系の数値計算って最終的に Ax=b を解く問題に帰着するケースが多いし
あ、そうとも限らないか^^;　訂正します
まあ何事も用途次第ということですかね

ttp://www.mm2d.net/c/c-15.shtml
Cでも同じような感じだね
ただしメモリ配置を明示できる点で若干のアドバンテージがあると解釈した

>>483
例えば，問題，各パラメータ，求解時間を格納する程度なら
構造体使っても差は出ないと思うけど，
行列やベクトルはメモリ上で連続になる配列を使うのがセオリーだね
CG法なら探査ベクトルや残差ベクトルもみんな1次元配列
行列も非零要素だけ抜き出した1次元配列とポインタを格納する配列があればいい

物理計算に構造体使ってるプログラムは見たことが無いｗ

>>484
うんにゃ，その認識で合ってる

ほとんどの物理現象は偏微分方程式で書き表わせるし，
偏微分方程式を離散化すれば線形方程式系になるよ

ここの人達は行列ソルバは自作したりします？
自分はLAPACKに頼りっきり…
ライブラリがないと何もできないような人間になっちゃいそうで怖いけどｗ

>>487
私は LAPACK 使ってます。自作できるほどのスキルが無いので。
それでも数値積分は泣く泣く自作した覚えがある。

>>488
自分もガウスルジャンドル積分くらいですかね、自作したのは
特殊関数とか自分で書くの無理だよなぁ…あれは

ソルバの研究中です
まだ駆け出しですが

＞＞４８０
解答ありがとうございます。

何回もすみません。


b = b * dble(i)
g = g + (a**i * x**i / b)
write(*,'(2f10.3,e11.3,i6,e11.3)') a, x, f, i, g
ここの意味がよく分からないので説明おねがいします。

あとgnuplotでのf(x)ーg(m)の表し方が分からないのでおしえてください。

>>492
ループ内1行目は>>477における m! の計算，
2行目は g(m) = g(m-1) + (a^m * x^m / m!) を計算してる
マクローリン展開を理解すれば自ずと分かるはず
3行目は問題どおり a x exp(a*x) m g(m) の形で標準出力してるだけ

gnuplotぐらい自分で何とか汁
ループ内に m と f-g をファイル出力するコードを追加して
それをgnuplotで読み込んで描画すればいい

この問題はmが増えるにつれてf(x)=exp(a*x)とg(m)の値が近づくことを確認しろって意図だろう
だからmが増えた時に f-g が0に収束していく様子がグラフに描ければいい

＞＞４９４
ありがとうございます

>>476ではmを手で与えて、データファイルを作成した。
今度は、誤差の絶対値eps（実数）を手で与え、そのeps以内に誤差が収まるためには
mがいくら以上であることが必要かを求めるプログラムを作成せよ。
この場合はどうすればよいですか?

ループ内で誤差を計算し，
ループを do while に書き換えるか if (〜) exit でループを抜ければいい
そのときにループが何回回ったかを数えろ

つーか質問の丸投げはよそうや
「〜って考え方で合てるか？」とか
「こう書いたが計算が合わない」みたいな質問の仕方にしようぜ

自分で考えた痕跡が無いと答える気が失せてくる

すいません
わかりました
ありがとうございます。

sinとdsinみたいな総称名と個別名はどちらを使うべきですか？
総称名はバグの元なので危険という意見もあれば、逆に見通しが良くなってバグも減るという考えもあると思います。

ただ、型変換のreal, float, dbleあたりは常に注意が必要だとは思いますが。

>>499
総称名は引数の型を問わないけど
個別名は引数の型が限定される

pi が単精度型だと
sin(pi) はOKでも dsin(pi) はエラーを起こす
ただし前者で得られる答えは単精度だから、
倍精度の答えが必要なら dsin(dble(pi)) とする必要がある

仮に pi が倍精度なら sin(pi) も dsin(pi) も同じ

この辺はコンパイラが勝手に辻褄合わせてくれるかもしれんけど

>>500
ありがとうございます。確かに Intel Fortran 11ではエラーになりました。

よく考えたら、piが単精度な時点で dsin(dble(pi)) とやっても得られる精度は sin(pi) と
全く変わらないわけですから、結局のところ個別名を使う理由なんてなさそうですね。

つまり、型の間違いによるバグは個別名を用いたところで解決できるものではなく、
型の間違い自体をしないよう代入・演算を気を付けるしかない、ということですね。

>>501
＞dsin(dble(pi)) とやっても得られる精度は sin(pi) と全く変わらない
まぁその通りだけど、sin(pi)を他の倍精度型の数値とかけ合わせたりする時が怖いね

型を明示する時以外は総称名で事足りるけど、
計算バグ除けには個別名が便利なこともあるって認識だな、俺は

>>502
確かに、いつも型を個別名にしておけば、書き間違えたときにコンパイラが教えてくれるというのは利点ですね。
やはり、少なくとも単精度・倍精度が入り交じるときは個別名を使ったほうが良さげですね。

内部副プログラムの関数（メインプログラムのcontains以下に書いた関数）の仮引数に
optional属性を付けると変なところで関数を抜けてしまうようで、なんだか挙動がおかしいです。

手元の本（数値計算のためのFortran90/95プログラミング入門、牛島）を見ますと、
「モジュール副プログラムではそのまま使えるが、外部副プログラムの場合はインターフェースが必要」
という記述があったのですが、もしかして内部副プログラムではoptional属性は使えないんですか？
というか、内部副プログラムの関数とモジュール副プログラムの関数は機能的に同じだと思っていたのですが、違うんですか？

>>504
>内部副プログラムではoptional属性は使えないんですか？

問題なく使える。

program main
implicit none

print *, f1(1.0)
print *, f1(3.0, c = 0.5)
print *, f1(1.0, 2.0, 3.0)

contains
real function f1(a, b, c)
real, intent(in) :: a
real, intent(in), optional :: b, c
! f1 = a + b * c
f1 = a
if (present(b)) f1 = f1 + b
if (present(c)) f1 = f1 * c
end function f1
end program main

optionalな引数をチェックなしで使ったり
誤って主プログラムで宣言した変数を内部副プログラム内で書き換えたり
してるんじゃないかな？

>>504
http://www.geocities.jp/eyeofeconomyandhealth/homepage/fuku7.html

ここら辺も参考にどうぞ

>>478
便乗質問ですが、定数って構造体の中では作れないんですか？
それってちょっと不便じゃないですかね？

>>504 さんの実例体験と書籍からすると、
使えるほうが例外なのかもね。
安全な回避策は今のメインを大きな副プログラムにして
それを呼ぶだけのメインを新たに作るとか、かなあ。

質問です。
読み込むデータが
１　0
２　4
３　5
：　：
24　0
25　6
：　：
48　3
：　：

というようになっており、1万行近く続きます。これを1列目が24の倍数になったらそれまで読んだ行を
すべて出力するようなプログラムを書きたいのですが、どう書けばよいのかわかりません。
1~24が1日分、25~48が1日分となっているので、分割して別ファイルに出力したいんです。
それぞれ別ファイルへというプログラムは書けるのですが、僕が思い浮かぶのは一列目に24の倍数がきたら、出力する
というプログラムしか思い浮かばず、24の倍数がきたら今まで読んだ行すべてを出力、次にまた24の倍数がきたら前の24の倍数
以後からの行すべてを出力というのはどう書いていいのか。。
どうかアドバイスいただけないでしょうか？

mod使って24で割った時の余りの値で制御すれば良いんじゃないかな

すまん。上は24で割って整数化しちゃえばのミス

俺だったら

m=1
do （ファイル終わるまで）
do i = m, m+23
（中略）
end do
m = m + 24
end do

と考えたがスマートじゃないなｗ

>>510,511
整数化というのはどういうことでしょうか？できれば、1列目の数字は変えたくないんですよね。
do i = 1 10000
read(50,*) number,sun
enddo
このiを使って24ずつにしゅつりょく分割できればいいんですが。。

>>512
10000行分書いていくのは、しんどいですよね。ただちょっと書き方参考にさせていただきます。

>>513
ループのiを24で割ると24行毎に値が変わる
1-24,25-48,...と組を作る為に、iをi-1に変えて
k=(i-1)/24として、装置番号に利用すると上手くいく

program kudasure
implicit none

integer :: i
integer :: number,sun
integer :: k
integer,parameter::num=11 !適当な数
character :: output*9,temp*5

do i=1,10000
k=(i-1)/24
write(temp,"(I5.5)") k
output=temp//".txt"

if(mod(i,24)==1) open(k+num, FILE=output, status="unknown")
write(k+num,*)i
if(mod(i,24)==0) close(k+num)
end do

endprogram kudasure

補足すると

write(temp,"(I5.5)") k

は5桁で空の部分に0を補って数字をtempに書き出す（ex 1→0001）やつだから

>>514
UNIT番号は機種依存だが、99までだったりすることが多いぞ。
ファイル名を動的に生成して開くほうが良いアルヨ。

>>514,515

ご丁寧にありがとうございました。
今日自分で書いてやってみます。
また、わからなかったら書き込むと思いますので、また見かけましたら
よろしくお願いします。

>>516
そうだったのか
知らなかった。マジサンクス

いえいえ

ふう

Windowsでg95を入れてプログラミングを始めたんですが、
アイコンへのドラッグ＆ドロップを実装したいんですがどうしたらいいんですか？
よろしくお願いします。

自己解決しました。ありがとうございました。

ttp://n.xxad.net/watch/sm13809685

【どうすればこの欲望の渦巻く堕落した世界に私の名を残せるだろうか】
※名を残したいとは思わないがｗ。
【狂乱のマシンを残せばいい】

【この狂乱の世界に相応しい野獣のようなマシンを】

【狼が来たからやってきた】

【お前の平和な街は今はもう安全ではない】

gfortranで可変FORMAT文を使えますか？

>>524
やってみたら良い。

< >つけてやってみたけど、できませんでした。
可変使わなくても対応できたので、問題はありませんでしたが。

質問です。


Fortran 2003でチェスのプログラムを書いています。

対戦相手の入力待ち時間内に並行してチェスのシミュレーションを行いたいのですがどうしたらいいでしょうか？

たとえばreadを使うと入力するまで計算が止まってしまいますが、計算を止めないまま予測計算を継続させられませんでしょうか？

>>527
F2003規格には、asynchronous I/O というものがある。
理屈の上ではそれを使えば良い。今のコンパイラで実装されているか不明。

野蛮になっていくが、OpenMPとかでOSのAPIでスレッドを投げる手もある。

＞OpenMPとか、MPIとか、コンパイラ並列拡張命令、OSのAPIでスレッドを投げる手もある。

>>528

サンクス！とりあえずIntel Fortran XE Compilerは対応してるっぽいわ。
ついでに関係ないけどMATLABもAsynchronous inputが2010aからできる。

長年の疑問が解決したわ。チョーありがと。

GPU Fortranで書こうかと思ってたけどCUDA Fortran高いからIntel FortranとMATLAB GPUでコードかくことにします。

>>530
おっおっ （ ＾ω＾）

MATLAB買う金あるなら、PGIのコンパイラをアメリカからダウンロードで買える気が・・

>>531

MATLABは海外の大学にいるんだけど工学部全体でまとめ買いしてるからタダなんだわ。
どっちにしろ6000円くらいだよ、アカデミックだと。ツールボックスついてないから無意味だけど。

毎学期技術料で2万くらいとられるけどな！

>>530
自分にスレだけどMATLABは完全にはAsynchronous I/Oに対応してなかった。。。
もうめんどくさいからC++で書こうかと思う・・・GPU CUDA C++サポートしてるし。
>>529
どちらにしろどうもありがと。

>>532
アメリカにいるならWindows用のPGI　Visual Fortranはただで使えるぞ。CUDAにも対応していたはず。
ttps://www.pgroup.com/purchase/freepvf.php
要件は自分で読んでくれｗ

>>534

パソコン一台でしか処理できない点を除けばオールオッケー！

>>534
注文した！やったね！たえちゃん！

おっおっ （ ＾ω＾）＜朕に感謝せよｗｗｗ　
日々、国民の皆様がすこやかに過ごされるよう切に願っておりますおっ！ｗｗｗ

fortran95を使ってます

今open、close文をやっているのですが、入力ファイルの作成方法がどんだけやっても、まったくわからないです

作り方もしくは参考になるサイトを教えてください

openしたら入力ファイルは作成されます。

fortran90使ってます

読み込んだファイルの列数を表示したいのですが、
なにかうまい方法はないですか？

>>540
列の区切り文字はなんだ？　というか、数えたい例を１行書いとくれ。

>>541

区切りはスペースあけてます
例えばこんな感じです

1 1 2 6 5 17 18 22 21
2 2 3 7 6 18 19 23 22
3 3 4 8 7 19 20 24 23
4 5 6 10 9 21 22 26 25

具体的に言うと、この様なテキストを動的割付の2次元配列に入れたいのですが、
その時の列数と行数を知りたいといったところです。
行数はカウントできたのですが、列数の仕方がわからなくて・・・

>>542
区切り文字が決まっているなら、割と簡単。
program vip
implicit none
character(136) :: text
integer :: i, n
read(9, '(a)') text
n = count( transfer(trim(text), ' ', size = len_trim(text)) == ' ' )
print *, n
end program vip

コンマとかが区切り記号なら。
n = count( transfer(trim(text), ' ', size = len_trim(text)) == ' ,' )

transferで有限長の文字列を１文字の配列に直して、そして区切り記号の数を数えている。
欲しい数は、区切り記号の数＋１個。

しか〜し、区切りが任意個のスペースの場合はめんどい。

>>543-544
とても参考になりました！
ありがとうございます

任意個空白区切りの場合（素朴法）　＆　>>543　をやや敷衍した物。
1 　1 2 6 5 17 18 　22 　21　データ例
module vip_m
implicit none
contains
integer function idelim(text)
　 character(*), intent(in) :: text
　 integer :: i
　 idelim = 0
　 do i = 1, len_trim(text) - 1
　　 if (text(i:i) == ' ' .and. text(i+1:i+1) /= ' ') then
　idelim = idelim + 1
　end if
　end do
　return
end function idelim
end module vip_m

program vip
use vip_m
implicit none
integer, parameter :: nmax = 136
character(nmax) :: text
character(1) :: tmp(nmax)
integer :: i, k, n, io

read(9, '(a)') text
tmp = transfer( trim(text), ' ', nmax )
n = count( tmp == ' ' )
print *, n, idelim(text)

end program vip

モジュール内で定義してある関数を、そのモジュール内のサブルーチンで使う時、
サブルーチンの頭で関数を宣言(real some_func)するとエラーになるようになｔった
のですが、これは仕様が変わったのでしょうか？

>>547
そんな事はない。
実例を上げてみてくれ。

>>548
こんな感じです。「integer :: some_func」をコメントアウトすれば通るのですが。
使う関数名は型宣言するものだと思っていたのですが、勘違いでしょうか。

module some_mod
contains
integer function some_func(var)
integer, intent(in) :: var
some_func = var + 1
return
end function some_func
subroutine some_sub1()
integer :: some_func
integer :: i
do i = 1, 10
write(6,*) i, some_func(i)
end do
end subroutine some_sub1
end module some_mod
program func_test
use some_mod
call some_sub1()
end program func_test

>>549
勘違いしていた。宣言側ではなく、使用する側での話か。

F77時代っぽい書き方だよな。
F90 では、副プログラムに関しては、INTERFACE BLOCK を書くのが普通になったし、
MODULE にいれておけば、INTERFACE 書かなくていいから、考えたことがなかったｗ

症状からすると、INTEGER some_func で型宣言すると、EXTERNAL func が暗黙でついてしまうようだが・・・
よほど規格書を詳しく見ないと本来の動作がなんなのか分からん。

こうやって無理やり代替名を与えれば通るようにはなるがｗ
integer function some_func(var) bind(c, name = 'SOME_FUNC')

MODULE内の副プログラムでは、この手の宣言しないのが一番普通だと思う。　

>>549 返信ありがとうございます。

とりあえずは、こういう宣言はしなくていいという事が分かって良かったです。
F90 の言語仕様をちゃんと勉強しようと思います。

allocateにラッパを被せたいのだが、Fortran90の仕様では無理？
配列を割り付けて、そのポインタを返し、
返り値statを見て、エラー処理をしたいのだが。

>>552
できるとおもうぞ。何を心配している？

Assertion failed: ( (last != 0) || (node == 0)) && (tree[last>>TREE_SHIFT][last&TREE_MASK].next == 0), file err_treelib.c, line 256
とエラーが出てしまいます。
原因わかる人いたら教えてください

LOGON TSS
ナツカシすなあ

メムバ、だっけ
ディレクトリにあたるの

フォートランをダウンロードしたいのですがどうすればいいのですか？
散々ググッたにもかかわらずよくわかりません。＞＜
オススメのサイトがあったら教えてください。

>>557
g95でも落としとけ。

g95は落としたのですがいざメモ帳にプログラムを書き込み
保存後、コマンドプロンプトで実行するとNo such file or directory
とエラーメッセージが出ます。どうすればいいのでしょうか？
ちなみにｇ95をタイプした後エンターキーを押すとNo input file
と出るのでフォートランのダウンロードには成功していると思います。

>>559
とりあえず、思いつくこと。
１．日本語の入ったディレクトリを避ける。
２．メモ帳でセーブするとき、拡張子がtxtにならないように” ”で括るなど注意すること。
３．ｇ９５　xxx.F90 のように拡張子も含めたファイル名を与えること。
こんくらい？

しょぼい統合環境でいいなら、フォースを使えルーク。
http://force.lepsch.com/

ありがとうございます。
とりあえずやってみます。

gfortran がタダだし
仕事場のマシンと自宅の私物Win機で
挙動がまあ同じでOpenMP使えるし
デバッグやチェックがIntelのより
親切だから使ってるけど・・・
g95ってイイの？

>>562
二つともフォークしたものだから、根本は変わらないのではないかと。
最近の様子は知らんが、g95はOOP機能を入れる方を優先していたが、gfortranはcoarrayとか並列化に熱心だったイメージ。

漏れのあやふやな記憶が正しければ、元々gfortranのオッサンがひとりでコツコツノロノロg77を90以降に拡張していたが、
gccとかとの統合の頃に、いろんな人が参入してスピードアップ。しかし、オッサンは方向性の違いからまた一人でやりたくなって
フォークした。

ガチ初心者です。

c chouwashindoushi psi v no kaisekikai

c bm:kansanshitsuryou bar:h-bar


implicit real*8(a-h, o-z)

c constants

bm=1.0d0 fk=1.0d0 ome=sqrt(fk/m) bar=1.0d0
…
と打ち込んで（値は適当ですが）コンパイルしようとしたら
エラーメッセージに
bm=1.0d0　…
　　　　 ^
Invalid form for assignment statement at (^)

と出てきました。
何がエラーなのかわかりません！
どう問題なんでしょうか？

一行一ステートメント　一期一会な
＞bm=1.0d0 fk=1.0d0 ome=sqrt(fk/m) bar=1.0d0
bm=1.0d0
fk=1.0d0
ome=sqrt(fk/m)
bar=1.0d0

564です。
返事が遅れましたが、、、
助かりました！感謝です！

>>566
--------------------
[[[[[[[ 返事が遅れましたが、、、 ]]]]]]]（きリッッッッ！！！！
----(きリきリッ
----
[[[[[[[[[[[[ 助かりました！感謝です！ ]]]]]]]]]]]]（きリッ！！キリッッ！キリッッ!!きリ!!!ッッッッ！！！
本当にバカだな

フォートランをやりはじめたばかりの超初心者です
サブルーチンに引数として渡すものと、
モジュールで渡すものはどういう根拠で区別すればいいんでしょうか？
全部引数で渡すこともできるとは思うのですが、
モジュールも組み合わせた方が良い場合もあるんだろうと思います
最初のうちに正しい渡し方を覚えておきたいので...

program quadeq

real a,b,c,d,x1,x2,rpart,ipart
read(*,*) a,b,c
d=b**2-4*a*c

if(d.ge.0) then
x1=(-b+sqrt(d))/(2*a)
x2=(-b-sqrt(d))/(2*a)
write(*,*) 'x1=',x1,'x2=',x2
else
rpart=(-b)/(2*a)
ipart=sqrt(-d)/(2*a)
write(*,*) 'real_part = ',rpart,' ',
'imaginary_part = ',ipart
end if
end

コンパイルできないとです；；
原因を教えてください

Intelのコンパイラ最新版　インストーラでエラーでるんだけど、同じ症状の人いない？
たぶん .net　Ver.4 絡みだと思うんだけど、インストールしてもうごかん・・・
英語版と日本語版の違い？
VS2010でEDITORの補助が増えたらしいので興味あるんだが・・

>>569
>'imaginary_part = ',ipart
この行を前の行び後ろに続けて１行とするか、６カラム目に継続記号を入れればおｋ。

>>568
個人的な考えで回答します。
基本的に、引数として渡す方がプログラムが見やすいと思います。
見やすいということが非常に大切です。
モジュールとしてまとめてしまうと、別の場所(そのサブルーチンとは離れた場所)に
宣言が書かれることになるため、変数が何を意味しているのか、そのサブルーチン独自の
変数なのか、などがぱっと見てわかりにくくなる恐れがあります。
特に、大規模なプログラム、サブルーチン毎に一つのファイルに分けているような
場合には、そのように感じます。
ただ、引数として並べたときに数行になってしまうくらいに引数が多い場合には、
モジュールを利用した方がいいかもしれませんが、その場合でも、サブルーチンを
呼び出すごとに変化するような変数は引数として渡した方がいいでしょう。

一方、全プログラムが一つのファイルに収まってしまうくらいに小規模なものは、
サブルーチンを一つのモジュールの中に入れてしまう方がいいでしょう。ただし、
その場合でも、なるべく引数を使って渡す方がいいと思います。

超絶初心者です。
学校のパソコンでfortran77を使っており、自宅のパソコンでもできるようにコンパイラをダウンロードしようとしたのですが、
Salford FTN77 Personal Edition Compilerをインストールする途中で止まってしまいます。
WinZip Self-Extractor-ftn77pe.exeというウインドウのSetupを押したっきり変化がありません。
エディタは学校でCpadを使っているので自宅のパソコンにもダウンロードしたのですが、コンパイラがうまくダウンロードできているのかもわからず
コンパイラのパスもわからず、どうしたらよいのかわかりません。
いろいろ調べてはみたのですが…

>>572 大変参考になりました、ありがとうごさいます

fortran 77使ってますが、C++からfortranの関数を呼び出すコードのコンパイルを試みてます。

f77 -c -o main.o main.f
f77 -c -o sub.o sub.f
g++ -c Jikkou.cxx
g++ main.o sub.o Jikkou.o -lg2c

するとここで、以下のようなエラー吐きます。

main.o: In function `main_':
main.f:(.text+0x12c): undefined reference to `ran_'

main.fの中で、乱数出す為にRAN()を使ってるんですが、g++はこの関数が分からないようです。
どうにかg++にRAN()を分からせる方法はありませんか？

>>575
おれにはわかRAN(笑)　なんつって（ ＾ω＾）

www

rand() で試すか・・・乱数生成くらい自分でつくれ。

>>575
CからのFortran呼び出しでは、サブルーチン名の規約が処理系によって異なるのでマニュアルを良く読むしかない。
頭にアンダースコアがついたり、全部大文字になったり、引数のバイト数が尻についたり千差万別。

コレをエレガントに避けるにはFortran2003の、C呼び出し規約互換の命令があるのでそれを使うべし。

ただいまの場合、Fortranのライブラリをリンクしていないだけもするので、まあマニュアル読んで頑張れ。
安直な手としては、Fortran90の乱数ルーチンを呼べば外部ライブラリ絡みの厄介は避けられると思う。

Intel Fortran Compilerを使っています。
ファイルから数値を読みとり、配列に値を代入したいのですが、NaNが含まれるためか
「forrtl: severe (24): end-of-file during read, unit 12」
というエラーが出て実行できません。うまく読み取らせる方法はないでしょうか。

real(8) a(100), b(100), c(100), d(100)
do i = 1, 100
read(10,*) a(i), b(i), c(i), d(i)
enddo

0.1000E+00 0.1234E+01 0.0000E+00 　　　 NaN
0.2000E+00 0.5678E+02 0.0000E+00 0.0000E+00
……
1.0000E+01 0.3456E+04 0.0000E+00 　　　 NaN

文字列として読めばいいんじゃない？

あるいは書式に'err=行番号'をつけるとか

>>580
文字列として読み込めばできました。
ありがとうございました。

>>579は本当にそれで解決できたのか？
読み込んだ数値をどう使うのか知らないけど、自分には、
なんで文字列として読み込んで満足できるのかワカラン

>>579, 582さん
エラーの内容は、基番12のファイルの最後以降を読み込もうとした、
ということ。基番10の'NaN'の読み込みのことではありません。
別の原因、別の場所でエラーが発生していると思います。

本当に解決したの？ 　583さんと同意見です。

4次元配列を2次元に配列にまとめたいのですが、
The shapes of the array expressions do not conform.と出てコンパイルできません
なにかうまい方法ってありませんか？

具体的には
do i = 1, 10
do j = 1, 10
do k = 1, 100
do l = 1, 100
y(i,j) = y(i,j) + x(k,l,i,j)
enddo
enddo
enddo
enddo

といった計算を行いたいのですが・・・

「4次元配列を2次元に配列にまとめたい」っていう目的と
実際に提示されているコードが自分には結びついていないけど、
配列変数の定義で配列の次元やサイズはちゃんと書いた？
real(8) :: x(100, 100, 10, 10)
real(8) :: y(10, 10)
みたいな記述が必要だと思うけど。

>>586
配列の定義はしてあります

x(k,l,1,1)のkとlの合計10000個の数をy(1,1)に足しこんで
次はx(k,l,1,2)のkとlの合計10000個の数をy(1,2)に足しこんで

といった具合にプログラムを書きたいです
説明下手ですいません

>>585 で何が不満なんだ。

f95バリバリならこうか？FORALLはうんこ機能と言われているがｗ
program Console1
implicit none
! Variables
real :: y(10, 10), x(100, 100, 10, 10)
integer :: i, j
! Body of Console1
y = 0.0
x = 1.0
forall (i = 1:10, j = 1:10) y(i, j) = y(i, j) + sum( x(:, :, i, j) )
print *, y(1, 1)
stop
end program Console1

forall ってウンコというか
使い間違いし易いのよね。
最初の4重doループグリグリでいいと思うけどね。
コンパイラが勝手に最適化してくれるだろうし。

PGI Visual Fortran 11.9 のコマンドラインでのコンパイルについて質問です。

以下のようなソースコード goodbye.f90 を作りました


program goodbye
write(*,*) 'Good-bye'
stop
end program goodbye



これを C:\tmp\に保存して


c:\tmp>pgf90 goodbye.f90


と打ったのですが，カーソルが次の行に行って点滅するだけで，
コンパイルはされませんでした。
何がまずかったのでしょうか？

>>590
それだけでは分からんｗ
まぁ、end program goodbye の後ろで改行しないまま、end of file に成っている可能性があるかな。

>>591
レスありがとうございます。
改行はしてありました。

なお，コンパイルを Ctrl + C で強制終了すると，

cleaning up after signal(2)

というメッセージが表示されます。これについても意味が分からないです。

>>592
それは強制終了の事後処理というメッセージだろう。

症状としては、end of file が見つからない感じだから、ｃｔｒｌ−Zないしｃｔｒｌ−Dあたりをうってみたらどうかな？

>>592
報告が遅れすみません。
再々インストールを行ったら，治りました。
ソフテックに訊いたら，ウイルス対策ソフトがブロックしていたらしいです。

>>594
そういえば一時期ノートン先生がIntelVisualFortranの実行ファイルを削除しまくって閉口したことがあったわ。

http://mixi.jp/view_bbs.pl?id=65775704&comment_count=0&comm_id=88243
mixiしてる方に質問させていただきます。
この拡散方程式を種々の初期条件、境界条件のもとに解くことによって、
横軸：時間、縦軸：ペイント全体の濃度合計
→ペイント全体の濃度の合計の時間変化を
fortranにてプログラム化できますでしょうか。。。

>>596
具体的にとこで困ってるんですか？
離散化とか、境界条件の設定とか、具体的に言ってください。

あとインプリシット、エクスプリシットどちらで解くつもりですか？

>>597さん
拡散方程式(２階偏微分方程式)の離散化は自分で行えます。
その拡散方程式を種々の初期条件、境界条件で解くプログラム
が作りたいのです。
私がわからないのは、拡散方程式や種々の条件をfortranにて示す方法が
まったくわかりません。。。
周りに相談できる方がいないので本当に困っています↓↓
この問題は、「ある気体を透過しない母体に、気体透過性のペイント
を貼った場合、その透過性のペイントの表面の気体濃度を正弦的に変化
させた場合、ペイント内部ではどのように気体濃度が変化していくか
ということを見るためのプログラム」作りたいのです。
言い換えるならば、縦軸：ある時間におけるペイント全体の濃度の合計
横軸：0<=t<=tm(tm→とりうる時間の最大値)
を作りたいのです。

>>597

とは違うけど、インプリシット(陰解法、連立方程式を反復的に解く)
とエクスプリシット(陽解法、中心差分のように、漸化式を反復的に繰り
返して解く？だったっけ。もう、忘れた)の意味わかってます？

プログラムを作成する前段階のアルゴリズムの定式化は本当にできて
いるんですか？

定式化ができた上で、fortaranでの表現法がわからないというのなら
fortranの勉強をするべきだと思うけど。

もうMATLABを使ったらと思うけど、取り組んでる問題をみるとMATLAB
ではヘタすると数日、計算がかかるかもしれませんね。

離散化ができているなら、流動方程式を解くプログラムに、拡散項を組み込むだけですから、そんなに難しくないと思います

流動方程式のプログラムなら適当な参考書にあると思います...

あと、境界条件が不安のようですけど、ディリクレ境界条件のようですから、ポイントセントラルでプログラムを組めば少しらくだと思います

const int N = 100;
const double q = 10.0, dt = 0.00001, Dm = 30.0, t0 = 2.0, K = 1.0, pi = 3.141592, f = 3.0;
double C[N], dC[N];
double dx = q/N;
for (int i = 0; i < N; ++i) C[i] = 0;　// 初期条件
for (double t = 0; t < t0; t += dt) {
　C[0] = 0;　// 境界条件1
　C[N - 1] = K*sin(2*pi*f*t);　// 境界条件2
　for (int i = 1; i < N - 1; ++i) dC[i] = (Dm*(C[i + 1] - 2*C[i] + C[i - 1])/(dx*dx))*dt;
　for (int i = 1; i < N - 1; ++i) C[i] += dC[i];
}
for (int i = 0; i < N; ++i) cout << C[i] << endl; // t = t0
を誰かfortranに変換できる人いますか？

返信ありがとうございます。
１ふぉｒｔｒ

返信ありがとうございます。
fortran始めたばかりなので正直知識不足は否めないです。
おそらくエクスプリシットだと思います。

アルゴリズムの定式化についてはまだ研究初期段階なので今後
修正する可能性はあります。
修正に応じてプログラムの方も修正する予定なのですが、
まずベースとなるプログラムが欲しくて、
このような投稿になってしまいました。
不快でなければ協力のほうをよろしく
お願いいたします。

http://soudan1.biglobe.ne.jp/qa7094085.html

program vip
implicit none
integer, parameter :: n = 100
real(kind = 8), parameter :: q = 10.0d0, dt = 0.00001d0, dm = 30.0d0, t0 = 2.0d0, k = 1.0d0, f = 3.0d0
real(kind = 8), parameter :: pi = 4.0d0 * atan(1.0d0)
real(kind = 8) :: c(n), dc(n), t, dx
integer :: i
dx = q / n
c = 0.0d0
t = 0.0d0
do
if (t > t0) exit
c(1) = 0.0d0
c(n) = k * sin(2.0d0 * pi * f * t)
forall (i = 2:n - 1) dc(i) = dt * dm * ( c(i + 1) - 2 * c(i) + c(i - 1) ) / dx**2
c = c + dc
t = t + dt
end do
write(*, *) c
stop
end program vip
正しいかどうかは知らんｗ

デフォルトの名無しさんへ
ありがとうございます！！！もとになるベースがないので本当に
困っていました↓↓
返信大変うれしいです！！！

>>604,607
二度と書き込むなよ、このクズ野郎

どこのスレも変なのがうろついてんな

困っちゃん

枯れ木も山のにぎわいってことでｗ

>>608

お前の方が余程クズで二度と書きこむなよ、と思うけど。







あと、死ね

http://qanda.rakuten.ne.jp/qa7112807.htmlにある質問なのですが
これをFORTRANでできる方はいらっしゃいますか？

あとこのプログラミングをC++でやってみたところ
const int N = 100;
const double q = 10.0, dt = 0.00001, Dm = 30.0, t0 = 2.0,
K = 1.0, pi = 3.141592, f = 3.0;
double C[N], dC[N];
double dx = q/N;
for (int i = 0; i < N; ++i) C[i] = 0;　// 初期条件
for (double t = 0; t < t0; t += dt) {
　C[0] = C[1];　// 境界条件1
　C[N - 1] = K*sin(2*pi*f*t);　// 境界条件2
　for (int i = 1; i < N - 1; ++i) dC[i] = (Dm*(C[i + 1] - 2*C[i] + C[i - 1])/(dx*dx))*dt;
　for (int i = 1; i < N - 1; ++i) C[i] += dC[i];
double s = 0;
double s_prev;
for (int i = 1; i < N - 1; ++i) s += C[i]*dx;
if (t) {
　double ss = (s - s_prev)/dt;
　cout << ss << endl;　// ∂s/∂t を出力
s_prev = s;
}
といちおやってみたのですが、間違いを指摘できる方はいますか？
あとできればこのプログラムをFORTRANに変換もしていただきたいです。。。
重複と長文同時に失礼いたしました。
図はhttp://qanda.rakuten.ne.jp/qa7112807.htmlにあります。
全開質問に答えていただいた方、本当にありがとうございました。

c++ はよくしらんが、｝ が複数個欠けてるような。
エディタはたいてい対になる(){}[] を色で表示してくれるので、それを助けに頑張れ。

Courant 条件は・・
(1/2) (dx*dx) / Dm= 0.5 *(0.1)^2 / 30.0 > dt
なので、よさげ。{}対をチェックしたら、ちゃんと走ると思うよ。

>>613
Fortran95への移植。内容は見ていないｗ　Sinカーブのように振動しているから、まぁこんなもんではないか？
program vip
implicit none
! Variables
integer, parameter :: n = 100
real(8), parameter :: q = 10.0d0, dt = 0.00001d0, dm = 30.0d0, t0 = 2.0d0, &
k = 1.0d0 , pi = 4.0d0 * atan(1.0d0), f = 3.0d0
real(8) :: c(n), dc(n)
real(8) :: dx, t, s, s_prev
integer :: i
! Body of Console3
dx = q / n
c = 0.0d0
t = 0.0d0
s_prev = 0.0d0
do while( t < t0 )
c(1) = c(2)
c(n) = k * sin(2.0d0 * pi * f * t)
forall(i = 2:n - 1) dc(i) = (dm * (c(i + 1) - 2 * c(i) + c(i - 1)) / dx**2) * dt
c = c + dc
s = sum(c) * dx
write(9,*), 't=', t, (s - s_prev) / dt
s_prev = s
t = t + dt
end do
end program vip

>>615 のバファリンのようなやさしさに感動した

>>615の人に質問
DO文ではなくFORALL文を使った理由は何ですか？

>>617
１行で書けるから。

f90とf95ってまったくコードの違うものなのですか？

>>619
f95 は f90 の小改訂版。Fortran は f90 以降、メジャーウップデートと
マイナーウップデートを交互にやることになってる。f95 は f90 を使ってみて
いや〜んな点を、かゆいとこに手が届くように修正してある。

あと変なの湧くからsageろよ。

>>620
交互にマイナーとメジャーのアップデートとのことですが、
それなら、今のJIS規格であるFortran2003はメジャーアップデート
みないなものなのでしょうか。

>>621
F2003はメジャーアップデートどえす。
目玉はオブジェクト指向の要素の導入ですが、その他色々便利な改良もあります。
でもあまりに欲張り過ぎて、まだ完全装備のコンパイラが殆ど無いという事実ｗ

ちなみにF2008はマイナーアップデートのはずなのに、coarrayのような並列化命令が入った規格が確定した。

coarray よさげだよね。
でも完全装備は先かな・・・intel はやる気あるみたいだけど。
>>615
ほどんど重箱隅レベルだけど、個人的には
c(:) = c(:) + dc(:)
s = sum(c(:)) * dx
と配列であることを明示したい。・・・ドツボ防止策。

intel は　coarray 結構実装してきてるんじゃない？　サンプル以外試してないから詳細は分からんが。

ｃ(:)の配列明示は理念としては理解できるし、ｆ７７から移行してきたばかりの昔はやっていたのだが、
f90初期の頃、処理系によってはコロン指定子がつくと無条件にコピーを作って作業して戻すという事をするので、
最適化の観点からなんとなく避けるようになってしまった。最近の処理系はそういうことはしなくなったようだが・・・

サブルーチンの引数でも、生で書くとそのまま渡すが、コロンがつくとcopy in / copy out　になって最適化が落ちるとか、
c(0:n) の場合、生で渡すとサブルーチン側でも配列は０から始まるが、コロンがつくとコピーが渡るので１から始まるとか、
その辺の規則がよく分からなっかったので、なんとなく惰性で配列明示しないで生で行ってるｗ

>>624
ためになった。

ついでに質問。coarray って分散メモリのシステムでもＯＫなのかな？
MPI とどちらがいいのかなとか、知ってたら書いて。

>>625
coarray は分散メモリーシステム用だな。
MPIより簡単で同等の性能が出るという主張をしている感じ。
分散メモリーマシン時代になってからのクレイが熱心にやっていて、それが規格に導入された？

うん

>615さん
ではf90に書き直すとそう変わるのですか？？

ほぼ>>615に等価だと思う。

>>626
coarrayて分散メモリー用なの？
http://www.xlsoft.com/jp/products/intel//compilers/fcl/12/comp-for-lin/comp/tutorial/tutorial-4/co_t4b.htm#co_t4b
のサンプルを見ると、OpenMP拡張とは混在できないとか書いてあるので、
てっきり、SMP用かと思っていた。

なんか前に
ノード内の共有メモリ空間でも OpenMP で書くより
ノード間も含めてもう全部 MPI でやり取りした方が、
速い(＋楽！)という意見があったような。

なんだろう、共存が面倒な割にはパフォーマンスがあんまりなので
もうしない、という事なんだろうか？

>>630
coarrayでググるとCrayのパワポとか出てくるから味噌。
現今のcoarrayは結局MPI上に実装されているので、MPI利用した簡便並列という感じ。

単純な仕事の分割の場合、MPI で書いた並列プログラムは
ループの前後に冗長に似た行が続くだけだから、その部分の追加を
ちゃっちゃとやってくれるのなら、うれしいね。

Open/MPはスレッド並列だが、MPIやCo-arrayはプロセス並列になってる。

gfortran, バージョン4:4.4.5-1 に関する質問です。

モジュールのバージョンが違うとかで下記のエラーメッセージが出てきます。
Fatal Error: Parse error when checking module version for file 'functions.mod'

これってどう解決すればいいですか？

コンパイルしなおせ

>>636
え？コンパイルしようとするとこのエラーが出るんですけど

gfortran hogehoge.f90
ってな具合です

functions.mod というのは、module functions がコンパイルされたときに生成されるインターフェース情報のファイルで
これのversion　が違うというのだから、これをもう一回コンパイルするべきじゃないかと思う。
最近gfortranをversion upしていたりするならよくある事。

そうでないならば、情報が足りなすぎてなんとも言えない。

>>638
あ、そのとおりです。
しかもフォートランを最近うぷしました。

これってuse mod名をメインプログラムに書き込む前に、一度コンパイルしてmodファイルを作成しなくてはいけないということですよね？

>>635
Debianスレでも聞いてるじゃねーか

>>640
え？聞いてないよ。
今確認したらコピペされててビビった。
信じて

関連するファイルを全部コンパイルしなおせ。

たかだか１行の英文エラーメッセージ。単語１１個くらい字引で引けや！インテリの端くれだろｗ
ピカチュウの鳴き声を解釈するより易しいだろ。
それで解決したのか？

>>643
解決致しました。アニキ

一度プログラム中のuse (mod名) 部をコメントアウトし、コンパイルすることでモジュールファイルを生成。
その後コメントアウトした部分を復活させて、再度コンパイルしたら無事できました

>>643
それにしてもググってもこれと同じ症例があまりなかったのに驚いた。
みんなどうしてるのかな〜

module functions　も、引用する program と一緒のファイルに入れてるのか？
それなら、program の前に来るようにして先にコンパイルされるようにすれば問題ないと思うぞ。
gfortran使ってないから分からんが常識的にはそれでおｋ。

それが駄目なら、ファイルを別々にして分割コンパイルすればよろし。

>>646
アニキありがとうっす！
それやってみます。

アニキとサブとファンクチオン

program moho
implicit none
logical:: asshappy
logical external:: playresult
asshappy=.false.
call meet_my_aniki()
do while (.NE. asshappy)
asshappy=playresult()
enddo
end program moho

REAL GAY
COMMON //ASS, HOLE
IF ( ANIKI .AND. SUB ) GOTO 69
69 CONTINUE
END

>>649
(.NOT. asshappy)
>>650
リアルホモｗ
でもその文だとなにも変化がないような気がする。

今だにGOTO文使ってるとかおっさんかよｗ

＼　　　/　.::::::::::::::::::::::::;;:;;::,ｯ､::::::　　 )　　く 　 ホ 　す
　 ＼　l　 ,ｯｨrj,rｆ'"'"'"　　　 lﾐ:::::::　く　　 れ　 モ　 ま
　　　　 Y　　　　　　　　　　 ,!ﾐ:::::::　ヽ　　な　 以　 な
｀ヽ､　　|　　　　　　　　　　 くミ::::::::　ﾉ 　 い　 外　 い
　　　　 |､__　 ｬー--＿ニゞ　｀i::::,rく　　　か　 は
｀`''ｰ-　ゝ､'ｌ　 ￣ﾞ´彑,ヾ　　　}::;! ,ﾍ.）　　！　 帰
　　　　　 ﾞｿ　　　""""´｀　　　　　〉 L_　 　 　 っ
　　　　　 /　　　　　　　　　　i　　,　　/|　　　　て　　　　r
≡＝-　〈´　,,.._　　　　　　　 i　　't-'ﾞ　|　,へ　　　　　,r┘
,､yx=''"　｀ｰ{ﾞ　＿, -､　　　　;　　l　　　ﾚ'　　ヽr､⌒ヽ'
　　　　　　　 ﾞ､`--─ﾞ　　　　　 /!　　　　　　　　 `､
　　_,,､-　　 　 ﾞ､ ー''　　　　／　;　　　　　　　　　　 `､
-''"_,,､-''"　　　 ﾞ､　　　　／;;'　,'　　/　　　　　　　　 、＼
-''"　　　 ／　　　`ｰ─''ぐ;;;;' ,'　　ノ　　　　　　
　　 ／／　　　 /　　　　　ヾ_､=ﾆﾞ

まさかfortranスレでこのAAを見ることがあるとは

GO TO文に怯えるとは、まだまだ坊やだな。
後藤おじさんが手ほどきしてあげよう。
まず全裸になれ。

なりました。
つっ次はどうすればいいのでしょうか?

>>656
次はションベンだ！
program o2chan

implicit none
character(10) :: siri

call sub( siri )
print *, siri
stop

contains

subroutine sub( anal )
character :: anal(*)
anal(1:10) = transfer('しょんべん', 'a', 10)
end subroutine sub

end program o2chan

　　　　　　　∧∧
　　　　　　 (д｀*　)
　　　　　　 （⊃⌒*⌒⊂)　　anal(*)
　　　　　 　　/__ノωヽ__）

みなさま、はじめまして。
趣味でとある連続計算をしています。
いままで、連立一次を解くサブルーチンをつかっていたのですが、
わけあって逆行列単体でシンプルなサブルーチンが必要となりました。
ネット上にN×Nの非対称行列を逆行列にして返してくれるサブルーチンってありますでしょうか？
ご存知なら教えてくださいませ。

LAPACK？

>>659
逆行列は、Ax=bのbに（0010..0)のたぐいの単位ベクトルをｎ本入れて計算するのが基本だぞ。
LU分解のサブルーチンでやればn回繰り返すのも楽勝だ。

>>660
ぐぐってみました。
ソースがHTMMLでUPされてないと、見る気になりませんでした。
すみません。

fortran 逆行列
でぐぐって一番上にでてきたサイト(とある大学の先生のページ)の
「逆行列を求めるサブルーチン」
を実行したんですが、配列異常のエラーが出ました。

連立一次を解くプログラムがあれば、それで逆行列くらい求まるだろ。
それが逆行列を求める定石だぞ。
単位行列のを、１列ごとに連立一次を解いて、それを行列の列が尽きるまで繰り返す。
出てきた答えを、列ベクトルとして並べたものが逆行列だ。

ネットで探して見つからないのは、そうやって計算するから、一次方程式のサブルーチンがあれば十分だからだ。
あるとしても、単に１次方程式を解くルーチンを呼び出すルーチンになってるだけ。

特にLU分解で連立一次方程式を解いてあると、毎回一から解くのではなく、前回の途中経過が利用できるので早く済む。
計算量はN^3に比例するので大きい時はLU分解使わないのは非常識。
問題が小さければ、好きな方式でやっておｋ。計算法にこだわるほうが非常識。

くりゅろふ法自作試みて結局投げだして
ガウス消去に逃げた事があるw

逆行列計算は奥深いよね。
いまでも計算科学ではメインのテーマの一つ。

で、やっぱりLapack とか MKL とか、かな

LU分解とか一列ごとに解くとか言ってるけど、直接法ならまとめて解けばいいと思うよ。AX=BのBを単位行列にすれば良いだけの話。

inplicitrial*8
とかまだやってるやついるの？

綴り間違えすぎｗ　それに　implicit real*8(a-h,o-z)　 だろｗ
耄碌しすぎだぞ爺さん。

今のヤングは implicit none で A'GoGo!!

さらにヤングｗだと
コンパイルオプションで暗黙の型宣言は消去！

ttp://ksirius.kj.yamagata-u.ac.jp/~shibata/kbg/invmx.f
このプログラムを参考にしています
implicit real*8(a-h,o-z)
double　plicision
で書きなおしたら、しっかり計算されません。
どうしましょ

>>669
全く見ていないがｗ　よくあることなので推測で書くｗｗｗ
常数のたぐいは手で1.0から1.0d0のように倍精度に直していないからだろう。
特に収束判定に絡んでくると、単精度のままだとゴミビットが入って精度が腐る。

その前にスペルミスを突っ込みたい。
コンパイルはできたんか？

すみません1.0d0と1.0_8はなにか違いますか

>>672
kind種はベンダー依存なので、今の機種でkind=8が倍精度を意味しているからといって
他の機種でそうなるとは限らない。
ただ世間の慣例としてバイト数を置いていることが多いので、大抵は同じ意味。

もっとも上品にはKIND種をパラメータ変数にしておく。日本NAGの書き方推奨ページに行くとグダグダ書いてある。