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/
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/
（ （￣ ￣） ／....:::;:;:;:;:;:;:;:;:;:;:;:;:;:;:;:;:;.....ヽ､／￣￣￣￣＼／
￣ ￣ 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 |／
きたねー言葉使うな糞野郎！
！
NetLibへ池。６６時代のものがごまんとあるぜ。
IMSLは有料パッケージソフトだぞ。
LAPACKに積分はふくまれておらんの？
Netlib行ったが数が多すぎて何がなんだかわからない。
これをどう使えばいいかすらわからない。
見た感じプログラム単体が落っこちているだけにしか見えないが
どうやってってこれを呼び出して使うんだろうか。
LAPACKは元々、LINPACK＋EISPACKで、それぞれ密行列の線形方程式と固有値問題の
ルーチン集なので積分は主題から外れていると思う。
LINPACK、EISPACKは６０年代の米英の国家プロジェクトだったように記憶している（記憶曖昧ｗ）。
LAPACKは基本的にそれを７７化してBLASの上に統一化し、一部アルゴリズムを改良・拡充した
ものだろう。その辺の経緯や順序はよくわかんね。
>>13
NetLibとはそういうところだｗ 埃の積もった土蔵見たいなもんだ。
結局、膨大すぎて何が欲しいか分かってないと使えないｗ
NumericalRecipesのページへ行ってみるのが吉かもしれない。
とりあえず専門家的にはベストで無いにしろ、素人が馬鹿な外れを引かされることは無いはず。
ただ本版のソースはしばしばBug入りなのでめくらめっぽう使うのはﾔｳﾞｧｽ
有料ソースは直っていたりするｗ
難易度も予算も
というか、フリーで有名どころのパッケージには積分ルーチンは無い予感
結果が
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
１時間も時間を食いました。
INTERFACE宣言してから使おう。
が正しい。
>>19はあまりに中途半端。きっちりと厳格なF90式で書くか、F77式に暗黙の型を守るべき。
昔はLinux用の一部ルーチンだけだったが、いまや全ルーチンを網羅したようだ。
漏れはAMDなので、誰かIntelChipでも使えるか試してくれ。
http://developer.amd.com/acml.jsp
fortranで作られたexeなんだけどソースくっつけるのがめんどくさいんで・・・ｗ
FORTRAN(90)の書籍で、定番の物ってありますか？
とても解りやすい入門書を教えていだたけませんか？
プログラミング自体、あまり経験がありません。
Ｃを少しやりました。
よろしくお願いします。
その常用対数、自然対数の値を計算して表示するプログラムを作れ
という課題が出たのですがサッパリ分かりません。
どうかお力添えください。
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
死ぬほど違う。
倍精度までは数値演算プロセッサの助けでハードウェア演算するが、
４倍精度はソフトウェアで浮動小数点演算をするので遅くなる。
６４BitCPUなら若干状況は改善するが本質においては変わらない。
DECのアルファーCPUは４倍精度も微妙に早かったが、特殊な命令セットがあったためらしく、
６４bitCPUだからといってよくなるわけでもないようだ。
DO文で書かないといけないのですが
どうすればいいんでしょうか・・
詳しいﾚｽｻﾝｸｽ！
残念ながらあまりいい本は無い。
今、手に入りやすいのはこれ。
http://www.amazon.co.jp/gp/product/4563014095
しかし、Ｆ７７の癖を引きずっていて、あまりいい書き方をしていない。
著者は京大のおっさんらしく、ＷＥＢで同等の内容のものを見られる。
アドレスはアマゾンのレビューに書いてある。
英語でよければ、この辺とか、
http://www.lahey.com/elfpage.htm
この辺のＥｓｓｅｎｔｉａｌ Ｆｏｒｔｒａｎあたりが学習用サブセットになっていて、純粋Ｆ９０を学べる。
http://www.fortran.com/F/books.html
しかし、現実世界ではＦ７７ソースがあふれているので微妙ではある。
実際上記のものは教育用としてもはやらなかった。
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
において，
c=(0.d0, 0.d0)
として初期化するのも，
c=0.d0
として初期化するのも全く同じことなのですか？
その常用対数および自然対数の値を 計算して表示するプログラムを作れ．
なお，計算結果は書式の指定をして，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
何処がおかしいんでしょうか?
read (*,*) x,log(x),log10(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,)
ここら辺が間違ってる。
どこが間違いなのかわからないときは、とりあえずコンパイラ様にソースを渡して
文句を言われたところを上から順に直していくのも１つの手だ。
だれかくだすれの次スレ立ててくれ。
>>35をもういっぺん見ろ！
>>39
重複スレをクダスレにするか？
タイトル一緒だから混乱しそうｗ
WRITE (6,*) が分かりません
>>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
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
こうですか?
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
おかげさまでここまで出来たんですが、どうしても*******が出てしまいます＞＜
「変数の型」を再勉強してから出直してこい。
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を整数表示に指定することはできない．
>>37の問題を読む感じでは、
write(*, '(i7, f15.7 f15.7)') int(x), log10(x), log(x)
を予期していると思われる。
>>46はもう少し貰った回答例をよく吟味しろｗ
http://www.moecity.jp/home/fortran/
ｗｉｋｉよく知らない。自由に使ってみてくれ。
英語か日本語でお願い
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
＞推奨する書き方 real:: x(n), y(n)
＞x(1:n) = y(1:n)
＞
＞
＞×推奨しない書き方 x(:) = y(:) または x = y
これはかなりうそ臭いと思う。
ｘ＝ｙ の場合と違って、コロン指定子を使うと、いったん暗にコンパイラが
テンポラリ配列にコピーするのでたいてい余分な時間がかかる。
しかし漏れが疑問に思うのは、F77では変数と配列に同じ名前をAとA(10)を指定できて
別物に扱われる様に規格で決まっていたと思うのだが、F90で矛盾が起きないのか・・・
それとも、これはJIS規格だけでANSIやMILでは違っていたのか？
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では（当然）
○ g95では（当然） A は配列全体と解釈された。
そうじゃなくて
REAL A, A(10)
が可能だったと思う。
hoge.f:1:
REAL A,A(10)
1 2
Invalid declaration of or reference to symbol `a' at (2) [initially seen at (1)]
実際のコンパイラでどうなっていたかは知らん。
integer a(1:100)=0
みたな初期化の仕方はよくやるんですか？
○みたいな
orz
よくやる。
しかしその初期化はリンク・ロード時に行われるものなので、
実行行でa(1:100)=0と書くのとは微妙に違う。
ﾚｽありがとう
↑これって通る？
integer:: じゃないと通らなかった気が。
参考プログラムというものが教科書にあり、まわしたのですが
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
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)
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
そうでした
integer :: a(1:100)=0
じゃないと通らんです．
くだすれFORTRAN（超初心者用）その2
http://pc8.2ch.net/test/read.cgi/tech/1164121236/
CPUはcore2duoとXeonのどちらが良いのだろうか
コンパイラはインテルの非商用フリーのやつを使おうと思っている
Xeon3000番台と比較するならCore2Duoかな
うちは
CPU：Pentium4 3.2GHz
メモリ：1GB
コンパイラ：Intel Fortran コンパイラ
メモリを2GBにしろと進言したが無理だった。
CPU: Pentium4 3.4GHz
RAM: 1.5GB (DDRPC3200ﾃﾞｭｱﾙ)
ｺﾝﾊﾟｲﾗ: Intel Fortran Compiler 9.1
OS: Windows XP pro SP2
友人のC2D (T7200) 搭載機使って計算したらノートPCのくせに俺の
ﾃﾞｽｸより2倍以上速くて鬱ｒｙ
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)
エラーはない（おめでと♪）けど、変なのが４つあったよ。
なんか、1122testっていうファイルのDLAX, DLAXRっていう名前が変な感じ。
ちょっと見た方がいいかも・・・（ｂｙふぉとらんこんばいら4.03号）。
今昔の先輩のプログラミングを解析中なだけで、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
わたしもよくわかんないけど、
DLAXRを呼んでるけど、このファイルではまだ定義されてないから
型・（引数の数とかの）形式があってるか不安・・・
よければファイルに書いといてね
たぶん他のファイルで定義されてるんだろうから
とりあえず適当にくっつけるけど、間違っても怒らないでね（ｂｙふぉとらんこんぱいら）
DLAX, DLAXR って名前から察すると、富士通の数値計算ライブラリSSL II を使ってるっぽい。
シンタックス・エラー： 解決されない助詞「で」
nice
それは>>78にあるように、SSLIIのライブラリを読んでいるので、それが無いと言っている。
SSLIIのライブラリがあるコンパイラ（Fujitsu、Lahey）の奴でやるか、
あるいは他のライブラリで置き換える必要がある。
LU分解っぽいので、LAPACKとか何でもあるべ。
よければ、教えてください＞＜
ちなみにfortran77を使っています＞＜/
言ってることがよく分からないが、２次元行列を列ごとにソートして、
上から下へ大きい順に並べなおしたいということか？
それなら簡単だ。
１．まず１次元ソートのサブルーチンを用意する。これを SORT(n, x) としよう。
２．列ごとにソートする。
DO １０ i = 1, m
CALL SORT(N, X(1, i))
10 CONTINUE
これでおｋ。
でもちょっと違うのです＞＜
たとえば、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
という風に並べたいのです
それならもっと簡単だ。FORTRAN７７は配列の先頭アドレスしか渡さない。
次元数とか大きさは気にしない。
また、FORTRANの配列はメモリー上では、列方向に隣り合って１次元に並んでいる。
ゆえに、N*Mの要素数で１次元のソートをすればいい。
１．１次元のソートを用意する。
２．CALL SORT（N*M、X）
これでOK
元々mixiに書き込んだんだから、解決したらmixiにも書くように。
つーか>>83のレベルでは>>86の回答は理解できないんじゃないか。
FORTRANはこういうイカサマ技によって、１次元整数０クリアをつっかって
多次元整数、実数などのクリアーに使えた。一種のPolymorphism ｗｗｗｗｗ
確かに速く知りたいという気がはやりすぎて、いろんな所で質問していますね
ちゃんと自分で理解できたら書いていきたいとは思っているのですけどね
でも、やはりちょっと理解できません＾＾；
よければ、何か例とかあるとありがたいのですけど＞＜
ありがとうございます！
>>78さんの言うとおり、富士通SSLⅡを使っている様子です。
>>82
LaheyやLAPACKとはどんなものですか？？
ttp://www.google.co.jp/search?q=lahey+fortran
ttp://www.google.co.jp/search?q=lapack+fortran
Laheyはコンパイラの自社独自開発をやめて、今では富士通のフロントエンドを使っている。
富士通はLinux版しか出していないが、LaheyはWindows版も出している。
富士通が絡んでいるので、プロフェッショナルバージョンを買うとSSLがついてくる。
あとSSLは名古屋大学のNUPACKと密接な関係があるので、netNUMPACが参考に
なるかもしれないし、ならないかもしれないｗｗｗ
漏れは富士通系はメインには使わなかったのでよく知らんｗ
LAPACKは固有地問題と線形問題用のサブルーチン集だ。これは国際標準と化してきているので
どこの機械でも動いている。
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でやれば、はじめ二次元にランダムだった数が、
列方向にソートされて戻ってくる。
＞確かに速く知りたいという気がはやりすぎて、いろんな所で質問していますね
それはマルチポストと言って、マナー違反だな。
ここをよーく読んどくように。
http://ja.wikipedia.org/wiki/%E3%83%9E%E3%83%AB%E3%83%81%E3%83%9D%E3%82%B9%E3%83%88
＞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になる)．なんで？？
配列の演算はすべて並列計算として処理されるので 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
では意味が違うのだ。
おお～よくわかりました！
ありがとう
REAL::A(N, N), B(N, N)
Ａ ＝ Ａ ＊ Ｂ
のような状況を考えれば、今のＦＯＲＴＲＡＮの仕様がもっともだと分かる。
行列操作がらみで、この仕様はとても便利。いちいち自分でテンポラリにコピーしなくていい。
あと一応、並列化されたときのことも考慮されてそういう仕様になっているらしい。
TOP カテ一覧 スレ一覧 2ch元 削除依頼 ▲
Google Maps API 質問箱 (305)
Coqスレ (138)
【ツール】 フレームワーク作成スレ 【ライブラリ】 (152)
D言語 Part29 (116)
CLDC+MIDP+携帯電話用Javaスレッド part 9 (919)
プログラミング言語 Scala 7冊目 (930)
--log9.info------------------
ぶっこたんFX◆163 (200)
【USD/JPY】ドル円専用スレ Part18948【$\】 (949)
【AUD/USD】豪ドルスト専用スレpart58 (124)
【沈没】無敵艦隊 １０６港目【寸前】 (426)
愉快なコテが冷静とのあいだで××したい6 (799)
【AUD】豪ドルのレートに一喜一憂するスレpart2224 (129)
【EUR/JPY】ﾕｰﾛﾄﾞﾙはどうなる?2600【ﾕﾛﾄﾞﾙ専用】 (549)
◆◆◆神様がどんな質問にも返答するぜ8◆◆◆ (190)
【NZD】キウイを語るスレ425匹目【kiwi】 (642)
【日経CNBC】総合実況スレ part133 (589)
【(・ω・) 】スワップ派のFX part7【ハイレバTRY】 (844)
（´・ω・｀）あのー そのー3 (250)
三空ADVackアフィ_鈴木隆一■リアルトレード24GP_11 (658)
【EUR/CHF】ユロスイすーいすい♪Part11【下限1.2】 (977)
【AUD】豪ドルのレートに一喜一憂するスレpart2223 (1001)
【EUR/USD】ﾕｰﾛﾄﾞﾙはどうなる?2659【ﾕﾛﾄﾞﾙ専用】 (1001)
--log55.com------------------
勘違い女が勘違い起こした話
失恋して後悔した事を書き込むスレッド
大好きな彼女が俺を振って不倫してるけど復縁したい
【都市伝説】天才天界の正体を探る【YouTube】
何の制約もなかったら、何したい？
てす
自分が情けなくてしょうがない
何に感謝してる？