くだすれ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/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 で書かれたパケなんですね！
これは貴重な情報です、ありがとうです！

**デフォルトの名無しさん** · 2014/08/31(日) 23:33:34.12

興味でお聞きしたいのですがfortranで株の発注のようなことはできたりしませんか？
例えば証券ホームページににログインし発注操作するなど・・
そういう言語じゃないのかもしれませんが・・

**デフォルトの名無しさん** · 2014/08/31(日) 23:49:50.53

fortranからC言語の関数呼び出せるので...
それじゃ邪道ですかw

まあHFTも金勘定も数値計算だしいいんじゃないでしょうかね

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

そう言えば、Fortranって少数第n位以下を切り捨てみたいな関数ないんだね…
Excelで言うrounddownみたいなやつ。
今、ファイル書き出すのに一々四捨五入されちまって困りマクリングでつ、はい(´･_･`)

**デフォルトの名無しさん** · 2014/09/01(月) 13:00:12.21

>>580
簡単だし自作したら良いじゃない？

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

>>579
ありがとうございます
なるほど他の言語を読み出せばいけそうなんですね
今後やってみたいです

**デフォルトの名無しさん** · 2014/09/02(火) 07:37:26.28

>>581
やはり自作しかないのね…
てかラウンドダウン程の関数が無い事に驚きだお
どっかでパッケージモジュールとして保管出来るようにするかね

**デフォルトの名無しさん** · 2014/09/02(火) 08:58:57.02

そんなに端数が気になるならもうひと桁増やして丸めればよかろうに

**デフォルトの名無しさん** · 2014/09/02(火) 20:48:24.61

むしろ、簡単に実装できるのにわざわざ関数がある言語のが珍しいかと

**デフォルトの名無しさん** · 2014/09/02(火) 21:36:00.72

>>584
0.4999999999を小数点以下3桁で切り捨て表示するために4桁で文字列化しても内部処理が丸めだと0.500になってしまうかもしれんぞ
素直に1000倍して切り捨てて1000で割るがよろしい

**デフォルトの名無しさん** · 2014/09/03(水) 02:46:13.39

>>580
なぬ？途中の計算の話じゃなくて、I/O 時の話か？
それならあるぞ。F2003もしくはF2008の機能だと思うが。

write(*, '(rd, f15.1)') 1.45 ! 1.4
write(*, '(ru, f15.1)') 1.45 ! 1.5

このほか rz, rn, rc もある。
それぞれ round up, round down, round to zero, round to nearest, round to compatible
最後の二つは最近接への丸めだが、等距離の時に、偶数への丸めかゼロから遠い方への丸めかだと思う
まぁ詳しくはマニュアルとIEEE754-2008の仕様をよく読んでくれ。

**デフォルトの名無しさん** · 2014/09/03(水) 13:00:08.93

昼休みに調べたら rd ru の類はF2003だった。
あと rp processor dependent という処理系依存型もあるようだ。
あとOPEN文で指定する方法もある。

しかし、正直いって>>580の基本的発想が間違っている。
他の人達のアドバイスにあるように、そもそも内部二進表現を十進表現に
直すときに丸めなどが入ってしまう。そのうえでさらに十進表現を
丸めようとするのが Fortran 的センスではおかしい。
だからみんな当惑した反応しか返さない。

とはいえ自然科学ではない商習慣等での丸め等もあるようだから
別にそれがいけないわけでなく、それに合った言語例えば
BCD演算できるCOBOLなどで計算すべきのみ。

**デフォルトの名無しさん** · 2014/09/03(水) 22:10:06.94

2015出たんだね

**デフォルトの名無しさん** · 2014/09/04(木) 02:18:06.97

Fortran コンパイラメーカ
Intel PGI Pathscale 　　（Intel NVidia AMD/ATI) アクセレレータ御三家

Absoft NAG 　　　　　　　　独自色　

Lahey/Fujitsu Silverfrost 　斜陽族

GFortran g95　　　　　　　　GNU Free系

こんなとこか？　

**デフォルトの名無しさん** · 2014/09/04(木) 02:51:19.67

Free系にOpen64と（一応）dragoneggがあるぞ
dragoneggはgcc4.8.3以降はサポートしてないしflangは開発止まってるし
この辺が活発化すれば若い人がたくさん来ると思うんだけどな

**デフォルトの名無しさん** · 2014/09/05(金) 04:17:28.12

GNUでFortranやAdaが熱心にサポートされているのは、米国防総省あたりからの
金が流れているからなのかな？

意識高いｗ進歩的ｗCS学科若年層には人気は全くないのに。
それともおっさん層しかやってないからなのか？
そのわりにGNU COBOLとかPL/IとかALGOL、PASCAL、Modula-2とかやらないし。

**デフォルトの名無しさん** · 2014/09/05(金) 08:06:03.87

>>587
でもfortran95でしょ、ほとんどサポしてるのは

**デフォルトの名無しさん** · 2014/09/07(日) 03:06:28.08

Cray, IBM, PGI, Intel が F2003 完全対応したな。
NAGはあと一歩。
GFortranは進捗どうですか？

スパコンだと　SXは90 Hitachiも90?　Fujitsuは95？　
それぞれつまみ食い的に拡張してあるが。

**デフォルトの名無しさん** · 2014/09/07(日) 04:28:58.48

日立のはコンパイルオプション-hf95で問題なく使えてたがあれ完全準拠じゃなかったのかな
まあ完全準拠だったらわざわざfortran90コンパイラとは名乗らないか

**デフォルトの名無しさん** · 2014/09/10(水) 14:57:01.18

そういや、fortranのデータの取り扱い方について詳しくかいせつしてる本てあまりないね。

**デフォルトの名無しさん** · 2014/09/13(土) 13:38:42.01

fortranのデータってなんじゃい？ファイルのことか？

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

>>596
format read writeのことなら、コンパイラーの取扱説明書が（独自仕様拡張部分も含めて）もっとも詳しいよね。

**デフォルトの名無しさん** · 2014/09/14(日) 15:27:09.67

intel fortranはIMSLも非商用は無償にしてくれないかな・・

**デフォルトの名無しさん** · 2014/09/15(月) 02:35:29.49

IMSLは会社が違うから無理だべ。
IntelはCPU売るためにソフトのコンパイラを撒き餌としてばら撒いてもペイするが、
IMSLの方としては特に得しないし。

富士通がLaheyにOEM供給した時にコンパイラにSSL2とかただでバンドルしていたが。
NECがASLをちょー高値で売っていたことを考えれば、ちょー太っ腹だった。

それはそうと日本の大学の大型計算機センターにたまっていたライブラリも、
NETLIBみたいにネットに公開して吐き出せばいいのに。

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

>>598
Gfortran なんやけど…(´･_･`)
もうやめようかな。

Gfortran で書いたプログラムをIntel fortran でもコンパイルすること出来ますか？

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

GForce思い出したわ

**538** · 2014/09/16(火) 03:14:29.14

>>538ですが、いまだに解決しません。
>>543の「数値が指定カラムからはみだしてないか」という意味がよくわかりません。
指定カラムというのはどこで見ることができるのでしょうか。
9.0と11.0ではビルドのUIが大きく変わっており、ビルドの条件が同じなのかどうか自信がありませんが、ほぼデフォルト状態でビルドしています。

**デフォルトの名無しさん** · 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の解説をしてもらえないでしょうか
まったく意味がわからない