FORTRAN �W

1 名前：デフォルトの名無しさん [2006/11/12(日) 17:13:35 ] FORTRAN総合スレッドです。 FORTRAN77/90/95/2003, HPF等、FORTRAN全般に関する話題を書き込んで下さい。 ・FORTRAN関連情報 お約束だが、まずは自分で調べよう www.google.co.jp/ x86で動作するFortranコンパイラ www.nminoru.jp/~nminoru/programming/x86-fortran.html The Fortran Company www.fortran.com/ Polyhedron Software www.polyhedron.co.uk/ Google Directory - Fortran directory.google.com/Top/Computers/Programming/Languages/Fortran/ ・前スレ FORTRAN III pc8.2ch.net/test/read.cgi/tech/1104724162/ FORTRAN II pc5.2ch.net/test/read.cgi/tech/1068351911/ Fortran pc2.2ch.net/test/read.cgi/tech/1003214738/ 2 名前：デフォルトの名無しさん [2006/11/12(日) 17:14:45 ] ・フリー／無料評価版のコンパイラ The G95 project （Fortran95） www.g95.org/ GCC (g77) （今のところ FORTRAN77） gcc.gnu.org/ Salford FTN77/95 Personal Edition （Fortran77/95, Windows, 要登録） www.salfordsoftware.co.uk/software/downloads/compilers.html FTN77/95のインストール方法 www.coastal-env.k.u-tokyo.ac.jp/koibuchi/fortran.htm the F Programming Language （Fortran90のサブセット, Win/Linux/Solaris） www.fortran.com/F/compilers.html Open Watcom Fortran Compiler （FORTRAN77） www.openwatcom.org/ Intel Fortran Compiler （Fortran95, 無料評価版(Win)／非商用バージョン(Linux), 要登録） www.xlsoft.com/jp/products/download/download.html Lahey/Fujitsu Fortran （期限付きトライアル版は終了した模様） www.lahey.com/downloads.htm ・ドキュメント Compaq Visual Fortran オンライン・ドキュメント www.xlsoft.com/jp/products/intel/cvf/docs/index.html Intel Fortran Compiler 日本語版マニュアルのダウンロード （要登録） www.xlsoft.com/jp/products/intel/download.html 富士通Fortranオンラインマニュアル www.cc.nao.ac.jp/fsunman/japanese/Fortran/index.html FORTRAN77 standard www.swcp.com/~walt/F77_std/f77_std.html www.fortran.com/fortran/F77_std/rjcnf0001.html JISデータベース www.jisc.go.jp/ 3 名前：デフォルトの名無しさん [2006/11/12(日) 17:48:19 ] 　　ヽ／l l　ニ|ﾆ　　　　　　　　　　 ,.､-''"..;:;:;:;:;:;:;:... `'ヽ､ 　　（　　　（￣　　￣）　　　　　 ／....:::;:;:;:;:;:;:;:;:;:;:;:;:;:;:;:;:;.....ヽ､／￣￣￣￣＼／ 　　　￣　　 ￣　　　　　　　　 i＿;;､:＿;､;＿;､;､;､､ｨｯ.;:;:;:;:;: /　　興　　F　　　　 君 　,.､-──-- ､.,_　　　　　,､　　|　　　　　　,,,,,,　　/ ;:;:;:;:;:;: |　　 味 　 O　　　 　 ： 　　　　　　　　　　｀`''--イ ,）,､,! '''''　　　　　　　　＼ ;;;;;;;;;_|　　 が 　 R　　　 　 ： 　　　　　　　　　　　　　 ヾｰ'ﾞ |ﾋﾆﾆュ　ｬﾆ,ニニ､＞ 〉;; / _|　　 あ　　T 　　　　 　 ﾉ　l　 ﾊ　 l　ヾ トイ ｀!ﾞl）_j　　 '　iﾘ__, ｀　 }ii l f'ﾄ〉 　 る　　R 　　　 _,,.ノ　_ﾉ / ﾉ ノ　ノﾉ!_丿　| 　 l　　　｀ " '''　　 }ii ﾘﾉ |　　 の 　 A　　|＼＿＿ 　　　　　 ノ　,､ｨ'-＝z=F [_　　 .l!　.{　　　、　　　　 ｨ!ii;}' ﾉ|　　 か　　 N 　 | 　　　-‐''ﾞ_ノ ,ノ　 'ﾞ (ｿ 　 ヽ　　 {!　ﾞｰ<⌒'　　　　 ,ﾐi;i;}ｰ'ﾞ |　　 ね　　 　 　 | 　　､ ィッ＞f｢　　_,,二-　　 ヽ.　 }i､　 -=＝=-'　　ﾘiii;ﾂ 　 |　　 ？　　 　 　 | 　 　｀〒T〔!|　　　　　　　r ,_ﾉ　_ﾉ}lli,　　-r=‐　　,i;llilili|　　 ＞ ＿＿＿＿_／｀ヽ､ 　　 　 ﾞ､ヽ`!　 l　　　＿ _｣　／/　'}llli,　,;i|i;,　,,ii;ilililll'ﾞﾘ ／ ￣ l l　　　　　　,､ ''⌒ﾞヽ、 　　　　 `ﾄ.､!　 lj　　（＿_l、/　 | 　 ﾞﾄ!llllllllllliillllllllﾂﾞ_､-ﾞ　 　 /,l l　　 　 　 / 　　　　　l　ﾞﾄ、　　　　　t'ﾞ | 　 |　　| |､'lﾄllﾄllﾄﾂ "´　　　 ／/ l l 　 　 　 / 　　 ,ｨ、化ｧ ',＼ 　 　 　 l　〉　 | 　 | | ﾞ､ /／∧　　　 ／ /　 l l　　　　 l （爪(（､`ｰ'′ ',　｀''t‐--'′〉ﾄ､ | 　 | |. Vﾊ彡 ∧　 ／　　　　|_L､　　i　|　／ ￣￣￣｀¨`''ー--ﾆL_　　　`!､ `!　l￣`''┴--┴'-'ﾞ-─…''"´_,, ｨ|　　l　|/ 　　　　　　　　　　　 �h　　 ヽ ｀ ', | Ｏ｀`''────…'''"´ 　Ｏ',ﾞ:､ 　l　|／ 4 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 18:34:32 ] 機種依存文字使うなボケ失せろ 5 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 19:07:04 ] >>4 きたねー言葉使うな糞野郎！ 死ね！ 6 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 19:30:21 ] 数値積分が使えるフリーの数値計算ライブラリありませんか？ 7 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 19:58:47 ] >>6 NetLibへ池。６６時代のものがごまんとあるぜ。 8 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 20:01:14 ] IMSLがあるんじゃね？ 9 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 20:02:08 ] >>8 IMSLは有料パッケージソフトだぞ。 10 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 20:08:40 ] IMSLは10万以上する予感 LAPACKに積分はふくまれておらんの？ 11 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 20:08:54 ] コンパイラより高い・・・ 12 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 20:10:38 ] このスレのエキスパートたちは有料パッケージとか使っておられるのですか？ 13 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 20:13:17 ] >>7 Netlib行ったが数が多すぎて何がなんだかわからない。 これをどう使えばいいかすらわからない。 見た感じプログラム単体が落っこちているだけにしか見えないが どうやってってこれを呼び出して使うんだろうか。 14 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 20:32:24 ] >>10 LAPACKは元々、LINPACK＋EISPACKで、それぞれ密行列の線形方程式と固有値問題の ルーチン集なので積分は主題から外れていると思う。 LINPACK、EISPACKは６０年代の米英の国家プロジェクトだったように記憶している（記憶曖昧ｗ）。 LAPACKは基本的にそれを７７化してBLASの上に統一化し、一部アルゴリズムを改良・拡充した ものだろう。その辺の経緯や順序はよくわかんね。 >>13 NetLibとはそういうところだｗ　埃の積もった土蔵見たいなもんだ。 結局、膨大すぎて何が欲しいか分かってないと使えないｗ NumericalRecipesのページへ行ってみるのが吉かもしれない。 とりあえず専門家的にはベストで無いにしろ、素人が馬鹿な外れを引かされることは無いはず。 ただ本版のソースはしばしばBug入りなのでめくらめっぽう使うのはﾔｳﾞｧｽ 有料ソースは直っていたりするｗ 15 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 21:19:12 ] 初心者は富士通のサブルーチン付きコンパイラ買うのがbestじゃね？ 難易度も予算も 16 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 21:25:25 ] Intel Math Kernel Libraryでok 17 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 21:36:09 ] MKLって積分ルーチン入ってるの？ 18 名前：デフォルトの名無しさん mailto:sage [2006/11/12(日) 21:39:47 ] 多分ない というか、フリーで有名どころのパッケージには積分ルーチンは無い予感 19 名前：デフォルトの名無しさん mailto:sage [2006/11/13(月) 19:48:52 ] 以下のプログラムを実行すると、 結果が 1016010618 となります。 なぜでしょうか？訳が分かりません 環境はg95@cygwinです。 program Interpolation print *,Lorentz(1.5) end program Interpolation real function Lorentz(x) real :: x Lorentz = 1/(1+25*x*x) return end 20 名前：デフォルトの名無しさん mailto:sage [2006/11/13(月) 20:30:25 ] メインプログラムにLorentzの型定義がないからinteger扱いされてんじゃねーの？ 21 名前：19 mailto:sage [2006/11/13(月) 21:30:09 ] implicit noneで解決しました。 １時間も時間を食いました。 22 名前：デフォルトの名無しさん [2006/11/14(火) 00:43:00 ] 自由書式で書くなら関数・サブルーチンは内部手続きにするか、MODULEにするか INTERFACE宣言してから使おう。 23 名前：デフォルトの名無しさん mailto:sage [2006/11/14(火) 01:07:20 ] >>20,22 が正しい。 >>19はあまりに中途半端。きっちりと厳格なF90式で書くか、F77式に暗黙の型を守るべき。 24 名前：デフォルトの名無しさん mailto:sage [2006/11/14(火) 02:41:24 ] INTEL　FORTRAN用のBLAS/LAPACKだが、AMDのcore-math-libraryにWindows用もうｐされている。 昔はLinux用の一部ルーチンだけだったが、いまや全ルーチンを網羅したようだ。 漏れはAMDなので、誰かIntelChipでも使えるか試してくれ。 developer.amd.com/acml.jsp 25 名前：デフォルトの名無しさん [2006/11/14(火) 15:36:00 ] fortran中で他のexeを実行できるようにするにはどうしたら良いですか？ fortranで作られたexeなんだけどソースくっつけるのがめんどくさいんで・・・ｗ 26 名前：デフォルトの名無しさん mailto:sage [2006/11/14(火) 16:32:26 ] >>23 FORTRAN(90)の書籍で、定番の物ってありますか？ 27 名前：デフォルトの名無しさん [2006/11/14(火) 16:53:09 ] これからフォートランを勉強するのですが、 とても解りやすい入門書を教えていだたけませんか？ プログラミング自体、あまり経験がありません。 Ｃを少しやりました。 よろしくお願いします。 28 名前：デフォルトの名無しさん [2006/11/14(火) 18:18:30 ] ｘの値を、1,10,100,1000,10000,100000,1000000の順に変えて その常用対数、自然対数の値を計算して表示するプログラムを作れ という課題が出たのですがサッパリ分かりません。 どうかお力添えください。 29 名前：デフォルトの名無しさん mailto:sage [2006/11/14(火) 20:51:20 ] >>24 Core Duoでもｲｺﾞｲﾀ 32bit ifort版はシングルプロセッサ用に最適化されているみたいだが SSE/SSE2が使えればIntel CPUでもOKのようだ。 >>25 SYSTEM とかそんな感じの関数なりサブルーチンなりが用意されている（ことが多い）。 >>28 X の値を 1 に変える（X に 1 を代入する） X = 1. X の常用対数を計算して Y に代入する Y = LOG10(X) X の自然対数を計算して Z に代入する Z = LOG(X) X, Y, Z の値を順に表示する WRITE(*,*) X, Y, Z 30 名前：デフォルトの名無しさん mailto:sage [2006/11/14(火) 23:23:10 ] 倍精度実数と4倍精度実数では実行速度は変わってしまいますか？ 31 名前：デフォルトの名無しさん [2006/11/14(火) 23:28:24 ] >>30 死ぬほど違う。　 倍精度までは数値演算プロセッサの助けでハードウェア演算するが、 ４倍精度はソフトウェアで浮動小数点演算をするので遅くなる。 ６４BitCPUなら若干状況は改善するが本質においては変わらない。 DECのアルファーCPUは４倍精度も微妙に早かったが、特殊な命令セットがあったためらしく、 ６４bitCPUだからといってよくなるわけでもないようだ。 32 名前：デフォルトの名無しさん [2006/11/14(火) 23:31:25 ] >>29 DO文で書かないといけないのですが どうすればいいんでしょうか・・ 33 名前：30 mailto:sage [2006/11/15(水) 00:03:03 ] >>31 詳しいﾚｽｻﾝｸｽ！ 34 名前：デフォルトの名無しさん mailto:sage [2006/11/15(水) 00:35:58 ] >>26 残念ながらあまりいい本は無い。 今、手に入りやすいのはこれ。 www.amazon.co.jp/gp/product/4563014095 しかし、Ｆ７７の癖を引きずっていて、あまりいい書き方をしていない。 著者は京大のおっさんらしく、ＷＥＢで同等の内容のものを見られる。 アドレスはアマゾンのレビューに書いてある。 英語でよければ、この辺とか、 www.lahey.com/elfpage.htm この辺のＥｓｓｅｎｔｉａｌ Ｆｏｒｔｒａｎあたりが学習用サブセットになっていて、純粋Ｆ９０を学べる。 www.fortran.com/F/books.html しかし、現実世界ではＦ７７ソースがあふれているので微妙ではある。 実際上記のものは教育用としてもはやらなかった。 35 名前：デフォルトの名無しさん mailto:sage [2006/11/15(水) 00:43:57 ] >>32 　　　PROGRAM chinpopo 　　　DO 10 i = 0, 6 　　　 x = 10.0**REAL(i) y = LOG10(x) z = LOG(x) WRITE(6, *) INT(x), y, z 10 CONTINUE 　　　END >>35 36 名前：デフォルトの名無しさん [2006/11/15(水) 14:47:06 ] 　　complex(8) c において， 　　c=(0.d0, 0.d0) として初期化するのも， 　　c=0.d0 として初期化するのも全く同じことなのですか？ 37 名前：デフォルトの名無しさん [2006/11/15(水) 17:59:05 ] DO文を使用して，xの値を，1,10,100,1000,10000,100000,1000000 の順に変えて その常用対数および自然対数の値を 計算して表示するプログラムを作れ． なお，計算結果は書式の指定をして，x, log10(x), log(x) の順に表形式で表示すること． また，計算結果が整数となるものは整数型(In)，実数となるものは実数型(Fn.m)または指数型(En.m)で出力すること． program report5 implicit none integer :: x real :: y,z write (*,*) 'x=' read (*,*) x,log(x),log10(x) do 10 i = 0.6 x=10.0**real(i) y=log(x) z=log10(x) write (6,) end do stop end program report5 何処がおかしいんでしょうか? 38 名前：デフォルトの名無しさん mailto:sage [2006/11/15(水) 18:31:55 ] >>37 read (*,*) x,log(x),log10(x) どこに入れたいの？ たぶん、その穴じゃないわよ 39 名前：デフォルトの名無しさん mailto:sage [2006/11/15(水) 19:20:43 ] >>37 >read (*,*) x,log(x),log10(x) >do 10 i = 0.6 >x=10.0**real(i) >y=log(x) >z=log10(x) >write (6,) ここら辺が間違ってる。 どこが間違いなのかわからないときは、とりあえずコンパイラ様にソースを渡して 文句を言われたところを上から順に直していくのも１つの手だ。 だれかくだすれの次スレ立ててくれ。 40 名前：デフォルトの名無しさん [2006/11/15(水) 23:55:01 ] >>37 >>35をもういっぺん見ろ！ >>39 重複スレをクダスレにするか？ 41 名前：デフォルトの名無しさん mailto:sage [2006/11/17(金) 00:48:27 ] >>40 タイトル一緒だから混乱しそうｗ 42 名前：37 [2006/11/17(金) 10:16:25 ] >>35 WRITE (6,*) が分かりません >>37の問題の書式の指定はこれでいいんでしょうか 43 名前：デフォルトの名無しさん [2006/11/17(金) 11:23:12 ] WRITE(6,　'(I10, F15.7, F15.7)') x, y, z もう少し教科書嫁！！ こんな感じでｗ 44 名前：デフォルトの名無しさん mailto:sage [2006/11/17(金) 12:00:23 ] >>37 がんばったけど書式がうまくいかなかった。ゴメン c DO文を使用して，xの値を， c 1,10,100,1000,10000,100000,1000000 の順に変えて c その常用対数および自然対数の値を 計算して表示する c 計算結果は書式の指定をして，x, log10(x), log(x) の順に c 計算結果が整数となるものは整数型(In)， c 実数となるものは実数型(Fn.m)または指数型(En.m)で出力． Program report5 Implicit none INTEGER i REAL x do i = 0, 6 x = 10.0**i write(*, *) aint(x), log10(x), log(x) end do stop end c 実行結果 c 1.00000 0.000000E+00 0.000000E+00 c 10.0000 1.00000 2.30259 c 100.000 2.00000 4.60517 c 1000.00 3.00000 6.90776 c 10000.0 4.00000 9.21034 c 100000. 5.00000 11.5129 c 0.100000E+07 6.00000 13.8155 45 名前：デフォルトの名無しさん [2006/11/17(金) 12:28:48 ] program report5 implicit none integer :: x real :: y,z write (*,*) 'x=' DO 10 i = 0, 6 x = 10.0**REAL(i) y = LOG10(x) z = LOG(x) WRITE(6,　'(I10, F15.7, F15.7)') x, y, z 10 CONTINUE END do stop end program report5 こうですか? 46 名前：デフォルトの名無しさん [2006/11/17(金) 17:12:34 ] Program report5 Implicit none real :: i,x do i=0, 6 x=10.0**i write(*, '(i7, f15.7 f15.7)') x, log10(x), log(x) end do stop end 結果 ******* 0.0000000 0.0000000 ******* 1.0000000 2.3025851 ******* 2.0000000 4.6051702 ******* 3.0000000 6.9077554 ******* 4.0000000 9.2103405 ******* 5.0000000 11.5129251 ******* 6.0000000 13.8155107 おかげさまでここまで出来たんですが、どうしても*******が出てしまいます＞＜ 47 名前：デフォルトの名無しさん mailto:sage [2006/11/17(金) 19:09:14 ] >>46 「変数の型」を再勉強してから出直してこい。 48 名前：デフォルトの名無しさん mailto:sage [2006/11/18(土) 00:05:43 ] >>46 write(*, '(i7, f15.7 f15.7)') x, log10(x), log(x) 　　　↓ write(*, '(f15.7, f15.7 f15.7)') x, log10(x), log(x) にしてみよう・・・ xはrealで宣言されてるからformatを整数表示に指定することはできない． 49 名前：デフォルトの名無しさん [2006/11/18(土) 00:44:07 ] >>48 >>37の問題を読む感じでは、 write(*, '(i7, f15.7 f15.7)') int(x), log10(x), log(x) を予期していると思われる。 >>46はもう少し貰った回答例をよく吟味しろｗ 50 名前：デフォルトの名無しさん [2006/11/18(土) 01:00:15 ] ＰｕｋｉＷｉｋｉ　つくってみた www.moecity.jp/home/fortran/ ｗｉｋｉよく知らない。自由に使ってみてくれ。 51 名前：デフォルトの名無しさん mailto:sage [2006/11/18(土) 19:42:09 ] f90スタイルで学べるいい文書はありませんかね。 英語か日本語でお願い 52 名前：デフォルトの名無しさん mailto:sage [2006/11/18(土) 20:03:11 ] >>51 Fortran 標準コーディングルール ttp://www.mri-jma.go.jp/Project/mrinpd/coderule.html European Standards For Writing and Documenting Exchangeable Fortran 90 Code ttp://www.meto.gov.uk/research/nwp/numerical/fortran90/f90_standards.html 上の和訳 ttp://www.mri-jma.go.jp/Project/mrinpd/eurostand.html 53 名前：デフォルトの名無しさん [2006/11/18(土) 21:18:26 ] >>52 ＞推奨する書き方 real:: x(n), y(n) ＞x(1:n) = y(1:n) ＞ ＞ ＞×推奨しない書き方 x(:) = y(:) または x = y これはかなりうそ臭いと思う。 ｘ＝ｙ　の場合と違って、コロン指定子を使うと、いったん暗にコンパイラが テンポラリ配列にコピーするのでたいてい余分な時間がかかる。 しかし漏れが疑問に思うのは、F77では変数と配列に同じ名前をAとA(10)を指定できて 別物に扱われる様に規格で決まっていたと思うのだが、F90で矛盾が起きないのか・・・ それとも、これはJIS規格だけでANSIやMILでは違っていたのか？ 54 名前：デフォルトの名無しさん mailto:sage [2006/11/18(土) 22:09:41 ] >>53 REAL A DIMENSION A(10) DO 10 I=1,10 A(I) = REAL(I) 10 CONTINUE A = 1.23 WRITE(*,*) A WRITE(*,*) (A(I), I=1,10) END 規格上どうなってたかは知らないが g77だと A = 1.23 の部分でコンパイルエラーになった。 g95では（当然） 55 名前：54 mailto:sage [2006/11/18(土) 22:12:21 ] × g95では（当然） ○ g95では（当然） A は配列全体と解釈された。 56 名前：デフォルトの名無しさん [2006/11/18(土) 23:02:07 ] >>54 そうじゃなくて REAL A, A(10) が可能だったと思う。 57 名前：デフォルトの名無しさん mailto:sage [2006/11/19(日) 00:16:02 ] hoge.f: In program `MAIN__': hoge.f:1: 　　　　 REAL A,A(10) 　　　　　　　　　1 2 Invalid declaration of or reference to symbol `a' at (2) [initially seen at (1)] 58 名前：デフォルトの名無しさん [2006/11/19(日) 01:39:12 ] たしか変数名と配列名は規格上の定義では別物だった記憶がある。 実際のコンパイラでどうなっていたかは知らん。 59 名前：デフォルトの名無しさん mailto:sage [2006/11/20(月) 08:30:40 ] Fortran90なんですが， integer a(1:100)=0 みたな初期化の仕方はよくやるんですか？ 60 名前：デフォルトの名無しさん mailto:sage [2006/11/20(月) 08:32:30 ] ×みたな ○みたいな orz 61 名前：デフォルトの名無しさん [2006/11/20(月) 09:01:16 ] >>60 よくやる。 しかしその初期化はリンク・ロード時に行われるものなので、 実行行でa(1:100)=0と書くのとは微妙に違う。 62 名前：59 mailto:sage [2006/11/20(月) 15:25:10 ] >>61 ﾚｽありがとう 63 名前：デフォルトの名無しさん mailto:sage [2006/11/20(月) 17:30:12 ] integer a(1:100)=0 ↑これって通る？ integer::　じゃないと通らなかった気が。 64 名前：質問1/3 [2006/11/20(月) 18:45:56 ] 拡散方程式の境界値問題をオイラー陽解法で解くための 参考プログラムというものが教科書にあり、まわしたのですが errorは出ないものの上手くまわりませんでした。 アドバイスを願います。 ちなみにMX=41,KM=30,DT=0.01としています。 PARAMETER(NX=51) DIMENSION U(NX),UU(NX) C**** INPUT & CALCULATE PARAMETERS WRITE(*,*) 'MX (<52: MESH POINTS) ? KM (NUMBER OF TIME STEP) ?' READ(*,*) MX,KM WRITE(*,*) 'DELTA T: DT ? ' READ(*,*) DT DX=1./FLOAT(MX-1) R=DT/DX**2 IH=(MX+1)/2 I5=.05/DT IF(I5.EQ.0) I5=1 CALL OUTPUT(U,NX,MX,1) C**** INITIAL CONDITION DO 10 I=1,MX X=FLOAT(I-1)/FLOAT(MX-1) IF(I.LE.IH) THEN U(I)=X ELSE U(I)=1.-X END IF 10 CONTINUE 65 名前：質問2/3 [2006/11/20(月) 18:46:56 ] C**** MAIN LOOP DO 20 K=1,KM U(1)=0. U(MX)=0. IF(MOD(K,I5).EQ.1) CALL OUTPUT(U,NX,MX,2) C DO 30 I=2,MX-1 UU(I)=R*U(I-1)+(1.-2*R)*U(I)+R*U(I+1) 30 CONTINUE DO 40 I=2,MX-1 U(I)=UU(I) 40 CONTINUE C IF(ABS(U(IH)).GE.10000.) THEN WRITE(*,*) 'DIVERGE!' STOP END IF C 20 CONTINUE C CALL OUTPUT(U,NX,MX,3) STOP END C SUBROUTINE OUTPUT(U,NX,MX,MM) CHARACTER Z(60,23) DIMENSION U(NX) 66 名前：質問3/3 [2006/11/20(月) 18:48:34 ] IF(MM.EQ.1) THEN DO 10 J=1,23 DO 10 I=2,MX-1 Z(I,J)=' ' 10 CONTINUE DO 20 J=1,23 Z(1,J)=':' Z(MX,J)=':' 20 CONTINUE DO 30 I=1,MX Z(I,1)='-' Z(I,23)='-' 30 CONTINUE END IF IF(MM.EQ.2) THEN DO 40 I=1,MX J=U(I)*40+.001 IF(J.GT.23) GO TO 40 Z(I,J)='*' 40 CONTINUE END IF IF(MM.EQ.3) THEN DO 50 J=23,1,-1 WRITE(*,600) (Z(I,J),I=1,60) 50 CONTINUE 600 FORMAT(1H ,60A1) END IF RETURN END 67 名前：59 mailto:sage [2006/11/20(月) 21:40:21 ] >>63 そうでした integer :: a(1:100)=0 じゃないと通らんです． 68 名前：デフォルトの名無しさん [2006/11/22(水) 00:03:15 ] くだすれFORTRAN（超初心者用）その2 pc8.2ch.net/test/read.cgi/tech/1164121236/ 69 名前：デフォルトの名無しさん mailto:sage [2006/11/22(水) 01:57:44 ] 今度、計算機を3台ばかし構築せねばならないんだが、 CPUはcore2duoとXeonのどちらが良いのだろうか コンパイラはインテルの非商用フリーのやつを使おうと思っている 70 名前：デフォルトの名無しさん [2006/11/22(水) 03:56:04 ] Xeon5000番台と比較するならXeon、 Xeon3000番台と比較するならCore2Duoかな 71 名前：デフォルトの名無しさん [2006/11/22(水) 11:20:02 ] キャッシュに載るならそれでもいいが、メモリーアクセス多いならItanium2で池 72 名前：デフォルトの名無しさん mailto:sage [2006/11/22(水) 14:12:45 ] お前らが使っている計算機のスペックきぼんぬ うちは CPU：Pentium4 3.2GHz メモリ：1GB コンパイラ：Intel Fortran コンパイラ メモリを2GBにしろと進言したが無理だった。 73 名前：デフォルトの名無しさん mailto:sage [2006/11/22(水) 14:54:28 ] 便乗しようじゃあないか． CPU: Pentium4 3.4GHz RAM: 1.5GB (DDRPC3200ﾃﾞｭｱﾙ) ｺﾝﾊﾟｲﾗ: Intel Fortran Compiler 9.1 OS: Windows XP pro SP2 友人のC2D (T7200) 搭載機使って計算したらノートPCのくせに俺の ﾃﾞｽｸより2倍以上速くて鬱ｒｙ 74 名前：デフォルトの名無しさん [2006/11/22(水) 15:35:21 ] これはなんですか？？？教えてください！！ NO ERRORS,4 WARNINGS [<TDOFS>FTN77 Ver 4.03] NO ERRORS [<FORVIB>FTN77 Ver 4.03] NO ERRORS [<COEF0>FTN77 Ver 4.03] NO ERRORS [<COEF1>FTN77 Ver 4.03] NO ERRORS [<SOLU>FTN77 Ver 4.03] NO ERRORS [<F_TAN>FTN77 Ver 4.03] NO ERRORS [<Z_GK>FTN77 Ver 4.03] NO ERRORS [<INTARG>FTN77 Ver 4.03] NO ERRORS [<INTARG3>FTN77 Ver 4.03] NO ERRORS [<INTARG4>FTN77 Ver 4.03] NO ERRORS [<INTARG5>FTN77 Ver 4.03] WARNING the following symbols are missing: DLAX D:\MyDo\Fortran\1122test.obj (D:\MYDO\FORTRAN\1122TEST.FOR) DLAXR D:\MyDo\Fortran\1122test.obj (D:\MYDO\FORTRAN\1122TEST.FOR) 75 名前：デフォルトの名無しさん mailto:sage [2006/11/22(水) 15:50:10 ] >>74 エラーはない（おめでと♪）けど、変なのが４つあったよ。 なんか、1122testっていうファイルのDLAX, DLAXRっていう名前が変な感じ。 ちょっと見た方がいいかも・・・（ｂｙふぉとらんこんばいら4.03号）。 76 名前：デフォルトの名無しさん [2006/11/22(水) 16:07:08 ] ありがとうございます☆ 今昔の先輩のプログラミングを解析中なだけで、FORTRAN自体を理解してませんが、恥を忍んでお伺いします。 DLAXRというのはたぶんこの部分ですが、どこか変でしょうか？ c **** DLAXR *************************************** call DLAXR(XX,A,k,hmn,ALU,b,IPP,VWW,icond) do 110 j=1,mn+1,1 jj=j-1 X(jj)=XX(j) 110 continue do 111 j=mn+2,hmn,1 jj=j-mn-1 Y(jj)=XX(j) 111 continue c *********************************************************** return end 77 名前：デフォルトの名無しさん mailto:sage [2006/11/22(水) 16:51:21 ] >>76 わたしもよくわかんないけど、 DLAXRを呼んでるけど、このファイルではまだ定義されてないから 型・（引数の数とかの）形式があってるか不安・・・ よければファイルに書いといてね たぶん他のファイルで定義されてるんだろうから とりあえず適当にくっつけるけど、間違っても怒らないでね（ｂｙふぉとらんこんぱいら） 78 名前：デフォルトの名無しさん mailto:sage [2006/11/22(水) 17:24:56 ] >>74 DLAX, DLAXR って名前から察すると、富士通の数値計算ライブラリSSL II を使ってるっぽい。 79 名前：デフォルトの名無しさん mailto:sage [2006/11/22(水) 17:45:18 ] すれ違いかも知れないが、CoreDuo、Core2DuoでLinuxで動くの？ 80 名前：デフォルトの名無しさん mailto:sage [2006/11/22(水) 18:00:54 ] >Core2DuoでLinuxで シンタックス・エラー： 解決されない助詞「で」 81 名前：デフォルトの名無しさん mailto:sage [2006/11/22(水) 20:45:13 ] >>80 nice 82 名前：デフォルトの名無しさん [2006/11/22(水) 23:35:23 ] >>74 それは>>78にあるように、SSLIIのライブラリを読んでいるので、それが無いと言っている。 SSLIIのライブラリがあるコンパイラ（Fujitsu、Lahey）の奴でやるか、 あるいは他のライブラリで置き換える必要がある。 LU分解っぽいので、LAPACKとか何でもあるべ。 83 名前：デフォルトの名無しさん [2006/11/23(木) 14:05:27 ] n行m列の行列があり、そこからn行m列のそれぞれの値を一列に直して、それぞれに値について大きい順に並べなおし、元のｎ行ｍ列に直してみたいのですけどどんなプログラムの方法がありますか？ よければ、教えてください＞＜ ちなみにfortran77を使っています＞＜/ 84 名前：デフォルトの名無しさん [2006/11/23(木) 16:00:07 ] >>83 言ってることがよく分からないが、２次元行列を列ごとにソートして、 上から下へ大きい順に並べなおしたいということか？ それなら簡単だ。 １．まず１次元ソートのサブルーチンを用意する。これを SORT(n, x) としよう。 ２．列ごとにソートする。 DO １０ i = 1, m 　CALL SORT(N, X(1, i)) 10 CONTINUE これでおｋ。 85 名前：デフォルトの名無しさん [2006/11/23(木) 18:03:09 ] ＞＞84さん指摘ありがとうございます でもちょっと違うのです＞＜　 たとえば、1　5　 9　13 　　　　　2　6　10　14 　　　　　3　7　11　15 　　　　　4　8　12　16 について、 1 2 3 … 16 といってから順番を変更して 16　12 8 4 15　11 7 3 14　10 6 2 13　 9　5 1 という風に並べたいのです 86 名前：デフォルトの名無しさん [2006/11/23(木) 19:56:32 ] >>85 それならもっと簡単だ。FORTRAN７７は配列の先頭アドレスしか渡さない。 次元数とか大きさは気にしない。 また、FORTRANの配列はメモリー上では、列方向に隣り合って１次元に並んでいる。 ゆえに、N*Mの要素数で１次元のソートをすればいい。 １．１次元のソートを用意する。 ２．CALL　SORT（N*M、X） これでOK 87 名前：デフォルトの名無しさん mailto:sage [2006/11/24(金) 11:12:38 ] >>83 元々mixiに書き込んだんだから、解決したらmixiにも書くように。 88 名前：デフォルトの名無しさん [2006/11/24(金) 11:44:03 ] >>87　 つーか>>83のレベルでは>>86の回答は理解できないんじゃないか。 FORTRANはこういうイカサマ技によって、１次元整数０クリアをつっかって 多次元整数、実数などのクリアーに使えた。一種のPolymorphism　ｗｗｗｗｗ 89 名前：デフォルトの名無しさん [2006/11/24(金) 14:36:03 ] ＞86さんありがとうございます 確かに速く知りたいという気がはやりすぎて、いろんな所で質問していますね ちゃんと自分で理解できたら書いていきたいとは思っているのですけどね でも、やはりちょっと理解できません＾＾； よければ、何か例とかあるとありがたいのですけど＞＜ 90 名前：デフォルトの名無しさん [2006/11/24(金) 14:36:17 ] >>77　 ありがとうございます！ >>78さんの言うとおり、富士通SSL�Uを使っている様子です。 >>82 LaheyやLAPACKとはどんなものですか？？ 91 名前：デフォルトの名無しさん mailto:sage [2006/11/24(金) 20:57:18 ] >>90 ttp://www.google.co.jp/search?q=lahey+fortran ttp://www.google.co.jp/search?q=lapack+fortran 92 名前：デフォルトの名無しさん [2006/11/24(金) 23:12:56 ] >>90 Laheyはコンパイラの自社独自開発をやめて、今では富士通のフロントエンドを使っている。 富士通はLinux版しか出していないが、LaheyはWindows版も出している。 富士通が絡んでいるので、プロフェッショナルバージョンを買うとSSLがついてくる。 あとSSLは名古屋大学のNUPACKと密接な関係があるので、netNUMPACが参考に なるかもしれないし、ならないかもしれないｗｗｗ 漏れは富士通系はメインには使わなかったのでよく知らんｗ LAPACKは固有地問題と線形問題用のサブルーチン集だ。これは国際標準と化してきているので どこの機械でも動いている。 93 名前：デフォルトの名無しさん mailto:sage [2006/11/24(金) 23:33:34 ] >>89 PROGRAM unko REAL a(5, 2) CALL kintama(5 * 2, a) DO 10 i = 1, 5 PRINT *, a(i, 1), a(i, 2) 10 CONTINUE END C SUBROUTINE kintama(n, x) REAL x(n) DO 10 i = 1, n x(i) = REAL(i) 10 CONTINUE RETURN END 実行結果 1.000000 6.000000 2.000000 7.000000 3.000000 8.000000 4.000000 9.000000 5.000000 10.00000 Press any key to continue サブルーチン　金玉　で、1次元配列ｘに１〜ｎまで順番に数を代入している。 一方、メインルーチンから呼び出すときは、２次元配列を引数（ひきすう）に与えている。 要素数は行×列＝１０個で渡している。　返ってくる結果を見ると、列方向に連続した値で戻って来ている。 同様のことをSORTでやれば、はじめ二次元にランダムだった数が、 列方向にソートされて戻ってくる。 94 名前：デフォルトの名無しさん mailto:sage [2006/11/25(土) 08:07:15 ] >>89 ＞確かに速く知りたいという気がはやりすぎて、いろんな所で質問していますね それはマルチポストと言って、マナー違反だな。 ここをよーく読んどくように。 ja.wikipedia.org/wiki/%E3%83%9E%E3%83%AB%E3%83%81%E3%83%9D%E3%82%B9%E3%83%88 95 名前：デフォルトの名無しさん mailto:sage [2006/11/26(日) 17:18:04 ] ＦＯＲＴＲＡＮのエロゲとかないの？ 96 名前：デフォルトの名無しさん [2006/11/27(月) 19:00:46 ] >>53 ＞x(1:n) = y(1:n) この書き方って例えば次のような漸化式に使うと ちゃんと計算してくれない． x(1)=0 x(2)=1 x(3:n) = 2*x(2:n-1) - 4*x(1:n-2) 自分に自分を代入するような 式だと駄目っぽい (途中から全部0になる)．なんで？？ 97 名前：デフォルトの名無しさん mailto:sage [2006/11/27(月) 21:01:20 ] >>96 配列の演算はすべて並列計算として処理されるので x の要素が 0 で初期化されていた場合、 x(3) = 2*x(2) - 4*x(1) = 2*1 - 4*0 = 2 x(4) = 2*x(3) - 4*x(2) = 2*0 - 4*1 = -4 x(5) = 2*x(4) - 4*x(3) = 2*0 - 4*0 = 0 x(6) = 2*x(5) - 4*x(4) = 2*0 - 4*0 = 0 ・・・ となる。 do i=2,n x(i) = x(i-1) + 1 end do と x(2:n) = x(1:n-1) + 1 では意味が違うのだ。 98 名前：96 mailto:sage [2006/11/27(月) 22:18:19 ] >>97 おお〜よくわかりました！ ありがとう 99 名前：デフォルトの名無しさん [2006/11/28(火) 00:17:59 ] >>96 REAL::A(N, N), B(N, N) Ａ　＝　Ａ　＊　Ｂ のような状況を考えれば、今のＦＯＲＴＲＡＮの仕様がもっともだと分かる。 行列操作がらみで、この仕様はとても便利。いちいち自分でテンポラリにコピーしなくていい。 あと一応、並列化されたときのことも考慮されてそういう仕様になっているらしい。 100 名前：デフォルトの名無しさん [2006/11/29(水) 13:24:49 ] >>74のプログラムにおいて、ファイルの定義とはどの部分にどのように書き込むものなのですか？ 全体のプログラムの中のどこかに隠れているのでしょうか？

---

---

[ 続きを読む ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef