くだすれ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/

**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かー,読みやすさを優先してくれないからなー
コンパイラのマニュアルは殆ど英語になるのが悲しい

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

>700,703
Fortran90の日本語の解説本で良いなら、
　bit別冊　詳解Fortran90　1993年12月号別冊　共立出版
という本があります。これは『Fortran90 Explained』の日本語訳です。
但し、雑誌の別冊ということもあり入手は極めて困難でしょう。
あるとすれば、理工系の大学図書館が一番可能性が高いと思います。
参考までにどうぞ。

**デフォルトの名無しさん** · 2014/12/29(月) 01:02:13.12

Fortranを使われている方は他にはどんな言語を扱っているのですか？
pythonとかでしょうか

**デフォルトの名無しさん** · 2014/12/29(月) 09:04:30.59

>>705
perl

**デフォルトの名無しさん** · 2015/01/04(日) 14:13:23.64

>>700
おいらもFORTRAN77使いなんで、小文字のプログラムに違和感ありまくり

**デフォルトの名無しさん** · 2015/01/04(日) 16:56:46.08

>>707
こっちはfortran90からの新参なんで、77のソース見ると大文字、goto文ありでなかなか新鮮です。

ところでfortran2003の機能すべてカバーしてるコンパイラは商用しかないのかな。
gfortran使いなんでカナピー(´･_･`)

**デフォルトの名無しさん** · 2015/01/04(日) 20:02:59.44

商用だけですな
しかし商用Fortranコンパイラは趣味で買うには高すぎるので
gfortranが対応してる範囲で書いてます

**デフォルトの名無しさん** · 2015/01/04(日) 23:43:36.28

>>708
Linuxならifortのnon-commercial板、MacならFreePGIがあるけど
Windowsは無いみたいだね。

**デフォルトの名無しさん** · 2015/01/04(日) 23:48:08.93

>>710
変換ミス
（誤）non-commercial板
（正）non-commercial版

gfortranの2003フルサポートは何時になるのやら。

**デフォルトの名無しさん** · 2015/01/05(月) 11:27:48.33

>>710
勉強用なら
VirtualBox(個人使用は無料だったハズ) + Ubuntu日本語Remix仮想ハードディスクイメージ
で、とりあえずLinux環境作っちゃうのも手なんかね

…もちろん、LinuxやBSDに詳しい人にはこの構成では不満が出てくるだろうが
(特にUbuntuの標準GUI、嫌いな人結構多いしな)
知らん人がサクッとUnix系の環境作るには楽だと思った

**デフォルトの名無しさん** · 2015/01/05(月) 12:04:08.73

>>712
>VirtualBox(個人使用は無料だったハズ) + Ubuntu日本語Remix仮想ハードディスクイメージ
>で、とりあえずLinux環境作っちゃうのも手なんかね

Fortranの勉強のためならこれが正解だと思う。わざわざLinux用の領域をハードディスク上に
確保して…というのは面倒だしね。自分も次は仮想マシンを使うと思う。
UbuntuのGUIが気に入らない場合は、VirtualBoxにインストールできるか知らないけれど、
RHELクローンやOpenSUSEもあることだし。

**デフォルトの名無しさん** · 2015/01/05(月) 12:13:27.17

VirtualBoxに入れるLinuxなんて、GUI抜きで容量減らしてインストールして
sshでもっと使い易い端末エミュレーターから操作するのが普通じゃないか

**デフォルトの名無しさん** · 2015/01/05(月) 12:35:46.63

>>714
もちろんそれが理想的だし、解ってりゃ大したことないんだけどね

**デフォルトの名無しさん** · 2015/01/07(水) 02:53:13.20

でも最近の Linux intel fortran も windows と同じ GUI でインストールできる
ようになってるし～楽ちんしたいのよ～

**デフォルトの名無しさん** · 2015/01/09(金) 00:18:45.90

ifortやgfortranのオプションについて教えてください．

Makefileを書くとき，コンパイルオプションやリンクオプションをCFLAGSとLDFLAGSに分けて書くのが一般的みたいですが，
Intel Fortranコンパイラのユーザー・リファレンス・ガイドやgfortran/gccのマニュアルを読んでも，
どのオプションがコンパイルオプションで，どれがリンクオプションか明記されていませんよね？

コンパイラとリンカに与えるべきオプションは，それぞれ，どうやって見分けたら良いんでしょうか？

**デフォルトの名無しさん** · 2015/01/16(金) 00:14:55.27

intel fortranの非商用版終了したっぽいな

**デフォルトの名無しさん** · 2015/01/16(金) 07:28:02.50

>>718
は？
ソースは？

**デフォルトの名無しさん** · 2015/01/16(金) 10:37:10.07

>>719
intelのdeveloper zoneで明言してるよ
https://software.intel.com/en-us/forums/topic/533638

**デフォルトの名無しさん** · 2015/01/16(金) 11:08:31.95

>>720
ありがとうm(__)m

**デフォルトの名無しさん** · 2015/01/20(火) 00:53:26.59

薄汚い乞食が多すぎたんだな

**デフォルトの名無しさん** · 2015/01/20(火) 02:58:09.86

>>722
確かに研究用途で無償版使われているって書いてるな

**デフォルトの名無しさん** · 2015/01/20(火) 16:42:37.43

すみません。本読んだりGoogle等で調べましたが、いまいちわからなかったので質問させてください。
result.datというデータの中にある本来１７行２１列で表示するデータが
#(A,B)

(1,0)
　・
　・
　・
(17,0)
(0,1)
　・
　・
　・
(17,1)
(0,2)
　・
　・
　・
(17,2)
　・
　・
　・
(17,20)
というように１列で表示されているデータがあります。
openを使ってこれを読み込んで本来の１７行２１列の形にして別のデータとして出力したいのですが、
どなたか知恵をお貸しいただけませんか？またできればデータの(A,B)のAの部分だけ、Bの部分だけを取り出したいのですが、
どうすればよいのでしょうか？

**デフォルトの名無しさん** · 2015/01/20(火) 17:17:10.54

>>724
過去レスヨメカス

**デフォルトの名無しさん** · 2015/01/20(火) 17:19:19.79

>>725
ソウイワナイデオシエロヨ

**デフォルトの名無しさん** · 2015/01/20(火) 17:39:47.89

>>725

すみません。過去すれにあるのですか？
ざっとしか見てなかったので気づいてませんでした。
もう一度探してみます。

**デフォルトの名無しさん** · 2015/01/24(土) 14:38:14.65

cosやexpとかの組み込み関数が学校のifortだと通るんだけど家のgfortranだと通らない助けて
no implicit typeって言われる

**デフォルトの名無しさん** · 2015/01/24(土) 16:14:31.36

>>728
その関数に渡す実引数の宣言で処理系依存の型パラメータを使っているとか？
ソースを晒せば適切な助言を受けられると思う。

**デフォルトの名無しさん** · 2015/01/24(土) 17:04:25.97

４倍精度から倍精度にしたら実行できました
gcc4.6以降は4倍もサポートしてるそうなんですがおかしいですね…
ありがとうございました

**デフォルトの名無しさん** · 2015/01/25(日) 14:22:03.88

>728

fortran 4.8.1
で
program test

real(16) :: D2R

D2R = 0.0174532925199432957692369076848861271344_16

write(*, *) '1/e= ',1/exp(1.0_16)
write(*,*) 'cos(30)= ',cos(30.0*D2R)

end program test
を実行すると
1/e= 0.367879441171442321595523770161460873
cos(30)= 0.866025403784438646763723170752936161
という答えが得られます
少なくとも30桁以上の精度がでているとおもいますが、

**デフォルトの名無しさん** · 2015/01/25(日) 14:24:04.46

>731

>fortran 4.8.1
は
gfortran 4.8.1

のtypoです。しつれいしました。

**デフォルトの名無しさん** · 2015/01/27(火) 15:08:58.38

>>730
もう解決しているかも知れないけれど、処理系がFortran2008と4倍精度をサポートしていれば
iso_fortran_envモジュールのreal128という定数を使う言語仕様に準拠した方法がある。

program qfloat

use iso_fortran_env
implicit none

real(kind = real128) :: q

q = 1 / exp(1.0_real128)
print *, q

end program qfloat

**デフォルトの名無しさん** · 2015/01/30(金) 19:29:21.36

超初心者です。質問させていただきます。
forrtl:severe(24): end-of-file during read,unit 110,file C:tkshmd\www01.dat
というエラーが出たのですがどうすれば解決できますか？

教授から卒業研究用に渡されたもので、自分で書いたものではありません
11年前に卒業した先輩が書いたものらしいのです

fortranに関する知識は一切ございません。すみません

**デフォルトの名無しさん** · 2015/01/30(金) 19:35:01.91

>>734
どんな環境なのか、OSとコンパイラくらい書こうよ。
エラーメッセージによるとwww01.datの途中に
ファイルの終端を表すend-of-fileが挿入されているっぽいが。

**デフォルトの名無しさん** · 2015/01/30(金) 20:02:06.39

>>734
データの読み込み中にEOFが来た
つまり、読み込む変数に対してファイルのデータ数が足りない

先ずは実行できるデータファイルを貰って
自分の環境でプログラムが動くことを確認したら？
自分の計算はその実行できるファイルをいじってする

何より教授に聞けよ

**デフォルトの名無しさん** · 2015/01/30(金) 20:49:39.27

>>734
736で解決すると思ったけど
フォルダアドレス名、もしかして名前の母音抜き？
違ったらごめんね
もしそうだったとしたら、今後フォルダ名の投稿の時とか気をつけてね
私、昔似たようなことやっちゃったので。

関係ないが気になったのでつい

**デフォルトの名無しさん** · 2015/02/04(水) 05:53:12.64

ちょいとみなさんのお知恵拝借させてもらえたら幸い。

二分法かなんか使ってイテレーションさせようかと思ってるんだけど、ある解に収束するまでの途中のデータが出力ファイルに書き込まれて困ってる。

収束した時のみの結果を出力させたいんだけどどうしたもんだろ？

すごく簡略化した例だけど

do i = 1, imax
....
! F(x(i),y(i),z(i)) は非線形で数値計算でその都度計算するしかない
F(x(i),y(i),z(i))= .........

write(15,*) x(i),y(i),z(i)

! イテレーションする部分
If(s(x(i),y(i),z(i)) < 1.0d-5) then
Exit
Else
....
Endif
Enddo

**738** · 2015/02/04(水) 05:54:22.46

連投すんません
S(x(i),y(i),z(i))は新しい関数ですm(__)m

よろしくお願い申し上げます。

**デフォルトの名無しさん** · 2015/02/04(水) 21:13:09.14

>>738
収束ループを出てから、書いたら？

**デフォルトの名無しさん** · 2015/02/04(水) 21:21:39.45

求めた解はどこかで使うんでしょ？なら，使う直前にwrite文を書けばいいのではないか？

というか，
>二分法かなんか
『なんか』ってなんだ？そんなきき方でどう答えろというんだろうか？
その例の書き方だと，そもそも二分法の使い方を理解しているのかと，そっちの方が心配になる．

**738** · 2015/02/10(火) 07:40:57.77

>>741
すまん！書き方悪かった。
「なんか」と書いたのは二分法等イテレーションする場合と言いたかった。

**デフォルトの名無しさん** · 2015/02/11(水) 01:54:42.49

debugをしたら
user breakpoint called from code at 0x7c94120e
と出てきたのですが、どういった意味なんでしょうか

**デフォルトの名無しさん** · 2015/02/11(水) 02:18:15.77

環境書こうよ……
まず、>>1の注意事項を読んでくれ

**デフォルトの名無しさん** · 2015/02/11(水) 02:35:40.04

>>743
回線切って逝ってよし

**デフォルトの名無しさん** · 2015/02/11(水) 10:40:19.88

>>743
デバッガがプログラムを呼び出して実行してブレークポイントに到達して止まった
つまり意図通り
ブレークポイントを誰がどうやって設置したかは環境次第

**デフォルトの名無しさん** · 2015/02/11(水) 22:51:10.13

コンパイラはどれを使っていますか？

**デフォルトの名無しさん** · 2015/02/11(水) 23:16:27.84

ｇ95

**デフォルトの名無しさん** · 2015/02/15(日) 01:10:05.40

Rの関数(cor(x,y)など）をFortranプログラム上で使えるようにできないでしょうか
Rに入っている関数のsoファイルを利用してなどできたりしないでしょうか

**デフォルトの名無しさん** · 2015/02/15(日) 13:38:58.18

>>749 あなたは Fortran ではなく R を使うべきなのではないのでしょうか

**デフォルトの名無しさん** · 2015/02/15(日) 14:45:03.59

Rがコマンドから使えるなら普通にSYSTEMすればいいんじゃないの

**デフォルトの名無しさん** · 2015/02/15(日) 22:02:11.88

>>750-751
事情がありFortranで組む必要があったので仕方なく・・
ありがとうざいました

**デフォルトの名無しさん** · 2015/02/15(日) 22:03:52.96

そもそもFORTRANにはRなんか手も足もでないような統計ライブラリが存在するはずだが

**デフォルトの名無しさん** · 2015/02/16(月) 00:11:19.19

>>753
IMSLかMKLの統計サマリーですかね？
使いたいのは比較的新しめの統計量で、
そこには用意されていないけどRにはあったので・・

**デフォルトの名無しさん** · 2015/02/16(月) 04:49:19.07

相関係数求めるだけだろ？
http://tama.green.gifu-u.ac.jp/~tama/soft/sample_prog/correlation/

**デフォルトの名無しさん** · 2015/02/19(木) 01:32:20.14

今更ながらなんだがフォートランとオクターブの組み合わせ最高だよな。
今まではfortran とgnuplot だったけど

**デフォルトの名無しさん** · 2015/02/19(木) 21:39:16.47

お絵かきできんの？

**デフォルトの名無しさん** · 2015/02/21(土) 13:41:52.59

オクターブってMATLABのパチモンだっけ？
最近もメンテされてるの？

**デフォルトの名無しさん** · 2015/02/21(土) 20:21:39.11

>>757
余裕でお絵描きできるよ。
お絵描き機能は gnuplot のを踏襲してるぽ

>>758
まだまだ現役

**デフォルトの名無しさん** · 2015/02/22(日) 06:07:07.89

お絵かきとは２ｃｈのお絵かき機能のこと。
この板はできないのか？

**デフォルトの名無しさん** · 2015/02/24(火) 12:10:45.71

3/3から2chのapi始動するけど移住しなくても大丈夫かね
とりあえずdat廃止されるまではここはNavi2chで巡回するけど

**デフォルトの名無しさん** · 2015/02/25(水) 18:07:19.77

>>761
移住先の候補があるなら教えて欲しい。

**デフォルトの名無しさん** · 2015/02/27(金) 19:13:57.57

とりあえずおーぷん2chにたてるだけたてといた
http://toro.open☆2ch.net/test/read.cgi/tech/1425031440/
（☆はNGワード対策）
何事もないことが一番だけどそんときは避難所代わりに使ってください現に今人多杉状態だけど
ここ質問スレだけど他にfortranスレないしfortran全般に関する喋り場ならおーぷんでいいよね？転載嫌いとかアフィリエイト嫌いとかのひとがいたらすまんが

純粋な質問に関してはstackoverflow日本語版も始まったしそっちも賑わってほしいなあと思ってる（現在3件しかない）
本家stackoverflowの方はfortranに関してもかなりお世話になってるし

**デフォルトの名無しさん** · 2015/02/28(土) 14:09:46.29

お若い方。
この世間に疎い爺に何が起きているのか３行くらいで教えてくれんかね？
カチュ～シャが使えなくなってから専ブラともお別れじゃよアッハッハ

２ｃｈがたらこのものでなくなったらしい話は風のうわさに聞いたども、
移住とは何のことじゃね？

**デフォルトの名無しさん** · 2015/02/28(土) 14:12:07.22

>>763
おーぷん2chは知らなかった。何かあればそちらを使わせてもらいます。
stackoverflowの方はgoogleでの検索がヒットした時ぐらいしか行かなかったけれど、
日本語版はありがたいね。こちらも覗きに行きます。

**デフォルトの名無しさん** · 2015/02/28(土) 14:20:37.08

>>764
私もよく知らないが、現在の専ブラが使っているdatが廃止される。2chを専ブラで利用するためには
2chが提供するAPIが必須となるが、このAPIは2chから使用許諾を受ける必要がある。
よって今のようなオープンソースの専ブラの開発、提供は難しいようだ。

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

>>492
私も見たいと思ったのですが…
該当章はなくなってしまったようですね．
（もう8か月以上だから仕方ないか）

**デフォルトの名無しさん** · 2015/03/09(月) 22:37:50.67

おいおいTip集に分離されただけで記事は残ってるぞ
nagのトップページから見てみ

**デフォルトの名無しさん** · 2015/03/10(火) 23:05:49.08

Hey!Yo!Brother!

Fortran2003でOOPを勉強してみたいんだが、NAGしか入門する門は無いのかい？
ジーザスは狭き門より入れと言ったが、おじちゃんは滅びの道でもいいから
広き門から入って太い道を通りたいんだＹｏ!

できればジーザスの言葉と同じアラム語かおじちゃんと同じ言葉の日本語で頼むよ！
敵性語はノンノン！

**デフォルトの名無しさん** · 2015/03/17(火) 02:57:32.45

>>784
ああ、それは腹立つな
Fortran Graphic、Fortran OpenGLで検索するとそこそこ出てくる
けど、まともにメンテナンスされてるのが無さげでもある……

**デフォルトの名無しさん** · 2015/03/21(土) 20:51:39.73

プログラム中のx行目?y行目のソースコードを出力させるwrite文の書き方はないでしょうか

**デフォルトの名無しさん** · 2015/03/22(日) 10:00:02.71

>>771
ねーよ

**デフォルトの名無しさん** · 2015/03/22(日) 14:18:33.06

>>771
プログラム単体では無理
単純に任意のファイルのx行目を表示するプログラムなら簡単
指定行まで読み飛ばして、指定行を十分長い文字列で読み込めば良い
それでソースファイルを読んで表示するならできるけど
ソースファイルは常に必要

**デフォルトの名無しさん** · 2015/03/22(日) 15:16:23.54

Quineが書きたいという話ではなくて？

**デフォルトの名無しさん** · 2015/03/22(日) 16:35:49.10

>>773
なるほどopenして目的行寸前まで空readしてread×行数すればいいということですね
解決できそうです。ありがとうございました
>>774
初めて聞きました。勉強になりました

**デフォルトの名無しさん** · 2015/03/24(火) 21:45:10.36

>>775
俺にもお礼言えよ、クソ野郎

**デフォルトの名無しさん** · 2015/03/26(木) 23:15:58.87

3次元のスカラーデータを吐く、x,y,zのdoループをompにしようとしたんだけど、スタックオーバーフローだなんて怒られた

どーすんべ？

**デフォルトの名無しさん** · 2015/03/27(金) 02:28:15.56

安直にはコンパイラ、リンカのオプションでスタック増やせ。

**デフォルトの名無しさん** · 2015/04/28(火) 17:28:05.66

MinGWのGCC4.9.2ののFORTRANで, -fcase-upperのようなオプションってありますかね?
昔のバージョンにはあったみたいなんですが, 現在のは削除されたようで.

**デフォルトの名無しさん** · 2015/05/10(日) 17:08:53.50

Intel fortran無料版のライセンスが改定されて
オープンソースへの貢献を証明しなきゃならなくなってた
https://software.intel.com/en-us/qualify-for-free-software/opensourcecontributor

**デフォルトの名無しさん** · 2015/05/11(月) 01:08:44.47

Gfortran とIntel fortran のそれぞれの利点って何かね( ´ ▽ ` )ﾉ

