くだすれ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/06/16(月) 10:00:44.78

cuda fortranについて初歩的な質問なのですが、
PGI workstationで、-cublasとオプションを付けてコンパイルしようとすると、
『fatal error LNK1104: cannnot open file 'libcublas.lib'』とエラーが出てしまいます。
この場合、どうしたら良いですかね?
どなたか教えてくださると助かります。

**デフォルトの名無しさん** · 2014/06/16(月) 10:13:18.23

>>477
パスが通ってないんだろ

そのエラーにあるファイルがあるか探す
ファイルがあるパスをLオプションで追加でとりあえず通るはず
例えば-L/usr/local/cuda/lib -cublas

**デフォルトの名無しさん** · 2014/06/16(月) 15:38:44.47

>>478
探して見たところ、'libcublas.ipl'、'libcublasemu.lib'という名前のファイルしか見つかりませでした。
そこで、これらのファイルが保存されているフォルダのパスを追加してコンパイルしたところ、
『error lnk2019: unresolved external symbol referenced in function』
というようなエラーが出てきてしまったのですが、
どうすれば良いか教えていただけませんか？

**デフォルトの名無しさん** · 2014/06/16(月) 15:50:37.89

>>479
試しにcublasをcublasemuに変えてみて
libcublas.libがなかったのならもしかしてパスはすでに通っていたかもしれません

**デフォルトの名無しさん** · 2014/06/16(月) 16:27:09.31

コード中の'cublas'を'cublasemu'に変更し、
'-lcublas'を'-lcublasemu'に変更してコンパイルしたところ、
『unable to open MODULE file cublasemu.mod
compilation aborted』
と出てしまいました。。。

**デフォルトの名無しさん** · 2014/06/16(月) 16:43:59.97

>>481
コードは元のままでいいです
申し訳ない

**デフォルトの名無しさん** · 2014/06/16(月) 16:54:43.67

>>482
元のコードでやってもやはり出てしまいますね。

**デフォルトの名無しさん** · 2014/06/16(月) 17:46:47.99

>>483
use cublasemuになっていませんか？

**デフォルトの名無しさん** · 2014/06/17(火) 09:43:20.11

>>484
いえ、ちゃんとuse cublasってなっています。

**デフォルトの名無しさん** · 2014/06/17(火) 13:19:25.73

use cudafor
use cublas
が指定されている

-Mcudaまたは-Mcuda:4.0あるいは-Mcuda=emu

unresolved external symbolと言われたら
libが足りないので-lで追加
*.libがないと言われたら
-Lなどでパス追加
*.modがないと言われたら
-Iでパス追加

参考
http://www.softek.co.jp/SPG/Pgi/TIPS/opt_cudaF.html
http://www.softek.co.jp/SPG/Pgi/TIPS/public/accel/cublas40.html

**デフォルトの名無しさん** · 2014/06/17(火) 16:58:12.46

>>486
-Lパスでディレクトリを指定して、-cublasemuをつけてコンパイルしたところ、
"libcublasemu.lib : error lnk2019: unresolved external symbol referenced in function"
が、ものすごい数のが出てきたのですが、どうしたらいいですかね？
一個ずつ-lで追加していくしか無いのでしょうか。

**デフォルトの名無しさん** · 2014/06/17(火) 17:24:37.11

cublas使わなければ普通に動くの？

**デフォルトの名無しさん** · 2014/06/18(水) 11:22:33.36

>>488
はい。cublasが関わるところを全てコメントアウトしたところ、動きました。

**デフォルトの名無しさん** · 2014/06/18(水) 18:11:03.14

>>489
可能な限り少しずつコメントアウトを外していき
その都度コンパイルが通るようにオプションを調整してみてください
最初はuse cublasだけなど
大量のエラーが少しずつ出ることを期待します

cublasは呼び出し方法が旧いものや新しいもので3種あるようですが
呼び出し方と CUDA toolkit のバージョン指定はあっていますか？

pgi cuda Fortranはデバッグ用にGPUを使わないエミュレートモードが用意されています
コンパイルオプションとリンクするライブラリの両方で切替えるようです
どちらを使用していますか？
cublasemuはエミュレートモード用で
オプションに-Mcuda=emuを併せて指定します

**デフォルトの名無しさん** · 2014/06/20(金) 16:00:17.10

Ｃのatoiとかatofのように，文字列としての数字の列を
integerなりrealの値に変換してくれる函数はあるでしょうか。
当方は，MinGWのgfortran95を使っています。

**デフォルトの名無しさん** · 2014/06/20(金) 19:25:02.85

