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

**デフォルトの名無しさん** · 2012/08/16(木) 20:58:00.31

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

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

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

●過去スレ
くだすれ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/

**デフォルトの名無しさん** · 2014/09/17(水) 01:15:32.64

>>603
当該 READ 文の FORMAT 文を見ると分かる。
たぶん昔のプログラムだから READ(99, 1000) N,B,C... みたいになっているだろう。
その場合行番号付で 1000 FORMAT(1H , 云々) となっている FORMAT文を見ればいい。
その行と引用しているREAD文を2chにコピペすればいい。

どの行かわからんというならDEBUGオプションでtracebackをオンにしてやれば
エラーで死んだところで行番号が出るはず。

**デフォルトの名無しさん** · 2014/09/18(木) 00:30:49.37

Fortranってどの程度配列メモリに記憶出来てるのかな？

例えば
do i = 1, 100
Read(20,*) (a(i,j),j=1,20
End do
で write(*,*) a(99,15)とかどこまで記憶出来てるもんなんやろ？

**デフォルトの名無しさん** · 2014/09/18(木) 21:41:43.06

eclipseである程度長いサブルーチンも何十個か呼び出しているプログラムで
call systemが通らくなるのですがなにか原因ありますか？
ブラウザソフト立ち上げながらでもダメだったのでメモリがらみと思うのですが

ざっくりしすぎですいませんがなにかヒントあればおねがいします・・

**デフォルトの名無しさん** · 2014/09/19(金) 00:43:16.70

>>605
質問の意味が分からんよ？？？
Fortranじゃ宣言分はがっつり確保されとるよ。

>>606
タスクマネージャーとか ps コマンドとかでメモリーの使用状況を見てみたら？
call system が通らないのはメモリー満杯くさい。

**デフォルトの名無しさん** · 2014/09/19(金) 01:31:53.81

>>607
ありがとうございました
ps alxで見るとrssがまさに満杯でいっぱいのようでした
cat /proc/meminfoで見るとフリーが１０Gぐらい残っていて
大丈夫と思っていたのですが違いが分かっておらず
違いを勉強しようと思います
もしご存知でしたらまた教えていただけると幸いです

**デフォルトの名無しさん** · 2014/09/19(金) 03:08:40.00

>>607
そうなの？

**デフォルトの名無しさん** · 2014/09/19(金) 22:55:20.41

>>608
http://tweeeety.hateblo.jp/entry/20140130/1391054452
色々まとめてあるみたい。

**デフォルトの名無しさん** · 2014/09/20(土) 20:11:45.61

>>610
ありがとうございます。べんきょうになりました。

メモリ節約のために可能な限りallocatableに変更したら余裕ができて解決しました。
個人利用のプログラムには必要ないと思っていましたが大事でした
ありがとうございました

**デフォルトの名無しさん** · 2014/09/22(月) 04:03:57.88

配列a=/2,3,4,5,8,5,4......./みないなのがあったとして, 最大値の8すなわち,a(5)までをファイルに書き込みたいんだけど…

do i = 1, maxloc(a)
write(20,*) a(i)
Enddo
ではエラーになるよね。おそらくmaxloc(a)がrank 1 でdo loop にはrank 0 (スカラー値)じゃなきゃいかんってことなんだろう
配列をスカラー値に戻す関数ってないのかね？

**デフォルトの名無しさん** · 2014/09/22(月) 08:08:38.00

配列は要素を指定すればスカラーになるだろ

**デフォルトの名無しさん** · 2014/09/23(火) 00:55:12.38

>>612
次元を指定すればおｋ。１次元でもあえて１を指定すればスカラーで返してくれる。
do i = 1, maxloc(a, 1)
１要素配列を返す奴らに使える、基本的なテクニックですぞい。

改行要らなければ
write(20,*) a(:maxloc(a, 1))
で１行で済む。

**デフォルトの名無しさん** · 2014/09/23(火) 02:10:29.02

>>614
いやぁ助かります。
ここは参考書に載ってないテクを知ってる方がたくさんいるからやめられない…

このようなテクって自分で見出すもの？それとも何か参考書みたいのってあるのですか？

あと、科学数値計算やってる方に聞きたいのですが、最適な刻み幅ってどう決定してますか？
対数刻みやらなんやらありますが、結局自己流やらなんやらでしっくりしたものに出会えず時間を無駄にしてる感たっぷりなんです

**デフォルトの名無しさん** · 2014/09/23(火) 12:59:19.57

>>615
maxloc とかに dim 引数が入ったのは F95 以降で、
F90 時代は要素1個の配列を返されて誰もが窮していた。

Numerical Recipes の F90 版では、SUM(MAXLOC(a)) みたいな形でこの問題を
避けていたが、俺的にはこれは無いなって感じだった。（元々あんまいい本じゃなかったし）
仕方ないから律儀に要素1個の配列を宣言して使っていて、もう解決策はないと諦めていた。
なので f95 が出てからも、この方法に全く気付かずにいた。

ある日、他人のソース眺めていて、１次元配列なのに maxloc の類で dim
引数指定していたのがあって、こいつ馬鹿じゃねーのｗｗｗと思ってよく見たら
このスカラー返しのテクニックだった。馬鹿なのは私でした。

たぶん１０年以上無駄な1要素配列のプログラム書いていたｗｗｗ

**デフォルトの名無しさん** · 2014/09/26(金) 01:28:02.38

誰か>>612 >>614の解説をしてもらえないでしょうか
まったく意味がわからない

**デフォルトの名無しさん** · 2014/09/26(金) 02:38:03.65

>>617
配列 a(7) = (/ 2,3,4,5,8,5,4 )/ があるとする。
最大値の位置を返す組み込み関数 maxloc に a を渡すと maxloc(a) は
１次元で１要素の配列 (/5/) を返す。

これを左辺で取るには integer :: m(1) が必要。スカラー変数では受け取れない。
m = maxloc( a ) !ok 　; integer :: n ; n = maxloc(a) !error

要するに do i = 1, maxloc(a) は、do i = 1, m(:) と同じ意味になり文法エラーが出る。
>>612 の質問は、do i = 1, m(1) に当たる事を実行したいので、
配列をスカラー値に戻す方法はないか？と質問している。

それに対する>>614の答えは、maxloc( a, dim = 1 ) を使えと言っている。
これは配列の要素を返す。普通は多次元配列に対して用いて、行や列の成分を
取りだすのに使う。

**デフォルトの名無しさん** · 2014/09/28(日) 02:23:58.55

>>618
1次元配列の特定の1つの値を返すのに、それは1次元の1要素の「配列」として返されるからdo文では使えないということですね。それが配列の次元(dimension)を1と言うと解決するあたりがしっくりきませんね・・・。

**デフォルトの名無しさん** · 2014/09/28(日) 03:11:58.13

>>619
多次元配列を maxloc に渡すと、最大値の座標を (/ 1, 2, 3 /) みたいな感じで
一次元配列で返す。この成分を取り出すのが dim 引数。
１要素配列でも成分取り出せる。
まぁしっくりくるまで文法書を舐めるように読んでくれｐｒｐｒ

**デフォルトの名無しさん** · 2014/10/08(水) 00:58:36.59

x(100)があるとして、計算しているのになぜか例えばx(40)から値が0になってしまいます
考えられる原因ないでしょうか
無茶ぶりですいません

**デフォルトの名無しさん** · 2014/10/08(水) 01:09:28.89

さすがに無理
もう少しヒントをくれ玉へ

**デフォルトの名無しさん** · 2014/10/08(水) 07:15:04.02

>>621
日本語になってない。
逝ってよし

**デフォルトの名無しさん** · 2014/10/08(水) 08:00:58.16

今時逝ってよしとかねえよ

**デフォルトの名無しさん** · 2014/10/08(水) 20:12:54.14

>>622-623
ありがとうございます
確かにひどかったです・・すいません

ただ、症状がなんとも表現できていないので
もう少しうまく表現できるよう考えてみようと思います

**デフォルトの名無しさん** · 2014/10/09(木) 00:38:44.90

>>624
そんなこと言ったら今時fortran なんかねーよと言われるぞ。
そんなことないだろ？つまり逝ってよしもそんなことないのだ。

**デフォルトの名無しさん** · 2014/10/09(木) 02:25:56.01

>>625
xが途中で変わってしまうなら

xの代入文や計算式が間違ってないか探す(特に範囲指定)
途中で変数を書き出し、変数の変化を追う
デバッガを使って変数に代入されるタイミングを調べる

くらいしか思い付かない
実は0になる計算をしていたなんてことはないよね

**デフォルトの名無しさん** · 2014/10/10(金) 03:19:01.38

まぁ大方書き込んだつもりになっているが、ループが３９までしか回ってないという
パターンだろう。
初期値に1e38とかNaNとか書いておけばいい。

**デフォルトの名無しさん** · 2014/10/11(土) 16:00:47.58

>逝ってよし
逝ってよしってフレーズ最後に見たの何年前だろう？１０年くらい前？
ちょっとクスッと来てしまった。

悪くないよ～

**デフォルトの名無しさん** · 2014/10/12(日) 11:19:06.43

>>627
>>628
ありがとうございました
基本的なミスだろう、というご教示を参考に見なおしたところ解決しました
x(40)以降を計算するのに地味にyを使っていたのですが、
yを計算するサブルーチンをcallする前にxを計算するサブルーチンをcallしてしまっていただけでした・・
ズコーって感じかもしれませんが私としてはレスいただき大変助かりました

**デフォルトの名無しさん** · 2014/10/12(日) 22:43:53.80

gfortranで
call system(cd ..)
とやっても上のディレクトリに移動できません。

他のコマンドは動くのですが…

これは仕様なのでしょうか？

**デフォルトの名無しさん** · 2014/10/13(月) 20:38:43.03

>>631
systemからのcd実行なら仕様。
別のshellが起動してcdを実行して処理を終了するので、呼び出しもとのshellとは無関係。
GNUならディレクトリを変更したい場合はCHDIRが使える。

**デフォルトの名無しさん** · 2014/10/13(月) 21:42:58.17

>>632
ありがとうございます。
if文でCHDIRに置き換えれば良さそうですね。
仕組みもわかりました。

**デフォルトの名無しさん** · 2014/10/14(火) 10:25:49.65

とか言って失敗しましたｗ　CHDIRをつかうとshからpremission deniedと怒られます。

if (IPTCMD(1:3).eq.'cd ') then
DSTDIR=IPTCMD(4:)
call CHDIR(DSTDIR)
end if

のようなコードです。SELinuxのせいかと思いましたが違いました。
ご教示いただければ幸いです。

**デフォルトの名無しさん** · 2014/10/14(火) 10:53:43.84

追加すると、cdコマンドは成功しています。でも、

sh: 1: ..: Permission denied

と怒られますｗ

**デフォルトの名無しさん** · 2014/10/17(金) 08:36:55.32

素朴な疑問なんですが
write( *,101) n
write(10,101) n

こういう式って一行にまとめられないんですか

**デフォルトの名無しさん** · 2014/10/17(金) 14:01:54.72

>>636
write ... ; write ,,,

**デフォルトの名無しさん** · 2014/10/18(土) 03:21:07.25

言い換えます
write( *,101) a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u,v,w,x,y,z
write( 10,101) a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u,v,w,x,y,z
write( 11,101) a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u,v,w,x,y,z
みたいな文を一行でまとめたいのです

write(10,101) n

**デフォルトの名無しさん** · 2014/10/18(土) 04:32:28.37

10,11だと外部ファイルか
それだと複数の装置番号に同時に出力する方法はない（0以上のスカラーでなければならない）

変数をまとめたいなら配列nにa-zを代入すればできる

**デフォルトの名無しさん** · 2014/10/18(土) 06:08:15.41

>>638
出力するユニットが同じならサブルーチンでまとめれば良いんじゃない？

**デフォルトの名無しさん** · 2014/10/18(土) 07:53:57.67

fortranからRを呼び出すようなことはできないでしょうか？
system call?からRを呼び出すぐらいしかないでしょうかね？

**デフォルトの名無しさん** · 2014/10/18(土) 08:55:47.73

gnuplotでもそうだが
Rのソースをfortranで生成するようにして
system callから実行させればいいだろ

**デフォルトの名無しさん** · 2014/10/18(土) 09:10:45.91

>>638
こんなのかな

iu(1:3) = (/6,10,11/)
do i=1,size(iu);
write(iu(i),101) ...
enddo

**デフォルトの名無しさん** · 2014/10/18(土) 09:18:33.82

考えてみれば四行に増えたなｗｗｗｗｗｗｗ

**デフォルトの名無しさん** · 2014/10/19(日) 01:11:43.08

たしか coarray を使うとプロセスが複数走って、それぞれのプロセスで
別々にI/Oできるから形式上は１行でできるぞ。
まぁこれがやりたいことではないのは重々承知だがｗｗｗ

**デフォルトの名無しさん** · 2014/10/19(日) 13:24:38.74

>>642
サンクスです
ライブラリみたいのはないのかもしれませんね

**デフォルトの名無しさん** · 2014/10/20(月) 00:32:03.39

配列内の最大最小を、区間指定でできませんか？
maxだとできなかったんですが

**デフォルトの名無しさん** · 2014/10/20(月) 01:10:20.17

区間指定ってなんだよｗ　
maxlocのことか？それともmaxval(a(20:20))のことか？

他人に理解できる表現で頼むぜｗ

**デフォルトの名無しさん** · 2014/10/21(火) 04:04:46.35

>>647
部分配列に代入して、maxを使うという意味？

**デフォルトの名無しさん** · 2014/10/21(火) 08:48:40.55

＞615
亀レスですが、最適な刻み幅については、

森口繁一『数値計算術』に
「題4章　刻み幅の自動調節」
という記述があります。
ご参考までに

**デフォルトの名無しさん** · 2014/10/22(水) 21:06:47.85

>>648
>>649
a(1:20)があったとして、10から20の間だけのaのmaxという意味で書きました
わかりづらくてすいません

**デフォルトの名無しさん** · 2014/10/23(木) 00:20:38.34

>>650
ありがとうございますm(__)m
さっそく買わせて頂きました。
まだ手元に届いてないですが、その章は結構ページ割かれてますか？

**デフォルトの名無しさん** · 2014/10/23(木) 00:27:54.57

ならば　MAXVAL(a(10:20)) だ。

**デフォルトの名無しさん** · 2014/10/23(木) 08:52:15.53

>652
p.195～206
なおプログラムはpascalで書かれています

**デフォルトの名無しさん** · 2014/10/23(木) 23:13:40.41

>>653
できました！ありがとうございました

**ホモおじさん** · 2014/10/24(金) 00:54:53.39

>>655
お安い御用さ。
おじさんに何でも聞いて呉れ玉へ！

>>652
ぼうず。
積分の自動刻み幅調節なら森正武の本に載ってるのがあるぞ。
岩波の７７本にプログラムもある。

**デフォルトの名無しさん** · 2014/10/24(金) 01:31:35.96

>>654
ありがとうございますm(__)m

>>656
有益な情報ありがとうございますm(__)m

お叱りを受けるかも知れませんが、指数関数が入った数値計算で手を焼いております

具体的には
exp(x)/(exp(x)-1)
といった関数を考え, 1から20000くらいまでの数値を出力したいのですが, 桁落ちしてNan やら Inf やら出て困っております.

定石としてはx を何かに規格化するのかな？とは思うのですが, x の範囲が1から20000と幅広くなかなか思いつきません.
何か妙案はございますでしょうか？

**ホモおじさん** · 2014/10/24(金) 03:18:53.64

>>657
xが大きい方で計算したいなら通分して 1/(1-exp(-x)) の形で計算するのが定石だが、
どっちにしろ x=20～30くらいでもう 1/1=1 でその先は無理に計算しなくてもいい。

この辺は数値計算術にいっぱい書いてある。
森口繁一の本は古いけど結構面白い。
数値計算術～計算数学夜話＜＜数値計算工学の感じで真面目になる。

それにつけても蔵王のお釜が白濁液まみれというニュースが興味深いよね。

**デフォルトの名無しさん** · 2014/10/28(火) 00:52:02.95

Fortran で realloc に相当することをやりたんだけど、どうしたらいいでしょうぁ。

**デフォルトの名無しさん** · 2014/10/28(火) 03:20:58.44

>>659
宣言でallocatable
使う前にallocateで割り当て
大きさ変える前にdeallocateしてallocateで再割り当て
中身を確実に保存したかったらdeallocate前に別の変数にコピー
これで一応配列の大きさを変えれる

**デフォルトの名無しさん** · 2014/10/28(火) 03:55:09.13

Fortran2003なら　move_alloc を使う。
alloctae(tmp(new_size))
new_size(1:size(array)) = array
call move_alloc(tmp, array)
move_alloc は番地（ポインタ）を写すだけなのでメモリー内容のコピーは１回で済む。

modern fortran explained 15.5.3 Transferring an allocation 参照

**デフォルトの名無しさん** · 2014/10/29(水) 00:46:52.02

ちょっとスレチかもしれませんが…
無限級数の和が1になるものうち、\Sigma^{\infty}_{0} (1/2)^n よりも収束速度がゆっくりなものって何かありますか？

刻み幅に使いたいんですが中々思いつかないですm(__)m

**デフォルトの名無しさん** · 2014/10/29(水) 12:33:02.25

質問がおかしいので、たぶん刻み幅っていうより問題設定が間違っている。

**デフォルトの名無しさん** · 2014/10/30(木) 00:25:40.34

>>663
つまりは足し上げて1に収束するような数列を知りたいんです。

(1/2)+(1/2)^2+(1/2)^3+(1/2)^4+.....+(1/2)^n->1 ですよね？

これだと収束スピードが速すぎるので、もう少し緩めに収束する数列を知りたいんです

**デフォルトの名無しさん** · 2014/10/30(木) 22:03:37.99

intel fortran compilerが、非商用版があると聞いたんですけど
ホームページに行ってもありませんでした
もう非商用版は終了したんでしょうか

**デフォルトの名無しさん** · 2014/10/30(木) 23:06:25.70

>>665
ありますよ

**デフォルトの名無しさん** · 2014/10/30(木) 23:08:34.49

>>666
自分の探し不足ですかね
よければURLおしえていただけませんか？

**デフォルトの名無しさん** · 2014/10/30(木) 23:46:30.56

>>667
ttp://pen.agbi.tsukuba.ac.jp/~RStiger/hiki2/?Fortran+(Intel+Parallel+Studio+XE)+%A4%CE%A5%A4%A5%F3%A5%B9%A5%C8

この辺参考にしましたよ(^^)

**デフォルトの名無しさん** · 2014/10/31(金) 06:45:08.15

非商用ページにたどり着く方法, 頻繁に変わって不便だよなあ
ライセンス登録して.licファイルなりシリアルコードなり持ってるなら直接DLするって手もあるが

ttp://registrationcenter-download.intel.com/akdlm/irc_nas/4584/parallel_studio_xe_2015.tgz

**デフォルトの名無しさん** · 2014/10/31(金) 18:03:58.39

>>668
>>669
無事ダウンロードできました。ご親切にありがとうございました。

**デフォルトの名無しさん** · 2014/11/01(土) 22:28:08.66

amazon ec2で計算してる人いる？

**デフォルトの名無しさん** · 2014/11/02(日) 10:01:51.45

>669
御呈示のアドレスはLinux版ですが、Windows版
について情報をおもちではないでしょうか

**デフォルトの名無しさん** · 2014/11/02(日) 10:23:39.60

>>672
無償版ってこと？Linuxのみ無償なんだよ
よく読めよ、カス

**デフォルトの名無しさん** · 2014/11/02(日) 20:58:20.47

乞食か

**デフォルトの名無しさん** · 2014/11/06(木) 01:07:08.08

Linux版以外も学生非商用は無料だろ。

**デフォルトの名無しさん** · 2014/11/06(木) 07:51:33.96

>>675
嘘つくなよ、ゴミ屑

**デフォルトの名無しさん** · 2014/11/16(日) 23:31:47.83

2つの実測値の相関係数を出したいのですが、
fortranでできないでしょうか
実測値は単純な直線傾向ではないです
エクセルで手動でやることも考えたのですが、
実測値が１００ケースほどあって、
今後も増えることを考えるとプログラムで作りたいと思いました
もしよければアドバイスいただけないでしょうか

**デフォルトの名無しさん** · 2014/11/17(月) 09:46:48.96

>>677
Rを使用するのが一番手っ取り早いと思います。
Rを調べてみてあなたのやりたいことができそうなら、あえてプログラムを作る必要もないのでは。

**デフォルトの名無しさん** · 2014/11/17(月) 13:53:58.07

>>678
Octave でもいいんじゃないの？

**デフォルトの名無しさん** · 2014/11/17(月) 18:52:41.93

>>678,679
いろいろソフトがあるのですね。
調べて試してみようと思います。ありがとうございました。

**デフォルトの名無しさん** · 2014/12/13(土) 01:43:27.96

Fortranで書かれている遺伝的アルゴリズムのライブラリが
手に入るところはないでしょうか

**名無しさん＠そうだ選挙に行こう** · 2014/12/13(土) 18:32:39.00

ちっとはググれw

http://www.cuaerospace.com/Technology/GeneticAlgorithm.aspx

**デフォルトの名無しさん** · 2014/12/15(月) 00:30:38.37

MPICH または OPENMPI で Fortran2008 bindings を使う方法を教えろ下さい。
ドキュメントにはサポートと書いてあるが、インストールしても MPI_F08.MOD が
生成されない＠Ubuntu

どうか助けろ下さい。こっちのスレへ行けという指示でもおｋ

**デフォルトの名無しさん** · 2014/12/15(月) 19:55:14.30

ランダム生成のコードのについて質問です。
ifortではコンパイルできたのですが、gfortranではエラーが出ました。
以下がエラーメッセージです。

""$ gfortran routine.f90
routine.f90:32.25:

call random_seed( put = (/c * 999/) )
Error: Size of 'put' argument of 'random_seed' intrinsic at (1) too small (1/12)""

putの大きさが本質的に小さすぎる、や（=c　とした時）本質的に配列である。といったメッセージだと思い、いじってみたのですが解決方法が分かりません。よろしくお願いします。

あと / ~ / のスラッシュの意味もよろしければお願いします。

**デフォルトの名無しさん** · 2014/12/15(月) 20:47:21.21

試しにgfortran4.9でrandom_seed(size=i)（iはinteger）を実行したらi=12だった
つまりgfortranの場合はputの引数が大きさ12以上の1次元整数配列でなければならない
(/c*999/)は大きさ1の1次元配列で、その値がc*999という意味だからエラーになる

**デフォルトの名無しさん** · 2014/12/15(月) 20:51:20.48

>>682
ありがとうございます
ググれと言いながら教えていただいた優しさに感謝します
使わせていただきます

**デフォルトの名無しさん** · 2014/12/15(月) 22:17:44.08

>>685
丁寧な説明ありがとうございます。

cを配列にすると上あるシステムクロックのcountの取得のところでエラーが出てと堂々めぐりをしてましたが、シードが12の配列でも要素一つ変えれば乱数にはなるということなんですね。

スラッシュが配列を表しているとは気づきませんでした。勉強になります。

**デフォルトの名無しさん** · 2014/12/17(水) 02:40:01.63

LU分解のサブルーチンが作れません(ToT)
サンプルがあったらください…

**デフォルトの名無しさん** · 2014/12/17(水) 02:41:44.43

>>688
検索して、それっぽいのが出てきたんですがメインプログラムに噛ませてもうまくいきませんでした…

**デフォルトの名無しさん** · 2014/12/18(木) 08:43:35.24

>>689
なら氏ね

**デフォルトの名無しさん** · 2014/12/18(木) 14:54:41.87

MacOSX10.10.1, gfortran4.9にてコンパイル後、実行時に発生する下記のエラーの解決策はあるのでしょうか？

At line 68 of file test.f90 (unit = 1, file = 'parameter.txt')
Fortran runtime error: Cannot match namelist object name 44

namelistの不一致は見当たらず、Windows7, g95では正常に動作しました。

**デフォルトの名無しさん** · 2014/12/18(木) 23:40:19.37

計算流体力学を習い初めの学生なんですが、Radeonでfortranの計算をする場合にはどのソフトをインストールすればいいのでしょうか？
MinGW Command PromptとTDM-GCC-64というのを使っています。
出来ればインストール方法が書かれたページも一緒に知りたいです。

**デフォルトの名無しさん** · 2014/12/19(金) 20:10:17.60

open文で開いたcsvファイルをread(51,*)～～などで
読み込む操作が何回かある場合、
並列化はできたりしないでしょうか
読み書きはやっぱり並列化は無理でしょうか

**デフォルトの名無しさん** · 2014/12/20(土) 11:31:57.54

『if … else if … end if』と『select case …』はどっちが速いのか，ふと思って計測してみた．
cygwin gcc 4.8.3@core-i5で繰り返し計測してみたけど，時間差はばらつきの範囲内で有意な差はなさそうな感じ．
二つの構文ともCPU時間同じ，機能も同じで，いまいち，使い分け方が分からないんだよねぇ．．．

使ったコードはこれ↓
module mod_util
implicit none
integer(kind=4) :: t1, t2, t_rate, t_max, t_diff
contains
real(8) function elapsedTime()
implicit none
call SYSTEM_CLOCK(t2, t_rate, t_max)
if ( t2 < t1 ) then
t_diff = t_max - t1 + t2
else
t_diff = t2 - t1
endif
elapsedTime = t_diff/real(t_rate)
end function elapsedTime
end module mod_util

**694** · 2014/12/20(土) 11:33:00.40

program testSelectIf
use mod_util, only : t1, elapsedTime
implicit none
real (kind=8) :: randX
integer(kind=4) :: i, j, randI
#ifdef ifelse
write(*,'("# if ... else if ... ")')
#else
write(*,'("# select case")')
#endif
call init_randomSeed(972384)
call SYSTEM_CLOCK(t1)
do i = 1, 100000000
call random_number(randX)
randI = int( randX * 10, kind=4 )
#ifdef ifelse
if ( randI == 1 ) then; j = j + 1
else if ( randI == 2 ) then; j = j - 2
else if ( randI == 3 ) then; j = j + 3
else if ( randI == 4 ) then; j = j - 4
else if ( randI == 5 ) then; j = j + 5
else if ( randI == 6 ) then; j = j - 6
else if ( randI == 7 ) then; j = j + 7
else if ( randI == 8 ) then; j = j - 8
else if ( randI == 9 ) then; j = j + 9
else if ( randI ==10 ) then; j = j - 10

**694** · 2014/12/20(土) 11:33:45.85

end if
#else
select case (randI)
case ( 1 ); j = j + 1
case ( 2 ); j = j - 2
case ( 3 ); j = j + 3
case ( 4 ); j = j - 4
case ( 5 ); j = j + 5
case ( 6 ); j = j - 6
case ( 7 ); j = j + 7
case ( 8 ); j = j - 8
case ( 9 ); j = j + 9
case ( 10 ); j = j - 10
end select
#endif
end do
write(*,'("j=",I0)') j
write(*,'("#time:", F10.3)') elapsedTime()
end program testSelectIf

**デフォルトの名無しさん** · 2014/12/21(日) 02:39:49.62

>>693
I/O中に計算とかしたいという意味での並列化なら、F2003でasynchronous I/Oが出来る
ようになった。同じファイルの別々の位置のデータを並列に読み出したいという意味なら
そもそもハードウェア的に無理だろう。

>>694
最適化の問題というより、計算型GOTOの後継機能といったものではないか？
構造化プログラミングがホットだった時代は、GOTO 文を追放するための基本構造
のひとつとされていた気がする。

CASE文の中に入るのが、整数型、文字型、論理型のように離散値をとる場合に
限られるので、IF..ELSE IFの連鎖との使い分けは、比較で＝が意味を持つ型か
或いは実数型などのように＝があったとしても基本的に＜、＞的な大小比較しか
意味を持たない場合を、意識して明示的に区別していますよというアピールではないか。

**694** · 2014/12/22(月) 19:17:39.41

>>697

明快なお答えありがとうございます．
これまでうっすらと使い分けルールの必要性を感じながらも，
「きっとこれで良いよね？良いんだよね？」と恐る恐る使い分けていました．少し安心しました．

**デフォルトの名無しさん** · 2014/12/22(月) 21:51:51.17

>>697
ありがとうございます
やりたかったことはcsvファイルがいくつかあり、
それを順番に読み込みreadするdoループを並列化できないかなということでした

教えていただいた後者の方にあたりそうな気がしまして、
勉強になりました。

**デフォルトの名無しさん** · 2014/12/23(火) 16:29:15.14

永らくFORTRAN77で飯を食ってたのですが,より新しいFORTRANへ移行してみたくなりました
なにか良書はありましょうか?

職業プログラマで社内で請負で流体や構造体のメッシャーを組んでます
まあ,ベースのプログラムも77,新規も77縛りの事が多いので特に困ってはないのです
しかし,昔,汎用機世代の6しか書けないような先輩方を古臭いと笑っていたのですが自分もそんな年になってしまい
後輩らに笑われたくないので,ちょっと勉強しておきたいと考えています

世の90の初心本(それ以降は書店でも見ないなー)等は基本的な制御構文等が厚く書かれており
ホントに身に着けたい新たな概念のintent等がおざなりになってしまっています
古臭い77書きに光を......
#C++もやってるのでクラスの概念なんかも問題はありません

**デフォルトの名無しさん** · 2014/12/23(火) 16:43:41.88

日本語の書籍に関しては仰る通りで、特に2003以降は全くないはず
まずはnagのサイトにあるFortran入門とFortran 2003入門が無難かと
ただ、これらは規格を網羅してるわけではないので
完全な仕様を知りたい場合はコンパイラのマニュアルや
JIS規格(JIS X 3001-1:2009)も参照する必要あり

**デフォルトの名無しさん** · 2014/12/23(火) 20:48:23.28

>694

selectは単一の変数でしか使えないよ
selectはどうかするとジャンプテーブルで一発の比較で処理してくれる可能性が高い
最近のコンパイラはelseifでもそうするとか聞くけど,どうだろうね?
elseifは良く引っかかるのを前に持って来て早めに引っ掛ける事ができるとかも聞くねー

個人的には,読みやすさを優先するなー

#Cのswitchの議論と同じなので探してみると良いよ

**700** · 2014/12/23(火) 20:51:07.24

ありがとうございます
まずはサイトのそれらを読んでみます

その後は......JISかー,読みやすさを優先してくれないからなー
コンパイラのマニュアルは殆ど英語になるのが悲しい