**デフォルトの名無しさん** · 2015/05/11(月) 01:47:31.11

え、gfortranに利点なんかあるの？
ああ、うん、gccツールチェインの一部なのは便利だよね

**デフォルトの名無しさん** · 2015/05/11(月) 06:54:00.10

ハードウェアで二倍の性能出す金払うくらいだったら、Intelのコンパイラ買えよ。
ソフトには金出せないってどんだけ唯物論のアカなんだよｗ

**デフォルトの名無しさん** · 2015/05/11(月) 10:33:22.97

>>782
なに、そんなレベル違いますのん？

**デフォルトの名無しさん** · 2015/05/13(水) 06:59:19.62

>>784
10年くらい前は差は歴然としてました。
最近はどうなのでしょうね？

**デフォルトの名無しさん** · 2015/05/15(金) 02:52:46.16

intelは自動並列化があるから中級者は使用を薦める

**デフォルトの名無しさん** · 2015/05/31(日) 00:46:49.58

このレス数からしてFortran使ってる人はほとんどいないのかな？

**デフォルトの名無しさん** · 2015/05/31(日) 01:15:28.56

ここでの質問は学生がほとんどのように見える
一般人は使わないだろうから少ないのは確かだろうけど
調べたり相談できる人が多いんじゃない？

**デフォルトの名無しさん** · 2015/05/31(日) 05:24:16.45

