トップページ⇒tech
991コメント434KB

FORTRAN Ⅳ

■ このスレッドは過去ログ倉庫に格納されています
0001デフォルトの名無しさん2006/11/12(日) 17:13:35
FORTRAN総合スレッドです。
FORTRAN77/90/95/2003, HPF等、FORTRAN全般に関する話題を書き込んで下さい。

・FORTRAN関連情報
お約束だが、まずは自分で調べよう
http://www.google.co.jp/

x86で動作するFortranコンパイラ
http://www.nminoru.jp/~nminoru/programming/x86-fortran.html

The Fortran Company
http://www.fortran.com/
Polyhedron Software
http://www.polyhedron.co.uk/
Google Directory - Fortran
http://directory.google.com/Top/Computers/Programming/Languages/Fortran/

・前スレ
FORTRAN III
http://pc8.2ch.net/test/read.cgi/tech/1104724162/
FORTRAN II
http://pc5.2ch.net/test/read.cgi/tech/1068351911/
Fortran
http://pc2.2ch.net/test/read.cgi/tech/1003214738/
0002デフォルトの名無しさん2006/11/12(日) 17:14:45
・フリー／無料評価版のコンパイラ
The G95 project （Fortran95）
http://www.g95.org/
GCC (g77) （今のところ FORTRAN77）
http://gcc.gnu.org/
Salford FTN77/95 Personal Edition （Fortran77/95, Windows, 要登録）
http://www.salfordsoftware.co.uk/software/downloads/compilers.html
FTN77/95のインストール方法
http://www.coastal-env.k.u-tokyo.ac.jp/koibuchi/fortran.htm
the F Programming Language （Fortran90のサブセット, Win/Linux/Solaris）
http://www.fortran.com/F/compilers.html
Open Watcom Fortran Compiler （FORTRAN77）
http://www.openwatcom.org/
Intel Fortran Compiler （Fortran95, 無料評価版(Win)／非商用バージョン(Linux), 要登録）
http://www.xlsoft.com/jp/products/download/download.html
Lahey/Fujitsu Fortran （期限付きトライアル版は終了した模様）
http://www.lahey.com/downloads.htm

・ドキュメント
Compaq Visual Fortran オンライン・ドキュメント
http://www.xlsoft.com/jp/products/intel/cvf/docs/index.html
Intel Fortran Compiler 日本語版マニュアルのダウンロード （要登録）
http://www.xlsoft.com/jp/products/intel/download.html
富士通Fortranオンラインマニュアル
http://www.cc.nao.ac.jp/fsunman/japanese/Fortran/index.html

FORTRAN77 standard
http://www.swcp.com/~walt/F77_std/f77_std.html
http://www.fortran.com/fortran/F77_std/rjcnf0001.html

JISデータベース
http://www.jisc.go.jp/
0003デフォルトの名無しさん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　|/
　　　　　　　　　　　 丨　　 ヽ ｀ ', | Ｏ｀`''────…'''"´ 　Ｏ',ﾞ:､ 　l　|／
0004デフォルトの名無しさん2006/11/12(日) 18:34:32
機種依存文字使うなボケ失せろ
0005デフォルトの名無しさん2006/11/12(日) 19:07:04
>>4
きたねー言葉使うな糞野郎！

死ね！
0006デフォルトの名無しさん2006/11/12(日) 19:30:21
数値積分が使えるフリーの数値計算ライブラリありませんか？
0007デフォルトの名無しさん2006/11/12(日) 19:58:47
>>6
NetLibへ池。６６時代のものがごまんとあるぜ。
0008デフォルトの名無しさん2006/11/12(日) 20:01:14
IMSLがあるんじゃね？
0009デフォルトの名無しさん2006/11/12(日) 20:02:08
>>8
IMSLは有料パッケージソフトだぞ。
0010デフォルトの名無しさん2006/11/12(日) 20:08:40
IMSLは10万以上する予感

LAPACKに積分はふくまれておらんの？
0011デフォルトの名無しさん2006/11/12(日) 20:08:54
コンパイラより高い・・・
0012デフォルトの名無しさん2006/11/12(日) 20:10:38
このスレのエキスパートたちは有料パッケージとか使っておられるのですか？
0013デフォルトの名無しさん2006/11/12(日) 20:13:17
>>7
Netlib行ったが数が多すぎて何がなんだかわからない。
これをどう使えばいいかすらわからない。
見た感じプログラム単体が落っこちているだけにしか見えないが
どうやってってこれを呼び出して使うんだろうか。
0014デフォルトの名無しさん2006/11/12(日) 20:32:24
>>10
LAPACKは元々、LINPACK＋EISPACKで、それぞれ密行列の線形方程式と固有値問題の
ルーチン集なので積分は主題から外れていると思う。

