FORTRAN Ⅳat TECH
FORTRAN Ⅳ - 暇つぶし2ch178:デフォルトの名無しさん
07/01/08 20:39:25
INTERFACE
LOGICAL FUNCTION PrintMessage()
!DEC$ ATTRIBUTES STDCALL :: PrintMessage
END FUNCTION
END INTERFACE

LOGICAL LX
LX = PrintMessage()

しかし、リンク時に_PrintMessageが見つかりません、と
リンクエラーが出ます。

解決方法はあるでしょうか?


なお、VCから上記VCを呼び出すケースもあるため、
両方Fortranで書けば?というのは無しでお願いしたいです。

一応環境としては、
Visual Studio .NET 2003 + Intel Visual Fortran v9.0
になります。

よろしくお願いします。

179:デフォルトの名無しさん
07/01/09 00:52:12
>>178
呼び出し規約がSTDCALLの場合、外部名はデフォルトですべて小文字と解釈されてしまうので
大文字・小文字が混在する場合は以下の様にする

!DEC$ ATTRIBUTES STDCALL, DECORATE, ALIAS: 'PrintMessage' :: PrintMessage

C側の PrintMessage() 関数も __stdcall でコンパイルされている必要がある。
C++の場合は extern "C" も必要。

また、IVFではLOGICALの内部表現(整数値)が奇数のとき真、偶数のとき偽となるので
他言語の関数とLOGICALを受け渡しする場合は注意すること。

180:178
07/01/09 21:34:40
>>179
レスありがとうございます。

なるほど・・。
F側の大文字小文字対策はマニュアルに書いてありましたが、
いろいろ調整しても結果が同じだったので関係ないと思ってました。
C側の対策は全く見えてませんでした^^
正直Cはド素人なので上記のようにいわれてすぐに対応できませんが
参考にして調べてみます!
ありがとうございました。




181:デフォルトの名無しさん
07/01/09 23:27:06
>>180
IntelFortranのマニュアル見るとちまちまと書いてあるよ。

DECORATE修飾子が新設されたので、DEC時代に場合分けして書いていたのが
すっきりしてうれしい。



182:デフォルトの名無しさん
07/01/10 19:42:57
>>180
ATTRIBUTES STDCALL, ... ではなく ATTRIBUTES C, ... を使えばC言語の方はいじらなくてもおk。

183:デフォルトの名無しさん
07/01/12 23:52:31
三つの数値データを与えてその三つの数値の最大値、最小値、中央値、
平均値、標準偏差、不偏分散を求められるプログラムを作成せよ。
わからなくて大変です。
お願いします

184:デフォルトの名無しさん
07/01/16 12:41:49
Fortressどう思います?
数式をそのまま書けるとか、並列化にコンパイラレベルで対応とか・・・
なかなかよさそげですが。
URLリンク(japan.zdnet.com)
URLリンク(slashdot.jp)

185:デフォルトの名無しさん
07/01/16 12:45:26
>>183
てっとりばやく答えが欲しいなら
スレリンク(tech板)
に行くといいのでは?

186:デフォルトの名無しさん
07/01/16 20:18:01
既に答えが出てる。 ていうかマルチ。

187:デフォルトの名無しさん
07/01/17 03:15:36
>>184
Fortranの資産使えないのにFortran風にする意味がいまいち良く分からん。

188:デフォルトの名無しさん
07/01/19 12:05:30
>>184
pseudocode をそのままコンパイルしてくれるならうれしい。


189:デフォルトの名無しさん
07/01/19 12:55:53
>>187
え、使えないの?

190:デフォルトの名無しさん
07/01/19 23:38:03
Salford ftn95でlapack95を使いたいねんけど、どーすればいいん??

191:デフォルトの名無しさん
07/01/20 00:43:01
>>190
1.NetlibからLAPACKを拾ってくる。
2.どうにかしてblas.libとlapack.lib

192:デフォルトの名無しさん
07/01/20 00:44:49
>>190
1.NetlibからLAPACKを拾ってくる。
2.どうにかしてblas.libとlapack.libをつくる。
3.NetlibからLAPACK95を拾ってくる。
4.どうにかしてlapack95.libとモジュールファイルをつくる。
5. > ftn95 hoge.f90
6. > slink hoge.obj lapack95.lib lapack.lib blas.lib

193:デフォルトの名無しさん
07/01/20 08:07:50
>>192
どもども!
どうにかしてってところが問題だが、どうにかしてみます。
どうにかする方法がわかったら、またよろしゅ~!

194:デフォルトの名無しさん
07/01/20 11:45:48
77で動的配列ってできないの?
x(i,j)
ってしといて
ファイルからi,j読み込んでサイズ決めるみたいな感じにしたいんだけど・・

195:デフォルトの名無しさん
07/01/20 13:15:50
>>194
言語規格としては出来ない。

ただ、コモンとか上のルーチンで大きなワーク配列を取っておいて、
いったん必要なサイズを知った上で下のサブルーチンで作業領域で整合配列を
切って使って似たようなことは出来る。

C言語もK&Rの付録を見るとmallocをそういう感じで実現する例を出している。


196:デフォルトの名無しさん
07/01/20 13:38:47
90ならポインター使える。
PGIとか77でもポインター使える処理系もある。
ポインターの使えない77でもCとの合わせ技で可能。

197:デフォルトの名無しさん
07/01/20 16:25:44
>>193
LAPACK(BLAS):
・netlibから lapack-3.1.0.tgz を取ってくる。

blas.lib:
・lapack-3.1.0\BLAS\SRC の中のソース(*.f)を全部コンパイルする。
・出来たオブジェクトファイルを slim で blas.lib にまとめる。
 slim *.obj /file:blas

lapack.lib:
・lapack-3.1.0\INSTALL の中の dlamch.f, dsecnd.f, ilaver.f, lsame.f, second.f, slamch.f を
lapack-3.1.0\SRC にコピーする。
但し、FTN95には second.f, dsecnd.f で使われている ETIME関数が無いので
ETIMEを自作するか CPU_CLOCK, SYSTEM_CLOCK辺りを使って書き換える。
・lapack-3.1.0\SRC の中のソースを全部コンパイルする。
この際、slaruv.f, dlaruv.f をコンパイルするときは最適化オプションを外す。
・出来たオブジェクトファイルを slim で lapack.lib にまとめる。

tmglib.lib:
・lapack-3.1.0\TESTING\MATGEN の中のソースを全部コンパイルする。
この際、slaran.f, dlaran.f をコンパイルするときは最適化オプションを外す。
・出来たオブジェクトファイルを slim で tmglib.lib にまとめる。


198:デフォルトの名無しさん
07/01/20 16:27:59
>>197の続き)

LAPACK95:
・netlibから lapack95.tgz を取ってきて展開する。
・lapack95\SRCに移動。
・la_auxmod.f90 をコンパイルする。 → la_auxmod.mod, la_precision.mod が出来る。
・f77_lapack_single_double_complex_dcomplex.f をコンパイルする。 → f77_lapack.mod が出来る。
・f95_lapack_single_double_complex_dcomplex.f をコンパイルする。 → f95_lapack.mod が出来る。
・la_csygst.f90, la_csygv.f90, la_csygvd.f90, la_csytrd.f90, la_zsygst.f90, la_zsygv.f90, la_zsygvd.f90, la_zsytrd.f90
以外の la_*.f90 をコンパイルする。
・出来たオブジェクトファイルを slim で lapack95.lib にまとめる。


以上で出来たライブラリファイル(*.lib)を環境変数 LIB で指定したフォルダへ、
モジュールファイル(*.mod)を環境変数 MOD_PATH で指定したフォルダへそれぞれ移す。


199:デフォルトの名無しさん
07/01/20 17:04:13
関数名を引数に使いたいのですが、EXTERNAL文を呼び出し側で使っても
うまくいきません。
下のプログラムなのですが、何処が間違っているんでしょう?

External func
real func,SUM
write(*,*) SUM(func)
end

real Function SUM(f)
External func
real f
SUM = f(1.0)+f(2.0)
return
end

real Function func(x)
real x
func = x
return
end



200:デフォルトの名無しさん
07/01/20 19:15:37
>>199
SUM() がF90の関数と混同されているんで内科医?
SUM0とかにすればおK

201:デフォルトの名無しさん
07/01/20 19:22:59
>>199
External func → external func, SUM

202:199
07/01/20 22:07:45
>>200>>201、ご指摘の方法で解決することが出来ました。
有難うございます。

203:デフォルトの名無しさん
07/01/20 22:37:50
>>197,198
ありがとう~、格闘の末何とかなりました。
オレもあんたくらい、pcを自由にあやつりたいもんですわ!

204:デフォルトの名無しさん
07/01/21 00:37:48
本来ならMakefileを書いて、テストプログラムも含めてmake(MK32)で処理させるべきなのだろうけど
FTN95の場合、他のコンパイラとオプションの指定方法が違ったりするから中々面倒くさい。

205:デフォルトの名無しさん
07/01/26 07:14:03
質問です。
OPEN (4,FILE=RESULT.OUT,FORM='FORMATTED')
このRESULT.OUTを、あるIと言う数値を読み込んで
例えばI=5000ならRESULT5000.OUTのようにしたいんですが、
出来るでしょうか?要するに計算のステップ数に応じて
自動的にファイル名をつけて出力してやりたいんですが、
OPENの中身をいじるだけじゃダメですかね?

206:デフォルトの名無しさん
07/01/26 10:04:30
>>205
その質問は千回くらい出ている。
過去ログ見れ



207:デフォルトの名無しさん
07/01/28 18:26:34
自分のマシンのフォートランコンパイラコマンドを汁にはどうしたら良いですか?

VINELINUX使ってます。

208:デフォルトの名無しさん
07/01/28 20:37:47
>>207
shell に g77 とか gfortran とか g95 とか ifort とか入れてみる。
使えるようなら no input files などと返ってくる。
なけりゃ自分でインストロールする。

209:デフォルトの名無しさん
07/01/28 21:47:58
全部使えませんでした。
ifortいれたいのですが、どこでダウンロードできますか?

210:デフォルトの名無しさん
07/01/28 21:54:39
>>209
つ google

211:209
07/01/28 22:35:38
何度もすみません。
インテルのコンパイラ入れようとしたら下記のメッセージで止まってしまいます。
解決方法を御存じのかたは教えてください。

./.././data/install_fc.sh: line 1459: ERROR: unable to find command "ls" !: command not found
./.././data/install_fc.sh: line 1463: Please add the location to the above commands to your PATH and re-run the script.: command not found
./.././data/install_fc.sh: line 1464: Please press Enter to continue...: command not found


212:デフォルトの名無しさん
07/01/28 22:39:51
>>211
> 解決方法を御存じのかたは教えてください。
夜間中学で英語を学びなおす。

213:デフォルトの名無しさん
07/01/28 22:48:29
つーか lsが見つからんて、どんな環境やねん

214:デフォルトの名無しさん
07/01/30 00:02:48
>>211
URLリンク(www.google.co.jp)


215:デフォルトの名無しさん
07/01/30 01:02:12
EclispのPhortranってどうよ?

とくにLinux上のIntelコンパイラとの組み合わせの使い心地 知りたい

216:デフォルトの名無しさん
07/02/03 07:25:11
Mac OS X 上のgfortran 4.0でコンパイルした時に
gfortran: spec failure: unrecognized spec option 'M'
というエラーがでるのですが何がまずのでしょう? ちなみに他のマシン(LINUX)のgfortranでは
エラーはでません。一応、走るのですが、writeでposition='append'を指定しても既存ファイルに
データを付け加えず、完全に上書きしてしまうので困っています。ちなみにLINUX上では
うまく行きます。

217:209
07/02/04 11:51:06
>>214
ありがとうございます。
感謝します。

218:デフォルトの名無しさん
07/02/15 00:11:54
NAGのFortranBuilderのNAGFortran95コンパイラのVersionって5.0?それとも5.1?

5.1なら買いたいんだけど、誰か知ってる人いる?


219:デフォルトの名無しさん
07/02/16 15:57:03
全角英数は気持ち悪いというセンスを身に着けてから出直してくれ。

220:デフォルトの名無しさん
07/02/16 19:36:46
>>219
おい、おまえ





