くだすれ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/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の解説をしてもらえないでしょうか
まったく意味がわからない

**デフォルトの名無しさん** · 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で書かれています