LINPACK、EISPACKは６０年代の米英の国家プロジェクトだったように記憶している（記憶曖昧ｗ）。
LAPACKは基本的にそれを７７化してBLASの上に統一化し、一部アルゴリズムを改良・拡充した
ものだろう。その辺の経緯や順序はよくわかんね。


>>13
NetLibとはそういうところだｗ　埃の積もった土蔵見たいなもんだ。
結局、膨大すぎて何が欲しいか分かってないと使えないｗ

NumericalRecipesのページへ行ってみるのが吉かもしれない。
とりあえず専門家的にはベストで無いにしろ、素人が馬鹿な外れを引かされることは無いはず。
ただ本版のソースはしばしばBug入りなのでめくらめっぽう使うのはﾔｳﾞｧｽ
有料ソースは直っていたりするｗ
0015デフォルトの名無しさん2006/11/12(日) 21:19:12
初心者は富士通のサブルーチン付きコンパイラ買うのがbestじゃね？
難易度も予算も
0016デフォルトの名無しさん2006/11/12(日) 21:25:25
Intel Math Kernel Libraryでok
0017デフォルトの名無しさん2006/11/12(日) 21:36:09
MKLって積分ルーチン入ってるの？
0018デフォルトの名無しさん2006/11/12(日) 21:39:47
多分ない

というか、フリーで有名どころのパッケージには積分ルーチンは無い予感
0019デフォルトの名無しさん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
0020デフォルトの名無しさん2006/11/13(月) 20:30:25
メインプログラムにLorentzの型定義がないからinteger扱いされてんじゃねーの？
0021192006/11/13(月) 21:30:09
implicit noneで解決しました。
１時間も時間を食いました。
0022デフォルトの名無しさん2006/11/14(火) 00:43:00
自由書式で書くなら関数・サブルーチンは内部手続きにするか、MODULEにするか
INTERFACE宣言してから使おう。
0023デフォルトの名無しさん2006/11/14(火) 01:07:20
>>20,22
が正しい。

>>19はあまりに中途半端。きっちりと厳格なF90式で書くか、F77式に暗黙の型を守るべき。
0024デフォルトの名無しさん2006/11/14(火) 02:41:24
INTEL　FORTRAN用のBLAS/LAPACKだが、AMDのcore-math-libraryにWindows用もうｐされている。
昔はLinux用の一部ルーチンだけだったが、いまや全ルーチンを網羅したようだ。

漏れはAMDなので、誰かIntelChipでも使えるか試してくれ。

http://developer.amd.com/acml.jsp
0025デフォルトの名無しさん2006/11/14(火) 15:36:00
fortran中で他のexeを実行できるようにするにはどうしたら良いですか？
fortranで作られたexeなんだけどソースくっつけるのがめんどくさいんで・・・ｗ
0026デフォルトの名無しさん2006/11/14(火) 16:32:26
>>23
FORTRAN(90)の書籍で、定番の物ってありますか？
0027デフォルトの名無しさん2006/11/14(火) 16:53:09
これからフォートランを勉強するのですが、
とても解りやすい入門書を教えていだたけませんか？
プログラミング自体、あまり経験がありません。
Ｃを少しやりました。
よろしくお願いします。
0028デフォルトの名無しさん2006/11/14(火) 18:18:30
ｘの値を、1,10,100,1000,10000,100000,1000000の順に変えて
その常用対数、自然対数の値を計算して表示するプログラムを作れ

という課題が出たのですがサッパリ分かりません。
どうかお力添えください。
0029デフォルトの名無しさん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
0030デフォルトの名無しさん2006/11/14(火) 23:23:10
倍精度実数と4倍精度実数では実行速度は変わってしまいますか？
0031デフォルトの名無しさん2006/11/14(火) 23:28:24
>>30
死ぬほど違う。　
倍精度までは数値演算プロセッサの助けでハードウェア演算するが、
４倍精度はソフトウェアで浮動小数点演算をするので遅くなる。