>>491
内部ファイルを使うのがセオリー
NAG Fortran入門 19.5 数値と文字の変換を行う方法
ttp://www.nag-j.co.jp/fortran/FI_18.html#ConvertingBetweenNumbersAndStrings

**デフォルトの名無しさん** · 2014/06/20(金) 21:53:29.14

>492
ありがとうございます。
こういう裏ワザがあったのですね。
さっそく使わせていただきます。

**デフォルトの名無しさん** · 2014/06/21(土) 17:36:05.91

裏ワザというかFORTRAN77の頃からあったし、
66時代にはENCODE/DECODE関数が用意されていたんだがな。

マニュアルをよく読めば、内部ファイルの項目が昔から書いてある。

**デフォルトの名無しさん** · 2014/06/21(土) 20:05:31.04

>464
そうでしたか。なにしろわたしが FORTRANを学生時代に始めた
ころはFORTRAN IVのころで80桁のパンチカードにパンチした
カードデックを窓口の届けると2日程で，棚にカードとラインプリ
ンタで打ち出した結果がおいてあるという時代(東大の計算センター
にあったHITAC　5020）のころで，FORTRAN77の出る前
ダイクストラの構造化プログラミングが騒がれ始めて，
カーニハンとリッチーの『ソフトウエア作法』をよんで，
構造化出来るRATFORを大学にあった小さなコンピュータで
実験しているころからの進歩を知りませんでした。
その後，CおよびC++に興味が移りFROTRANをご無沙汰してい
たのですが，最近のFORTRANで128bit処理ができることを
知り，まいもどってきた次第です。
市販のfortran90，95の入門書を見てみたのですが，見当ら
なかったので，質問したしだいです。

**デフォルトの名無しさん** · 2014/06/21(土) 22:49:19.55

はろけくも来つるものかなフォートラン

あっ季語がない

**デフォルトの名無しさん** · 2014/06/24(火) 20:49:58.43

Gfortran 使ってるんだが、データを読み込ませる際、2行ごとにデータを読み込ませることって可能？

**デフォルトの名無しさん** · 2014/06/25(水) 01:47:13.51

素直に全部読みなさい
変数に入った値は、必ず使えとも使うなとも言われてないだろう

**デフォルトの名無しさん** · 2014/06/26(木) 00:45:30.86

まぁ空のread文で１行飛ばせるけどな。
Fortmatの/でも行けたっけかな？曖昧みー

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

awkなりsedなりで一行飛ばしの中間ファイル作って、それをデータファイルにするとかじゃだめなの？
入力関係でつまづくなら、fortranにあわせたファイルを作成する方が早いのではと思ってしまう。

バイナリデータだとテキストみたいにいかないのか。バイナリデータだったらごめん。

**デフォルトの名無しさん** · 2014/06/26(木) 08:16:37.95

>>498->>500
レスどもです。直接ｎ行ごとに読み込むという作業はできないみたいですね。
私も中間ファイルを作って…と思いましたが、読み込みデータをいじってはいけないらしく断念しました。
解決案は以下の通り、(n行ごとに読みたい場合)
If(mod(i,n) .eq. 0) then
Read(20,*) hoge(i)
Else
Read(20, *)
Endif

**デフォルトの名無しさん** · 2014/06/26(木) 21:08:20.99

詰めないと使いにくくない？

i=1
do m=1, linenum
if(mod(m,n) .eq. 1) then
read(20,*) hoge(i)
i = i +1
else
read(20,*)
endif
enddo

**デフォルトの名無しさん** · 2014/06/27(金) 00:10:25.56

>>499の言うとおり、formatを'(i,5/)'とかにすれば良いのでは？

**デフォルトの名無しさん** · 2014/06/27(金) 00:38:41.62

>>502
あっ、もしかして>>501のやり方って途中空白行入ってることになるのかな？？？

>>503
それについて少し調べる必要があるかも。
具体的にどゆこと？

**デフォルトの名無しさん** · 2014/06/27(金) 00:42:32.62

formatのスラッシュでskipできるぞ。

program test
implicit none
character(len = 132) :: text
open(10, file = 'test.f90')
1 read(10, '(a, /)', end = 999) text
print *, trim(text)
goto 1
999 stop
end program test

**デフォルトの名無しさん** · 2014/06/27(金) 01:05:20.13

>>505
おお！スラがスキップすると言う指示内容なんですね。

**デフォルトの名無しさん** · 2014/06/27(金) 02:12:24.88

>>506
んだ。
writeの時は出力１行改行、readのときも入力改行で結局１行読み飛ばすことになる。

**デフォルトの名無しさん** · 2014/06/27(金) 09:18:57.76

久しぶりにfortranのプログラム見たが、FORTRAN77の時代と違って気持ち悪い。
これならC++の方がまだいい。テンプレート関係を除けば。

でも、fortranは４倍精度計算ができるからなあ。
並列処理やGPGPUに対応したライブラリも出ているし。

まあ、年寄はfortranからMatlabに移行しているみたいだが。

**デフォルトの名無しさん** · 2014/06/28(土) 00:46:11.64

おじいちゃん。
あんま結論を急がない方がいい。
配列演算気持ちいい。PUREでELEMENTALなら最高。

４倍精度は精度は4倍だが計算時間は１０倍超だから、
１０倍の時間の合間にアルゴリズム見直した方が吉。

**デフォルトの名無しさん** · 2014/06/28(土) 02:59:17.00

gccにもquadmathがあるんだからあまりFORTRANを高く評価しない方がいいよ

**デフォルトの名無しさん** · 2014/06/28(土) 13:41:05.28

並列計算させてdoループの中でwrite文でファイル書き込みしたら中身の順番めちゃくちゃで泣きました。
どーにもなんないんですかね？

**デフォルトの名無しさん** · 2014/06/28(土) 15:42:05.65

>>510
そもそもlibquadmathはFortranのために追加されたんだが

**デフォルトの名無しさん** · 2014/06/28(土) 16:56:46.35

>>511
write文だけ逐次処理させるか、ループの外に出すしかないと思う。

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

>>509

はぁ｡｡｡

配列演算気持ちいい、って何？　それ。
C**なら行列演算クラスライブラリ作って、演算子オーバーロード使えば終わりじゃん。

まさか、行列積が「・」でできると言うんじゃないだろな？

>４倍精度は精度は4倍だが計算時間は１０倍超だから、
>１０倍の時間の合間にアルゴリズム見直した方が吉。

何言ってんだ？この馬鹿。
物理計算によっては４倍精度の計算が必要なの。計算時間やアルゴリズム
の問題じゃない。人の言ってることをそのまま受け売りするな。阿呆

**デフォルトの名無しさん** · 2014/06/28(土) 22:07:05.28

>>514

C** → C++

**デフォルトの名無しさん** · 2014/06/29(日) 00:26:55.37

懐かしい単語を見かけたからって
いい大人が質問スレに何しに来てんだろね

それとも>>508は質問なのか？

**デフォルトの名無しさん** · 2014/06/29(日) 10:07:26.16

libquadmathの中を覗いてみたら、x86版はFPU使ってるけどx64版はSSE2使ってるな

**デフォルトの名無しさん** · 2014/06/29(日) 12:17:10.06

ポインタが表に出ない（もちろん POINTER は除く）、
配列長と次元が暗黙的に持てることは評価できると思うんよ
matlab/octave に勝るかというと微妙だけど（移植性くらい？

>>511
ソースコード> このループは順番ぐちゃっても大丈夫です
処理系> なら並列実行しても同じだよね！
という発想なので、そういう副作用も考慮しなきゃダメ

>>517
確か x86_64 は SSE の存在を仮定していいはず
まぁ x86 も SSE/AVX も quad precision floating-point の直接演算は無い

**デフォルトの名無しさん** · 2014/06/29(日) 12:46:18.15

> C++なら行列演算クラスライブラリ作って、演算子オーバーロード使えば終わりじゃん。
> まさか、行列積が「・」でできると言うんじゃないだろな？

冗長だけど、～すれば終わり、程度には書けますぜセンセー。

function dotp_r8(A, x) result(y)
double precision, intent(in) :: A(:,:), x(:,:)
double precision :: y( size(A,1), size(A,2) )

y = matmul(A, x)
end function
!----====----====----====----====----====----====----====
program test

interface operator(.dot.)
function dotp_r8(A, x) result(y)
double precision, intent(in) :: A(:,:), x(:,:)
double precision :: y( size(A,1), size(A,2) )
end function
end interface
double precision :: A(2,2), x(2,2), y(2,2)

A = reshape( (/1.D0, 1.D0, 1.D0, 1.D0/), (/2,2/) )
x = reshape( (/1.D0, 3.D0, 2.D0, 4.D0/), (/2,2/) )

y = A .dot. x

write (*, '(2F8.3)'), transpose(A); print *
write (*, '(2F8.3)'), transpose(x); print *
write (*, '(2F8.3)'), transpose(y)
end

**デフォルトの名無しさん** · 2014/06/29(日) 13:09:54.23

クラスライブラリで演算子オーバーロードって・・・遅そう。
素人が書いたマクロじゃ最適化効かないでしょ。
ソースの見かけで書けたって意味ないの。

大体、高速行列演算ライブラリを自分で作れるくらいなら、四倍精度のライブラリ
も自作すればいい。昔の爺さん世代はよく倍精度・四倍精度自作ライブラリを自慢していたｗ

まぁC++のBlitzなんて１０年以上前から夢を語ってノロノロ開発して
ようやく最近なんですけどね。

**デフォルトの名無しさん** · 2014/06/29(日) 14:46:16.10

>>513 >>518
なるほど、順番がどうなっても良いような時にしか使っちゃいけないポイですね。

どういった場面で並列化使えるのかな…
自分はある物理の数値計算に使ってるんだけども、順番が変わったらダメ。

それにしてもここにいる住人は凄い人が多くて、泣ける

**デフォルトの名無しさん** · 2014/06/29(日) 15:59:56.49

並列化は普通に物理の数値計算につかうけど、、、

**デフォルトの名無しさん** · 2014/06/29(日) 17:26:15.54

>>519

あ、そう。matlabみたいに糞な「：」が好きなのね。もういいから、死んでね
お爺ちゃんからこんな言葉を言われたお前はすでに終わっている。

>>520

その通り。面倒くさかったから書かなかったけど。
とにかく、一時オブジェクトの生成と.消滅、データのコピーが頻繁におきる。

それを避けるために右辺値参照が使えるかどうか。でも行列乗算や線形連立代数方程式
を解く場合には使えないかも。まだ、やってみてないからわからんが。

作業用行列のポインタ値を管理するグローバルオブジェクトを使うと
上の問題は回避できるが、残念ながらマルチスレッドには対応できてない。

これをどのようにグローバルスレッド対応にしたらいいか、戦略が
立たない状態です。

**デフォルトの名無しさん** · 2014/06/29(日) 17:29:17.08

>>519

しかし、ほんとコードが糞だよな。こんなの覚えるくらいなら
C++の方がましだわ

**デフォルトの名無しさん** · 2014/06/29(日) 17:34:36.24

>>522
いや、そうでしょうけど、出力の順番とか気にする場合が私の研究分野だとほとんどなので…
なにか勉強になるものがあればいいのですが

**デフォルトの名無しさん** · 2014/06/29(日) 17:57:50.26

David BaileyがFortran用とC++用の4倍、8倍、多倍精度ルーチン公開してた気がする。
http://crd-legacy.lbl.gov/~dhbailey/

**デフォルトの名無しさん** · 2014/06/29(日) 17:59:34.52

間違えた

グローバルスレッド対応→スレッド対応

**デフォルトの名無しさん** · 2014/06/29(日) 21:18:38.95

>>525
プロセス(OpenMPならスレッド)ごとに別ファイルに出力して、
後で適当なスクリプトでマージするとか。

でも普通はそんなことしないで、>>513 が書いてるように、
doループを計算部と出力部の2つのループに分けて、
計算部のみ並列化して出力部分は逐次処理だと思う。

**デフォルトの名無しさん** · 2014/06/29(日) 21:46:47.56

はいはいソーデスネーでっかい一時領域作ってなにしてんでしょーねorz

pure で elemental とゆーと配列を引数に放り込むとその全要素に
関数定義された処理が適用される、各処理は独立なのでまとめるなり
並列化なりをコンパイラーに期待、かな。演算子全然関係ないしorz

過去のシサンガーとかもう辛いだけ。
どこまでも未来に生きてください、Fortran スレになんか戻らないでください。

>>521
出力の順序だけ重要ならバッファして後から並べ替えるとか
浮動小数点演算的に結果を変えられないなら仕方ないかもだが

だって人口少ないし交流少ないし機会もないし意味すらないかもしれない、
たまに発言すれば馬鹿にされるから凄い人しか残らないわけで。

**デフォルトの名無しさん** · 2014/06/30(月) 00:33:26.64

>>518
それは分かる
でもエミュレーションしたら1/10程度じゃ済まないだろうから一応ハードウェアが
有効に働いてるんだな
double-double演算でもないみたいだし

>>526
面白そうだね

4倍精度が必要になるのは例えば反復法などの収束を加速する時などで、
単なる倍精度だと桁落ちのために誤差が生じてなかなか収束しない時に
有効だったりするね

**デフォルトの名無しさん** · 2014/06/30(月) 02:26:19.67

Fortranの強みは常に時代の最先端ハードウェアに対応していることじゃないかな。

最近でもIntelはDECのFortran買ったし、NvidiaはPGIを買収したし。
両社ともGPUアクセレレータを含む自分らの最新ハードとセットで、
同じタイミングで対応したコンパイラを出してきている。

本来ハードウェアとソフトウェアは表裏一体、精神と肉体のような
不可分の関係にあるはず。そしてFortranは未だその密接さを保っている。

ソフトウェアの人たちは、ソフトウェア単体での独立価値を信じたいようだけど
現実は違っている。

この現実認識の歪が過剰なFortranディスカウントキャンペーンとして噴出して
いるのではないだろうか・・・と新聞並に論理を飛躍させて終わるｗ

**デフォルトの名無しさん** · 2014/06/30(月) 13:21:06.69

>>529

そういうのに向かない計算アルゴリズム(隣のセルの計算結果を使う。つまり独立ではない。
そうなると同期の問題も発生する)もあるんだ。

>過去のシサンガーとかもう辛いだけ。
>どこまでも未来に生きてください、Fortran スレになんか戻らないでください。

で、お前何が言いたいんや？日本語大丈夫か？

**デフォルトの名無しさん** · 2014/07/01(火) 01:36:55.46

1 ：デフォルトの名無しさん：2012/08/16(木) 20:58:00.31
このスレッドは、他のスレッドでは書き込めない超低レベル、
もしくは質問者自身何が何だか分からない質問を勇気を持って書き込むスレッドです。
FORTRAN使いが優しくコメントを返しますが、
お礼はFORTRANの布教と初心者の救済と次期Fortran2008規格でのCOMEFROM文採用をお願いします。

**デフォルトの名無しさん** · 2014/07/03(木) 20:06:28.12

つかぬことをお聞きしますが
皆さんのfortranの本の勉強履歴はどんな流れですか？
私は冨田さんfortran90の本→辞書的にfortran90入門をかいつまむぐらいしか
やっておらずとりあえず拙いなりに読み書きはするのですが
なにか良い本あればおしえていただけないでしょうか

**デフォルトの名無しさん** · 2014/07/04(金) 01:38:26.65

故岡村先生の授業テキスト
やーっふぉ、やっふぉ、ふぉーとらんらんらん

**デフォルトの名無しさん** · 2014/07/04(金) 01:45:11.81

Fortran95/2003 Explained --> Modern Fortran Explained

**デフォルトの名無しさん** · 2014/07/05(土) 07:04:30.09

emacsでコメントアウトに使う関数fortran-comment-regionのデフォルトの動作が
行頭にCを挿入するのではなく,C$$（自由形式では!!$）を挿入するのは何か理由があるのでしょうか？
$つけたほうがメリットとかあるの？

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

absoft 9.0でf77fileをコンパイル成功、プログラム自体も上手く行く
absoft11.0でコンパイル成功、プログラム自体もインプットファイルのある部分の記述次第で上手く行く
ある記述部というのはアウトプットするデータ個数を指定するところです。２０個程度なら上手くいくのですが数百個と増やすとうまくいきません。
エラーメッセージは
illegal character in numeric input
で、インプットファイルにおかしい部分があると言っています。しかし9.0でコンパイルしたものはきちんとうまくいくのです
何かアドバイスいただけないでしょうか。プログラム自体は３万行程度あるというのと内容だけに書き写すのは難しいです。

**デフォルトの名無しさん** · 2014/07/07(月) 01:25:39.47

>>538
よく分からないけど、正常に動くなら9.0のままじゃ駄目なの？
出力するときフォーマットは指定してる？
してるなら*で指定なしにしてみたら

**デフォルトの名無しさん** · 2014/07/07(月) 09:24:44.41

>>539
formatを指定しています。そこを当たってみる事にします
9.0のままで良いのですが、引き継ぐものなので後の人のために11.0で動くよう乗り越えておきたいんです。

**デフォルトの名無しさん** · 2014/07/08(火) 00:03:54.23

Gfortran で定義してない変数があってもたまにエラーはかない時あって焦った。
デバッグで小一時間汗流したわ。
gfortranってたまにアホな子になるの？

**デフォルトの名無しさん** · 2014/07/08(火) 23:44:01.37

>>535
>>536
ありがとうございます
まだ私が理解できるレベルにないかもですがいつか読んでみたいと思います

**デフォルトの名無しさん** · 2014/07/09(水) 03:16:15.11

>>538
出ているエラーメッセージは、入力ファイルの問題で、数字の代わりに
（たぶん)文字が来ているということ。

