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コアでは確かにスピードダウンが
確認できましたが、４コアではスピードアップが確認できました、と回答。
まるでリフォーム詐欺のようだ・・・
単純な行列×ベクトルのみの計算でも最大３０％アップ程度なのだが
並列計算ってそんなものですか？
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を入れれば単精度は４倍速、倍精度は２倍速とか、素直に早くなる。
コンパイル時のメッセージにベクトル化云々と出れば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
バージョンアップ込みで年間数万円（いくらか忘れた）のサポート料を払っているので、
使ってみたのですが、本当にがっかりです。
サポート料金を１０倍にしてもいいので、まともな回答を期待したいです。
ググってますがなかなかいい情報が出てこないです。
私の最高は３割増が最高で、実務的な計算ではほとんど役に立ちません。
あれこれ考えている間に、時間をかけて計算したほうがいいのかと思えてきました。
366:デフォルトの名無しさん
07/10/24 12:59:13
>>363
いや詐欺ではない。
既存プログラムを自動並列化して場合、スピードダウンするのはごく普通の現象。
はじめから並列化を前提に設計して、かつ指示行を入れまくってもいまいちなのは、実は常識だｗ
367:デフォルトの名無しさん
07/10/24 15:53:27
>>365
スレ違いな話だけど、本気のガチンコサポートをするなら商売する側としては10倍でもペイしないと思うけどね・・・
そこはしょうがないっすよ
368:358
07/10/24 16:19:19
日本語がおかしかった。失礼しました。
>>367
あるソフトベンダーは年間１４０万円のサポートでいたれりつくせり。
（何が原因でしょうね？と問い合わせたら修正されてかえってくる）
そこは最優先でサポートを継続している。
別の４０万円のところは、それなりだけど解決のヒントはくれる。
（これらは言語ではなくてアプリケーションですが）
VisualFORTRANは「開発に問い合わせます」のまま放置・・・。
ただ、今回の件については決して特殊な質問ではなく
広告にうたってある性能が出ないよに始まり、
そちらのもっている理想的なプログラムでいいから出してくれといっても
出てこない。いったいどんな質問だったら答えることができるのだろう。
369:デフォルトの名無しさん
07/10/24 19:20:01
>>368
言語のサポートの場合、コンパイラのバグに対する対応の早さだと思う。
そういう点でＩｎｔｅｌはパッチの出方は早い。
あとサポート掲示板には、様々なレベルの質問が出ているのでこっちも見てみたらいいかも。
URLﾘﾝｸ(softwarecommunity.intel.com)
370:デフォルトの名無しさん
07/10/25 20:49:16
ＮＥＣ：世界最高性能のスパコン「ＳＸ－９」を発売
ＮＥＣは２５日、世界最高レベルの計算能力を持つスーパーコンピューター「ＳＸ－９」を開発し、
発売を始めたと発表した。１秒１６０億回の計算能力を持ち、最大限の５１２台を連結して運用し
た場合、１秒８３９兆回の計算が可能なシステムが完成する。
新開発の高性能ＣＰＵ（中央演算処理装置）を搭載した。国内外の大学や研究機関などが、
気象予測やロケットの設計、核融合のシミュレーションなどに利用することを想定している。
価格は１億５０００万円。
現在の世界最速のコンピューターシステムは米ＩＢＭが開発し、秒間２８０兆回の計算ができる。
ＮＥＣは、１０年の完成を目指す世界最速スパコンの開発計画を、日立製作所や富士通と進めて
いる。
URLﾘﾝｸ(mainichi.jp)
371:デフォルトの名無しさん
07/11/02 00:54:14
プレステ３用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ﾘﾝｸ(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を利用するプログラムを１行も変更することなく、性能を
更に高めることが可能です。
379:デフォルトの名無しさん
07/12/06 03:09:49
初心者で最小二乗法のプログラムが組めないんですがどうすればいいんでしょうか？
380:デフォルトの名無しさん
07/12/06 10:03:56
>>379
図書館へ行って、データ処理の本を読むんだ！
381:デフォルトの名無しさん
07/12/06 19:03:50
データ読み込んでΣ求めて公式通り計算すればおｋ
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
最近の Ｉntel 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
㌧クス
確認・把握した。
漏れは近年デバッグ実行をまったく使っていなかったのと、エラーの出る奴は
コマンドラインから実行したことがなかったので気づかなかったｗ
多分、デバッグなし実行でもトレースバックは呼び出し履歴に強制リダイレクトされているが、
（error出力かなんかで）
この場合呼び出し履歴を見るメニュー項目が無いというバグなのだろうｗ 空想だがｗ
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を試しています。ただ、解いている問題の係数行列が良くないのか
PETSｃだと収束しないことがあって、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とかフリーコンパイラも転がっていることだし、
７７でよければWatcomとか昔なつかしいのもあるし、g77のWin用もあるし、
Fortranで素直にやらせてやれｗ
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
正直質問の意味がよく分からん。
コンパイラのバグでなければ、ありえない症状なので
ここに書いてない部分がおかしいとしか言いようが無い。
ｆ77の範囲ならg77とか他のコンパイラもあるのでそれで試してみるとか、
salfordのｆ95で試してみるとかして見るのも一つの案。
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
エラーでないぞｗ
エラーメッセージをコピーしてみせてくれ。
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
お前ひどい奴だなｗｗｗ
>>420
最初ッからそう言ってくれれば、すぐに答えられる。
それはCVF時代のFAQで、日本語のパスやファイル名が入るとコンパイラが対応していない
のでエラーを吐く。２バイト文字を使わないディレクトリでやればおｋ。
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は使っちゃいかんｗ
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入門（FORTRAN７７入門)が普通において有って驚いた
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の続きが出来るというもので、当時の俺の想定外の悶絶機能で
何がおきているのか全く理解できなかったｗ
FORTRAN77で廃止になったはずなんだが、たいがいのコンパイラは密かに対応していた。
まぁ、怖いので書き直したがｗ
そういう事もあるので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 専用）つきで。
気になる人はチェキラ！ｗ
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が来そうだな。
ベクトル長を長くする方向だから８０年代のプログラミングモデルが復活するし、
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
長すぎて考える気にならん。
関数のパラメータも多すぎて、条件も分からん。
ある特定のパラメータで発散するのか？ルンゲ喰ったギル法だけで異常が出るのか？
どういう状況なのかも分からん。
答えようも考えようも無い。
もう少し親切心を持って聞けｗ
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
古本屋での相場は美本で３５０円くらいだったと思う。
漏れは６６版を１５０円で買ったような気がするｗ
>>470
>T"(t)+T'(t)+ωT(t)+B1(t)T(t)+B2(t)T^2(t)+B3(t)T^3(t)=B4(t)
よく分からんが、第三項はω^2じゃ無いか？
それでBｘ(ｔ）が微小量だとすれば、なんとなく意図が分からなくも無い。
減衰振動の単振り子に非線形項と強制振動項が加わったものなるので、
なんとなく動きのイメージがつかめるのでは。
パラメータによっては、当然発散も大いにありうるだろう。
そもそも正しく計算できているのか自信が無いのなら、Bx(ｔ）を全部０に置けば、
へたれて行く減衰振動になるので、そういう解析解のあるところからチェックしてゆけば
自分のプログラムの挙動に自信がもてるのでは？
詳しく知りたければランダウ=リフシッツの力学の非線形振動における共鳴の節を読めｗｗｗｗｗ
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は１と見間違いやすいからあまり使わないが、貴重な整数型なのでもったいない。
479:デフォルトの名無しさん
08/11/15 05:22:43
>>478
Lは大文字で使う
480:デフォルトの名無しさん
08/11/16 00:36:10
>>479
ボーランド記法？
昔PASCALに散々いじめられたから嫌ｗ
481:デフォルトの名無しさん
08/11/16 12:49:36
>>477
c を複素数に振るときには character 型は s で始めるようにしてるな。
>>478
l は先頭でだけ使うようにしてる。 先頭が数字の変数ははじかれるから
混同しなくてすむ。
482:デフォルトの名無しさん
08/11/26 01:24:43
東京で先週開かれたＩＳＯのＦｏｒｔｒａｎ部会でＣｏＡｒｒａｙＦｏｒｔｒａｎに関する扱いで
何か進展はあったのでしょうか？
つまりＣＡＦを排除できたのでしょうか？
483:デフォルトの名無しさん
08/11/29 16:57:56
>>482
東京会議の結論がうｐされたようだ。
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
インラインされるようにしろｗ
コンパイラが安心してグローバル最適化できるように、付加的な情報を与えてやれ。
グローバルな構造が人間が見て分かりやすいならコンパイラも見て分かりやすい。
ＩＮＴＥＮＴ指定はつける。
ＦＵＮＣＴＩＯＮは、ＰＵＲＥ （可能ならＥＬＥＭＥＮＴＡＬ）にする
とかすればDO..LOOP内で関数呼び出しとかしても並列化やベクトル化してくれるはず。
コンパイラによるが。
実際にベンチマークを取って、自分が勝利すれば、後は文句いわれないべｗ
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を満たさねばならない条件はいくつかあるので、ｋｗｓｋはマニュアル見てくれ）
ELEMENTALはPUREの中でさらにいい性質を持ったもの。
本来は１変数で定義してあっても、配列を引数にとって、配列を返すように使ってもいいもの。
たとえば、
REAL :: a(10), b(10)
b = SIN(a)
こういう使い方が組み込み関数では出来るわけだが、自作の関数でもこうやって使っても
いい場合にELEMENTAL （要素毎の意味） をつける。
これもコンパイラが利口なら、最適化をかけてくれるはず。
普通、DO..LOOP内に関数やサブルーチンの呼び出しがあると、呼び出し先で何が
あるかわからないので、最適化が抑止されるが、PUREやELEMENTALがあれば、
安心してコンパイラが最適化できる。（原理的にはｗ）
494:デフォルトの名無しさん
09/02/13 08:29:24
>>493
ありがとう。pure 属性は使い易そうだ。
495:494
09/02/13 11:01:13
loop内で頻繁に呼ばれる４個つくらいの 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を
あからさまに書かずにすむ。時々ソースがすごくすっきりして楽になることがある。
でも、ついでだから古い人間を擁護してみるかｗ
FORTRAN77 にも、むきだしの良さがあるｗ
Fortran90 以降だとソース上は簡潔になっても、影で動的にテンポラリ配列とか取るので
気持ち悪いところもある。長時間並んだJOB がStackOverflow とかで止まっていたり
すると殺してやりたくなる。古い時代の人はその辺が嫌なのかもしれないｗ
COMMON もグローバル変数だと思うとむかつくだろうが、 自分でガベージコレクション
しながら、変数領域を自由自在に再利用して操っていると思うと、結構気持ち良い。
メモリーが貴重だった時代には、利用価値が高かった。
IMPLICIT に関していえば IMPLICIT NONE をつけた上で、頭文字を守ればいい。
C と違ってローカルスコープの捨て変数が使えないから、頭文字を守らないと
いちいち宣言部に戻って型を確かめないといけないのでめんどい。
Fortran では数値の精度が致命的に重要で、倍精度計算中に１個単精度が混じるだけでも
すべての結果が単精度のゴミになるのだから、数値精度に無頓着の他のおしゃれ言語の
連中の言説を信じると痛い目に会わされるｗ
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
おｋ事故解決
-auto_scalar オプションでいけた
501:デフォルトの名無しさん
09/02/26 16:49:38
>>500
だたの
-auto
や
-static
でうまく行くこともあるよ。あとは実行時の環境変数の
setenv KMP_STACKSIZE 16M
とかデフォより大きめに取ると幸せなことが多い気がする。
502:デフォルトの名無しさん
09/02/26 16:50:50
>>500
ていうか、
-auto_scalar
ってデフォルトになってない？
503:デフォルトの名無しさん
09/02/26 17:23:54
-openmpだとデフォで動的（-auto）になるんだっけか
>>499は単にスタック不足なだけでは？
ところで、KMP_STACKSIZEってunlimitedに設定できないの？
504:デフォルトの名無しさん
09/02/26 17:58:56
動的っていってもヒープとスタックがあるよね。
Fotranの場合ってどういう風に割り当てられるんだろう・・・そういやあまり意識したことなかった。
505:デフォルトの名無しさん
09/02/27 02:43:33
>>504
元々はAUTOMATIC変数の類はスタックに、ALLOCATE命令のときはヒープだった気がする。
ただ歴史的経緯のせいでサブルーチン側で巨大な配列をとっているプログラムが多くて
それがSTACKOVERFLOWを引き起こすから、IntelFortranなんかはデフォールトでは
スカラー変数はスタックで、配列の類はヒープで確保するようになった気がする。
インテルの掲示板に、デフォールト・オプション切り替えの話が出ていた旗もするのだが
昔のことなので記憶がはっきりしない。
OpenMPを指定すると、オプションが切り替わってしまう話はこのスレに出ている。
URLﾘﾝｸ(software.intel.com)
FortranはFORTRAN66時代がSTATICに変数をとっていた関係で、
FORTRAN77で規格の上ではサブルーチン類の変数はAUTOMATICでありうるように
なったにもかかわらず、互換性の観点からデフォールトでSTATICに変数を
とっていることが多くて、それはFortran90でもしばらく続いていたという歴史的経緯がある。
506:デフォルトの名無しさん
09/02/27 11:54:19
>>505
そのスレのSteveさんが言ってるlocal scalar variablesって、配列じゃない（次元１の）フツーの変数のことだよね？
だとしたら、
If you say /Qopenmp (or /Qopenmp-stubs), this changes to /automatic - ALL local scalar variables are allocated on the stack.
のALL local scalar variablesはALL local variables（つまりローカルスカラ変数だけじゃなく、ローカル配列変数も含む）の間違いじゃないかな？
勘違いしてたらすいません
507:デフォルトの名無しさん
09/02/27 13:12:01
>>506
そうなんじゃないかな？
配列までスタックにとるからoverflowしていると考えるとつじつまが合う？
508:506
09/02/27 13:52:41
>>507
そうだよね。SteveさんIntelのスタッフなんだし書き間違わないでくれよぉ・・・
OpenMPの場合にデフォルトで全てのローカル変数（スカラも配列も）がスタックになるのは
スレッドセーフのためってあるけど、逆に言えばそれらをstaticに扱ったらスレッドセーフに
ならないってことだよね。それってやっちゃまずいことなんだろうか？
OpenMPまともに使った経験がないもんで・・・識者の方コメントおながいします
509:デフォルトの名無しさん
09/02/27 14:48:49
-openmp
つけてコンパイルすると、仮にautomaticは嫌いだっと
-noauto
つけても消されるはず。なのでopenMPを読んだ状態だと
有無をいわさずにぜんぶautomaticになる・・・と理解しているけど。
shared属性をつけてかつcommonブロックに貼り付けるとかして
やや無理なことをしないとstatic領域には配置されないと思ふ。
-staticはその名のイメージとは違って共有ライブラリを使わない
(というか実行ファイルに埋め込んじゃう)オプションだし。
spec.orgにあるOMP2001の情報のところに試した人たちが
実際に使ったオプションが載ってるので、自分のところのシステムに
近い(or 同じ)マシンがあったらそれと同じモノをだいたい使っているけど、
まあ大間違いではないと思う。お薦めできるのかは分らないけど・・・。
510:デフォルトの名無しさん
09/02/27 15:30:30
え、OpenMPだとサブルーチン内ではstatic配置ってできないわけ？
数値計算では馬鹿でかい配列は可能な限りstaticにするのが基本じゃない？
サブルーチン渡る度に配列をコピーとかやってられんでしょ？
その代わりメモリのアロケーションはコーディングする側がきちんと考えてあげないといけないけど。
この辺がリソースをある意味富豪的に使うことが前提とされるJavaやC#などのオブジェクト系言語との違いだね。
511:デフォルトの名無しさん
09/02/27 16:56:15
510はちょっと尋ね方が変だな。
副プログラムが並列のloop内にあるのか、
副プログラム内で並列のloopがあるのかで答は変わる。
512:510
09/02/27 19:38:52
的外れな質問してたのかな、ごめん・・・。
>>511の両方のケースとも知りたいので、よかったら教えてほしい。
流体とかの数値計算では行列の反復法とかはサブルーチン内でやるけど、
そういうケースでメモリのアロケーションはどうすべきなのか、という視点で興味があります。
513:デフォルトの名無しさん
09/02/28 00:58:18
>>512
引数に載ってれば、Fortranは参照呼出しなんだから配列のコピーが生じるはずもない。
CやPASCALとは違う。（まぁ最近はVALUE属性も出来るようになったが）
ここでは副プログラム内で新たに割り付けられる配列の話だろう。
それが動的に取られるとして、スタックに取られるか、ヒープに取られるかの違い。
514:510
09/02/28 01:28:49
>>513
＞副プログラム内で新たに割り付けられる配列
あ、そういう話か・・・ありがとう。
自分は一時作業用の配列も全部メイン側で用意してやって引数渡ししてるんだけど、
そういうのはサブルーチン内でスタックやヒープで取ってやる方がいいのかな？
皆さんどうしてますか？
515:デフォルトの名無しさん
09/03/03 06:16:03
allocatableを副プログラム内で作ってそれ(先頭アドレスやらサイズその他)を
メイン階層や他の副プログラムに行き渡らせるのが割りと
面倒臭いのでメイン階層でやってる。
メイン階層が長くなるのは正直あまり好みではない、けど
まあ好みの問題なだけなので気にしないようにしている。
allocatable使うと場合によってはプログラム内で明示するサイズより大きいメモリを扱えたり
(コンパイラはハンドルできるメモリサイズ ＜ OSがハンドル出来るメモリサイズ、なので)、
格子数(データ数)の違う処理を同じ実行イメージで出来る(コンパイルしなおさなくて良い)のが
いいけど、最近ではサイズをプログラム内で固定変数で明示していちいちコンパイルしている。
その方がミスは少ない気がする。気のせいかもしれないけど。
516:デフォルトの名無しさん
09/03/03 12:11:34
何がスタックやヒープに割り当てられるのか分からんかったのでググった。
ブログだけどｗ
URLﾘﾝｸ(monologuemidnight.blogspot.com)
module中の大域的な変数、save属性付きの変数→.bss
allocatable属性な配列→ヒープ
普通にサイズが宣言されている自動配列→スタック
手続き中の変数（mainでも同様）→スタック
ってことらしいが…。
この辺は規格ではないっぽいので、コンパイラの仕様をちゃんと確認したほうがいいな。
517:デフォルトの名無しさん
09/03/11 18:46:48
ifortで複数の数字を出力すると勝手に改行されるんだけど、改行を抑制するオプションってある？
gforntranでは改行されないんだけど。
real(8) a(5)
write(*,*) a(:)
とするとa(3)とa(4)の間で改行されちゃう。
518:デフォルトの名無しさん
09/03/12 04:57:59
>>517
気持ち悪いけど、下記のようにしてもだめかな？
クラシックな書き方だけど。
write(*, '(100f10.5)') a(:)
519:デフォルトの名無しさん
09/03/12 09:37:42
>518
ありがとん。
write(*,'(100f)') a(:)
でいけた。
520:デフォルトの名無しさん
09/03/15 01:01:35
改行抑制なら、write文ならadvance='no'を付ければできるはず
出力するものが決まってるなら、改行してしまう書式の箇所に\（バックスラッシュ）を入れてもいい
521:デフォルトの名無しさん
09/03/15 11:34:42
>520
advanceはフォーマット文がいるから、結局519と同じでは？
522:デフォルトの名無しさん
09/03/17 03:12:33
>>521
'(100f)' なんて書くよりはスマートな気がするけど？
ま、どうでもいいことだな
523:デフォルトの名無しさん
09/03/18 04:24:12
do i = 1, 5
write(*, '(f10.5,$)') a(i)
enddo
write(*,*) ''
はダメかな？
改行するなよ、の＄は非標準だけどDEC以来ほぼ標準と思っていいのではないだろうか
524:デフォルトの名無しさん
09/03/18 21:51:40
ちょっとテストしてみた＠Intel Fortran 10.1
integer::i
integer,parameter::imax=5
real(8)::f(imax)
do i=1,imax ; f(i)=i ; enddo
print *, "print *,f"
print *, f
print *, '---'
print *, "write(*,'(f)') f"
write(*, '(f)') f
print *, '---'
print *, "write(*,'(f)') (f(i),i=1,imax)"
write(*, '(f)') (f(i),i=1,imax)
print *, '---'
print *, "write(*,'(f\)') f"
write(*, '(f\)') f
print *, '---'
print *, "write(*,'(f)',advance='no') f"
write(*, '(f)',advance='no') f
print *, '---'
print *, "write(*,'(f)',advance='no') (f(i),i=1,imax)"
write(*, '(f)',advance='no') (f(i),i=1,imax)
print *, '---'
end
525:>>524の結果
09/03/18 21:52:34
print *,f
1.00000000000000 2.00000000000000 3.00000000000000
4.00000000000000 5.00000000000000
---
write(*,'(f)') f
1.0000000000000000
2.0000000000000000
3.0000000000000000
4.0000000000000000
5.0000000000000000
---
write(*,'(f)') (f(i),i=1,imax)
1.0000000000000000
2.0000000000000000
3.0000000000000000
4.0000000000000000
5.0000000000000000
---
write(*,'(f\)') f
1.0000000000000000 2.0000000000000000 3.0000000000000000 4.0000000000000000 5.0000000000000000 ---
write(*,'(f)',advance='no') f
1.0000000000000000
2.0000000000000000
3.0000000000000000
4.0000000000000000
5.0000000000000000 ---
write(*,'(f)',advance='no') (f(i),i=1,imax)
1.0000000000000000
2.0000000000000000
3.0000000000000000
4.0000000000000000
5.0000000000000000 ---
526:デフォルトの名無しさん
09/03/18 22:07:03
ごめん。
print文かwrite文か、また出力並びをfとするか(f(i),i=1,imax)とするかで、
出力結果が違うのかと思ってやってみたが、違いはないっぽい。
というわけで、>>524をもっとすっきりさせた。
integer::i
real(8)::f(3)
do i=1,3 ; f(i)=i ; enddo
print *, "write(*,*) f"
write(*,*) f
print *, '---'
print *, "write(*,'(f)') f"
write(*,'(f)') f
print *, '---'
print *, "write(*,'(f\)') f"
write(*,'(f\)') f
print *, '---'
print *, "write(*,'(f)',advance='no') f"
write(*,'(f)',advance='no') f
print *, '---'
end
527:デフォルトの名無しさん
09/03/18 22:19:08
↑の結果
write(*,*) f
1.00000000000000 2.00000000000000 3.00000000000000
4.00000000000000 5.00000000000000
---
write(*,'(f)') f
1.0000000000000000
2.0000000000000000
3.0000000000000000
4.0000000000000000
5.0000000000000000
---
write(*,'(f\)') f
1.0000000000000000 2.0000000000000000 3.0000000000000000 4.0000000000000000 5.0000000000000000 ---
write(*,'(f)',advance='no') f
1.0000000000000000
2.0000000000000000
3.0000000000000000
4.0000000000000000
5.0000000000000000 ---
結論
①書式省略すると、適当な位置で改行してくれる（規格で決まってるんかな？）
②書式指定すると、書式の中身が終わる度に改行
③書式に\を付けると改行抑制
④書式にadvance='no'（デフォルトはyes）を付けるとwrite文が終了したときに行う改行を抑制（②との違いに注意）
って感じかな。
改行抑制の\と改行指定の/が混ざるとどんな出力になるんだろう。
528:デフォルトの名無しさん
09/03/18 22:20:24
ごめん、>>526は
real(8)::f(3)
do i=1,3 ; f(i)=i ; enddo
は↓に修正。
real(8)::f(5)
do i=1,5 ; f(i)=i ; enddo
529:デフォルトの名無しさん
09/03/22 01:30:43
＞①書式省略すると、適当な位置で改行してくれる（規格で決まってるんかな？）
コマンドプロンプトの右端で改行だと思う
違うかなー
530:デフォルトの名無しさん
09/03/22 10:23:20
>>529
他のコンパイラでは改行しないので、コンパイラのせいだと思います。
>>524-528
ありがとうございます。
参考になりました。
write(*,'(f\)') f
が良さそうですね。
531:デフォルトの名無しさん
09/03/23 15:14:13
>>529
Fortranに限らず、プロンプトのサイズは関係ないよ
532:デフォルトの名無しさん
09/03/27 05:59:41
亀だけど、Ifortなら<>がお勧め。
多次元配列の最初の数を入れることが多いです
例
program main
implicit none
integer,parameter :: num = 9
integer :: ii,jj
real :: arry(num,num)
do ii=1,num
do jj = 1,num
arry(ii,jj) = ii*jj
enddo
enddo
write(6,'(<num>F)') arry
end program
533:デフォルトの名無しさん
09/03/27 09:10:33
>532
勉強になった。
ありがとん
534:デフォルトの名無しさん
09/03/28 00:52:37
>>531
WindowsはDOS窓のサイズで改行位置が変わるぞ。
535:デフォルトの名無しさん
09/04/03 21:43:51
>>532
<>で整数くくるとその数が使えるんだね・・・知らなかった。昔からあった機能？
>>534
そうなんだ。>>531はLinuxの場合ルだけかな
536:デフォルトの名無しさん
09/04/03 22:09:05
>>535追加
<実数変数>でも勝手に整数にキャストしてくれるね。
あと他のところにも、例えば
num1=5
num2=1
write(6,'(<num>F<num1>.<num2>)') arry
みたいにして使えるね。
537:デフォルトの名無しさん
09/04/04 02:01:36
>>536
Fortran2008 では、 *( i3, '-' ) 見たいな感じで ＊ が無限反復の指定子になるらしい。
538:デフォルトの名無しさん
09/04/08 09:50:52
*( i3, '-' )
を、スキー帽がぶったバカボンのパパっぽい大人と
すました子供に見えるようになると、やばいね。
539:デフォルトの名無しさん
09/04/08 12:40:58
やばいな
540:デフォルトの名無しさん
09/04/17 23:47:22
Fotranはじめようと思ってインストールして、
簡単なプログラムを実行しようとしたんだけど、
ld: crt2.o: No such file: No such file or directory
ってエラーが・・・。ググってもよくわからず・・・。
助けてください。
541:540
09/04/18 00:12:28
アンインストールして、はじめからやり直したらでけた。
542:デフォルトの名無しさん
09/04/18 20:15:44
>>542
おまえ、プログラムのセンスあるから、今後も缶ばれ！
543:デフォルトの名無しさん
09/04/19 12:16:56
>>542
自画自賛?
544:デフォルトの名無しさん
09/04/23 16:31:22
windowsXPでフォートラン90をフリーで使いたいのですけど、
cygwin以外でどうにか出来ますか？
545:デフォルトの名無しさん
09/04/23 22:10:01
free Fortran でぐぐれ
546:デフォルトの名無しさん
09/04/24 00:32:48
>>544
g95だな・・・
547:デフォルトの名無しさん
09/04/24 23:53:11
windowsでVMWareを立ち上げてUbuntuいれてgfortran
548:デフォルトの名無しさん
09/05/13 00:24:34
n次元配列をなめるルーチンを作成できる?
reshepe使って一次元に落とし込んで処理しようとしたら、segmentation faultになってしまった。
コードはこんな感じ。
-----------------------------
program main
integer :: arr(3,4,5)
integer :: i,j,k
do i=1,3
do j=1,4
do k=1,5
arr(i,j,k) = i+j+k
enddo
enddo
enddo
call check(arr)
end program main
subroutine check(arr)
integer,intent(in) :: arr(:,:,:)
integer,allocatable :: tmp(:)
integer :: i,elem=1
do i=1,size(ubound(arr))
elem = elem * (ubound(arr,i)-lbound(arr,i)+1)
enddo
allocate(tmp(elem))
tmp = reshape(arr,shape(tmp))
do i=lbound(tmp,1),ubound(tmp,1)
print *,tmp(i)
enddo
end subroutine check
-----------------------------
549:デフォルトの名無しさん
09/05/13 01:22:27
>>548
PRINT *, arr
でだめなのか？
550:デフォルトの名無しさん
09/05/13 20:37:56
>>548
ほらよ。つーかこれテンプレにしてちょ
・形状引継ぎ配列を引数に使う
・割り付け配列を引数に使う
・ポインタを引数に使う
・引数キーワードやoptional属性を使う
・関数の返値を配列にする
↑をしたいなら、副プログラムは明示的引用仕様にしなきゃならない。
■明示的引用仕様
・副プログラムを内部副プログラムにする（呼び出し側でcontains）
・副プログラムをモジュール副プログラムにする（モジュールでcontains）
・副プログラムが外部副プログラムならば、それを呼び出す側でinterfaceを用意する
■暗黙的引用仕様（FORTRAN77的なやり方）
・副プログラムが外部副プログラムで、それを呼び出す側でinterfaceを用意しない
なお、副プログラムとは関数やサブルーチンを指す。
551:デフォルトの名無しさん
09/05/13 23:25:40
>>548
多次元配列を一次元にするには PACK UNPACK を使うのが一般的だと思う。
Fortran90 で導入されている。
例題はマニュアルをみてちょｗ
552:デフォルトの名無しさん
09/05/14 23:40:08
FORTRAN90で引数の整合性をチェックする方法ってないかな？
よくそこでバグがおきるんだ。
553:548
09/05/14 23:50:44
>>549
print の部分は、本番仕様ではガチャガチャ変更したいところなのでございます。
554:デフォルトの名無しさん
09/05/14 23:51:21
>>552
>>550
555:548
09/05/14 23:51:43
>>550>>551
って～ことは、呼び出し側でinterfaceを用意して、packで一次元化するということか。
------------------------------------------
program main
integer :: arr(3,4,5)
integer :: i,j,k
interface
subroutine check(arr)
integer,intent(in) :: arr(:,:,:)
end subroutine check
end interface
省略
call check(arr)
end program main
subroutine check(arr)
integer,intent(in) :: arr(:,:,:)
integer,allocatable :: tmp(:)
integer :: i
tmp=pack(arr,mask=(arr .ne. 0))
do i=lbound(tmp,1),ubound(tmp,1)
print *,tmp(i)
enddo
end subroutine check
------------------------------------------
これでできました。ありがとうございます。
余談だけど、gfortranでは期待通りに処理されるのに、intel fortranだとダメみたい。
pack()の返り値が壊れてしまうっぽい。バグか?
556:デフォルトの名無しさん
09/05/15 00:45:08
>>555
tmpをALLOCATEしてないからかな？
ALLOCATE(tmp(SIZE([arr .ne. 0])))
これを入れておくと動く。
配列を返すものを、未割り当ての割りつけ配列に入れたときの挙動が、
F90の範囲ではプログラム側に問題があって悪いとおもうが、
新しい規格ではかなり自動で調整してくれるようになったので、
これがBugなのかどうか文法書を読み直さないと分からんｗ
[ ] はF2003で、f90なら(/ /) に対応。
557:556
09/05/15 00:49:27
ああ、Fortran95/2003 Explained 17.5.2 に そのものずばり、PACKの結果を
未割り当てALLOCATABLEに突っ込む例題が載っていた。
Fortran95までは出来なくて、F2003で出来るようになるものだそうだ。
あと、>>556で [ ] は無くても良かったｗ
558:デフォルトの名無しさん
09/05/15 01:22:11
>>556.
＞新しい規格ではかなり自動で調整してくれる
へーそうなんだ
いちいちallocateしなくてすむのは嬉しい
559:548
09/05/17 22:24:17
>>556>>557
フォロー、ありがとうございます。
allocateでintel fortranでも無事動くようになりました。
てか、2003いいね。期待大。
560:デフォルトの名無しさん
09/05/21 00:38:15
スレ違いだったら失礼します
FORTRANで他からデータを読み込むのに、そのデータの最初にデータ概要が書いてあるので
読み飛ばす方法があったら教えてください
初心者なので分かりにくくてすみません
561:デフォルトの名無しさん
09/05/21 01:30:02
>>560
適当な変数を定義して、読み込みだけさせて、その中身を
使わなければいいぢゃない。嘘データで良いから実例だし
てくれ。
# リアルなのを晒すとマズいかもしれんし。
562:デフォルトの名無しさん
09/05/21 01:54:07
>>561数字を読み込んで並べるだけなのですが
読み込むデータが
○月×日
気象データ
東京都日黒区
0.0001
0.0025
0.0025
0.0105
0.1568
1.5385
1.4657
12.368
こんな感じで数字の部分のみ読み込みをしたいです
563:デフォルトの名無しさん
09/05/21 02:05:29
データ構造が既知ならば読みたいところまで回していけばいいし、
レコード長が分かってれば指定のレコード部を直接読みに行ける、はず
もっと良い方法あったらﾌｫﾛｰよろ
564:563
09/05/21 02:08:22
あ、sequentialで書いたファイルをdirectで読む（あるいはその逆）って、できないんだっけ？
自信なくなってきた＞＜；
565:デフォルトの名無しさん
09/05/21 02:15:57
>>563すみません、もう少し詳しくお願いできますか？
566:デフォルトの名無しさん
09/05/21 09:45:54
>>562
いくつか方法はある。
最初のデータ概要の行数が決まっていて一定ならば簡単で、
CHARACTER TEXT*256
DO i = 1, 4
READ(7, *) TEXT
END DO
的に読み飛ばせばいい。
手抜きな方法としては、数値変数に文字列を読み込むとエラーが出ることを利用して、
10 CONTINUE
READ(7, *, ERR= 10) dat
で無理に読み飛ばすことも出来る。
本物のエラーのときや、データ概要が数値として読めてしまうと困るｗ
やや高度な方法としては、いったん文字列に読み込んで、その文字列の中身をしらべて
文字列を内部ファイルとして入力に使って数値を読むという方法もある。
データ開始地点まで適当にすすんで、ＢＡＣＫＳＰＡＣＥで適切なところまで戻る手もある。
567:デフォルトの名無しさん
09/05/21 12:20:07
>>566
２つ目の手抜き法はオイラもよく使う。
念のため、
icount = 1
20 continue
read(7,fmt=*) string
read(string,fmt=*,err=10,end=999) dat(icount)
icount = icount + 1
goto 20
10 write(*,*) string
goto 20
continue 999
みたいに文字列を一度かまして
読み込みは内部ファイルで処理して、
飛ばされたケースは表示させて目でチェックするとか。
あまりほめられないスパゲッチ～だけど。
568:デフォルトの名無しさん
09/05/23 20:59:00
>>562
本気でスピード求めるならdirectアクセスがベストだが
多分あなたの使い方じゃそうじゃなさそうだから、普通に次々読んでいって
使わないデータは破棄しちゃえばいいじゃん
>>562がASCIIデータだとすると、5行目からの数字が欲しいなら、
do i=1,5
read(unit) tmp
enddo
とやって捨てちゃえばいい
569:デフォルトの名無しさん
09/05/23 21:16:21
読み捨てるなら
DO I = 1, 5
READ(NUNIT, *)
END DO
でおｋ
570:デフォルトの名無しさん
09/05/23 23:01:20
そっか、捨てるんだから何もいらないか
571:デフォルトの名無しさん
09/06/09 12:40:21
FORTRANがCOLUMN-MAJORなことを調べてみたいのですが、
変数のポインタ値を知る方法はありますか？
572:デフォルトの名無しさん
09/06/09 12:49:09
>>571
無いが、サブルーチンに送ればわかる。
つまり、２次元配列を、サブルーチンの１次元形式引数の配列に渡せばよい。
573:デフォルトの名無しさん
09/06/09 13:18:41
ありがとうございます。確かにWRITEなんかにわたすことで
順番を知ることはできるのですが、メモリ上の位置を知ることは
できないのですね
574:デフォルトの名無しさん
09/06/09 13:32:54
>>573
標準文法ではないが、LOC 関数などがベンダー拡張として用意されていることが多い。
またリンカーでマップを出力させると静的な変数ならばアドレスが分からなくも無い。
配列のならびに関しては、そのままでは分からないが、さらに EQUIVALENCE でも
かければ分かるのかもしれない。
575:デフォルトの名無しさん
09/06/20 14:16:58
IFort10.1+マルチコア環境で使ってる人いるかな
OPENMPで並列化してるんだけど、
デュアルコア環境でトータルCPU使用率が50%から上がらない理由わかる？
１コアしか使ってなくて2コア目がアイドルになっちゃう
576:デフォルトの名無しさん
09/06/20 14:21:55
OSによるだろうけど、
setenv なんたら、
をちゃんとしてる？
577:575
09/06/21 16:49:20
setenvはよくわからないんだけど、OPENMP+自動並列化で5スレッド走ってる
環境はOSがXPProSP3で、Core2DuoE8500
578:デフォルトの名無しさん
09/06/21 18:23:29
OpenMPの場合は、環境変数でCPUの位置をbindしてやらないといけないんだっけ？
自動並列も同じだっけかな？
誰か詳しい人の降臨を待つしかないなｗ
579:デフォルトの名無しさん
09/06/22 09:45:02
>>577
なぜ最初からWinXPだとかけないのだろうか？
自分が使ってるlinux＋ifort9.1では、
setenv OMP_NUM_THREADS 8 (最後の８は適当なCPU数に置き換え)
で環境変数で並列数を指定することになってるけどな。Linuxでもシェルごとに違うはずだが。
実機がないのでWinXPだとどの程度違うのかわからんが、Intel提供の例では
setenv を set に換えるだけのようだ。
全体の監視＆同期とり用のあまり忙しくないスレッドが１個走るはずだけど、
全部で５つはなんか多い気もするがそんなものかもしれない気もする。
>>578
CPUの位置のバインドは自動並列やらOpenMPの場合、Linuxでは dplace や taskset を使うな。
指定しないとふらふら移動することがあるからややこしい。
まあ詳しいWin機でいじくりまくってる神の降臨を待とう。Winでの並列に興味あるし。
580:575
09/06/22 17:58:29
自己解決したので報告
ProcessExplorerでスレッドの状況を調べてたら、メインスレッド以外無負荷だた
並列化を見直して負荷をかけるようにしたら100%まで掛けられるようになった
迷惑掛けてスマソ
581:デフォルトの名無しさん
09/06/22 23:05:58
>>579
WinでのCPUバインドは環境変数でできた。
マニュアルに載っているはず。
ｲ7でハイパースレッディングの８コア運用していると、ちゃんとバインドしないと
性能が全く出ない。デフォで８スレッド運行してバッティングするし、４スレッドにしても
物理コアに割り当たるとは限らないようで放置プレイでは性能でない。
582:デフォルトの名無しさん
09/06/23 16:54:00
>>581
環境変数でCPUバインドっていうのは、別にWindowsに限らず他も同じだと思う。
Intel FortranだとKMP_AFFINITY、PGI FortranだとMP_BLISTかな。
というか、いい加減環境変数の設定とかをユーザにやらせるんじゃないくて、
コンパイラがシステム構成を見ながら勝手に使うCPUを指定してくれればいいのに…。
どういうスレッドの配置が効率がいいかなんて、予想できそうなもんなのに。
583:556
09/07/09 23:41:32
>>555
最近になって、Intel掲示板を見ていてようやく分かったが、オプションをつけないと
駄目らしい。
/assume:realloc_lhs
これをつけるとFortran2003解釈、つけないとFortran90/95解釈になるらしい。
584:デフォルトの名無しさん
09/09/07 13:50:42
age
585:デフォルトの名無しさん
09/09/08 07:43:50
>>583
お、参考にしよう。
おまけ・・・になるのか知らないけどw
メインがCでFortranコードをサブルーチンとして呼び出してリンクするときは
締めの実行ファイル生成時には icc でリンクして実行ファイルつくってもいいけど、ifort で
ifort -nofor_main -no-ipo
ってやると振る舞いが礼儀正しい事が多い、気がする。
-no-ipo はリンクするときだけね。これがないとifortは無茶をするようだw
586:デフォルトの名無しさん
09/09/27 01:29:32
>>nag の fortran 入門でおｋ
587:デフォルトの名無しさん
09/09/27 05:35:56
CERNLIBのmathlibのラッパをモジュールで作ろうと思っているのだけれど、
gfortranで例えば自己完結している共有オブジェクト(.soファイル)って
作れるもんなの？
一々$(cernlib -G Motif mathlib)をつけてコンパイルするの
面倒くさい。
ライブラリ作ったこと無いんでアレゲな文章だが、もし誰かわかるのなら。
ちなみに、OSはUbuntu8.04な。
588:デフォルトの名無しさん
09/09/27 07:43:50
追加で質問なんですが、
URLﾘﾝｸ(wwwasdoc.web.cern.ch)
のラッパーとして
URLﾘﾝｸ(gist.github.com)
というのを書いて
program main
use FunctionBinom
implicit none
integer::x,y
x=5
y=2
print *,binom(x,y)
end program
でコンパイルしてみると（もち、$(cernlib mathlib)はつけて）
/tmp/ccdsgC3z.o: In function `__functionbinom__binom_integer':
binom.f95:(.text+0x14): undefined reference to `kbinom_'
/tmp/ccdsgC3z.o: In function `__functionbinom__binom_double':
binom.f95:(.text+0x34): undefined reference to `dbinom_'
/tmp/ccdsgC3z.o: In function `__functionbinom__binom_real':
binom.f95:(.text+0x54): undefined reference to `rbinom_'
collect2: ld returned 1 exit status
というエラー出る。何で参照先が無いのかさっぱりわからない。
589:デフォルトの名無しさん
09/09/29 00:00:44
>>588
よく知らんが、関数呼び出しの名前規約が食い違っているとか。
エラーメッセージを見る限り、小文字でアンダースコア末尾付与の形式で
オブジェクトを探しに行っているが、確かに存在しているのかね？
名前一覧とか出力させるとか、バイナリエディターで直接のぞいてい見るとか。
590:デフォルトの名無しさん
09/09/29 03:11:28
external binom
つけとく、とか
interface 文で前もって引数部を教えてあげるとか、
_binom みたく下付線つけとくとか、かなあ・・・違うな。う～ん。
591:デフォルトの名無しさん
09/09/29 22:56:52
>>589
>関数呼び出しの名前規約が食い違っているとか。
他の関数（例えば、票にもあるATGという関数）だとエラー無くできるし、
自分がマニュアル読む限り名前規約の話もないのよね。
>名前一覧とか出力させるとか、バイナリエディターで直接のぞいてい見るとか。
バイナリエディタは避けたい。名前一覧どうやって出力させればいい？
>>590
モジュラー性上げたいのでexternalを使いたくないのよ。
引用仕様記述は書いてあるよ。
592:デフォルトの名無しさん
09/09/29 23:04:23
スパコンの人じゃないのでわからないけれど、スパコンの人たちって
CERNLIBとかは使わず、何のライブラリ使っているの？NAG?
フリーであるならそのライブラリのラッパに鞍替えするわ。
593:デフォルトの名無しさん
09/09/30 10:40:07
自作に決まってるだろ
他人のライブラリなんか遅くて使えないし
594:デフォルトの名無しさん
09/09/30 12:41:03
>>593
なぬ。それスパコン業界の一般論？>>593が自作好きというわけじゃなく？
新しいアルゴリズムで作る、というのならわかるのだけれど、
それ以外の速度的にもう上げられないだろう、上げる必要も無いだろう的
関数も自作かい。
595:デフォルトの名無しさん
09/09/30 15:42:43
行列ソルバに関しては、やっぱまだまだ自作の余地があるんじゃないの？
スパコンのアーキテクチャによって最適なアルゴリズムは違うだろうし。
x86なスパコン（クラスタというべきか？ｗ）なら適当な余所様のライブラリでいいかもしれないけど。
596:デフォルトの名無しさん
09/10/01 02:38:44
>>591
＞バイナリエディタは避けたい。名前一覧どうやって出力させればいい？
ライブラリアン用のコマンドを使えば、一覧出せると思う。
メインフレームなら出来たｗ
>>594
ライブラリはその性質からどうしても汎用だから、自分の問題に合わせた
特殊条件を使って早くすることはあるな。
597:デフォルトの名無しさん
09/10/01 12:47:41
なるほど。個人的な話だけど、いろんなことが繋がったわ。
ライブラリは確かに汎用的じゃないといけないわな。
行列ソルバ、つーか固有値計算は興味深い。
アルゴリズム色々調べているけど、最新のアルゴリズムってなによ？
まぁ特許とられて普通に使えないんだろうけど。
よく知らないけれど、アルゴリズムに特許はかけられるよね？
アルゴリズムを実装したアプリケーションに対してのみ特許かけられるんだけっけ？
598:デフォルトの名無しさん
09/10/01 12:48:45
>>596
>ライブラリアン用のコマンドを使えば、一覧出せると思う。
お礼言うの忘れてた。ﾄﾝｸｽ！
599:デフォルトの名無しさん
09/10/01 13:13:45
アルゴリズムなんかで特許とんのは浅ましい・・・
そんなふうに考えていた時期が（ｒｙ
600:デフォルトの名無しさん
09/10/06 10:10:52
do 600 i = 1, 100
write(*,*) '600 get'
600 continue
601:デフォルトの名無しさん
09/10/10 13:11:05
module mod
interface iface
integer function func(i, j)
integer :: i, j
end function func
end interface iface
end module mod
integer function func(i, j)
use mod
integer :: i, j
func = i * j
return
end function func
program main
use mod
print *, iface(10, 20)
end program main
-------------------------------------------------------
func = i * j
1
Error: Name 'func' at (1) is an ambiguous reference to 'func' from current program unit
上記のようなエラーがでてコンパイルできません
moduleを使わない、あるいはfunctionではなくsubroutineならコンパイルできました
ひょっとして総称名はfunctionには使えないんでしょうか？
602:デフォルトの名無しさん
09/10/10 14:29:25
>>601
integer function func(i, j)
の中の
use mod
が不要なのでは
603:デフォルトの名無しさん
09/10/10 15:44:58
>>602
できました。ありがとうございます
参考にしたサイトが間違っていたようですね
604:デフォルトの名無しさん
09/10/10 16:32:59
....その言い草はねぇだろうww
605:デフォルトの名無しさん
09/10/10 22:01:30
>>604
だったらそのサイトのURLを教えて欲しい
606:デフォルトの名無しさん
09/10/10 22:03:42
安価間違いました>>604じゃなく>>603です。すみません。
607:デフォルトの名無しさん
09/10/10 22:14:02
さらしていいものなのかな？
URLﾘﾝｸ(www.geocities.jp)
ここなんだけど副プログラムのほうにもuseしろって書いてありました
608:デフォルトの名無しさん
09/10/13 03:37:16
まずはgeocitiesごときにあるサンプルを鵜呑みにした自分を責めるところから始めるんだ
609:デフォルトの名無しさん
09/10/14 00:33:47
>>609 (インデント用)
module mod
contains
real function call_func(func, arg)
interface
real function func(arg)
real, intent(in) :: arg(2)
end function func
end interface
real, intent(in) :: arg(2)
call_func = func(arg)
return
end function call_func
end module mod
real function func(arg)
real, intent(in) :: arg(2)
func = arg(1) * arg(2)
return
end function func
program main
use mod
external func
real :: arg(2) = (/ 2.0, 3.0 /)
! dummy = func(arg)
print *, call_func(func, arg)
stop
end program main
610:デフォルトの名無しさん
09/10/14 00:36:31
>>609のコメントアウトした部分のコメントをはずすとコンパイルも通り、正常に動作するんですが
コメントアウトするとなぜかコンパイルすらできません
理由がさっぱりわからないんですが、なにが問題なんでしょうか？
コンパイラはgcc 4.2.4
エラーメッセージは↓です
print *, call_func(func, arg)
1
Error: Type/rank mismatch in argument 'func' at (1)
611:デフォルトの名無しさん
09/10/14 00:43:12
>>610
Intel VF 11.046 なら通る。
612:デフォルトの名無しさん
09/10/14 01:54:59
回答有難う御座います
gccコンパイラのバグということでしょうかね
汎関数を引数にもってる関数の総称名をつけようとするとコンパイルできなかったり(これは仕様？)
613:デフォルトの名無しさん
09/10/14 03:35:23
>>612
PRINT文のバグ臭いな。
FUNCTIONの返り値が、VALUEなのかADDRESSなのか今一不明なので
７７時代からFUNCTIONの返り値を代入を経ずに使うと、機種依存の振る舞いを
することが良くあった。
最近は大丈夫なことが多いが。
関数引数の総称名は、IVFでは出来た。多分、仕様上では出来ておかしくないと思う。
MODULE m_unko
INTEGER, PARAMETER :: kd = SELECTED_REAL_KIND(14)
INTERFACE gen
MODULE PROCEDURE sub1, sub2
END INTERFACE
ABSTRACT INTERFACE
REAL FUNCTION x_f(a)
REAL, INTENT(IN) :: a
END FUNCTION x_f
REAL(kd) FUNCTION d_f(a)
IMPORT
REAL(kd), INTENT(IN) :: a
END FUNCTION d_f
END INTERFACE
614:556
09/10/14 03:36:04
CONTAINS
REAL FUNCTION sub1(f, x)
REAL, INTENT(IN) :: x
PROCEDURE (x_f) :: f
sub1 = f(x)
RETURN
END FUNCTION sub1
REAL(kd) FUNCTION sub2(f, d)
REAL(kd), INTENT(IN) :: d
PROCEDURE(d_f) :: f
sub2 = f(d)
RETURN
END FUNCTION sub2
END MODULE m_unko
!
PROGRAM unko
USE m_unko
REAL :: x
REAL(kd) :: d
PROCEDURE (x_f), POINTER :: f
PROCEDURE (d_f), POINTER :: g
x = 4.0 * ATAN(1.0)
d = 4.0_kd * ATAN(1.0_kd)
f => SIN
g => DSIN
PRINT *, gen(f, x)
PRINT *, gen(g, d)
STOP
END PROGRAM unko
615:デフォルトの名無しさん
09/10/14 21:01:13
>>613,612
サンプルコードありがたいんですが、コンパイルできませんでした
abstract interfaceのところで最初のエラーがでるので原因を調べたんですが
hpとintelのマニュアルを検索してもabstract interfaceの項目が見つからりませんでした・・・
616:デフォルトの名無しさん
09/10/15 01:08:20
>>615
ABSTRACT INTERFACE はF2003の機能で、まだIntelのHelpに載ってないが、実装されている。
Crayとか（多分）IBMのマニュアルなら載っていると思う。まぁ素直にFortran95/2003 Explainedを
見るのがいいと思うが。GoogleBooksでもページが飛ばされてなければ見られる。
しかし、これは本質ではない。二度Interfaceを書くのが面倒だったから使ったまで。
これならどうだべ？
MODULE m_unko
INTERFACE gen
MODULE PROCEDURE sub1, sub2
END INTERFACE
CONTAINS
REAL FUNCTION sub1(f, x)
REAL, INTENT(IN) :: x
INTERFACE
REAL FUNCTION f(x)
REAL, INTENT(IN) :: x
END FUNCTION f
END INTERFACE
sub1 = f(x)
RETURN
END FUNCTION sub1
617:デフォルトの名無しさん
09/10/15 01:09:15
INTEGER FUNCTION sub2(f, x)
REAL, INTENT(IN) :: x
INTERFACE
INTEGER FUNCTION f(x)
REAL, INTENT(IN) :: x
END FUNCTION f
END INTERFACE
sub2 = f(x)
RETURN
END FUNCTION sub2
!
REAL FUNCTION rmult2(x)
REAL, INTENT(IN) :: x
rmult2 = 2.0 * x
RETURN
END FUNCTION rmult2
!
INTEGER FUNCTION imult2(x)
REAL, INTENT(IN) :: x
imult2 = INT( 2.0 * x )
RETURN
END FUNCTION imult2
END MODULE m_unko
!
PROGRAM unko
USE m_unko
REAL :: x
x = 4.0
PRINT *, gen(rmult2, x)
PRINT *, gen(imult2, x)
STOP
END PROGRAM unko
618:デフォルトの名無しさん
09/10/15 18:46:55
>>616,617
有難うございます。コンパイルできました。完璧です
主プログラムにinterfaceを書く方法は見つけたんですが
moduleにくくってしまえばinterfaceを書かなくてもできるんですね
619:デフォルトの名無しさん
09/10/16 08:07:09
FORTRANやめますか？人間やめますか？
620:デフォルトの名無しさん
09/10/16 08:50:16
つ DarkBASIC
621:デフォルトの名無しさん
09/10/16 23:00:56
Windows版NetBeansの環境でgFortranを使いたいんだが、
実行ファイルの設定をg77.exeからgfortran.exeに変えるだけでok？
もし、そのほかに設定事項等あれば教えて欲しい
622:デフォルトの名無しさん
09/10/17 20:44:18
moduleでオブジェクト指向っぽいことができるのはわかったんですが
c++とかでいうインスタンスを作ることはできないんでしょうか？
623:デフォルトの名無しさん
09/10/17 22:43:57
>>622
F2003なら出来る。F95以前では出来ない。
624:デフォルトの名無しさん
09/10/17 23:13:19
こんな感じかな？
とりあえず、IVF 11.1 では通った
誰か Fortran 95/2003 Explained の日本語訳作って
module mod1
implicit none
type Circle
real :: d
contains
procedure :: area => func1
end type Circle
contains
real function func1(this)
class(Circle) :: this
func1 = (this%d ** 2) * 3.14 / 4.0
end function func1
end module mod1
program main
use mod1
implicit none
type(Circle) :: c
c%d = 10.0
print *, c%area()
end program main
625:デフォルトの名無しさん
09/10/18 16:52:26
>>624
ありがとうございます
残念ながらgfortranでは2003が使えないようなので試せませんが・・・
626:デフォルトの名無しさん
09/10/19 14:40:01
実際、FortranでOOPはどうなの？
イラネって意見も散見するが・・・
627:デフォルトの名無しさん
09/10/19 16:11:22
乱数とか行列？
628:デフォルトの名無しさん
09/10/19 17:25:58
>>626
それで速度が落ちるならイラネ。落ちそうな気しかしないけど。
中田育男『基礎FORTRAN』（岩波書店）のコラムに
FORTRANでもCOMMON文とENTRY文を使って、オブジェクト指向設計に
従ったプログラムを書くことができる。
とあったけど、具体的なことは何も書いていなかった。
629:デフォルトの名無しさん
09/10/19 17:54:29
FORTRANで無理矢理それっぽいことをやったらどうなるかという例が
「オブジェクト指向入門」34章にあるようだ
URLﾘﾝｸ(www.seshop.com)
俺は初版しか読んでないけど、今でも事情は変わらんでしょ
630:デフォルトの名無しさん
09/10/19 20:01:44
>>629
読んでみたいけど、ちょっとお値段が・・・
そういやまともなオブジェクト指向の本を読んだことは一度もないや
良い機会だから買ちゃおうかしら
631:デフォルトの名無しさん
09/10/19 20:18:28
>>628
> COMMON文とENTRY文を使って、オブジェクト指向設計に
なぜかはわかりませんが、イヤな予感しかしません・・・
632:デフォルトの名無しさん
09/10/19 20:42:47
というか、OOP用に作られてない言語で安全なOOPをするのって無理だし・・・
安全でないOOPなんてOOPのメリットないに等しいし・・・
素人は大人しくJavaやC#使うほうが幸せになれる
633:デフォルトの名無しさん
09/10/19 21:13:30
>>629
初版のはFORTRAN77でトリッキーなカプセル化をする話だけだったから
たぶん>>631>>632の言う通り
今実際どうなのか
俺は2003とか不勉強なんでわからん
634:デフォルトの名無しさん
09/10/20 00:37:24
>>633
そういわず、２００３を勉強しようぜ！
ばりばりOOPでなくても
TYPE の拡張とかが、入れ子じゃなくて、同じ深さで構造体のメンバーを増やせるし。
OOPでない部分でも結構拡張があって便利。
日本語 パワポ
URLﾘﾝｸ(www.itscj.ipsj.or.jp)
英語
fURLﾘﾝｸ(ftp.nag.co.uk)
635:デフォルトの名無しさん
09/10/20 01:07:43
正直、OOPとそうでないものの違いがよくわかんない。
暴力的な書き方だけど、「因果関係が一意」という点では同じだし。
C++しか扱ったことないしあまり深く触らなかったからだけど、
OOPの深淵に辿り着かないと判らない機微もあるのだろうけどな。
636:デフォルトの名無しさん
09/10/20 02:20:34
>>629
どんな感じのことを言ってんすか？
簡単に教えてくだしあ
637:デフォルトの名無しさん
09/10/20 05:21:22
>>636
たいした話じゃなくて
ENTRY文をメソッドの代わりにしてカプセル化をまねるだけ。
継承とかそういう話はないし、動的にオブジェクトを作れるわけでもない。
本来の使い方でもないから危なっかしい。
ただ、それでも素のPascalよりは気がきいてるってことなんだなｗ
今じゃ2003もあるし
第二版はどうなってるか知らないよ。
638:デフォルトの名無しさん
09/10/20 05:43:47
>>635
自分で書く分にはOOPは邪魔なだけだよ
OOPが真価を発揮するのは他人のライブラリを使うとき。
639:デフォルトの名無しさん
09/10/20 12:58:28
ENTRY文はバグでやすいな。
COMMONの代わりに、共有変数を実現するのに、７７でENTRY使ったが、かえってめんどかったｗ
640:デフォルトの名無しさん
09/10/21 10:41:26
>>637
thx
2003でOOPを謳ってるってことは、ある程度まともには作られてんだろうね多分。
アドレス指定で偶然にそのインスタンスを破壊しちゃう、なんてリスクはないように
出来てるんだろう・・・（どうやって実現してんのかさっぱり想像できないけど）
641:デフォルトの名無しさん
09/10/21 21:31:44
すいません。最近始めた超初心者です。
簡単なシミュレートで
open(10,file='testout.out',status='unknown',form='formatted')
do n=0,nmax
計算略
write(10,30) t, a(1),a(2),a(3)
30 format(1H ,4(1pe12.4))
enddo
こんな感じで配列の値を表示しようと思うのですが
なぜか出力ファイルが改行されません＞＜
g77でコンパイルしています。コンパイルはできるのですが・・・
どなたかファイルに出力するときの改行の仕方をご教授していただきたいのです
お願いします
642:デフォルトの名無しさん
09/10/22 18:39:19
>>641ですが
出力ファイルをdosに直してませんでした
自己解決しますた
フォートランはフォーマット文がわかりにくいなぁ
643:デフォルトの名無しさん
09/10/24 00:54:54
>>642
FORMATの書き方が古い。
まぁｇ７７だから仕方ないのか。
ｆ９０で学びなおせばFORMATもそれほどつらくない。
644:デフォルトの名無しさん
09/10/24 03:17:58
Fortran77では動的メモリ割付け(allocate)は使えないの？
645:デフォルトの名無しさん
09/10/24 03:48:06
>>644
使えない。
646:デフォルトの名無しさん
09/10/25 00:06:12
>>644
大きくCOMMONを取って、必要に応じて自分で切り出して使うんだ。
647:デフォルトの名無しさん
09/10/25 00:56:39
>>644
大きくCOMMONがとれない環境なのです。
allocateしないと最大メモリまで使えない。
648:デフォルトの名無しさん
09/10/25 16:16:54
>>647
むかし過渡期の技術としてダイナミックコモンというのがあった。
ひょっとしたらコンパイラのオプションにあるかもしれない。
その場合、普通のスタティックなコモンより大きく領域が取れるかもしれない。
649:デフォルトの名無しさん
09/10/27 14:20:53
>>646
そんなこと教えちゃらめ～
650:デフォルトの名無しさん
09/10/28 20:57:22
allocateを使うとパフォーマンスが落ちるって本当ですか？
実際に使われるサイズに関わらず？
651:デフォルトの名無しさん
09/10/29 00:56:45
>>650
多分誤解。
静的配列と比べるならば、それは実行時にアロケート処理するのだから余計な時間はかかるだろう。
しかし、それは計算時間に比べれば瑣末なはず。
POINTERで配列をとった場合は、メモリーの割付が連続したメモリー領域になるとは限らないので
最適化がバリバリなされず、ALLLOCATEでとった場合（連続にメモリーを取る）とくらべて
パフォーマンスが落ちる可能性があるとは言われている。