>>788
なるほど
じゃあ私遠慮無く質問しまくっちゃっおうかな

**デフォルトの名無しさん** · 2015/06/05(金) 00:40:30.28

fortran の組み込み関数でminval ってのがありますよね．
ある配列a(要素は一つだけ)に関して
minval(a,dim=1,mask=a>0.0d0)すると値が1e-324みたく値がおかしくなるんですが原因として考えられるものは何でしょう？
ちなみに配列aの中身は
a(1)= 4.2e-2です。
なので期待する値は
minval(a,dim=1,mask=a>0.0d0)=4.2e-2になるはずなんですが。。。

**デフォルトの名無しさん** · 2015/06/05(金) 01:52:18.65

>>790
なんとなく初期化してない変数が引っかかっている気がする
minvalの動作不良を疑うならコンパイラとか環境を書こう
違うならソースをできるだけ書こう
長くてソース貼れないならこんな感じに晒すと検証が楽
https://ideone.com/JstGnQ

**デフォルトの名無しさん** · 2015/06/05(金) 23:13:40.04

>>790
単精度を入れているのに倍精度で値が返ってきているのが気になる。
なんか倍精度と単精度をまぜこぜにしているか、コンパイラが間違っていてまぜこぜにしている。

ｇ９５はデフォの変数長を変えたせいで混乱に陥った。その辺かな？
FortranはRealとIntegerの長さが等しくないといけないのにそれを崩した。
まぁ６４Bit時代には８バイト整数にしたいのはわかるけど、
そうなると倍精度がRealになって、倍精度が４倍精度になってハードウェア対応せずで憤死する。

https://ideone.com/m11vQN

こんな感じかな？
program test2ch
implicit none
real(8) :: a(1)
real :: b(2)
equivalence(a, b)
b(1) = 4.2e-2
b(2) = 0.0e0
print *, minval(a, dim = 1, mask = a > 0.0)
end program test2ch
実行結果
5.070575254129006E-315

**デフォルトの名無しさん** · 2015/06/06(土) 00:34:55.86

>>792
ideoneってfortranいけるのかｗ

**デフォルトの名無しさん** · 2015/06/06(土) 01:24:14.49

変数に値が入ってるかどうかの判断ってどうしてる？ポインタとか使わずになんかNull かどうか？っての判断できんのかね？
ちなみにgfortran 4.2 でDebian 8.0 64bit版での話なんだが。

**デフォルトの名無しさん** · 2015/06/06(土) 02:41:06.93

>>794　？？
ASSOCIATED　とか　ALLOCATED　じゃダメなのか？
普通の変数の場合はどうしようもないよな。