F77時代（かそれ以前）のプログラムの用だから、入力ファイルの数値の
フォーマットも指定されていると推測される。数値が指定カラムから
はみ出していないか見て味噌。

以前
　　300 400 200 0 0
これで右から２番目でおかしくなると言っていたが、たぶん４カラムづつ
入力にとってある。300 400 200 000 000 ならたぶん動くだろう。
今は 0 0 を読んで、数値Formatがおかしいと文句を言われているのだと思う。

**デフォルトの名無しさん** · 2014/07/09(水) 03:17:43.97

ではなぜ Absoft9 ではよくて、11では駄目なのかと問うかもしれない。
考えられることとしては、入力データの空白の扱いのデフォルトが変わった
ことが考えられる。FORTRAN最初期からの伝統で、入力データの空白は
数値としてはゼロを意味している。文法書のFORMAT分のところをよく読むと、
BZとBNというオプションがある。BZはBLANK ZEROでBNはBLANK NULL(NONE)
を意味していて、解釈を変えられる。昔のデフォルトはBZだった。
しかしF95かF2003くらいでたしかデフォルトをBNにしていた気がする。
今、手元に資料がないので確認できないがｗ
というわけで、absoft9は空白を０として読んで、0 0 を000と解釈したが、
Absoft11は0 0を数値として解釈不能とみなして、illegal character in numeric input
の実行時エラーをはいているのだろう。
またそうだとするとAbsoft9では０入力の実行回数が１回少ないはずである。