同感

221:デフォルトの名無しさん
07/02/17 00:00:33
>>219
ごめん 訂正させてもらう

NAGのFortranBuilderのNAGFortran九五コンパイラのVersionってⅤ.0?それともⅤ.Ⅰ?

Ⅴ.Ⅰなら買いたいんだけど、誰か知ってる人いる?

全角漢数字かローマ数字をつかえばいいわけだろ?

222:デフォルトの名無しさん
07/02/17 00:07:15
つまらん

223:デフォルトの名無しさん
07/02/17 00:22:25
フォートランって今でもリアルに使ってるとこあるのか?w
本物のプログラマは~、の時代でもすでに時代遅れのダメ言語だったのに

224:デフォルトの名無しさん
07/02/17 01:16:58
>>223
おまえ、ちょっとそこに正座しろ。
この世間知らずがw


いいか、電子計算機の本流はFortranだぞ!
スパコン・ランキングの機械で動いているプログラムは基本的にFortranだ。

大体マルチコアで並列化といっているが、MPIもOpenMPも対応コンパイラはFortranか
Cしかなくて、Cはあくまでおまけだ。大多数の人間はFortranでプログラム書いてる。

これから並列化が主流になるにつれ、Fortranの天下が再びやってくるのは間違いない!
悔い改めよ!さもなくば地獄の業火に焼かれるべし!
詳しくはマタヨシイエスのサイトを熟知すべし。


225:デフォルトの名無しさん
07/02/17 01:22:52
唯一神Fortran

226:デフォルトの名無しさん
07/02/17 01:36:08
>>221
「全角英数」の意味がわかるようになってから出直してくれ。

227:デフォルトの名無しさん
07/02/17 01:43:01
>>216
gfortran どやっていれたんですか?
俺のmacにははいってないんだてど

228:デフォルトの名無しさん
07/02/17 01:53:16
>>227
君のMacはハズレだったんだよ。

229:デフォルトの名無しさん
07/02/17 01:54:03
>>226
おまいら,うぜーから他所でやれ.

230:デフォルトの名無しさん
07/02/17 10:54:58
>>229
うるせー、>>218>>221のようなグロは見たくねーんだよ。

231:デフォルトの名無しさん
07/02/17 23:09:57
最近くだスレのほうがレベルが高いなw


232:デフォルトの名無しさん
07/02/21 15:40:26
>>218 5.1を搭載

233:デフォルトの名無しさん
07/02/21 22:54:23
>>232
㌧クス!


nagのサイトに調べにいったんだが、Windows用はまだ5.0だったので、だめかと思ったが。

安いやつ買って見っかー でもライセンスがウゼーんだよな。

234:デフォルトの名無しさん
07/02/25 00:31:39
Intel Visual Fortran 9.1に付属してくるデバッガ(idb)の詳しい日本語の資料や書籍ってないんでしょうか?

235:デフォルトの名無しさん
07/02/25 17:39:49
エクセルソフトがコンパイラマニュアルの日本語訳を出していたが、デバッガもあったかな?

236:デフォルトの名無しさん
07/02/25 23:30:36
>>235
XLsoftから日本語マニュアルを落としてみましたが、
詳しく載っているのはVisual StudioでのデバッグだけでIDBについては主要なコマンドの説明のみでした。
>IDB に関する詳細は、[スタート] メニューから利用可能なインテル・デバッガーのオンラインマニュアル (英語) を参照してください。
とのことで・・・このマニュアル、完全に日本語訳されてる訳じゃないんですね。
Language Referenceとか英文のままだし。

237:デフォルトの名無しさん
07/03/11 12:16:05
>>236
URLリンク(www.xlsoft.com)
日本語版スペシャルエディションでもなんか中途半端だし、しょうがないんじゃない?

コンパイルする時のオプションてなに使ってる?/O3とか普通に使って大丈夫なのかな・・・


238:デフォルトの名無しさん
07/03/11 12:20:26
Interl Visual Fortranのデバッガは糞だと思う。HPやCompaqの頃までのが一番良かった。

239:デフォルトの名無しさん
07/03/11 15:45:54
Compaq Visual Fortranの頃はVisual Studioのデバッガだけじゃなかったっけ

240:デフォルトの名無しさん
07/03/11 17:16:03
idbは概ねUNIX/Linuxのdbx/gdb互換(のコマンド体系)なので
これらについて書かれた資料を探せば参考になるかも知れん。

>>237
>/O3とか普通に使って大丈夫なのかな・・・
同時にデバッグ・オプションも付けとけばとりあえずおk、だが
最適化による命令の並べ替えやループのアンロールなどの影響で
ソースファイル通りにトレースできなくなる所が出てくる。

241:デフォルトの名無しさん
07/03/11 19:28:32
>>240
インテル・コンパイラはデバッグオプションつけても自動的に最適化は抑止されないのか?


242:デフォルトの名無しさん
07/03/11 20:07:27
>>241
ifort /debug:full hoge.f90
のようにデバッグ・オプション(/debug:full)だけ付けて最適化の指示を省略すると
すべての最適化が無効になる。

ifort /debug:full /O3 hoge.f90
のように最適化を明示(/On)すると、最適化した状態でデバッグ情報の生成を行う。

243:デフォルトの名無しさん
07/03/11 22:03:02
>>242
そういった資料はどこかにあるものだろうか?
必死こいてググって見たがイマイチよく理解できなくて・・・
何かよいものがあれば教えてください

244:デフォルトの名無しさん
07/03/11 22:43:38
>>243
XLsoftからDL出来る日本語マニュアルから。
>>242の内容なら、「アプリケーションのビルド」-「Debugging」-「デバッグと最適化」の所に載っています。

245:デフォルトの名無しさん
07/03/11 23:44:47
>>244
まさかそんなところに!!
ありがとうございます、熟読してきます。

246:デフォルトの名無しさん
07/03/20 19:03:09
ジョン・バッカス追悼age

John W. Backus, 82, Fortran Developer, Dies

John W. Backus, who assembled and led the I.B.M. team that created Fortran,
the first widely used programming language, which helped open the door to
modern computing, died on Saturday at his home in Ashland, Ore. He was 82.

URLリンク(www.nytimes.com)

247:デフォルトの名無しさん
07/03/21 11:13:38
【訃報】コンピュータ言語「FORTRAN」開発者、ジョン・バッカスさん死去…米・オレゴン州
スレリンク(newsplus板)

248:デフォルトの名無しさん
07/03/21 13:16:33
ジョンさん、懐かしい

249:デフォルトの名無しさん
07/03/21 14:02:59
【訃報】ジョン・バッカスさん死去 コンピューター言語「FORTRAN」開発
スレリンク(news板)

250:デフォルトの名無しさん
07/03/21 14:50:05
バッカスって酒飲みの

251:デフォルトの名無しさん
07/03/23 19:54:29
すいません、初心者です。
↓のようなデータを画面に出したいのですが
参考プログラムを教えていただけませんか?
1964/1/1 115 89 80 81 104 154 209 256 287 298 285 260 229 199 176 170 189
1964/1/2 155 111 89 81 91 124 171 227 266 293 292 274 252 221 190 174 181
1964/1/3 191 147 112 92 94 113 151 194 241 271 290 283 267 241 214 189 176
1964/1/4 222 185 147 117 101 107 134 168 208 248 274 284 274 256 233 208 185
1964/1/5 238 213 184 156 134 127 139 164 193 227 259 276 274 267 249 226 204
1964/1/6 248 234 214 190 165 151 149 161 181 206 234 257 267 266 257 243 222
1964/1/7 234 236 229 215 199 184 175 175 184 203 224 241 256 264 259 248 233
1964/1/8 216 228 234 229 221 213 203 192 190 197 210 225 239 252 258 254 243
1964/1/9 193 212 227 237 238 231 220 209 199 194 195 205 217 232 243 250 246
1964/1/10 163 181 203 226 238 243 240 229 217 205 195 196 202 215 232 245 251


252:デフォルトの名無しさん
07/03/24 07:34:06
開発者死んだのにお前ら流石にスルーし過ぎw

253:デフォルトの名無しさん
07/03/25 10:32:06
>252
手は合わせたが、2chで語る話題でもなかろう。

254:デフォルトの名無しさん
07/03/25 22:54:12
>>252
バッカスのヤオイ同人でも作ればいいのか?


255:デフォルトの名無しさん
07/03/25 23:23:03
まあ、バッカスさん自身はあまりFORTRANを好きじゃなかったみたいだしな。

256:デフォルトの名無しさん
07/04/10 08:28:32
writeをつかってバイナリデータを書き出しているのですが、
マシンによって出力されるバイナリが微妙に違い、
互換性のないバイナリファイルができてしまいます。
どのマシンでも共通のバイナリファイルを読み書きできるようにするには
どのようにしたらよいのでしょうか?


257:デフォルトの名無しさん
07/04/10 10:59:24
>>256
無理言うなw
どうしてもやりたければ、自分でビット列を生成して、それを1文字にして書け。
それでもあわせるのは不可能。

テキスト文字で出してIBM単密度形式の磁気テープでやり取りするのが一番確実だ。

大体、マシン毎にワード数(バイト数)が異なる場合がある。
昔のACOSは9ビット、CDCも妙なワード長だった。

文字コードにしても、IBM互換のEBCDIC系と(ただしメーカー毎に微妙に違う)、
DECやACOSのASCII系などが入り乱れている。

数値のフォーマットは会社毎にマシン毎にコンパイラ毎に異なるのが普通だ。

さらにはIBM、モトローラ系列BigEndianビット列LittleEndian、DEC、Intel系のLittleEndian系に
分かれている。これらもビット並びでさらに分裂するらしい。


とはいえ、明るい話もある。数値計算の世界では、もはやACOSもCDCもほとんど使われていない。
EBCDICもめったに見なくなった。ワード長は32ビットが普通だし、1バイトは8ビットに統一された。
数値のフォーマットもIEEE754にほぼ統一された。

残る問題はエンディアンだけだ。そのエンディアンすらコンパイラのオプションで大抵は、マシンの
自然体に関わらず、Big/Littleのどちらでも読み書きできるのが普通だ。

ゆえに、『エンディアン』をマニュアルで調べれば、チミの問題は解決する。

258:デフォルトの名無しさん
07/04/11 19:09:47
インテルのフォートランコンパイラifortを使ってコンパイルしています。
CPUはOpteronなのですが、-xWオプションをつけないとSSEやSSE2は全く使用されないのでしょうか?

259:デフォルトの名無しさん
07/04/12 00:11:39
>>258
デフォでは使用されなかったと思う。詳しくはHelpで確かめよ。

SSEは厳格にはIEEE754の規約を満たしていないせいか、
あるいは全バイナリー対応のためか?w


260:デフォルトの名無しさん
07/04/12 00:22:30
Mac用は-xPがデフォルトの様だから互換性のためなんでしょうね。

261:デフォルトの名無しさん
07/04/12 13:06:13
FORTRANで書き出したバイナリファイルの読み出し方は分かるのですが、
実験機器がはき出したバイナリファイルを読み込むにはどうしたらよいのでしょうか?
なお、FORTRANで書き出すと、今使っているマシンでは先頭に4バイトほどの識別データが付加されますが、
実験機器がはき出したバイナリファイルにはこれがついていません。


262:デフォルトの名無しさん
07/04/12 17:43:14
>>261
ファイルをOPENするときに form="binary" としたり (処理系依存)
access="stream" としたり (Fortran 2003)

263:デフォルトの名無しさん
07/04/13 02:58:02
>>261
DIRECT形式を使って書くという裏技がある。昔のコンパイラでやるにはこれか。

ただ今は>>262の書くように、F2003拡張の先取りでSTREAMモードが大抵あるのでそれを使うがよろしい。


264:デフォルトの名無しさん
07/04/13 21:39:29
面倒なのでCで読み込む。

265:デフォルトの名無しさん
07/04/25 16:57:37
STDCALL の関数を g95 から呼びたいのですが、どうするのがよいでしょうか。

266:デフォルトの名無しさん
07/04/25 21:35:56
>>257
>1バイトは8ビットに統一された。

以前は統一されてなかったの?

267:デフォルトの名無しさん
07/04/26 22:37:47
>>266
あらゆるものが統一されていなかった。


268:デフォルトの名無しさん
07/05/03 14:36:07
今は統一されているといっていいのか?

269:デフォルトの名無しさん
07/05/03 16:02:16
まあ昔に比べれば

270:デフォルトの名無しさん
07/05/07 06:10:19
fortranで
webの結果を得ることできる?

get URLリンク(www.hoge.com)
を行い
結果をファイルに残す

あるいは、戻り値 403 not found
かどうかの判定とか。
Linux *g77
です

271:デフォルトの名無しさん
07/05/07 06:15:37
できなくはないだろうけど、
別の言語を使った方がいいと思う。

272:デフォルトの名無しさん
07/05/07 08:01:53
うわ、返事はや!
ありがとう。

自由自在に書ける言語がFortranとVB系以外のbasicなんです。
検討してみます

273:デフォルトの名無しさん
07/05/07 10:14:30
>>270
表示するだけならこんなので動くと思うけど

web2console.py
----------------------------
#!/usr/bin/env python
import urllib

conn = urllib.urlopen("URLリンク(www.yahoo.co.jp))
data = conn.read()
print data
----------------------------

F2PY - Numerical Analysis with Python
URLリンク(scipy.dip.jp)
URLリンク(nihongo.homeip.net)

じぶんは fortran は使ったことないからよくわからないけど
fortran ともそれなりに相性も良いみたいだね…(正直よくわからん)

274:デフォルトの名無しさん
07/05/07 11:18:23
>>270
外人がg77でのやり方を解説しているページがある。
環境変数を使って変数のやり取りをしていた。

LINKがみつかんねーw

fortran cgi でググルと色々引っかかる。


275:274
07/05/07 11:24:19
>>270
ああごめん。
FORTRANでサーバー側のCGIを書くのかと思ったが逆か。
クライアントを書きたいのか。

Windows なら WinInet.Lib+Intel Visual Fortranでなら出来る。
Linuxでもその類のDLLをFortranから呼び出せば何とかなるのではないか?



276:デフォルトの名無しさん
07/05/08 13:33:07
【科学技術計算】C言語とFortran どっちがいい?
スレリンク(sim板)

629 名前:名無しさん@5周年[] 投稿日:2007/05/08(火) 07:49:36
FortranかCかなんて議論は馬鹿げているよ。
Fortranなんてもはや化石言語。
だから消去法でCを選ぶ他ない。
しかし君たち。
実はもっとも身近で簡単な科学計算用言語があるのだ。
それがExcelVBA。
答えは案外足元にあったわけだ。
いわゆる青い鳥ってやつだね。




(#^ω^)



277:デフォルトの名無しさん
07/05/08 21:42:56
fortranプログラムの作成依頼はこのスレでできますか?

278:デフォルトの名無しさん
07/05/08 23:37:48
>>277
くだすれFORTRAN(超初心者用)その2
スレリンク(tech板)

どちらかというと、こっちの方がやさしくしてくれるお

279:270
07/05/09 23:38:08
皆さんありがとう。

>>273
始め、pythonで「????」でしたが見ると便利ですね。
でもpythonのファイルIO知らないので調べてみます。

>>274
ええ、見る方です。でもcgi検索して見るとfortranでもcgi出来ますね。

FORTRAN(超初心者用)その2
でも聞いてみて良いのかな?

280:デフォルトの名無しさん
07/05/18 23:30:44
Visual Fortran 9 を使用されている方いますか?
購入を検討しているのですが,アクティベーションあるんかな?

281:デフォルトの名無しさん
07/05/19 00:48:26
>>280
ネット経由の認証も出来るし、ローカルでのライセンスファイルでの認証も出来る。
どちらかというとネット経由の認証へ移行したいと思っている節がある。

あと1年でライセンスは切れる。切れた後はインストールできない。再インストールも出来ない。
実行には問題ない。ただしバグフィックスは再インストールになるので出来ない。
うっかり消してしまうと、再インストールできない。

282:デフォルトの名無しさん
07/05/19 01:03:43
ライセンスって単年契約なんですか?恒久権利だと思ってた・・・やっぱり今のVisual Fortran 6使い続けようかなあ.

283:デフォルトの名無しさん
07/05/19 08:28:45
> 切れた後はインストールできない。再インストールも出来ない。
そんなことはない。
購入後、1年間の間にリリースされたマイナーバージョンアップまでは自由に使えるよ。
使えないのは、それ以降にリリースされた最新版とか、プレミアサポートとか。


284:デフォルトの名無しさん
07/05/19 13:01:56
ネット経由の認証がいるってことは,複数台のPCにインスコして同時に計算走らせるとかできないわけですよね?
(認証なくてもやっちゃだめなんだろうが)
結局この理由でofficeも2000,visual studioも6のまま・・・
Fortranも6でいいかー

285:デフォルトの名無しさん
07/05/19 14:30:29
>>284
契約したライセンス数を超えて(または契約した使用者以外が)「コンパイラ」を使うのはいけないが
コンパイルして出来た実行プログラムの方は配布が可能。
(評価版や非商用ライセンスの場合は別だが)

またそのプログラムの実行にIntelが定めた再配布可能ファイルが必要なら、それも含めて配布可能。

ちなみにシングル・ユーザー・ライセンスの場合、契約者個人が同時に1つしか使用しないのであれば
コンパイラを複数のPCにインストールすることも可能。

URLリンク(www.intel.com)

286:デフォルトの名無しさん
07/05/19 23:06:03
誰かeclipseを使ってフォートランのコンパイルしている人いませんか?
photranの使い方がわからなくて..

287:デフォルトの名無しさん
07/05/20 13:33:53
>>286
それをやろうとして、諦めたw
設定いろいろいじろうと思うとeclipseまで勉強しなきゃならなくて時間の無駄だと感じた。

設定できたら報告よろ!

288:デフォルトの名無しさん
07/06/09 00:09:26
空気読まずに言うけどフリーのリンク先にある

FTN77/95のインストール方法
URLリンク(www.coastal-env.k.u-tokyo.ac.jp)

95のインストール先とマニュアルの画面が全然違うのは仕様?

289:デフォルトの名無しさん
07/06/09 02:12:36
通りすがりの者だが・・・・
>>288
提供先がSalfortから関連会社のSilverFrostに移行したときに変わってしまったのさ。


290:デフォルトの名無しさん
07/06/09 17:02:26
>>286
fortranでeclipse使ってるとこってあるんかな?
fortranって科学技術系が多いし、他の一般業界の言語やコードに比べたらシンプルだからね・・
中にはviでガシガシの人もいるし

291:デフォルトの名無しさん
07/06/09 19:33:49
>>289
サンクス。
じゃあ、あの状態でFortran95をダウンロードするにはどうすれば……?

本来なら説明しているサイトに報告すべきだと思うけど、
SilverFrostでのダウンロード方法を知っている人がいたら教えてください。

292:デフォルトの名無しさん
07/06/09 19:55:34
>>291
>>141

293:デフォルトの名無しさん
07/06/11 10:43:46
>>292
サンクス。
上段のページからダウンロードして後、VisualStugio2003をインストールしてたから同じく2003でインストールしました。

あとは複数あるシステムの中で、Fortran教科書に載ってるプログラムを書いて実行するのに選ぶべきシステムとその方法さえわかれば……

294:デフォルトの名無しさん
07/06/13 22:07:13
Fortressについてどう思うよ?

295:デフォルトの名無しさん
07/06/13 23:23:23
オワタ

296:デフォルトの名無しさん
07/06/14 00:49:28
すいませんが、FORTRAN77で、二階非線形常微分方程式をルンゲ・クッタ法を使ってプログラムを組みたいのですが、
一階なら分かったのですが、高階になるとよくわかりませんorz


297:デフォルトの名無しさん
07/06/14 00:54:41
>>296
それはプログラムがわからないの?
それとも数学がわからないの?

298:デフォルトの名無しさん
07/06/14 23:34:11
>>296
ルンゲ喰ったなんて、数値計算の本の、頻出事項じゃんか。
図書館に行ったほうが2chに来るより効率的。


299:デフォルトの名無しさん
07/06/15 00:09:50
ルンゲ警部を喰っちゃうなんて・・・

300:デフォルトの名無しさん
07/06/15 10:39:29
ペッパーランチでペッパー警部に食われました(正笑

301:デフォルトの名無しさん
07/06/15 13:00:52
>>300
アッー?

302:デフォルトの名無しさん
07/07/07 05:25:58
FortranからLinuxのsystem関数を利用して出力したモノを
変数に納めたいのですが、どうすればよいのでしょうか?

例えば

call system('echo 2')
でa.outを実行すると正常に2が出力されるのですが、
その2を変数に代入しようと。

integer var
call system('echo 2',var)

のようにしてみたのですが、エラーになりました。
character型も試してみたのですが、やはりエラーになりました。

303:デフォルトの名無しさん
07/07/07 12:43:32
>>302
コンパイラなによ。

ま、マニュアル読んだほうが早いと思うけど。

304:デフォルトの名無しさん
07/07/07 13:22:12
>>302

call system('echo 2 > hoge.txt')
open(10, file='hoge.txt', status='OLD')
read(10,*) var
close(10,status='DELETE')
write(*,*) var


305:デフォルトの名無しさん
07/07/07 14:54:26
>>303
pgiコンパイラだそうです。
ちょっとマニュアル探してみます。

>>304
コンパイラの問題なのかDELETEのところがちょっとエラーでできなかったのですが。
close(10)としたら変数に2を入れることができました。
ありがとうございました

306:デフォルトの名無しさん
07/07/07 14:57:11
あ、すみません。
>>304でできました。

307:デフォルトの名無しさん
07/07/07 15:06:41
あと、System関数を利用せずに外部ファイルの行数を取得する方法ってありますか?
例えばhoge.txtの中身が
----------
apple
123 456
pie
----------
であれば、そのファイルの行数を読み取りinteger型の変数に3を
代入するという操作を行いたいのです。

もちろんSystem関数の利用でできたのですが、それだとプラットフォームが
Linuxに限られてしまうので、Fortranでファイルの処理を行うようにすれば
各種プラットフォーム上でも使えるソースになって便利だなと思いました。
この操作に限らず、Fortranでファイルを編集するための勉強はどのような
ところで勉強すればいいのでしょうか?

308:デフォルトの名無しさん
07/07/07 17:41:38
>>307
普通に空readでend of file が出るまでまわして、変数に一足していけばいい。


309:デフォルトの名無しさん
07/07/09 08:19:26
>>308
できました。
[Fortran]$ cat num_of_tmp
sfejoiwa
seflkew
sfepo pweko
wsefk

sefpowe posef

sdf s;ef
[Fortran]$ ./a.out
x = 8

ありがとうございました(><)

310:TESTちゅう
07/07/12 17:03:55
フィボナッチのプログラムが急ぎで必要です。だれか教えてください!!!

311:デフォルトの名無しさん
07/07/12 21:18:17
      WRITE(*,*) 'N='
      READ (*,*) N
      DO 10 I=1,N
      IF0=IF1
      IF1=IF2
      IF ( I.GT.2 ) THEN
            IF2=IF1+IF0
      ELSE
            IF2 = 1
      END IF
   10 WRITE(*,*) I,IF2
      STOP
      END

312:デフォルトの名無しさん
07/07/17 22:54:57
Numerical Recipes の第3版が出るらしいな。


313:デフォルトの名無しさん
07/07/18 03:14:53
確かに本家の看板に出てるね
URLリンク(www.nr.com)
情報さんくす


314:デフォルトの名無しさん
07/07/18 03:17:27
↑ってC++がメインになってないか(笑)
微妙にスレちがいかも

315:デフォルトの名無しさん
07/07/18 23:37:05
>>314
ホントだ
裏切られた!www
怪しからん奈!


316:デフォルトの名無しさん
07/07/20 03:17:16
>>314
FORTRAN版は出す気あるのかな?

初版の非Fortran版はトランスレータにかけただけのような
gotoバリバリコードで出版していたのだから、無茶なFortranコードでも
著者は別に気にしないよなw


317:デフォルトの名無しさん
07/07/20 04:09:06
>>316
どうだろ?C++の機能をつかって実装しているようだから、
以前の FORTRAN -> それ以外の言語 の移植とは訳が違いそう

NRの2nd(1992)から3rd(2007)への以降に15年掛かっているね
たぶん3rdに取り組み始めた頃は、まだC++が輝いてたんだろ

複素数の取り扱いが、クソややこしいテンプレートになっていたり、
整合配列に問題があったりと、C++は数値計算に向いてない
c99の方がずっとマシ(ただし、こっちは対応コンパイラがPGIぐらいしか無い)
Fortran90かFORTRAN77のまま、アルゴリズムだけ進歩させた方が良かったと思うよ 

NRのメリットは安い値段とソースコード付きってとこだったが、
C++のみになると、あんまりメリット感じないね C言語だとGSLとかあるしね
どうせ自由に配布できないんだったら NAG, IMSLの方がいいんじゃないか?
少なくとも漏れはNR 3rdイラネ

318:デフォルトの名無しさん
07/07/20 13:01:52
>>317
あれれ2版からも15年たってるのかよw

しかしアルゴリズムの進歩版とかのってはいないのかね?
第2版では、初版に無かったnon-numerical algorisms の章が一番役に立ったw

CRC16とかHuffman符号化とか、情報工学的なアルゴリズムが、数値計算との垣根を乗り越えてきた
タイミングだったのか。


319:デフォルトの名無しさん
07/07/20 23:48:31
今日本屋に行ったらこんな本があった。

Fortress言語―マルチコア時代の並列化プログラミング言語

もう出すのかよw 
ざっと見た感じ糞本だったがw



320:デフォルトの名無しさん
07/07/22 12:21:50
まったく話を聞かないしなぁ>Fortress
スラドでも議論があったけど、将来も普及する見込みは厳しそうね

321:デフォルトの名無しさん
07/07/23 01:55:57
>>320
よく知らんが、SUNのマルチコアCPUに大量のスレッドをばら撒いて並列計算させるという
コンセプトの言語なんだろ?

アイデアとしては面白いんで内科医?

Haskellっぽい仕様も入れているらしいし、次世代の踏み台にwww


322:デフォルトの名無しさん
07/07/23 10:15:40
>>321
オープンにしたからSUNの石専用ってわけじゃないのでは?

確かに踏み台ぽいなw

323:デフォルトの名無しさん
07/07/23 23:44:19
Thoughts on Language Design
New challenges require new solutions
Guy L. Steele Jr.
URLリンク(www.ddj.com)

まぁ、それなりに考えがあってFortressを作っているようだ。
Sunがアメリカのスパコン助成から脱落したので、さいさきは暗いだろう。

日本(富士通)で拾う振りをして、アメリカに疑心暗鬼を起こさせるのが吉だと思うwwww

324:デフォルトの名無しさん
07/07/24 09:34:34
>日本(富士通)で拾う振りをして、アメリカに疑心暗鬼を起こさせる
スゲーwwww
策士っすねw

325:デフォルトの名無しさん
07/07/25 00:35:12
>>324
最近の最速スパコン競争は、なんか土建屋の既得権争いっぽいから、
釣りに走って遊ぶほうが、たしかにいい気がするw


326:デフォルトの名無しさん
07/08/12 22:19:24
「スカラ型とベクトル型を両輪に」,次世代スパコン構成の狙いを語る
URLリンク(techon.nikkeibp.co.jp)



327:デフォルトの名無しさん
07/08/13 21:56:30
>>326
53. 次世代、ペタスケール、T2K (2007/8/12)
URLリンク(grape.mtk.nao.ac.jp)

すっぱいぶどうの更新w

328:デフォルトの名無しさん
07/08/19 02:10:13
Cygwin上で動くfortran90コンパイラってありませんか?
すみませんが どなたか情報をお持ちでしたら教えて下さい。

329:デフォルトの名無しさん
07/08/19 18:45:37
>>328
以前
URLリンク(g95.sourceforge.net)
を試したら動いた

330:デフォルトの名無しさん
07/08/19 19:08:37
>>328
URLリンク(g95.sourceforge.net)

でも、なんでパッケージに入ってないんだろう?

331:デフォルトの名無しさん
07/08/19 20:03:47
CygwinのGCCが4.x.xになったらこっちが入るからじゃないかな
URLリンク(gcc.gnu.org)


332:デフォルトの名無しさん
07/08/19 20:17:52
>>329 - 330
ありがとうございまs。試してみます!


333:デフォルトの名無しさん
07/08/23 23:01:43
Steve Lionel のコラムで、FORTRANで書かれた(多分)世界初のアドベンチャーゲームネタが出ている。
URLリンク(softwareblogs.intel.com)

334:デフォルトの名無しさん
07/08/23 23:11:05
マジカヨ

335:デフォルトの名無しさん
07/08/30 23:10:18
Intel fortranて、バージョン10になってから、MKLがおまけで付いてくるようになったな。
MKL買った人カワイソス=おれ

336:デフォルトの名無しさん
07/08/30 23:28:08
おまけと言ってもMKL有りとなしじゃ、けっこう価格差が…

337:デフォルトの名無しさん
07/08/31 00:31:24
>>335
忘れたがMKLのライセンスからをコンパイラのライセンスに転換できて2個分になったんじゃ中田ッけ?
ガセだったら酢マンコw


338:335
07/08/31 01:05:11
9→10へのバージョンアップのときに、スタンダード版からプロフェッショナル版への変更が
安く(1万円くらい)できただけだった。いずれにしても、MKL単体で買って損した気分だ。

339:デフォルトの名無しさん
07/09/07 13:02:07
URLリンク(www.nag-j.co.jp)

NAGのcompiler IntelMac版登場 トライアルができる

340:デフォルトの名無しさん
07/09/08 21:10:14
intel fortran compiler for windows 9.0 を使用していて

URLリンク(icl.cs.utk.edu)

の手順に従いlapackとblasを導入しようとしてますが

1. C:\Program Files\LAPACK_3.1.1_for_Windows\LAPACK 3.1.1 for Windows\ に移動
2. set_environment.bat x86 Intel
 と入力(Blasはデフォルトを使用)
4. nmake /nologo all

とすると

xilib.exe /OUT:"C:\Program Files\LAPACK_3.1.1_for_Windows\LAPACK 3.1.1 for Windows\blas.lib" *.obj
が実行されているときに次の旨の表示がでてうまくいきません。
"lib"を実行しています
"lib"がみつかりません

>>197
の slim に相当する形でxilibが使われているものだと思うんだけど、どう使用すればうまくいくのかわかりません。
ちなみにxilib に代えて xilink /lib を使っても同じです。

blas.lib lapack.lib を作るにはどうすればよいでしょうか?

341:デフォルトの名無しさん
07/09/09 13:51:43
質問です。
fortran95を使っているのですが

do i=1, size(A)
A(i) = hoge + i * fuga
end do

みたいのを

A = hoge + [なにか] * fuga

のようにして添字を配列演算に用いて一行にまとめられないのでしょうか?
それと

do i=1, size(A)
hoge = 複雑な式
A(i) = hoge * fuga
B(i) = hoge * piyo
end do

のような場合、

A = 複雑な式 * fuga
B = 複雑な式 * piyo

と書く事で二度計算されてしまうのは無駄ですし
式が長くなって読みづらいのですが、配列演算に対する一時変数のようなものは無いのでしょうか。

342:デフォルトの名無しさん
07/09/09 21:33:40
>>340
Visual Studioに含まれているlib.exe/link.exeへのパスが通ってない。

通常はset_environment.batから呼ばれている
"%IFORT_COMPILER91%\IA32\Bin\ifortVars.bat"
の中でVSへの設定も行われるはずだが
環境変数 IFORT_COMPILER91 はIVF9.1用の設定なので
IVF9.0なら IFORT_COMPILER90、10.0なら IFORT_COMPILER10
に変更する必要がある。

ifortVars.bat を実行してもlib.exe/link.exeが使えない場合は
ifortVars.bat を修正するか、予めVS中の VSVars32.bat か
VCVars32.bat を実行しておく必要がある。

343:デフォルトの名無しさん
07/09/09 23:48:01
>>341
A = hoge + (/(i,i=1,size(A))/) * fuga

344:デフォルトの名無しさん
07/09/10 17:01:44
>>343
ありがとうございます。こんな素直な方法があったんですね。
後者については今のところ難しいでしょうか。

あと配列演算に使える関数が作れるといいのですが、これも難しいですかね。
配列をもらって配列を返すようなタイプだと一次的に戻り値用の配列を確保しなくてはいけないので現実的ではないんです。

345:デフォルトの名無しさん
07/09/10 19:36:29
>>344
しかし>>343のようにやると、いったん暗黙のテンポラリ配列が生成されて
それから積が取られることになるので、DO LOOPの展開式とは等価ではない。

>あと配列演算に使える関数が作れるといいのですが、これも難しいですかね。
elementalだったかpureオプションつければ可能だった気がする。
調べてくれw



346:345
07/09/11 01:37:16
>>344
MODULE m_elemental
CONTAINS
ELEMENTAL REAL FUNCTION watson(x)
IMPLICIT NONE
REAL, INTENT(IN) :: x
watson = x**2
RETURN
END FUNCTION watson
END MODULE m_elemental
!
PROGRAM holmes
USE m_elemental
IMPLICIT NONE
INTEGER :: i
REAL :: x(10), y(10)
DO i = 1, 10
x(i) = REAL(i)
END DO
y = watson(x)
PRINT '(10f7.1)', x
PRINT '(10f7.1)', y
STOP
END PROGRAM holmes


1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0
1.0 4.0 9.0 16.0 25.0 36.0 49.0 64.0 81.0 100.0
Press any key to continue . . .

elemental ⊂ pure だったw


347:デフォルトの名無しさん
07/09/13 12:33:44
>>345
試してみましたが確かにテンポラリ配列作られてますね。
elementalは便利ですね。ありがとうございました。

348:デフォルトの名無しさん
07/09/13 19:16:35
>>347
a(1:n) とか部分配列の類もテンポラリを作ることが多いので、
ソースコードが短くなるからと調子に乗って、DO LOOPの代わりに乱用すると、
思わぬところでスタックオーバーフローとかで憤死する。

ベクトル機の場合は一旦テンポラリ配列に展開したほうが早くなるだろうが、
パイプラインのスカラー機の場合はDO LOOPのままの方が早いのではないかと思う。

まぁ勝手に最適化するかも知れんし、トレードオフが色々あるので、なんとも胃炎が。


349:デフォルトの名無しさん
07/09/18 15:26:34
age

350:デフォルトの名無しさん
07/09/18 22:17:04
ageも満足にできんのか!!

351:デフォルトの名無しさん
07/09/21 00:26:58
intel Visual fortran 10.0 日本語版 発売

352:デフォルトの名無しさん
07/09/21 01:08:09
「日本語版」つーてもFortranの場合、コンパイラのメッセージは英語のままだからありがたみがないなぁ。
日本語マニュアルは別途ダウンロードできるし。

353:デフォルトの名無しさん
07/09/21 23:56:06
質問です。
fortran77を使っているのですが、コンパイラ時にはノーエラーで、実行のときランタイムエラーと表示されるのですが、

WARNING - Label 170 has not been referenced
NO ERRORS,1 WARNINGS [<MAIN@>FTN77 Ver 4.03]
Creating executable: sinchouryuu.exe
*** Unable to create executable file: sinchouryuu.exe

実行可能なファイルの作り方を教えて下さい。

あと、実行時のエラーの
00401000 MAIN# [+174f]

て書いてあるのですが、どういう意味でしょうか?


因みに170行は
166:* cnbの転置行列
167:
168: do 70 j=1,10
169: do 80 i=1,6
170:
171: cbn(j,i)=cnb(i,j)
172:
173: 80 continue
174: 70 continue

こんなんなんです。

354:デフォルトの名無しさん
07/09/22 00:01:51
ヤコビアン行列がつくれません!!


どういった風に作るのか、作り方をどなたかおしえてください!

355:デフォルトの名無しさん
07/09/22 21:55:29
>>353
>fortran77を使っているのですが
それは「規格の種類」であってこの場合、それよりも使用しているコンパイラや
OSの名前を書いてくれた方がありがたい。

> WARNING - Label 170 has not been referenced
「170行」ではなく文番号の170が参照されていない。

> *** Unable to create executable file: sinchouryuu.exe
> 実行可能なファイルの作り方を教えて下さい。
おそらくsinchouryuu.exeを実行中で書き換えられないのだろう。
すべてのsinchouryuu.exeを終了させてからコンパイルするのだ。

> 00401000 MAIN# [+174f]
> て書いてあるのですが、どういう意味でしょうか?
MAIN#(プログラムの開始位置)から 5919バイト目の位置でエラーが発生した。
エラーの種類はFTN77の場合、ダイアログの上の方に書いてある。

>>354
つ google

356:デフォルトの名無しさん
07/10/16 01:37:18
FTN95をインストールしました。
plato3(日本語不使用)でコンパイル実行可能でした

fcpadを利用しては出来ないですか?
コンパイラパス
C:\Program Files\Silverfrost\FTN95\ftn95.exe
と入力
*** Invalid executable file suffix - only .EXE is permitted
とでてコンパイルできない。

秀丸でマクロを使おうかと思いきや、
PATHが通らず失敗。

誰か、ご指導願います・・・。

357:デフォルトの名無しさん
07/10/16 02:17:02
>>356

くだすれFORTRAN(超初心者用)その2
スレリンク(tech板:355番)

358:デフォルトの名無しさん
07/10/23 23:51:02
Intel Visual Fortranで/Qparallelを使って並列化したところ、
あまり高速化しませんでした。プログラムによってはむしろ遅くなることも。
サポートに問い合わせたら、こちらでも2コアでは確かにスピードダウンが
確認できましたが、4コアではスピードアップが確認できました、と回答。
まるでリフォーム詐欺のようだ・・・
単純な行列×ベクトルのみの計算でも最大30%アップ程度なのだが
並列計算ってそんなものですか?


359:デフォルトの名無しさん
07/10/24 01:28:55
マシンの構成によるし・・・
QuadコアなCPU?

360:デフォルトの名無しさん
07/10/24 01:48:48
>>358
自動並列化は、そんなもんだよ。スパコンでも同じ。

まぁ並列化よりSSEの方が簡単に簡単にスピードが出るようになる。

361:デフォルトの名無しさん
07/10/24 03:04:44
SSEは適当に-O3とかするだけでも有効になるはず(多分

Qparallelって要するに自動でOpenMPディレクティブ入れてるだけだよね?
最近のOpenMPはどうなんだろう・・・古い古いとよく言われるけど、今後のマルチコアCPUで通用しなくなるのだろうか

362:デフォルトの名無しさん
07/10/24 09:34:52
>>360
行列積とか素朴な演算ならSSEを入れれば単精度は4倍速、倍精度は2倍速とか、素直に早くなる。
コンパイル時のメッセージにベクトル化云々と出ればSSE使ってる。


363:358
07/10/24 11:02:26
ありがとうございます。
現在はOpteron1210(DualCore)を使っています。
インテルからはまるでQuadCoreのCPUを買えというような回答だったので、
詐欺商法のように聞こえました。
SSEは使っています。
CPUのクロック数があがらなくなった現在では並列計算にしないと
スピードアップが望めなくなると思い、並列を検討しています。
自分の勉強不足も大きいのですが、誇大広告にがっかりです。


364:デフォルトの名無しさん
07/10/24 11:20:56
Intelのサポートに聞いたの?
こういうのって大抵まともな(役に立つ)回答は期待できないものでは?
MSよりはマシだけど・・・

結局ググるなり2chなりどっかのコミュニティなりで探すしかないんだよね

365:358
07/10/24 12:58:59
>>364
バージョンアップ込みで年間数万円(いくらか忘れた)のサポート料を払っているので、
使ってみたのですが、本当にがっかりです。
サポート料金を10倍にしてもいいので、まともな回答を期待したいです。
ググってますがなかなかいい情報が出てこないです。
私の最高は3割増が最高で、実務的な計算ではほとんど役に立ちません。
あれこれ考えている間に、時間をかけて計算したほうがいいのかと思えてきました。


366:デフォルトの名無しさん
07/10/24 12:59:13
>>363
いや詐欺ではない。
既存プログラムを自動並列化して場合、スピードダウンするのはごく普通の現象。
はじめから並列化を前提に設計して、かつ指示行を入れまくってもいまいちなのは、実は常識だw


367:デフォルトの名無しさん
07/10/24 15:53:27
>>365
スレ違いな話だけど、本気のガチンコサポートをするなら商売する側としては10倍でもペイしないと思うけどね・・・
そこはしょうがないっすよ

368:358
07/10/24 16:19:19
日本語がおかしかった。失礼しました。

>>367
あるソフトベンダーは年間140万円のサポートでいたれりつくせり。
(何が原因でしょうね?と問い合わせたら修正されてかえってくる)
そこは最優先でサポートを継続している。
別の40万円のところは、それなりだけど解決のヒントはくれる。
(これらは言語ではなくてアプリケーションですが)
VisualFORTRANは「開発に問い合わせます」のまま放置・・・。
ただ、今回の件については決して特殊な質問ではなく
広告にうたってある性能が出ないよに始まり、
そちらのもっている理想的なプログラムでいいから出してくれといっても
出てこない。いったいどんな質問だったら答えることができるのだろう。


369:デフォルトの名無しさん
07/10/24 19:20:01
>>368
言語のサポートの場合、コンパイラのバグに対する対応の早さだと思う。
そういう点でIntelはパッチの出方は早い。

あとサポート掲示板には、様々なレベルの質問が出ているのでこっちも見てみたらいいかも。
URLリンク(softwarecommunity.intel.com)

370:デフォルトの名無しさん
07/10/25 20:49:16
NEC:世界最高性能のスパコン「SX-9」を発売
 NECは25日、世界最高レベルの計算能力を持つスーパーコンピューター「SX-9」を開発し、
発売を始めたと発表した。1秒160億回の計算能力を持ち、最大限の512台を連結して運用し
た場合、1秒839兆回の計算が可能なシステムが完成する。

 新開発の高性能CPU(中央演算処理装置)を搭載した。国内外の大学や研究機関などが、
気象予測やロケットの設計、核融合のシミュレーションなどに利用することを想定している。
価格は1億5000万円。

 現在の世界最速のコンピューターシステムは米IBMが開発し、秒間280兆回の計算ができる。
NECは、10年の完成を目指す世界最速スパコンの開発計画を、日立製作所や富士通と進めて
いる。

URLリンク(mainichi.jp)


371:デフォルトの名無しさん
07/11/02 00:54:14
プレステ3用FORTRANが動くようになった?
URLリンク(www.ibm.com)

URLリンク(cell.fixstars.com)



372:デフォルトの名無しさん
07/11/02 13:20:39
このままPS2クラスタのようにお遊びで終わるのかと思ったが、そうでもないのか・・?
でも低価格の普及ゲーム機を使えるってのがそもそものウリなんだから、今のXboxやwiiに苦戦してる状況じゃ厳しいんじゃないか?

373:デフォルトの名無しさん
07/11/02 23:44:01
東北大学が NEC の世界最速スパコン「SX-9」16台導入へ

著者: japan.internet.com 編集部 プリンター用 記事を転送

NEC は1日、東北大学情報シナジー機構情報シナジーセンターから、世界最速の
ベクトル型スーパーコンピュータ「SX シリーズ モデル SX-9」16台からなる大規模
スーパーコンピュータシステムを受注したことを発表した。新システムは、最大性能
26.2テラフロップス(テラフロップス:1秒間に1兆回の浮動小数点演算能力)に達し、
SX シリーズでは国内最大規模のシステムになるという。

同センターは、全国の大学や研究機関の研究者が利用する共同施設として、常に
最高・最新鋭の計算機システムを設置し、航空宇宙、地球環境、IT、ナノテク分野など
最先端の学術研究を支援。他では実行できない大規模シミュレーションの実行環境の
整備、専門的立場からの利用者プログラムの高速化推進を積極的に行い、長年研究者
にとって使い勝手の良いシステムが構築されている。現在、センターの利用者は学内外で1,500人を超えるという。

同センターは、その前身である大型計算機センターが1986年に SX シリーズを導入して以来、
ベクトル型スーパーコンピュータを採用している。
URLリンク(japan.internet.com)


374:デフォルトの名無しさん
07/11/06 14:00:01
阪大は10ノードらしい。
URLリンク(opentechpress.jp)

375:デフォルトの名無しさん
07/11/06 23:51:09
>>374
阪大センターはハッカーにやられてから、セキュリティーが強化されすぎて
使いにくくなったって知り合いが言ってたなー。
数年前のことだが。



376:デフォルトの名無しさん
07/11/23 16:54:32
休みだし、何か新しいネタを出してくれ。


377:デフォルトの名無しさん
07/11/23 17:13:22
暇ならFortran2003の処理系作ってくれ。

378:デフォルトの名無しさん
07/11/29 13:22:06
━━━━━━━━━━━━━━━━━━
     日本ニューメリカルアルゴリズムズグループ株式会社
      LAPACKサンプルプログラム集を公開開始
        LAPACKの利用がより身近に!
━━━━━━━━━━━━━━━━━━
日本ニューメリカルアルゴリズムズグループ株式会社(東京都中央区、以下
日本NAG)は、netlibにより公開されている線形代数パッケージ「LAPACK」の
サンプルプログラム集を公開いたします。これによりLAPACKを利用するユーザ
プログラムの開発が、よりスピーディーに行えるようになります。

公開URL:URLリンク(www.nag-j.co.jp)

また、日本ニューメリカルアルゴリズムズグループ株式会社ではより高性能な
LAPACKを、NAGライブラリ製品(NAG Fortran SMP Library)と共に提供してい
ます。

■ LAPACKについて
LAPACKはLinear Algebra PACKageの略で、netlib(URLリンク(www.netlib.org)
により公開されている高性能線形代数パッケージです。(日本語ではエルエー
パック、もしくはレイパックと発音されます)連立一次方程式、最少二乗問題
固有値問題、特異値問題を解くためのサブルーチンが含まれています。

■ NAG Fortran SMP Libraryについて
1500以上のサブルーチンから構成される並列科学技術、統計計算ライブラリ
です。1500のサブルーチンの中には並列化済みのLAPACKルーチンも数多く含ま
れており、LAPACKを利用するプログラムを1行も変更することなく、性能を
更に高めることが可能です。


379:デフォルトの名無しさん
07/12/06 03:09:49
初心者で最小二乗法のプログラムが組めないんですがどうすればいいんでしょうか?

380:デフォルトの名無しさん
07/12/06 10:03:56
>>379
図書館へ行って、データ処理の本を読むんだ!


381:デフォルトの名無しさん
07/12/06 19:03:50
データ読み込んでΣ求めて公式通り計算すればおk

382:デフォルトの名無しさん
07/12/07 22:27:07
>>379
いますぐ俺の口座に3000円振り込むんだ!

383:デフォルトの名無しさん
07/12/09 12:57:37
最小二乗法くらいならいくらでも転がってるだろ・・・

384:デフォルトの名無しさん
07/12/12 20:46:38
コンパイルはできたのですが、こんなメッセージが出てきて
途中で止まってしまいます。

Runtime error from program:c:\users\○○○○\documents\fortran\traje_0point.exe
Unknown floating point exception 91a1
Floating point co-processor fault at address 00401ef8


00401ad0 GOENTR [+0428]

00401000 main [+059c]


eax=00000b77 ebx=00000096 ecx=00000000
edx=0365fc38 esi=0365fcdc edi=0365fcd8
ebp=0365fc94 esp=0365fbd8 IOPL=0
ds=0023 es=0023 fs=003b
gs=0000 cs=001b ss=0023
flgs=00010206 [NC EP NZ SN DN NV]

00401ef8 fsqrt
00401efa fmulp st(1)
00401efc fmul [0044a2c0]

誰かこのエラーの意味ならびに解決策をご教授していただけないでしょうか?
ちなみに○○○○のところは私の名前が入ってるので、伏せておきました。


385:デフォルトの名無しさん
07/12/12 22:23:51
>>384
コプロの機械語の命令に対応物が無いのだから、コンパイル時に対象機種以外の命令が入るような
最適化をしたのではないかと思われる。
Pen4SSE2を無印ペンチアムで動かすとか、AMD用SSE命令をIntelで動かすとか。

そうだとすれば、最適化を落とすか、対象機種を明示的に指定するか、全機種対象にオプションを変えればいい。



386:デフォルトの名無しさん
07/12/13 01:46:25
>>384
ちょうどfsqrtのところ(00401ef8)で例外が起きているので
GOENTRの中でSQRT()に負の数を渡す等の不正な計算を
しているように見えるが。

387:デフォルトの名無しさん
07/12/13 16:42:46
386に一票


388:デフォルトの名無しさん
07/12/13 23:13:02
コンパイラは何だろう?


389:デフォルトの名無しさん
07/12/14 18:44:44
384です。
自己解決しました。
ありがとうございました。


390:デフォルトの名無しさん
07/12/15 00:00:18
>>389
今後、他の方の問題解決に役立つこともありますので
出来ればコンパイラの種類や解決方法を公開できる範囲で良いので
書き込んでもらえると有難いです。

391:デフォルトの名無しさん
07/12/22 10:35:08
CONTINUE

392:デフォルトの名無しさん
07/12/23 10:22:03
stop
end


393:デフォルトの名無しさん
07/12/25 23:51:38
最近の Intel Fortran って Traceback 機能が効いてなくね?


394:デフォルトの名無しさん
07/12/26 21:28:28
そうか?
-traceback オプションででちゃんとソース名と行番号まで表示されるぞ。

395:デフォルトの名無しさん
07/12/26 21:50:42
年明けでよかったらv10.1で試してみるよ

396:デフォルトの名無しさん
07/12/26 23:18:42
>>395
頼む。やってみてくれ。

漏れの環境だと、Windows+VisualStudioのビルドオプションでtracebackをつけているのだが、
9.1だとTracebackを吐くが、Ver10.0やVer.10.1だと普通の実行時エラーしか出ない。

ツール+オプションでコンパイラだけ変えているはずなので、オプション指定は共通だと思うのだが・・

397:デフォルトの名無しさん
07/12/27 01:48:12
手元のPCに10.1しか入ってないので他のバージョンと比較できないが、10.1の場合:

・VS上で[デバッグ] -> [デバッグ開始] (F5) で実行
 VSのデバッガ上でブレイクがかかる(「呼び出し履歴」にスタックトレースを表示)

・VS上で[デバッグ] -> [デバッグなしで開始] (Ctrl + F5) で実行
 トレースバック情報を表示せずに終了

・別途コマンドプロンプトを開いて実行
 トレースバック情報を表示して終了

となった。

398:デフォルトの名無しさん
07/12/28 01:15:49
>>397
㌧クス
確認・把握した。

漏れは近年デバッグ実行をまったく使っていなかったのと、エラーの出る奴は
コマンドラインから実行したことがなかったので気づかなかったw

多分、デバッグなし実行でもトレースバックは呼び出し履歴に強制リダイレクトされているが、
(error出力かなんかで)
この場合呼び出し履歴を見るメニュー項目が無いというバグなのだろうw 空想だがw




399:デフォルトの名無しさん
07/12/29 19:18:05
せっかくVisual Studioを使っているんだからデバッガも使おうぜ。

400:デフォルトの名無しさん
07/12/30 11:51:04
>>399
いや、昔は使っていたけど、なんというか、Fortran向きでない感じがして。
巨大配列のはみだしとかの検出に向いてないせいかな?


401:デフォルトの名無しさん
08/02/21 23:04:14
Fortran2003の処理系はよ出せ

402:デフォルトの名無しさん
08/02/21 23:57:48
>>401
IBMのコンパイラが一番進んでるって聞いたんだけど。
どうなのかね?
NAGも最近おとなしいし。


403:デフォルトの名無しさん
08/03/16 17:47:06
中規模疎行列連立方程式を解くのにFortran (Intel Win & Linux)を使おうとしてます。
小規模のものは、MATLABで、不完全三角分解前処理+BiCGSTABとかGMRESで解いてきました。
しかし方程式の大きさが1万を超えるあたりから、MATLABの組み込みソルバーのメモリ効率が悪くて
Fortranに戻ってきた次第です。

今、とりあえずMKLを入れてみて、その中でおすすめ(?)ぽかったPARDISOを使ったコードを組んでみて
サイズ5万くらいだとそれなりの速さで解いてくれることがわかりました。MKLのマニュアルには、普通の
疎行列用反復ソルバもあると書いてありますが、サイズが10万とか超えるようになると、PARDISOのような
疎行列直接法より、反復法の方が良いのでしょうか。

ほかのソルバライブラリでは、並列用PETScを試しています。ただ、解いている問題の係数行列が良くないのか
PETScだと収束しないことがあって、PARDISOやMKLを使おうかと思っています。

並列計算に移行するちょっと前くらいの大きさの疎行列連立方程式を解くのにおすすめのライブラリあったら
教えてください。

404:デフォルトの名無しさん
08/03/16 18:23:17
>>403

LAPACK か ATLAS ではダメ?そんな大規模な計算をやった事が
無いので、使い物になるかどうかは分かりません。

Linux ならパッケージでインストールできるはず。
ウチの Debian には aptitude で入れた。

405:403
08/03/16 18:44:46
>>404さん
レスサンクスです。

LAPACKとATLASも使ってます!MKLもLAPACKとかLAPACK95を使います。素晴らしい
ライブラリだと思います。
漏れがやりたいのは、疎行列連立方程式なので、LAPACKとかATLASをベースにして
そのラッパーが欲しいという感じです。LAPACKがすでにBLAS/ATLASのラッパーなわけ
ですが。。。



406:デフォルトの名無しさん
08/03/16 19:57:43
>>405
一般論からすれば、反復法か直説法かの選択は行列の内容によるんでね?
たちの悪い行列でなければ、反復法の方が有利だと思うけど。
直説法なら確実に有限回で終了することが分かっているが、反復法は直交性が数値誤差で腐るから
数学的にはおかしいが、現実には有限回で終わらなかったりするし。

行列が、ほんとにスカスカの疎か、singularに近いか、固有値が縮退しまくっているかとか、
その辺は使用者しか分からんから、なんともいえないような。

中身の分からん奴に機械的に使わせるから安全第一で行きたいとか、運用面にも依存するし。


407:405
08/03/16 21:00:50
>>406さん
仰るとおりっす。
とりあえず、自分一人で使う(研究)ので、ブラックボックス的で安全なものでなくとも
高速で、導入やコーディングが楽なのが希望です。

行列の状態は悪いです。4x4のブロックに分かれてる感じで、離散化を細かくやると
主対角が卓越して、対角ブロックのみに不完全分解系の前処理掛けると、
それなりに収束しますが、係数行列のサイズが増大するというジレンマです。

疎行列性は、5万x5万で、非零要素が0.1%切ってるくらいです。

ご指摘のように、問題依存性が高いので、なかなかご意見いただくのは難しいかなとは
おもっております。

とりあえず、みなさん、10万x10万くらいの疎行列連立方程式(差分法や有限要素法
なら頻出だと思います)を解くときって、どんな言語で、どんな実装してるんですかね。。。
反復法ソルバもSaadのtemplateの擬似コードから自分で書いてるのか、とりあえず
いろいろソルバライブラリ試してるのか。。。たぶん、疎行列の格納形式から考える必要が
あると思うので、そのあたり、Fortran使いの皆さんがどうやってるのか知りたいです。

408:デフォルトの名無しさん
08/04/05 16:41:03
Steve Lionel (Intel) のコラムで、A=B+CとA(:)=B(:)+C(:)での後者の問題点について
触れられていた。最近のインテルコンパイラは少し賢くなって、この程度の単純な場合なら
後者で無駄なテンポラリへのコピーをしないようになったようだ。

とはいえ、まだ問題ありと。
URLリンク(softwareblogs.intel.com)


409:デフォルトの名無しさん
08/05/17 02:00:17
EXCEL VBAである計算をしてたのですが、時間がかかって仕方がないので
PCにFORTRANをインストールしてfortranでプログラミングやろうかと
思っています。学生時代にFORTRAN 77を使った経験がある程度。
たしかMS FORTRAN ぱわーすてーしょんつうのも使ったことがあります。

専門家の多そうなこのスレで恐縮なのですが、現行のwinXPが入ってるPCで
動かすとよさげなFORTRANコンパイラは何がいいでしょうか・・。
計算内容は簡単なモンテカルロ法です。



410:デフォルトの名無しさん
08/05/17 10:28:08
>>409
URLリンク(www.nminoru.jp)
見れ。無料でとりあえず試すならcygwin+g77 かminGW+g77でしょう。
コストパフォーマンスならインテルコンパイラか。

どっちにせよ今EXCEL VBAで組めているのなら、今からわざわざFortranなんぞやらずに、
コンパイラ買う金で速いCPU買ったほうがマシだと思うけどね。

411:デフォルトの名無しさん
08/05/17 22:05:46
>>どっちにせよ今EXCEL VBAで組めているのなら、今からわざわざFortranなんぞやらずに、
>>コンパイラ買う金で速いCPU買ったほうがマシだと思うけどね。

レスありがとうです。
VBAって異常に遅いなと思っていたのですが、CPU速くすれば
FORTRANに追いつける程度なのか・・・。
そこらへんがよくわかってないので。


412:デフォルトの名無しさん
08/05/17 22:13:50
VB2008Expressでも使えば?
VBAよりはだいぶ速いしタダだよ

413:デフォルトの名無しさん
08/05/17 22:41:52
どのくらいスピードを要するのか知らんがVBはインタープリタだろう?
SilverFrostのFTN95とかフリーコンパイラも転がっていることだし、
77でよければWatcomとか昔なつかしいのもあるし、g77のWin用もあるし、
Fortranで素直にやらせてやれw


414:デフォルトの名無しさん
08/05/19 22:54:50
モンテカルロ法とかだとCPU周波数にほぼ比例して
速くなると思っておいていいのですか・・・


415:デフォルトの名無しさん
08/05/19 23:35:44
>>411
その速いCPUで動かせばFortranはもっと速くwwww
マジレスするとパラメータサーベイとかなら今の4コアとか使えばかなり効率はあがる。
さすがにBASICがFORTRANを超えることは無理。
自分のプログラミングスキルと作るコードの複雑さから、コーディングやデバッグの
労力を見積もって自分で判断すべし。別にfortranとかVBAに限らんけどな。

VBAはクソ遅いと聞くが使ったことがないのでなんともだが、
Free BASIC VS intel Fortran ではFORTRANが速いといっても高々数倍程度~コードによってはBASICが勝っている。
URLリンク(shootout.alioth.debian.org)
さすがにIntel Fortranはgnu fortranやJava, Mono Ruby なんかには楽々勝利。

>>414
だいだいYES。と言いたいがモンテカルロ法だけでは「乱数を使った確率的解法」ぐらいしか
わからんからなんとも。ランダムウォークのシミュレーションなんかだとCPU周波数に
比例してくれるようなコーディングができる場合が多い。

416:デフォルトの名無しさん
08/05/22 09:28:10
質問です。
salfordのf77を使っているのですがgoto文のループ内で配列を毎回初期化してたら
ループが途中で終わってしまうのですがどうしたら良いでしょうか?
何でか分からないのですが、初期化しないと繰り返してくれます。あと、以下のプログラムの文
で初期化がi,j共に31の繰り返しだと少し長く繰り返します。
goto文の中にはガウスの消去法のサブルーチンがあります。

100 continue
do i=1,32
do j=1,32
kakunou(i,j)=0
enddo
enddo

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
中略(以下このkakunou(32,32)に計算値を代入してます)
   途中でガウスの消去法のサブルーチンを呼び出してます
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
goto 100

すいませんが、どなたか教えてください。

417:デフォルトの名無しさん
08/05/22 11:44:17
>>416
正直質問の意味がよく分からん。
コンパイラのバグでなければ、ありえない症状なので
ここに書いてない部分がおかしいとしか言いようが無い。

f77の範囲ならg77とか他のコンパイラもあるのでそれで試してみるとか、
salfordのf95で試してみるとかして見るのも一つの案。


418:初心者
08/05/23 15:40:54
失礼します.
Compaq Visual fortran 6.1を使用しているものですが,
練習で以下のプログラムを書いたのですが,cpmpileを押し,
「.F」に変換しても「1errors」が出てしまいます.
なぜ出るのか教えていただけないでしょうか?
初歩中の初歩で申し訳ないです.

INTEGER WIDTH,HEIGHT,AREA
READ(5,*) WIDTH,HEIGHT
AREA=WIDTH*HEIGHT
WRITE(6,*) WIDTH,HEIGHT,AREA
STOP
END


419:デフォルトの名無しさん
08/05/23 23:58:39
>>418
エラーでないぞw

エラーメッセージをコピーしてみせてくれ。

420:初心者
08/05/26 12:58:26
デフォルトの名無しさん
早速の返信ありがとうございます.
エラーは以下のように出ます

--------------------Configuration: Text1 - Win32 Debug--------------------
Compiling Fortran...
C:\Documents and Settings\Admin\デスクトップ\Text1.F
f90: Severe: Invalid argument
... file is 'C:\Documents and Settings\Admin\デスクトップ\Text1.F'
Error executing df.exe.

Text1.obj - 1 error(s), 0 warning(s)


自分も以前OSが2000を使用していたときは何のエラーも出なかったのですが,
XPに変えてからエラーが出るようになったように感じます.

お手数をかけて申し訳ないです.
よろしくお願いします.

421:デフォルトの名無しさん
08/05/26 17:55:54
そりゃおめぇ,”デスクトップ”で実行しようたって無理な話よ
マイドキュメントに移動させな

422:デフォルトの名無しさん
08/05/26 23:07:33
>>421
お前ひどい奴だなwww

>>420
最初ッからそう言ってくれれば、すぐに答えられる。
それはCVF時代のFAQで、日本語のパスやファイル名が入るとコンパイラが対応していない
のでエラーを吐く。2バイト文字を使わないディレクトリでやればおk。



423:初心者
08/05/27 13:41:47
デフォルトの名無しさん
解決しました!!!
簡単な質問に親切に答えていただいてありがとうございました!!!

424:初心者
08/05/30 15:38:10
重ね重ね質問すみません.

fortranをまわした結果をexcelかテキストファイルに出力したいのですが,
どのようなコマンドを打てばよいのでしょうか?
また,長い結果(10万桁以上)をファイルに出力したいのですが,
御教授お願いします.

425:デフォルトの名無しさん
08/05/30 16:32:44
本人はできるだけ鼻につかない文章で丁寧に質問してるつもりなんだろうな。

426:デフォルトの名無しさん
08/05/30 17:59:33
>>424
それ、先生に聞いた方がいいよ。お金払ってるんでしょ。

427:デフォルトの名無しさん
08/05/30 19:26:16
program main
integer :: i = 1

open(unit=10,file='hoge.txt',status='unknown')

do
write(10,'(i4)') i
i = i + 1
end do

close(10)

end program main


「ご教示」を「ご教授」って書いちゃう男の人って。。。

428:デフォルトの名無しさん
08/05/31 00:19:23
>>427
hogeはUNIX系の連中の使う言葉だ。
メインフレームの誇りを持つFortrannerは使っちゃいかんw



429:デフォルトの名無しさん
08/05/31 11:14:03
foo bar baz は?

430:デフォルトの名無しさん
08/05/31 16:02:34
鋼鉄の撃墜王

431:デフォルトの名無しさん
08/06/01 18:28:25
>>415
これみるとIntel C/C++使ったほうがいいみたいだな

432:デフォルトの名無しさん
08/06/01 19:41:18
>>431
単にFortranの書き方知らないだけだべ。


433:デフォルトの名無しさん
08/06/01 19:44:49
つまりFortranを知ってる人が少ないと。
技術者を集めにくいと。

434:デフォルトの名無しさん
08/06/04 01:08:45
>>431
INTEL Fortran と C++ はバックエンド部は共通だろうから、そもそもそれほど大きな違いが
出るはずが無い。


435:デフォルトの名無しさん
08/06/07 20:13:52
慶大の生協運営の書店に
FORTRAN入門(FORTRAN77入門)が普通において有って驚いた

436:デフォルトの名無しさん
08/06/07 20:25:26
>>435
浦昭二君、名誉教授だから当然でしょ

437:デフォルトの名無しさん
08/06/07 20:33:09
もうf77は教えないほうがいいだろ。常識的に考えて。

438:デフォルトの名無しさん
08/06/07 20:50:45
>>437
講義の教科書ではなくて、研究室用だと思うよ。
研究室にはf77の遺産があるから、参考図書として必須。

439:デフォルトの名無しさん
08/06/07 23:48:11
>>435
大学の生協なら売ってるのが普通だと思ってた。
俺は慶応じゃないけどね

440:デフォルトの名無しさん
08/06/08 00:29:15
阪大にも置いてあるよ。
どこでも置いてあると思う。

441:デフォルトの名無しさん
08/06/08 00:51:54
浦昭二のFortran66の方はさすがに絶版になってしまったな。



442:デフォルトの名無しさん
08/06/08 01:30:36
>>441
遺産プログラムがFortran66(ブロックIFなし)で書かれていて驚愕した orz

443:デフォルトの名無しさん
08/06/08 02:26:31
>>442
俺のもらった遺産プログラムは、拡張DO LOOPを使っていたせいでコンパイラを通らず
すごく困った。

拡張DO LOOPとは、DO LOOP の中からGOTOで飛び出して、なんかやった後DO LOOPの
中に帰ると、またさっきのLOOPの続きが出来るというもので、当時の俺の想定外の悶絶機能で
何がおきているのか全く理解できなかったw

FORTRAN77で廃止になったはずなんだが、たいがいのコンパイラは密かに対応していた。
まぁ、怖いので書き直したがw

そういう事もあるのでFORTRAN66の教科書も無くしてはいけないと思うですよ。

444:デフォルトの名無しさん
08/06/08 11:39:03
そうか、塾生は君付けなんだなw

445:デフォルトの名無しさん
08/06/18 02:23:23
今、本屋で買える範囲で、おすすめのFortranの教科書は何ですか?

絶版は困ります。アマゾンの古本で買えるならいいですが。


446:デフォルトの名無しさん
08/06/18 04:13:42
Michael Metcalf, John Ker Reid, "Fortran 90/95 Explained"

447:デフォルトの名無しさん
08/06/18 10:04:59
>>445
Fortranを勉強する目的などを明確にせよ。
 規格:わからない、f77、f90以上
 目的:講義対応、研究室対応、趣味、業務対応

以下はとりあえずの推薦
規格わからない 講義対応
 Fortran 90プログラミング 冨田博之
 URLリンク(www.amazon.co.jp)

f77規格 講義対応
 FORTRAN77入門 浦昭二
 URLリンク(www.amazon.co.jp)

f90規格 研究室対応
 数値計算のためのFortran90/95プログラミング入門 牛島省
 URLリンク(www.amazon.co.jp)

448:デフォルトの名無しさん
08/06/18 10:10:18
テンプレに追加きぼん

449:デフォルトの名無しさん
08/06/19 00:50:25
>>446-447
さんくす


450:デフォルトの名無しさん
08/06/29 02:12:38
>>434
バックエンド共通ってだけで速度が同じになるなんてありえない。
字句解析や最適化は言語仕様とコンパイラ実装に強く依存する。

Fortranのシンプルさはその点強みだったけど、
最近のFortranはその他言語にどんどん近づいているからねえ。

451:デフォルトの名無しさん
08/06/30 03:11:59
>>450
>>415のベンチマークの話だから。

452:デフォルトの名無しさん
08/07/03 00:26:47
>>450
ベンチマークサイトのような、汎用アルゴリズムならfortranは速くなくて不思議は無いですよ。
C/C++とfortranのユーザー数考えたら、intelだってどちらにリソース割くかは自明でしょう。
ナンバークランチング用途ならfortranは多分まだ速いと思いますけど。


453:デフォルトの名無しさん
08/09/06 00:53:21
Intel Fortran の次期バージョンベータテストが始まっている。
V.11が来年一月までの期限付きだが、無条件に配布されている。
Windows版はVisualStudio(Fortran 専用)つきで。

気になる人はチェキラ!w


454:デフォルトの名無しさん
08/09/21 14:52:15
もともとIntel FortranとC++は非商用目的なら無料で使えるんじゃない?

455:デフォルトの名無しさん
08/09/21 23:57:43
>>454
無期限はLinuxだけじゃないか?


456:デフォルトの名無しさん
08/09/22 15:10:37
Linux版だけだな
Linux版ならIntel FortranとIntel C++が無期限で無償で使える

457:デフォルトの名無しさん
08/09/22 15:12:37
非商用目的限定の話だけどね

458:デフォルトの名無しさん
08/10/07 00:50:36
1秒間プロセスを寝かせるのはsleepという関数があるけど,
10msとか100msとか,ミリセカンドオーダーで寝かせるようにするには
どうすればいいですか?
Intelコンパイラを使っているので,Intelコンパイラ依存の関数でも
いいです.
今は,無駄な計算ループを作って待たせているのですが,
それだとCPU負荷になるので,できればCPU負荷にならないで
10msとか100msとか待たせられると嬉しいです.

459:デフォルトの名無しさん
08/10/07 01:03:04
>>458
Windowsなら、USE IFWIN して Win32 の API の Sleep 関数でやれた気がする。
Linux は知らないが、同じようにやれると思う。


460:デフォルトの名無しさん
08/10/07 01:30:30
nanosleep(2)

461:デフォルトの名無しさん
08/10/07 03:06:27
>>458
ifort依存ならWin/Linux共に SLEEPQQ で
USE IFPORT
CALL SLEEPQQ(3000)
とミリ秒で指定できる。

但しあくまでも指定できるだけであって実際の休止時間をミリ秒単位の精度で
制御できるかどうかは別だが。

462:458
08/10/07 10:32:59
>>459-461
皆さん,ありがとうございます.
461さんのやり方でできました.
それほど正確性を求めていないので,このやり方で十分です.
ありがとうございました.

463:デフォルトの名無しさん
08/10/24 15:08:54
おまいら、CUDAでFORTRANがサポされたらどうする?逝く?

464:デフォルトの名無しさん
08/10/24 17:01:26
キュダってなに?

465:デフォルトの名無しさん
08/10/24 20:17:32
>>464
GPGPUの開発環境 >CUDA
GPUを計算に転用して10倍以上の高速化を実現するらす
URLリンク(journal.mycom.co.jp)

466:デフォルトの名無しさん
08/10/25 10:45:52
倍精度に対応したのか

467:デフォルトの名無しさん
08/11/11 02:36:04
CUDAよりインテルのLarrabeeが来そうだな。
ベクトル長を長くする方向だから80年代のプログラミングモデルが復活するし、
Intelのコンパイラは必ずサポートするだろうし。



468:467
08/11/11 14:06:01
ルンゲ・クッタ・ギル法を用いて連立の方程式を解きたいのですが,
時間刻みを小さくしても発散してしまうのですが,何かサブルーチン内に欠陥があるのでしょうか?
QA0=QA
CKA0=H*FNCT1(X0,YA0,YB0,W0,W1,A,A0,A1,B,B1,H0,L,PI)
RA1=(CKA0-2.0*QA0)/2.0
YA1=YA0+RA1
QA1=QA0+RA1*3.0-CKA0/2.0
XX=X0+H/2.0
CKA1=H*FNCT1(XX,YA1,YB1,W0,W1,A,A0,A1,B,B1,H0,L,PI)
RA2=(1.0-1.0/sqrt(2.0))*(CKA1-QA1)
YA2=YA1+RA2
QA2=QA1+RA2*3.0-(1.0-1.0/sqrt(2.0))*CKA1
CKA2=H*FNCT1(XX,YA2,YB2,W0,W1,A,A0,A1,B,B1,H0,L,PI)
RA3=(1.0+1.0/sqrt(2.0))*(CKA2-QA2)
YA3=YA2+RA3
QA3=QA2+3.0*RA3-(1.0+1.0/sqrt(2.0))*CKA2
XX=X0+H
CKA3=H*FNCT1(XX,YA3,YB3,W0,W1,A,A0,A1,B,B1,H0,L,PI)
RA4=(CKA3-2.0*QA3)/6.0
YA=YA3+RA4
X=X0+H
QA=QA3+3.0*RA4-CKA3/2.0
return
end
function FNCT1(X0,YA0,YB0,W0,W1,A,A0,A1,B,B1,H0,L,PI)
FNCT1=-2.0*W1*H0*YA0-(W1**2.0+(W0**2.0)*(A/A0+(B**2.0)/(2.0*L*A0)-(L*A1*B)/(2.0*A0)))*YB0+((W0**2.0)*(3.0*L*A1)*(YB0**2.0))/(PI*A0)
* -((W0**2.0)*(PI**2.0)*(YB0**3.0))/(4.0*L*A0)+(W0**2.0)*((-2.0*B1)/(PI*(W0**2.0))+(4.0*(L**2.0)*A1*A)/((PI**3.0)*A0)+(2.0*L*A1*(A**2.0))/((PI**3.0)*A0)
* -(2.0*(L**3.0)*(A1**2.0)*B)/((PI**3.0)*A0))
return
end
長々と申し訳ないです.お願いします

469:デフォルトの名無しさん
08/11/12 00:29:13
>>468
長すぎて考える気にならん。
関数のパラメータも多すぎて、条件も分からん。

ある特定のパラメータで発散するのか?ルンゲ喰ったギル法だけで異常が出るのか?
どういう状況なのかも分からん。

答えようも考えようも無い。
もう少し親切心を持って聞けw


470:468
08/11/12 03:21:51
>>469さん

すみません.どう書いたらよく伝わるのか分からず,書いてしまいました.
反省しております.

方程式は,T"(t)+T'(t)+ωT(t)+B1(t)T(t)+B2(t)T^2(t)+B3(t)T^3(t)=B4(t)
のような方程式で,
B1~B4の中にも変数が存在してtが変化した際には変数も変化させたいのです.

パラメータを変えた際には異常が現れなかったりします.

わかりにくくて申し訳ありません.よろしくお願いします

471:デフォルトの名無しさん
08/11/12 15:42:56
ブックオフで培風館のFORTRAN77入門が100円だった
さすがブックオフ本の価値が分かってない

472:デフォルトの名無しさん
08/11/13 00:52:36
>>471
古本屋での相場は美本で350円くらいだったと思う。
漏れは66版を150円で買ったような気がするw

>>470
>T"(t)+T'(t)+ωT(t)+B1(t)T(t)+B2(t)T^2(t)+B3(t)T^3(t)=B4(t)
よく分からんが、第三項はω^2じゃ無いか?
それでBx(t)が微小量だとすれば、なんとなく意図が分からなくも無い。

減衰振動の単振り子に非線形項と強制振動項が加わったものなるので、
なんとなく動きのイメージがつかめるのでは。
パラメータによっては、当然発散も大いにありうるだろう。

そもそも正しく計算できているのか自信が無いのなら、Bx(t)を全部0に置けば、
へたれて行く減衰振動になるので、そういう解析解のあるところからチェックしてゆけば
自分のプログラムの挙動に自信がもてるのでは?

詳しく知りたければランダウ=リフシッツの力学の非線形振動における共鳴の節を読めwwwww


473:デフォルトの名無しさん
08/11/13 01:00:09
>>471
教科書で使ってたりするから、古本の流通量は結構多い気がする。

474:デフォルトの名無しさん
08/11/13 02:46:31
カバーの焼けはあるけどおそらく新品
改訂版じゃなくて初版の方だったので返品在庫の放出かしら

配列のとこやってるけど、配列への入力とかあってすごいね
宣言なしで唐突に出てくるループカウンタにもびっくりするが

475:デフォルトの名無しさん
08/11/13 10:19:44
>>474
えぇ~?
暗黙に型宣言されてるぜ!

最近のスクリプト言語で、これに近い形に先祖がえりしているのがワロス


476:デフォルトの名無しさん
08/11/13 14:07:56
implicit none しても結局は暗黙のやつに則ることは多いよね。
character,logicalは無いので c と l で始めるようにしたり。
暗黙の型宣言はあれはあれで良いものだ。

477:デフォルトの名無しさん
08/11/14 21:11:06
いや、Cが複素数でZが倍精度複素数だろ

478:デフォルトの名無しさん
08/11/15 00:55:13
>>476
lは1と見間違いやすいからあまり使わないが、貴重な整数型なのでもったいない。

479:デフォルトの名無しさん
08/11/15 05:22:43
>>478
Lは大文字で使う

480:デフォルトの名無しさん
08/11/16 00:36:10
>>479
ボーランド記法?
昔PASCALに散々いじめられたから嫌w

481:デフォルトの名無しさん
08/11/16 12:49:36
>>477
c を複素数に振るときには character 型は s で始めるようにしてるな。
>>478
l は先頭でだけ使うようにしてる。 先頭が数字の変数ははじかれるから
混同しなくてすむ。

482:デフォルトの名無しさん
08/11/26 01:24:43
東京で先週開かれたISOのFortran部会でCoArrayFortranに関する扱いで
何か進展はあったのでしょうか?

つまりCAFを排除できたのでしょうか?

483:デフォルトの名無しさん
08/11/29 16:57:56
>>482
東京会議の結論がうpされたようだ。
URLリンク(www.nag.co.uk)




484:デフォルトの名無しさん
08/12/06 18:55:27
The shapes of the array expressions do not conform.
って出てくるのですが、どういうことですか?

485:デフォルトの名無しさん
08/12/11 16:23:19
配列間違ってないか?

486:デフォルトの名無しさん
08/12/13 02:11:30
C Pad for Salford FTN77
でCPU TIME を出力できるようにする方法はありますか?

487:デフォルトの名無しさん
08/12/13 18:02:53
>>486
HIGH_RES_CLOCK@

使い方は Library Reference に載ってる。

488:デフォルトの名無しさん
08/12/23 09:11:38
>>487
亀ですが
ありがとうございます。
出来ました!

489:デフォルトの名無しさん
09/01/24 19:01:34
関連スレが新しくなりました。

くだすれFORTRAN(超初心者用)その4
スレリンク(tech板)l50


490:デフォルトの名無しさん
09/02/12 00:23:26
うちの研究室の人たちがFortranでプログラミングしているけれど
お前らgoto使いすぎ・・・
というか自分の提出したプログラムがgoto文に書き換えられている・・・

subroutine、functionを使いすぎるのはスパコン計算に向かないと言っていたけれど
勉強してきた構造化プログラミングの手法自体がfortran向きでないということなの?

491:デフォルトの名無しさん
09/02/12 02:32:22
>>490
インラインされるようにしろw

コンパイラが安心してグローバル最適化できるように、付加的な情報を与えてやれ。
グローバルな構造が人間が見て分かりやすいならコンパイラも見て分かりやすい。

INTENT指定はつける。
FUNCTIONは、PURE (可能ならELEMENTAL)にする
とかすればDO..LOOP内で関数呼び出しとかしても並列化やベクトル化してくれるはず。
コンパイラによるが。

実際にベンチマークを取って、自分が勝利すれば、後は文句いわれないべw


492:デフォルトの名無しさん
09/02/12 12:52:18
>>491
intent属性つけると微妙に速くなるよね。
特にloop内で頻繁に呼ばれるやつだと全体で5から20%くらい違ってくる。
最近のコンパイラはF77で書いててもintentをわかってくれるから楽だ。

で・・・・pure,elementalの事はよくわからないので詳しく!

493:デフォルトの名無しさん
09/02/13 00:04:13
>>492
PUREもELEMENTALもFortran95で導入されたキーワードで、関数やサブルーチンの
副プログラムの性質をコンパイラに教える修飾語。

PURE FUNCTION unko(x) などのように宣言の頭につけて使う。

PUREというのは、大雑把には副作用がないこと。それゆえ、並列に動作させても
値が変わらないので、コンパイラは並列化など最適化をやれる。
(PUREを満たさねばならない条件はいくつかあるので、kwskはマニュアル見てくれ)

ELEMENTALはPUREの中でさらにいい性質を持ったもの。
本来は1変数で定義してあっても、配列を引数にとって、配列を返すように使ってもいいもの。

たとえば、
REAL :: a(10), b(10)
b = SIN(a)
こういう使い方が組み込み関数では出来るわけだが、自作の関数でもこうやって使っても
いい場合にELEMENTAL (要素毎の意味) をつける。 

これもコンパイラが利口なら、最適化をかけてくれるはず。

普通、DO..LOOP内に関数やサブルーチンの呼び出しがあると、呼び出し先で何が
あるかわからないので、最適化が抑止されるが、PUREやELEMENTALがあれば、
安心してコンパイラが最適化できる。(原理的にはw)



494:デフォルトの名無しさん
09/02/13 08:29:24
>>493
ありがとう。pure 属性は使い易そうだ。

495:494
09/02/13 11:01:13
loop内で頻繁に呼ばれる4個つくらいの function や subroutine に pure 属性つけただけで
並列(OpenMP)の時にちょびっと速くなった。わーい。
もともとdo-loopで一つずつ処理しているので elemental を使える状況は少ないけど、
もう少し理解してから使ってみよう。
inlineを少し深くするようになったからかコンパイル時間がすこし長くなったけどまあ気にしない。

496:デフォルトの名無しさん
09/02/17 22:12:50
>>491
なるほど、ためになったわ。サンキュー。PUREを付けるだけで簡単に速くなったわ。

担当のは基本的に古い人間だから、あまりコンパイラ任せにするなとか言うんよな…。
Fortran90っぽい77形式の文も、「完成した後で直します」とか言ってそのままだったり。
implicit realもcommon文もどうにかしたいけれど、きっと文句言うんだろうな。

497:デフォルトの名無しさん
09/02/18 02:22:54
>>496
PURE 属性はもともと HPF で導入されたものなので、それなりのメーカならコンパイラ
にも最適化のノウハウの蓄積があるはず?
ELEMENTAL も一変数用にサブルーチン書いて、配列を渡せばいいので、DO..LOOPを
あからさまに書かずにすむ。時々ソースがすごくすっきりして楽になることがある。

でも、ついでだから古い人間を擁護してみるかw

FORTRAN77 にも、むきだしの良さがあるw
Fortran90 以降だとソース上は簡潔になっても、影で動的にテンポラリ配列とか取るので
気持ち悪いところもある。長時間並んだJOB がStackOverflow とかで止まっていたり
すると殺してやりたくなる。古い時代の人はその辺が嫌なのかもしれないw

COMMON もグローバル変数だと思うとむかつくだろうが、 自分でガベージコレクション
しながら、変数領域を自由自在に再利用して操っていると思うと、結構気持ち良い。
メモリーが貴重だった時代には、利用価値が高かった。

IMPLICIT に関していえば IMPLICIT NONE をつけた上で、頭文字を守ればいい。
C と違ってローカルスコープの捨て変数が使えないから、頭文字を守らないと
いちいち宣言部に戻って型を確かめないといけないのでめんどい。
Fortran では数値の精度が致命的に重要で、倍精度計算中に1個単精度が混じるだけでも
すべての結果が単精度のゴミになるのだから、数値精度に無頓着の他のおしゃれ言語の
連中の言説を信じると痛い目に会わされるw



498:デフォルトの名無しさん
09/02/18 13:19:26
pure属性付けたサブルーチン内で
pure属性付けた他のサブルーチン呼ぶ時には
interfaceで引数部分の性質を教えて置かないといけないんだね。
f77書式で書いてるのに interface intent が
現れるのでムズムズするけどムズムズするだけなのでそのまま。

499:デフォルトの名無しさん
09/02/26 16:08:37
OpenMPを試してみてるんだけど、いまいちうまくいかないんだが・・・
ifort -QParallel でコンパイル+実行するとうまくいくんだけど、
同じソースをifort -Qopenmp でコンパイル+実行すると
Program Exception - stack overflow
っていうエラーが出る
ソース中にOpenMP構文は use omp_lib だけ
これってどうしたらいいのかな?

500:499
09/02/26 16:22:10
おk事故解決
-auto_scalar オプションでいけた


次ページ
最新レス表示
レスジャンプ
類似スレ一覧
スレッドの検索
話題のニュース
おまかせリスト
オプション
しおりを挟む
スレッドに書込
スレッドの一覧
暇つぶし2ch