６４BitCPUなら若干状況は改善するが本質においては変わらない。
DECのアルファーCPUは４倍精度も微妙に早かったが、特殊な命令セットがあったためらしく、
６４bitCPUだからといってよくなるわけでもないようだ。
0032デフォルトの名無しさん2006/11/14(火) 23:31:25
>>29
DO文で書かないといけないのですが
どうすればいいんでしょうか・・
0033302006/11/15(水) 00:03:03
>>31
詳しいﾚｽｻﾝｸｽ！
0034デフォルトの名無しさん2006/11/15(水) 00:35:58
>>26
残念ながらあまりいい本は無い。

今、手に入りやすいのはこれ。
http://www.amazon.co.jp/gp/product/4563014095
しかし、Ｆ７７の癖を引きずっていて、あまりいい書き方をしていない。
著者は京大のおっさんらしく、ＷＥＢで同等の内容のものを見られる。
アドレスはアマゾンのレビューに書いてある。

英語でよければ、この辺とか、
http://www.lahey.com/elfpage.htm
この辺のＥｓｓｅｎｔｉａｌ Ｆｏｒｔｒａｎあたりが学習用サブセットになっていて、純粋Ｆ９０を学べる。
http://www.fortran.com/F/books.html

しかし、現実世界ではＦ７７ソースがあふれているので微妙ではある。
実際上記のものは教育用としてもはやらなかった。
0035デフォルトの名無しさん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
0036デフォルトの名無しさん2006/11/15(水) 14:47:06
　　complex(8) c

において，

　　c=(0.d0, 0.d0)

として初期化するのも，

　　c=0.d0

として初期化するのも全く同じことなのですか？
0037デフォルトの名無しさん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

何処がおかしいんでしょうか?
0038デフォルトの名無しさん2006/11/15(水) 18:31:55
>>37
read (*,*) x,log(x),log10(x)
どこに入れたいの？
たぶん、その穴じゃないわよ
0039デフォルトの名無しさん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,)
ここら辺が間違ってる。
どこが間違いなのかわからないときは、とりあえずコンパイラ様にソースを渡して
文句を言われたところを上から順に直していくのも１つの手だ。


だれかくだすれの次スレ立ててくれ。
0040デフォルトの名無しさん2006/11/15(水) 23:55:01
>>37
>>35をもういっぺん見ろ！

>>39
重複スレをクダスレにするか？
0041デフォルトの名無しさん2006/11/17(金) 00:48:27
>>40
タイトル一緒だから混乱しそうｗ
0042372006/11/17(金) 10:16:25
>>35
WRITE (6,*) が分かりません
>>37の問題の書式の指定はこれでいいんでしょうか
0043デフォルトの名無しさん2006/11/17(金) 11:23:12
WRITE(6,　'(I10, F15.7, F15.7)') x, y, z

もう少し教科書嫁！！

こんな感じでｗ
0044デフォルトの名無しさん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
0045デフォルトの名無しさん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
こうですか?
0046デフォルトの名無しさん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


おかげさまでここまで出来たんですが、どうしても*******が出てしまいます＞＜
0047デフォルトの名無しさん2006/11/17(金) 19:09:14
>>46
「変数の型」を再勉強してから出直してこい。
0048デフォルトの名無しさん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を整数表示に指定することはできない．
0049デフォルトの名無しさん2006/11/18(土) 00:44:07
>>48
>>37の問題を読む感じでは、
write(*, '(i7, f15.7 f15.7)') int(x), log10(x), log(x)
を予期していると思われる。

>>46はもう少し貰った回答例をよく吟味しろｗ
0050デフォルトの名無しさん2006/11/18(土) 01:00:15
ＰｕｋｉＷｉｋｉ　つくってみた
http://www.moecity.jp/home/fortran/

ｗｉｋｉよく知らない。自由に使ってみてくれ。
0051デフォルトの名無しさん2006/11/18(土) 19:42:09
f90スタイルで学べるいい文書はありませんかね。
英語か日本語でお願い
0052デフォルトの名無しさん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
0053デフォルトの名無しさん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では違っていたのか？
0054デフォルトの名無しさん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では（当然）
0055542006/11/18(土) 22:12:21
× g95では（当然）
○ g95では（当然） A は配列全体と解釈された。
0056デフォルトの名無しさん2006/11/18(土) 23:02:07
>>54
そうじゃなくて
REAL A, A(10)
が可能だったと思う。
0057デフォルトの名無しさん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)]
0058デフォルトの名無しさん2006/11/19(日) 01:39:12
たしか変数名と配列名は規格上の定義では別物だった記憶がある。