ちなみに、なぜ空白がゼロなのか？これは昔は、紙テープ、パンチカードの
穴あき不良や読み取り不良が多くて、パンチカードのスペース（穴なし）
が出やすかったからと聞いたことがある。ホントかは知らない。

あとメール欄にsage と書くように。
スレが上がると、冷やかしが紛れ込んできてスレが荒れる原因になる。

**デフォルトの名無しさん** · 2014/07/09(水) 04:42:58.29

>>538
>>543,544 に加えてOPEN文のBLANK=も確認するといいかも。
手元の言語仕様のドラフト板ではFORTRAN77/90/95/2003共に
NULLがデフォルトになっているけれど。

それとabsoft9.0と11.0とでコンパイルオプションが違うというオチはない？
コンパイラは上記の様なデフォルトの解釈を変更するオプションを持っていたりするよ。
そのオプションが使われていた場合、ソースをいくら読んでも原因はそこにはないから注意ね。

**デフォルトの名無しさん** · 2014/07/10(木) 18:28:18.82

FORTRAN77のCOMMON文とEQUIVALENCE文について質問です.
!---
DIMENSION A(4)

EQUIVALENCE (A, V1)

COMMON /VS/ V1, V2, V3, V4

のようなコードをfortran90風に書き直したいのですがAをCOMMONブロックの変数で置換する以外にいい方法は無いでしょうか？
できればSubroutine等を増やすのではなく,Moduleの変数宣言部分だけで解決したいのですが...

