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

**デフォルトの名無しさん** · 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にはあったので・・