実際のコンパイラでどうなっていたかは知らん。



0059デフォルトの名無しさん2006/11/20(月) 08:30:40
Fortran90なんですが，

integer a(1:100)=0

みたな初期化の仕方はよくやるんですか？
0060デフォルトの名無しさん2006/11/20(月) 08:32:30
×みたな
○みたいな

orz
0061デフォルトの名無しさん2006/11/20(月) 09:01:16
>>60
よくやる。
しかしその初期化はリンク・ロード時に行われるものなので、
実行行でa(1:100)=0と書くのとは微妙に違う。

0062592006/11/20(月) 15:25:10
>>61
ﾚｽありがとう
0063デフォルトの名無しさん2006/11/20(月) 17:30:12
integer a(1:100)=0
↑これって通る？

integer::　じゃないと通らなかった気が。
0064質問1/32006/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
0065質問2/32006/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)
0066質問3/32006/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

0067592006/11/20(月) 21:40:21
>>63
そうでした

integer :: a(1:100)=0

じゃないと通らんです．
0068デフォルトの名無しさん2006/11/22(水) 00:03:15

くだすれFORTRAN（超初心者用）その2
http://pc8.2ch.net/test/read.cgi/tech/1164121236/
0069デフォルトの名無しさん2006/11/22(水) 01:57:44
今度、計算機を3台ばかし構築せねばならないんだが、
CPUはcore2duoとXeonのどちらが良いのだろうか
コンパイラはインテルの非商用フリーのやつを使おうと思っている
0070デフォルトの名無しさん2006/11/22(水) 03:56:04
Xeon5000番台と比較するならXeon、
Xeon3000番台と比較するならCore2Duoかな
0071デフォルトの名無しさん2006/11/22(水) 11:20:02
キャッシュに載るならそれでもいいが、メモリーアクセス多いならItanium2で池



0072デフォルトの名無しさん2006/11/22(水) 14:12:45
お前らが使っている計算機のスペックきぼんぬ

うちは
CPU：Pentium4 3.2GHz
メモリ：1GB
コンパイラ：Intel Fortran コンパイラ

メモリを2GBにしろと進言したが無理だった。
0073デフォルトの名無しさん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倍以上速くて鬱ｒｙ
0074デフォルトの名無しさん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)
0075デフォルトの名無しさん2006/11/22(水) 15:50:10
>>74
エラーはない（おめでと♪）けど、変なのが４つあったよ。
なんか、1122testっていうファイルのDLAX, DLAXRっていう名前が変な感じ。
ちょっと見た方がいいかも・・・（ｂｙふぉとらんこんばいら4.03号）。
0076デフォルトの名無しさん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
0077デフォルトの名無しさん2006/11/22(水) 16:51:21
>>76
わたしもよくわかんないけど、
DLAXRを呼んでるけど、このファイルではまだ定義されてないから
型・（引数の数とかの）形式があってるか不安・・・
よければファイルに書いといてね

たぶん他のファイルで定義されてるんだろうから
とりあえず適当にくっつけるけど、間違っても怒らないでね（ｂｙふぉとらんこんぱいら）
0078デフォルトの名無しさん2006/11/22(水) 17:24:56
>>74
DLAX, DLAXR って名前から察すると、富士通の数値計算ライブラリSSL II を使ってるっぽい。
0079デフォルトの名無しさん2006/11/22(水) 17:45:18
すれ違いかも知れないが、CoreDuo、Core2DuoでLinuxで動くの？
0080デフォルトの名無しさん2006/11/22(水) 18:00:54
>Core2DuoでLinuxで
シンタックス・エラー： 解決されない助詞「で」
0081デフォルトの名無しさん2006/11/22(水) 20:45:13
>>80
nice
0082デフォルトの名無しさん2006/11/22(水) 23:35:23
>>74
それは>>78にあるように、SSLIIのライブラリを読んでいるので、それが無いと言っている。

SSLIIのライブラリがあるコンパイラ（Fujitsu、Lahey）の奴でやるか、
あるいは他のライブラリで置き換える必要がある。

LU分解っぽいので、LAPACKとか何でもあるべ。

■ このスレッドは過去ログ倉庫に格納されています