**デフォルトの名無しさん** · 2014/07/10(木) 20:24:35.98

>> 546
↓使えるかな
http://www.nag-j.co.jp/fortran/FI_16.html
16.1.2 モジュール内の名前を別の名前で参照する方法

**デフォルトの名無しさん** · 2014/07/10(木) 22:05:54.56

>>547
ありがとうございます
早速以下のように使ってみたんですが
module VS
real :: V1, V2, V3, V4
end module VS
program
use VS, A => V1
...
A(1) = 1.0d0
...
end program

↑コードのA(1)に代入する部分で
error #6410: この名前は配列または関数として宣言されていません。
とgfortranのコンパイルエラーになります
dimension A(4)などと変数宣言すると
error #6401: この名前の属性が USE 文によってアクセス可能な属性と矛盾しています。
でコンパイルエラーになります
あと
use VS, A(1) => V1
のような書き方もだめでした
nag-j.co.jpのコード例だとAは変数とは扱われていないように思うのですがAを配列として扱うことはできないでしょうか？

**デフォルトの名無しさん** · 2014/07/10(木) 22:50:41.50

>>546
無理なんじゃないかな。
それはそれで最小にまとまっている。
最適化に問題があるとはいえ、無理してCOMMONやEQUIVALENCEを直す必要もない。

そのコードが意味しているのは、V1,V2,V3,V4を４つの実数変数として共有しているが、
これをエイリアスとして配列に結び付けて使いたいということになっている。

ところが、COMMONでないとV1,V2,V3,V4のメモリー上の番地が連続することを保証できない。
だから配列へのポインタを充てることができない。
（TYPEの中に入れれば、sequence属性で番地の連続性を保証できるが…
f2003のC言語との互換性機能も使えるかもしれないが、COMMONよりひどい気もする。）

ばらばらの番地いる変数を個々にFortranにはポインタの配列はない。
（配列へのポインタはある。これを使うにはV1~V4が連続した番地にいる必要がある。）
実変数へのポインタを要素として持つ構造体を定義して、その配列を作れば
いいが、それでは大幅な書き換えで意味がない。

あとは、Fortranの規格を外れるが、どの処理系もほぼ必ず持っている
Cray Pointer を使う手もあるが、これもModuleの中では納まらない。
実行文が必要。

COMMONとかEQUIVALENCEはメモリーの番地共有だから、便利なことができる。
最適化を抑止してしまうから、主要部分には使わない方がいいが、
技巧的小技としては、うまく使うと短く書けるので、使えばいいんでね？

**549** · 2014/07/10(木) 22:53:27.62

修正ミスｗｗｗ

＞ばらばらの番地いる変数を個々にFortranにはポインタの配列はない。
ばらばらの番地いる変数を個々に、ポインタの配列の要素に結び付けようにも、
Fortranにはポインタの配列はない。

**デフォルトの名無しさん** · 2014/07/11(金) 03:55:49.85

ちょいと手間だがこれじゃダメかい

module VS
real :: A(4)
real :: V1, V2, V3, V4
equivalence (V1, A(1)), (V2, A(2)), (V3, A(3)), (V4, A(4))
end module VS

**デフォルトの名無しさん** · 2014/07/11(金) 04:36:13.83

>>549
大変詳細な解説ありがとうございます
てっきりEQUIVALENCEは廃止予定事項だと思い込んでましたが便利ですね

>>551の書き方でいこうと思います
equivacenle(A, V1)
でも -chack all -warn allでエラーは出ませんでしたが実行時にメモリ破壊起こるのかな

**デフォルトの名無しさん** · 2014/07/11(金) 09:33:00.25

>>543,>>544,>>545ありがとうございます。
入力データの書式を指定する部分がを見た事がないのでそこを探してみます。ありがとうございます。

**デフォルトの名無しさん** · 2014/07/11(金) 15:22:52.88

>552
equivacenle(A, V1)
A(1)==V1 にはなるが、V2,V3,V4がV1の直後に連続的に配列されていなければ
A(2),A(3),A(4)とV2,V3,V4は全くの無関係となってしまいます。
MAPファイルを書かせてみれば、位置関係は分かります。

最近のコンパイラはかなりソース上の順番とは無関係に自由に
並べ替えるので、うまくいかない可能性大かと。

EQUIVALENCEはDeprecated Listに上がっていました。
この先廃止の見込みが高いので避けた方がいいリストかな？
deprecated -> obsolescent -> deleted の三弾活用の１段目。

**デフォルトの名無しさん** · 2014/07/12(土) 03:13:21.15

５５５

**デフォルトの名無しさん** · 2014/07/12(土) 11:23:27.92

>>555
ダミアンもどき

**デフォルトの名無しさん** · 2014/07/13(日) 00:10:07.50

サブルーチンの内部ではうまく計算されているのに、
メインでwriteするとよくわからない大きな数値になってしまっているのですが、
なにか解決のヒントないでしょうか？
断片的で申し訳ありません

**デフォルトの名無しさん** · 2014/07/13(日) 00:31:08.73

すいません自己解決しました
配列宣言の量が少なかったボンミスです
allocatable使わないとこういうことになるんでしょうか
お騒がせしました

**デフォルトの名無しさん** · 2014/07/19(土) 14:16:39.37

2014年になってもFORTRANが選ばれる理由
ttp://developers.slashdot.jp/story/14/05/12/0410259/2014%E5%B9%B4%E3%81%AB%E3%81%AA%E3%81%A3%E3%81%A6%E3%82%82FORTRAN%E3%81%8C%E9%81%B8%E3%81%B0%E3%82%8C%E3%82%8B%E7%90%86%E7%94%B1

まだ残っていきそうな感じだね

**デフォルトの名無しさん** · 2014/07/20(日) 02:33:05.33

FORTRANで書かれたBLAS/LAPACKをライブラリとして呼び出しているLLが
FORTRANと同等のスピードが出ることを大々的に喧伝すると、
それに釣られて、もうFORTRANは要らないのだ！とか言い出す人が、
情強を誇るネットの有名人にも多いので草不可避ｗ

**デフォルトの名無しさん** · 2014/07/27(日) 19:42:35.57

>>553の問題は解決したのかいな？

**デフォルトの名無しさん** · 2014/08/04(月) 11:41:31.29

バイナリ出力をするときにOPEN文に指定する
form='binary'はFortranの標準ではなく、独自拡張みたいですが、
標準のコーディングだけで'binary'と完全に等価なバイナリ出力をする方法はありますか？

form='unformatted'を指定すると、write文を書くたびに
データの最初と最後に4バイトのレコード長がつき、
form='binary'を指定すると前後4バイトがつかないみたいです。

**デフォルトの名無しさん** · 2014/08/04(月) 17:53:28.34

Fortran2003以降ならストリームが使える
http://www.nag-j.co.jp/fortran/fortran2003/Fortran2003_7_1.html#AUTOTOC_7_1

**562** · 2014/08/04(月) 22:17:03.55

>>563
ありがとうございます。
試してみます。

**デフォルトの名無しさん** · 2014/08/05(火) 02:23:56.59

directファイルを使うテクニックもある。
確実ではないが、大抵の処理系では、バイトまたはワード単位でべた書き
しているので、binaryの代わりに使える。

バイトかワードかは処理系によるが、コンパイラのオプションでどっちにも
できることが多い。

まぁ今の時代は素直に stream だろうが。

**562** · 2014/08/07(木) 21:48:00.99

>>565
またまたありがとうございます。今日、会社でstreamを試してみたら、
コンパイラが古いためか（Windows版Intel Fortran 9.1）、
コンパイル時に「streamなんて知らないよ」って言われました。
directで試してみることにします。

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

構造型にSave属性つけられないんですけど
構造型をModuleに入れてグローバル変数として使いまわすのは危険（未定義になる可能性があるん）でしょうか

あとFortranの質問じゃないかもしれませんが
Intel MKLのBLAS/Lapack95はpure procedureであるとInterface定義に明記されてるけど
BLAS/Lapack77はPureなんでしょうか

**デフォルトの名無しさん** · 2014/08/13(水) 11:45:31.24

↑すみません多分自己解決しました
構造体の定義のあと変数宣言したものにはsave属性つけられますね（定義しかしてなかった）

BLAS77は間違った引数渡すと標準出力（か標準エラー出力）になんか表示されるのでpureじゃないってことでしょうかね

**デフォルトの名無しさん** · 2014/08/13(水) 12:55:08.27

>>568
そもそもFORTRAN77には言語仕様上pureは存在しない。
Lapack95はFortran95以降のプログラムからLAPACKのprocedureを
利用し易くするためのラッパーライブラリに過ぎない。

**デフォルトの名無しさん** · 2014/08/13(水) 22:38:10.66

（コンパイル済みの）外部サブルーチンが真にpureかどうかは開発者/提供者の責任かと
例外処理は pure 要件から外していいと判断し、
pure と明記しておくことで最適化を妨げないようにしたのだろう

**デフォルトの名無しさん** · 2014/08/14(木) 02:43:06.91

Fortranの組み込みサブルーチンで、乱数がサブルーチンで関数になっていないから
他言語に比べて使いにくいなと思っていたが、FORTRAN は PURE じゃないものは
関数にしない方針で一貫していて、乱数はPURE（参照透明）じゃないから
関数にしなかったんだと分かって、なんか納得した。時刻サブルーチンとかも。

最近関数型、関数型うるさいけど、他言語は意外にその辺一貫してないよな。

まぁFORTRANの場合、依存性があると最適化の邪魔になるから、
依存性を嫌っているだけだけどｗ　結果的にはより徹底している。

**デフォルトの名無しさん** · 2014/08/21(木) 12:25:07.46

マットラボてどうなん？フォートランと比べて何か秀でた点ある？

**デフォルトの名無しさん** · 2014/08/21(木) 22:20:31.76

そんなん解決したい問題,実装したい理論によりけりやろ
まあ計算規模が小さくて使いたいライブラリがMATLABにあるならそっち使ったほうが幸せになれるぞきっと

**デフォルトの名無しさん** · 2014/08/23(土) 12:29:18.39

>>573
ま、それもそうだよな。
グラフのフィッティングソフトでいいのなかったので悩んでた。
無償でいいのないもんかねぇ。
Fittkだかもあんま美味だし…

**デフォルトの名無しさん** · 2014/08/24(日) 20:56:38.40

MINUIT使えば

**デフォルトの名無しさん** · 2014/08/28(木) 13:13:31.88

MATLABって個人使用なら五千円に値下げになってなかったけ？

**デフォルトの名無しさん** · 2014/08/30(土) 13:35:28.03

>>575
ググってもなかなか情報ないんだね。
でもfortrN で書かれたパケなんですね！
これは貴重な情報です、ありがとうです！