くだすれFORTRAN（超初心者用）その2

このスレッドは、他のスレッドでは書き込めない超低レベル、
もしくは質問者自身何が何だか分からない質問を勇気を持って書き込むスレッドです。
FORTRAN使いが優しくコメントを返しますが、
お礼はFORTRANの布教と初心者の救済と次期Fortran2008規格でのCOMEFROM文採用をお願いします。

●注意事項
・質問する前にGoogle等の検索サイトで検索しましょう。
・回答者にわかりやすい様に、質問内容はできる限り詳しく書きましょう。
・エラーの場合は起きた状況、環境（OS・コンパイラ）、エラーメッセージも詳しく書きましょう。


●前スレ
くだすれFORTRAN（超初心者用）
http://pc8.2ch.net/test/read.cgi/tech/1138063703/

●関連スレ
FORTRAN �W
http://pc8.2ch.net/test/read.cgi/tech/1163319215/

質問募集中age

関連スレFORTRAN �Wの>>64,65,66お願いします。

>>3
情報が足り無すぎる。
うまく回らないとは、どのようにうまく回らないのだ？

ランタイムエラーが出るのか？
エラーも出ずに実行は終了するが結果がおかしいというのか？

>4
すいません。
エラーも出ずに実行は終了するが結果がおかしいという感じですね。
ディバック機能でしょうか？使ったら3/3のDO 40〜40CONTINUEあたりが
悪いらしいのですが

>>5
正直その質問はスレの趣旨をはみ出している。
それは結果の意味まで理解していないと駄目なので、ここで聞くより周りの人に聞いたほうがいい。

f(x)＝sinx-sin2x/2+sin3x/3-sin4x/4・・・+(-1)^(n-1)sinnx/n
の値をx＝0,0.1,0.2,0.3,0.4,0.5,0.6・・・・0.9,2.0に対して計算し表示するプログラムを作成する
なお，n=8 として計算すること

program report6
implicit none
integer :: n,i
real :: s,a,x

write (*,*) ' n= '
read (*,*) n
do x=0,2.0,0.1
a=(-1)**(n-1)*sin(n*x)
s=0.0
do i=1,n
a=a/i
s=s+a
end do
end do
write (*,*) n,s
stop
end program report6
計算が上手くいきません

program report6
implicit none
integer :: i,n
real :: s,a,x

write (*,*) ' n= '
read (*,*) n

s=0.0

do x=0.0, 2.1, 0.1

do i=1,n
a=(-1)**(i-1)*sin(i*x)
s=s+a/real(i)
end do
write (*,*) x, s

end do
stop
end program report6

改良しましたが、駄目でした計算結果は

0. 0.
0.100000001 0.00579053164
0.200000003 0.0474120378
0.300000012 0.164783984
0.400000006 0.380745053
0.5 0.685271919
0.600000024 1.04163218
0.700000048 1.41077828
0.800000072 1.77801299
0.900000095 2.16398764
1.00000012 2.61202049
1.10000014 3.15931892
1.20000017 3.81087661
1.30000019 4.53395605
1.40000021 5.27820969
1.50000024 6.00946999
1.60000026 6.73514318
1.70000029 7.50308847
1.80000031 8.37220001
1.90000033 9.37193203
2.00000024 10.4772806
になりました・・・・

ちなみに正しい答は
X F
0.00000000 0.00000000
0.10000000 0.00579053
0.20000000 0.04162151
0.30000001 0.11737194
0.40000001 0.21596113
0.50000000 0.30452693
0.60000002 0.35636017
0.70000005 0.36914620
0.80000007 0.36723456
0.90000010 0.38597476
1.00000012 0.44803283
1.10000014 0.54729855
1.20000017 0.65155733
1.30000019 0.72307980
1.40000021 0.74425435
1.50000024 0.73125988
1.60000026 0.72567332
1.70000029 0.76794499
1.80000031 0.86911219
1.90000033 0.99973196
2.00000024 1.10534859

これです

自己解決しました
program report6
implicit none
integer :: i,n
real :: s,a,x

write (*,*) ' n= '
read (*,*) n

do x=0.0, 2.1, 0.1

s=0.0
do i=1,n
a=(-1)**(i-1)*sin(i*x)/real(i)
s=s+a
end do
write (*,*) x, s

end do
stop
end program report6

>>11
おめ♪

お願いします。
三次元の法線ベクトルを求めたいのですが、上手くいかないのです。
どうすればいいでしょうか？？
とりあえず、今僕が作ったのを貼り付けるので検討してください。

C HOUSEN VECTOR
DO K=1,KMX
DO J=2,JMX-1
　　　DO I=2,IMX-1

U=(Z(I+1,J+1,K)-Z(I,J,K))*(Y(I+1,J,K)-Y(I,J+1,K))-
& (Y(I+1,J+1,K)-Y(I,J,K))*(Z(I+1,J,K)-Z(I,J+1,K))

V=(X(I+1,J+1,K)-X(I,J,K))*(Z(I+1,J,K)-Z(I,J+1,K))-
& (Z(I+1,J+1,K)-Z(I,J,K))*(X(I+1,J,K)-X(I,J+1,K))

W=(Y(I+1,J+1,K)-Y(I,J,K))*(X(I+1,J,K)-X(I,J+1,K))-
& (X(I+1,J+1,K)-X(I,J,K))*(Y(I+1,J,K)-Y(I,J+1,K))

AK(I,J,K)=U/DSQRT(U**2.0D0+V**2.0D0+W**2.0D0)

BK(I,J,K)=V/DSQRT(U**2.0D0+V**2.0D0+W**2.0D0)

CK(I,J,K)=W/DSQRT(U**2.0D0+V**2.0D0+W**2.0D0)


END DO
END DO
END DO

DO K=2,KMX-1
DO J=1,JMX
DO I=2,IMX-1
　　　　U=(Z(I+1,J,K+1)-Z(I,J,K))*(Y(I+1,J,K)-Y(I,J,K+1))-
& (Y(I+1,J,K+1)-Y(I,J,K))*(Z(I+1,J,K)-Z(I,J,K+1))
　　　　V=(X(I+1,J,K+1)-X(I,J,K))*(Z(I+1,J,K)-Z(I,J,K+1))-
& (Z(I+1,J,K+1)-Z(I,J,K))*(X(I+1,J,K)-X(I,J,K+1))
　　　　W=(Y(I+1,J,K+1)-Y(I,J,K))*(X(I+1,J,K)-X(I,J,K+1))-
& (X(I+1,J,K+1)-X(I,J,K))*(Y(I+1,J,K)-Y(I,J,K+1))

AJ(I,J,K)=U/DSQRT(U**2.0D0+V**2.0D0+W**2.0D0)
　　　　　　BJ(I,J,K)=V/DSQRT(U**2.0D0+V**2.0D0+W**2.0D0)
　　　　　　CJ(I,J,K)=W/DSQRT(U**2.0D0+V**2.0D0+W**2.0D0)
　　　　END DO
END DO
END DO

DO K=2,KMX-1
DO J=2,JMX-1
DO I=1,IMX

U=(Z(I,J+1,K+1)-Z(I,J,K))*(Y(I,J+1,K)-Y(I,J,K+1))-
& (Y(I,J+1,K+1)-Y(I,J,K))*(Z(I,J+1,K)-Z(I,J,K+1))
　　　　　V=(X(I,J+1,K+1)-X(I,J,K))*(Z(I,J+1,K)-Z(I,J,K+1))-
& (Z(I,J+1,K+1)-Z(I,J,K))*(X(I,J+1,K)-X(I,J,K+1))
　　　　　W=(Y(I,J+1,K+1)-Y(I,J,K))*(X(I,J+1,K)-X(I,J,K+1))-
& (X(I,J+1,K+1)-X(I,J,K))*(Y(I,J+1,K)-Y(I,J,K+1))

AI(I,J,K)=U/DSQRT(U**2.0D0+V**2.0D0+W**2.0D0)
　　　　　　BI(I,J,K)=V/DSQRT(U**2.0D0+V**2.0D0+W**2.0D0)
　　　　　　CI(I,J,K)=W/DSQRT(U**2.0D0+V**2.0D0+W**2.0D0)

WRITE (*,*) I,J,K,DSQRT(U*U+V*V+W*W)*0.5D0

END DO
END DO
END DO

反応工学の課題を解くためにシンプソン法を使って積分する必要があるのですが
プログラミングがサッパリわからないのでプログラムが作れません！
プログラム例を貼り付けていただければそれを改変しあとは自分で何とかします
誰か私めに救いの手を差し伸べていただけないでしょうか？
よろしくお願いします！

fortran77でお願いします

上の人と同じですが、
三次元の法線ベクトルの書き方ってどんなのですか？？
正立方体だとどんな風に書けばいいですか？

fortran77でお願いします

>>16
program integral
implicit none
integer :: N = 36 ! 分割数
integer :: i
real :: x,h,a,b,sum ! x,区間幅、下限、上限、積分値
real,external :: f !被積分関数
real :: PI = 3.1415926
a = 0
b = 1
sum = 0
h = (b-a)/N
x = a
sum = sum + 1.0/3 * f(a)
do i = 1,N-1
x = a + h * i
if(mod(i,2)== 0) then
sum = sum + 2.0/3 * f(x)
else
sum = sum + 4.0/3 * f(x)
end if
end do
sum = sum + 1.0/3 * f(b)
sum = sum * h
print *,sum
end program integral

function f(x)
real :: x
f = x**2
return
end function f

>>16
fの中を変えて、求めたい積分にしてください。

>18
f90で書いたのは何か意味があるの?

>>6
すいません、スレチ失礼しました。

以下のようなプログラムでは、サブルーチン(sub) で配列Aに対して
segmentation errorが出るはずなのですが出ません。
このようなサブルーチンでもちゃんとエラーを出してくれるようなコンパイルオプションってあるんですか？

program main
implicit none
integer,parameter :: N=10
integer :: A(1:N)
call sub(N,A)
end program


subroutine sub(N,A)
implicit none
integer :: i,N,A(1:N)
do i=1,2*N ! segmentation error??
A(i)=i
enddo
return
end subroutine sub

>>21
あやまる必要はないべ

大学なら周りに聞いたほうが早いと思う。

>>22
その程度のは見出しではセグメンテーションエラーまでは行かない。
もっと大きく、１０＊＊３０くらいに書き込んで味噌。


配列はみだしを実行時に検出するオプションは、Subscript Check(添え字チェック） subchk
という名前で知られている。

これはマトモなコンパイラなら必ず付いている。しかし、個々のコンパイラによって違うので
マニュアルやヘルプをみて調べるしかない。

これは、実行が極端に遅くなるので普通debugの時だけ使う。

>>24
ﾚｽｻﾝｸｽｺ!

ifort (ver9.1) 使ってるんですが，
とりあえず/check:allとつけとけば色々検出してくれるようです．
このはみ出しもチェックしてくれました．

波数kの関数phik(k)から、xの関数phi(x)へと逆フーリエ変換をする
プログラムを作るのにてこずっています。
phi(x) = sum phik(k)*exp(i*k*x)だと思うのですが、どうもうまく作れません。
xは位置座標、kは波数で、それぞれ-5〜5、-1〜1の間でとっています。
試しに、phi(x)=2cos(0.1x)の波を作るため、k=-0.1と0.1以外は重み
(amp)を0にとりました。
何かがおかしいとは思って考えているのですが、わかりません。
どうかアドバイスを頂けたら幸いです。よろしくお願いします。

integer nx, imax, nk, jmax, i, j
real*8 xmax, dx, kmax, dk, x, k
complex*16 phi(500), phik(500), sumphik
parameter (xmax = 5.0d0, nx = 50, kmax = 1.0d0, nk = 10)
dx = xmax/nx ; imax = 2*nx + 1
dk = kmax/nk ; jmax = 2*nk + 1
do i = 1, imax
x = -xmax + (i-1)*dx
do j = 1, jmax
k = -kmax + (j-1)*dk
if (((k-0.1d0).lt.1.0d-3).or.((k+0.1d0).lt.1.0d-3)) then
amp = 1.0d0
else
amp = 0.0d0
end if
phik(j) = amp*cdexp((0.0d0, 1.0d0)*k*x)
sumphik = sumphik + phik(j)
end do
phi(i) = sumphik
write(*, *) x, dble(phi(i))
sumphik = 0.0d0
end do
end

>>26
if (((k-0.1d0).lt.1.0d-3).or.((k+0.1d0).lt.1.0d-3)) then
これが原因だ。ABSをつければ解決する。

IF ( (ABS(k-0.1d0).lt.1.0d-3) .OR.
& (ABS(k+0.1d0).lt.1.0d-3) ) THEN
ABSが無いと、ｋが0.1ないし-0.1以下のときは必ず条件を満たしてしまう。

あとsumphik = 0.0d0　の位置が悪い。入れるべき位置はここだ。
あと複素変数なんだから複素でクリアすべし。
do i = 1, imax
sumphik = (0.0d0, 0.0d0)
x = -xmax + (i-1)*dx

FORTRANでは変数の初期値は不定であって０になっているとは限らない。
変数は誰が使ったか分からないコップや茶碗だと思って、必ず洗ってから使うべし。

＝＝＝＝＝＝＝＝
ところで一般的なアドバイスだが、FORTRAN77で行くのなら、暗黙の型を絶対守るべきだ。
そうすることによりいちいち宣言部を見なくても型が分かる。

FORTRAN７７で暗黙の型を守っていない場合、もうそれは素人が書いたものとして
普通相手にされない。

普通のFORTRAN使いが、このプログラムを見た場合ｋ，ｋｍａｘがreal*8になっているのを見た時点で、
もう問題外としてその先を見るのをやめる。

変数を全部宣言すればいいじゃないかと思うかもしれないが、FORTRAN77は暗黙の型があるため
宣言忘れがエラーにならないのでヤヴァイ。たとえば、このプログラムではａｍｐの宣言を忘れている。

Ｆｏｒｔｒａｎ９０にはＩＭＰＬＩＣＩＴ　ＮＯＮＥがあるので暗黙の型宣言を無効化できるが、
それでも暗黙の型を守ったほうがいい。

26様、お答え下さりありがとうございました。
今までわからず、まだうなっておりました。絶対値を取っていないとは、
本当に馬鹿です。

FORTRANの一般的なアドバイスまでご教授頂き、本当に恐縮です。
変数はなるべく暗黙の型を使うと言うことは、今後は肝に銘じようと思います。
ampの宣言を忘れていることは、全く気づきませんでした。エラーの原因になる
可能性があり、大変危険ですね。

このたびはご助言頂き、本当にありがとう存じました。

1行目は27様の間違いです。

27様、返す返すこのたびは有り難うございました。

windows上で ifort 使ってるんですが
makefileを自動でサクッと作ってくれるツールないでしょか．．
csmkmkっていうよさげなｱﾌﾟﾘを見つけたんですが挙動が変です．

皆さんは手動で書いてるんでしょうか？

>>30
インテルのサポート掲示板にも似たような要望が出ていた。
automakeのようなツールを用意することに前向きな姿勢を見せていた。

漏れはモジュールの依存関係が複雑になってから面倒だから統合環境を使っている。

fortranにはCのようなインライン関数は無いんですか？

>>32
あるが、コンパイラによる自動最適化として行われる。
ベンダー依存のコンパイル指示行をいじれば、制御できる。

>>32
あるが、コンパイラによる自動最適化として行われる。
ベンダー依存のコンパイル指示行をいじれば、制御できる。

subroutineで戻り値もたせるにはどうしたらいいの？
うまいやり方あるんだよね？おしえて☆

PROGRAM REIDAI4_3
IMPLICIT NONE
INTEGER :: I,N
REAL :: S,SS,HEIKIN,SIGMA,BUNSAN
REAL :: A(100)
WRITE(*,*) 'N='
READ(*,*) N
DO I=1,N
WRITE(*,'(1X,A2,I2,A2)') 'A(',I,')='
READ (*,*) A(I)
END DO
S=0.0
DO I=1,N
S=S+A(I)
END DO
HEIKIN=S/N
SS=0.0
DO I=1,N
SS=SS+(A(I)-HEIKIN)**2
END DO
BUNSAN=SS/REAL(N)

SIGMA=SQRT(BUNSAN)
WRITE(*,*) 'ヘイキンチ ハ',HEIKIN
WRITE(*,*) 'ブンサン ハ',BUNSAN
WRITE(*,*) 'ヒョウジュンヘンサ ハ',SIGMA
STOP
END PROGRAM REIDAI4_3

何故かエラーが出ます・・・・・

>>36,37
超能力者を呼べ

Fortran90です。
mainルーチンでallocateした配列をサブルーチンに渡して、その中で
サイズを変えて再アロケート(deallocate->allocate)することは出来る
のでしょうか？

randomという数値ルーチンはwindows上で動かすと
問題なく動作しているようなのですが、
unix環境になるとまったく動かなくなります。
どうすればいいでしょうか？
Fortran90です。

>>40
マニュアル嫁！！

Fortran90固有の乱数サブルーチンはRANDOM_NUMBER（）だろ。
RANDOMはベンダー拡張のルーチンだ。

エラーが出たら文法書を嫁！　法律書と同じだ。トラブッたら文法書を開け！

>>37
どういうエラーが出るか書け。答えようがない。

ざっと見たところプログラムはおかしくない。
コンパイルして動かしてみたが、ちゃんと動いているように見える。
データ数３個でしか試してないがｗｗｗ

>>39
確か出来たはず・・・・・　最悪でもポインターで送ればできる。

コンパイラでやってみたらできた。一応F90厳格オプションを掛けたのでおｋだろ。

>>35
サブルーチンに戻り値を持たせたければ、引数を１こ増やせばいいだろ。

それが嫌ならFUNCTIONを使え。

しかしFORTRANの伝統では、FUNCTIONは返す値が１個で、引数はINTENT(IN)で
副作用なし、と取ることに不文律でなっている。数学のFunctionを連想するので、
これが自然に思える。

Cのような何でもFunctionにするのは愚かな浅知恵だと思う。

>>41
マニュアルって、どこにあるんでしょうか？

>>45
見つかるまで探せ！　幸せの青い鳥は足元にいるんだ！

下記のfortran90のプログラムをfortran77にするとどういうふうになるのか
教えていただけませんか？

PROGRAM vipper
IMPLICIT NONE
REAL :: theta1, theta2, phi1, phi2, alpha, anu, pi
REAL :: rot_z(3, 3), rot_y(3, 3), vec1(3), vec2(3)
CHARACTER :: ns, ew
pi = 4.0 * ATAN(1.0)
alpha = 6242.2 / 6370.0
anu = ( 180.0 - 86.6 ) * pi / 180.0
theta1 = ( 90.0 - ( 35.0 + 40.0 / 60.0 ) ) * pi / 180.0
phi1 = ( 139.0 + 12.0 / 60.0 ) * pi / 180.0
!
rot_z(1, :) = (/ COS(phi1), -SIN(phi1), 0.0 /)
rot_z(2, :) = (/ SIN(phi1), COS(phi1), 0.0 /)
rot_z(3, :) = (/ 0.0, 0.0, 1.0 /)
!
rot_y(1, :) = (/ COS(theta1), 0.0, SIN(theta1) /)
rot_y(2, :) = (/ 0.0, 1.0, 0.0 /)
rot_y(3, :) = (/ -SIN(theta1), 0.0, COS(theta1) /)
!
vec1(1) = SIN(alpha) * COS(anu)
vec1(2) = SIN(alpha) * SIN(anu)
vec1(3) = COS(alpha)
!
vec2 = MATMUL(rot_z, MATMUL(rot_y, vec1))
!

theta2 = ASIN(vec2(3)) * 180.0 / pi
phi2 = ATAN2(vec2(2), vec2(1)) * 180.0 / pi
!
PRINT *, theta2, phi2
ns = 'N'
IF (theta2 < 0.0) ns = 'S'
ew = 'E'
IF ( phi2 < 0.0) ew = 'W'
PRINT '(i4, a, i2, 1x, a)', INT(ABS(theta2)), ' ', NINT(( ABS(theta2) - INT(ABS(theta2)) ) * 60.0), ns
PRINT '(i4, a, i2, 1x, a)', INT(ABS( phi2)), ' ', NINT(( ABS( phi2) - INT(ABS( phi2)) ) * 60.0), ew
STOP
END PROGRAM vipper

よろしくおねがいします。

>>47
(1) ソースを自由形式→固定形式に変更する
(2) 変数の宣言をF77式に変更する（REAL :: hoge → REAL hoge）
(3) コメント( ! ) をF77式に変更する
(4) rot_z(1, :) = (/ ... /) 等配列に代入している部分は各要素毎に１つずつ代入する
(5) MATMULに相当する手続きを自作する
(6) 関係演算子( < ) をF77式( .LT. )に変更する

まだあるかも知れんが、とりあえずこんなところか。

>>43
どうも。
出来るのですか。とするとやり方が間違っていたのですね。
fortran90のポインタは使ったことありませんが調べてみます。

>>50
こんな感じでできた。IntelFortran
MODULE m_test
IMPLICIT NONE
CONTAINS
!
SUBROUTINE chinpo(kk)
IMPLICIT NONE
INTEGER, ALLOCATABLE, INTENT(IN OUT) :: kk(:)
ALLOCATE (kk(10))
kk = 999
RETURN
END SUBROUTINE chinpo
!
END MODULE m_test
!=============================================================
PROGRAM unko
USE m_test
IMPLICIT NONE
INTEGER, ALLOCATABLE :: n(:)
CALL chinpo(n)
PRINT *, n
STOP
END PROGRAM unko

[7] g77 -o r7 r7.f
/tmp/cc9A4yij.s: Assembler messages:
/tmp/cc9A4yij.s:136: Error: symbol `MAIN__' is already defined
/tmp/cc9A4yij.s:299: Error: symbol `MAIN__' is already defined


こういうエラーが出ます・・・・

>>52は>>37です。

windows上で使えるフリーのコンパイラはありますか？

>>52
はじめっからエラーメッセージを書くべし。

それはメインルーチンが複数あるということだから、メインルーチンと勘違いされるような
行があるということだ。

ｇ７７についてはよく知らんが、１３６や２９９がエラー行を表しているんじゃないか？
その辺の行を見て味噌。
一見空行でもみょうな制御コードが混じっていたりするとこういう事が起こる。

たとえば２バイトの日本語文字とか。
そもそもｇ７７は２バイト文字に対応しているのかね？

>>54
本スレの>>1をみろ！
FORTRAN �W
http://pc8.2ch.net/test/read.cgi/tech/1163319215/

三次元のポアソン方程式を解きたいのですが、どんなプログラムにしていったらいいでしょうか？？
fortran７７で書いているのですが、なかなか上手くいかなくて…
お願いします。教えてください。

もーーどーーしてもうまくいかないので質問させてください。

自宅のwinXPでfortran95を使いたくて、本スレにリンクしてある「FTN77/95のインストール方法」を参照していたのですが・・

まず、リンク先から「Salford FTN95 Personal Edition」をデスクトップにDLしました。

次に、インストーラを起動させてC:\Program Files\Salford Software\FTN95にinstallしました。

で、参照ページを見ると.NETを持ってない場合はDOSでテストしてみると書いてあったので、
DOSでcd C:\Program Files〜に移動して、ftn95とコマンド入力。
そしたら「Ver4.9 copyright 云々」と出たので多分OK。

次に、installしたフォルダに一緒に入っていたfcpad.exeを起動すると書いてあるので、起動。
するとコンパイラのパスを求められたので、C:\Program Files\Salford Software\FTN95と入れてOKを押しました。

そしたら、「コンパイラのパスが正しくありません」と言われてしまいます。
コンパイラftn95.exeを右クリック→プロパティで「場所」をコピーしてやってみてもやっぱりだめでした。

長文ごめんなさい、自分のやり方はどこが間違ってるのでしょうか？

>>57
>するとコンパイラのパスを求められたので、
>プロパティで「場所」をコピーしてやってみても


C:\Program Files\Salford Software\FTN95\ftn95.exe
でもダメ？

>>58

で　き　ま　し　た　。

なんか単純常識的なことを知らなかったようで・・
即レスありがとうございました。

>>59
いいのよ〜♪

>>57,59

CPad for Salford FTN77もいいけど FTN95付属のPlato3 IDEも使ってやれよー。
日本語が使えなかったりプラトン先生の顔が怖かったりするけど
ちゃんと自由形式で色分けしてくれてデバッガもついとるでよー。

>>47　世話が焼けるなー
PROGRAM vip77
REAL rotz(3, 3), roty(3, 3), vec1(3), vec2(3), tmp(3)
pi = 4.0 * ATAN(1.0)
degrad = pi / 180.0
alpha = 6242.2 / 6370.0
anu = ( 180.0 - 86.6 ) * degrad
theta1 = ( 90.0 - ( 35.0 + 40.0 / 60.0 ) ) * degrad
phi1 = ( 139.0 + 12.0 / 60.0 ) * degrad
C
rotz(1, 1) = COS(phi1)
rotz(1, 2) = -SIN(phi1)
rotz(1, 3) = 0.0
rotz(2, 1) = SIN(phi1)
rotz(2, 2) = COS(phi1)
rotz(2, 3) = 0.0
rotz(3, 1) = 0.0
rotz(3, 2) = 0.0
rotz(3, 3) = 1.0
C
roty(1, 1) = COS(theta1)
roty(1, 2) = 0.0
roty(1, 3) = SIN(theta1)
roty(2, 1) = 0.0
roty(2, 2) = 1.0
roty(2, 3) = 0.0
roty(3, 1) = -SIN(theta1)
roty(3, 2) = 0.0
roty(3, 3) = COS(theta1)

C
vec1(1) = SIN(alpha) * COS(anu)
vec1(2) = SIN(alpha) * SIN(anu)
vec1(3) = COS(alpha)
C
CALL mulmat(3, roty, vec1, tmp)
CALL mulmat(3, rotz, tmp, vec2)
theta2 = ASIN(vec2(3)) / degrad
phi2 = ATAN2(vec2(2), vec2(1)) / degrad
PRINT *, theta2, phi2
STOP
END
!
SUBROUTINE mulmat(n, a, b, c)
REAL a(n, n), b(n), c(n)
DO 10 i = 1, n
c(i) = 0.0
DO 20 j = 1, n
c(i) = c(i) + a(i, j) * b(j)
20 CONTINUE
10 CONTINUE
RETURN
END

使っているのは高校でも習う座標の回転だけだ。
高校では２次元、ここでは３次元だが、ｙ軸とｚ軸に関して廻すだけだから
高校のものと本質的には変わりはない。ＣＯＳ、ＳＩＮが湧いて来るのはそれが理由だ。

あと、ＦＯＲＴＲＡＮの三角関数はラジアンで与えることになっているから、
度表示をラジアンに直さなければならない。π/１８０はそこから湧いてくる。

また、角度表示は古代メソポタミア文明以来の６０進法が使われているので、
度、分、秒を度の小数点表示にしてから、さらにラジアンにしている。
それが６０で割ったりしている部分だ。

あと、座標としては極座標をとるが、北緯、東経などを極座標に直す必要がある。
座標系は右手系を取るのがふつう。（右手系とは右手のお父さん指がｘ軸、
お母さん指がｙ軸、のっぽのお兄さん指がｚ軸だ。）


考え方はこうだ。まず八王子が北極にあるとする。北極からミサイルをぶっ放した
到着点をしるのはたやすい。その地点をｖｅｃ１に入れている。
ここで北極点を本来の八王子の位置に持っていくように座標軸を回転させる。
このときｖｅｃ１にも同じ回転をかければ求める地点が分かる。それがvec2だ。

その他、問題文にもあるが、弧の長さと半径が分かれば角度が出る。
方位角ν（読み方はニュー速のニューだぞ）は、本来の八王子から見て真北から
時計回りにとると考える。（問題文からは判然としないが、常識的にはこう取るのが普通）
北極点から考える場合は南極向きで時計回りになる。

訂正

＞あと、座標としては極座標をとるが、北緯、東経などを極座標に直す必要がある。

角度の取り方は極座標の取り方をとって、ｘ、ｙ、ｚの直交座標系を右手系でとった。

>>62
解説乙。
本当は出題者が地球半径を与えている筈ですね。
もし地球を楕円体として考えた場合、かなり面倒そうですね。

>>66
それは強烈に難しいｗ

>>51
ありがとうございます。出来ました。
subroutineをmoduleで包まないといけないのですね。
最近77から移ってきたのでどうもmoduleの使い方がいまいち
分かりません。

subroutineのなかにsubroutineとかfunction定義できるの？
自分で実験したほうが早いって？ｗｗｗ

>>69
できる・・・・が、それは普通のサブルーチンやファンクションとは違っている。
親ルーチンの変数をプチ・グローバルに共有している。

うっかりiとかｊとかローカルに宣言しわすれたまま使ってしまうと、
親ルーチンのほうのiやjが書き換わってとんでもないことになる。

好みの問題だが、漏れはあまり好きくない。

CPad for Salford FTN77
のコンパイラのパスが分かりません。
教えて下さい

>>71
>>57のへん見れ

>>72
C:\Program Files\Salford Software\FTN95\ftn95.exe
と入力しても
”コンパイラのパスが正しくありません”と表示されます
ちなみにOSはXPです

>>73
エクスプローラで"ftn95.exe"を検索したら？

>>73
つーかお前の使いたいコンパイラはFTN77かFTN95かどっちなんだ？
FTN77しか入れてないのにftn95.exeを探しても見つかりっこないぞｗ

FTN77です

もぅ、ぉにぃちゃんしっかりしてよ〜

>>71
なんとかできました

>>78
どうやったら「なんとかでき」たの？

FTN77peをダウンロードして,コンパイラのパスに
FTN77pe〜exeと入力してら、できました

FTN95ってどうよ？

>>81
コンパイル時/実行時のエラーチェックがけっこう強力で
デフォルトで未使用変数の警告なんかも出してくれるので、
入門用やエラーチェック用にイイ感じ。
Personal版がフリーで使えるし。

その他の特徴としては
・Visual Studio 2003/2005上で使える（もちろんVSは別途購入する必要がある）
・VSがなくても独自のIDEが付属する（日本語は使えないが）
・（今時のFortranに必要かどうかはともかく）インラインアセンブラが使える
・（これも実用性のほどはともかく）.NETなプログラムも作れる
・GET_COMMAND_ARGUMENT()、[...] 等Fortran2003の一部に対応（他のコンパイラでも大抵使えるレベルだが）

>>82
�dクス。
意外に高機能だな。
サポートとかはどうなんだろう？パッチとか結構でてるんだろうか？

intel fortran ver 9.0 for winを使っております。
自分でプログラムを組んでコンパイルできるところまでは到達したのですか実行すると
forrtl: severe (168):Program Exception - illegal instruction
Image PC Routine Line Source

Stack trace terminated adnormally.

と実行してくれません。なにがおかしいのでしょうか?
ソースをのせるべきなんでしょうがnetnumpacから主に引っ張ってきている分権利がややこしそうです。

よろしくお願いします。

program main
common/ numbb / nnn,lll,dz,ameson
integer :: nnn, lll
real*8 :: dz, ameson, dr(251),scr,hx
open (unit=10, file='functest0.txt', status='unknown')
open (unit=11, file='functest1.txt', status='unknown')
open (unit=12, file='functest2.txt', status='unknown')
open (unit=13, file='functest3.txt', status='unknown')
open (unit=14, file='functest4.txt', status='unknown')
open (unit=15, file='functest5.txt', status='unknown')
nnn = 6
lll = 5
dz = 60.0d00
ameson = 413.536d00
dr(1)=exp(-8.8)*dz
hx=5.d-2

do l=2, 251
dr(l)=dr(1)* exp(hx*(l-1))
end do
zero = 0.0
y = 0.0
do i = 1, 251
call GASND (zero, dr(i), screen, 40, Y, ICON)
scr = y
write(10,'("i="i2," dr(i)="e15.8," scr="e20.10)')i,dr(i),scr
end do
close (unit = 10)
close (unit = 11)
close (unit = 12)
close (unit = 13)
close (unit = 14)
close (unit = 15)
end

function screen(ddr)
! it's the Dnl(r)
common/ numbb / nnn,lll,dz,ameson
real*8 :: a1, a2, a3, al1, arn
integer :: nnn, lll, l5
real, external :: DPLAGG
integer, external :: factorial
real, intent(in) :: ddr
a0 = ameson
arn = - a0 * ddr
aa1 = 2.0d00*dz/(nnn*a0)
nl = nnn-lll - 1
nll = nnn + lll
a2 = factorial(nl)
a3 = factorial(nll)
a4 = 2.0d00*dble(nnn)*a3*a3*a3
a5 = aa1*(a2/a4)
l5 = 2*lll + 1
l6 = 2*lll
a6 = (arn)**(l6)
al6 = dble(l5)
al7 = DPLAGG(nll,al6, arn)
al8 = al7*al7
ss = a5*exp(arn)*al8*a6
screen = ss
write(11, '("nnn="i2," ss="e20.10)') nnn, ss
end function screen

function factorial(n)
integer :: fact
integer, intent(in):: n
fact = 1
do i = 1, n ,1
fact = fact*n
end do
factorial = fact
write (12, '("n="i1," fact="i10)') n, fact
end function factorial

FUNCTION DPLAGG (N, A, X)
IMPLICIT REAL (8)(A - H, O - Z)

END FUNCTION DPLAGG
SUBROUTINE GASND (A, B, FUN, N, Y, ICON)
IMPLICIT REAL (8)(A - H, O - Z)
EXTERNAL FUN
DIMENSION X (50), W (50)
SUBROUTINE WGLEGD (NP, PT, WT, ICON)
! LEGENDRE-GAUSS FORMULA
IMPLICIT REAL (8)(A - H, O - Z)
DIMENSION X (25, 50), W (25, 50), PT (NP), WT (NP)

netnumpacからのソースはプログラム名だけにしました
変数の定義があったりなかったりで非常にみづらいとは思いますがよろしくお願いします。

>>84 これは試したかい？
Re: forrtl: severe (168) Program Exception - Illegal instructiom
Reply Quote
Turns out that the problem was caused by an older-generation processor not understanding
newer instructions. The application had been compiled with the "Generate most optimized code"
(/fast) setting, which implies /arch:host. Unfortunately, the project settings box display doesn't
reveal the implied switches, leading to this sort of problem.

Steve

http://softwareforums.intel.com/ISN/Community/en-US/forums/thread/106764.aspx

>>83
SilverfrostのホームページからForumを見に行くと
ここのAnnouncementsにパッチの情報も載っている。

>>90
英語をよみくだくと /arch:hostのコンパイラオプションをしろ　ですよね
やってみましたが変化なしでした。
環境はpen4 2.4BGHz 1024MB 845Echipset winxp sp2なんですがね

>>92
debugモードでtracebackもonにしてどのソース行でおかしくなってるか調べてみたら？

>>92
Pentium 4の2.4BGHzだと拡張命令はSSE2までだから、
どこかにSSE3等を使ったコードが紛れ込んでるんですかね。

別のPCでコンパイルしたライブラリをリンクしているのなら
そのライブラリがSSE3命令を使っている可能性も考えられます。

>>93
debugモードは使ったことがないけどやってみます
>>94
ライブラリはもってきていないはずなので
SSE3をつかったコードが紛れ込んでるんですかね

>>95
つーかVisualStudio使わないでコマンドラインでやってるの？デフォはdebugのはずだが。

debugモードを使わないという神経が分からんぞい！
まずは警告オプションも全部ｏｎにしてやりんしゃい！

fortran77でファイルに書き込みをするときに
open(2,file='filename.txt')
とかけばfilename.txtに結果が書けるのですが、

ひとつのプログラムでa=1〜100まで変化したときに
（do a=1,100 〜continueを利用）
a=1のときの結果はfilename1.txtに記録
　　　・・・
a=100のときの結果はfilename100.txtに記録
するにはどのようにすればいいのでしょうか？

可能であればa=i,b=j(i,jに自然数が入る)のときに
filename_a=i_b=j.txt
に書き込めるように、お願いします。

>>97
前スレでも出てたと思うが、FORTRAN77の場合だと
fmt, fname は適当な長さの文字列変数として

write(fmt,100) int(log10(real(i)))+1, int(log10(real(j)))+1
100 format('(''fname_'',I',I1,',''_'',I',I1,',''.txt'')')
write(fname, fmt) i, j
open(10, file=fname, ...

書式に 'I0' を指定すると「その数値を表現するのに最小の欄幅」を取ってくれる処理系なら

write(fname,'(''filename_'',I0,''_'',I0,''.txt'')') i,j
open(10, file=fname, ...

個人的には数値は適当な桁数にそろえて
filename_001_001.txt, ..., filename_010_010.txt, ..., filename_100_100.txt, ...
としたほうが後でデータファイルをいじったりする時に都合がよいように思う。

write(fname,'(''filename_'',I3.3,''_'',I3.3,''.txt'')') i,j
open(10, file=fname, ...

>>96
コンパイラマシンは混んでいていつもコンパイルのみに使っていたのでdebugモードは未経験でした

debugしてみたとこと
dr(1)=exp(-8.8)*dz がひっかかっており
sub = -8.8
dr(1) = dz *exp(sub) としたら動くようになりました。

exp(-8.8)がSSE3なりを使用しているんでしょうか

最適化オプションを無効にすればSSE/2/3は使われないはずだが。

>>99
ReleaseでSSE抑止命令（CPUを旧型に指定）すれば、SSEが原因なのかどうか分かる。

定数式だから変な最適化でバグったのかね？

求めるＰ３の値がちょうど３（限りなく近い）値になるときのdaの値
を小数点６桁ぐらいまで求めたいのですが、どうすれば次のプログラムを
効率よくdaの値が一回で出せるプログラムにできるか知恵をお貸しください。

DIMENSION E22(2),oo(2)
REAL k,pi,a,o,z,Kf,kf1
pi=3.1415926
DO 30 da= 1.000 , 0.800 , -0.002
DO 30 n4= 8 , 8

k=5.319
a=k/0.529177
o=(a**3)/2.0*da
oo(1)=o
z=1.0

rs=(3.0*o/(z*4.0*pi))**(1.0/3.0)
Kf=(9.0*pi/4.0)**(1.0/3.0)/rs
Ef=Kf**2
c=2.0*1.0/(1.0+0.0155*pi/Kf)
Ei=1.79186*z**(5.0/3.0)/(3.0*o*rs)
P=-Ei*1.4703*10**4
E0=-4.42*z/(3.0*o*rs**2)+0.916/(3.0*o*rs)+0.031/(3.0*o)
P0=-E0*1.4703*10**4
Rmk=3.182
q0=0.89*2*0.3947587
uk=q0*Rmk*cos(q0*Rmk)/(sin(q0*Rmk)-q0*Rmk*cos(q0*Rmk))
E1=-2.0*pi*z*2.0/o**2*(z*Rmk**2*(1.0+2.0/3.0*uk))
P1=-E1*1.4703*10**4
DO 20 L= 1 , 2
E2=0.0

IF(L.EQ.1) GO TO 100
o=o*0.996
oo(L)=o
a=(o*2.0)**(1.0/3.0)
rs=(3.0*o/(z*4.0*pi))**(1.0/3.0)
Kf=(9.0*pi/4.0)**(1.0/3.0)/rs
Ef=Kf**2
c=2.0*1.0/(1.0+(0.0155*pi)/Kf)

100 DO 10 n1= -n4 , n4
DO 10 n2= -n4 , n4
DO 10 n3= -n4 , n4
IF(n1==0.and.n2==0.and.n3==0)GO TO 10
IF(MOD(n1+n2+n3,2)/=0)GO TO 10

G=2.0*pi/a*SQRT(FLOAT(n1*n1+n2*n2+n3*n3))
f=G**2/(2.0*(G**2+c*Kf**2))
ae=1.0+(1.0-f)*4.0*pi*z*2/(o*G**2)*3.0/(2.0*Ef)*
&(1.0/2.0+((4.0*Kf**2-G**2)/(8.0*Kf*G))*ALOG(ABS((2.0*Kf+G)/(2*Kf-G))))
V=(-4.0*pi*Z*2.0/(o*G**2))*((1.0+uk)*cos(G*Rmk)-uk*sin(G*Rmk)/(G*Rmk))

E3=1.0/(4.0*pi*2.0/(o*G**2))*V**2*(1.0/(1.0-f))*(1.0/ae-1.0)
E2 = E2 + E3

10 CONTINUE
E22(L)=E2/2.0
20 CONTINUE

E8=(E22(2)-E22(1))/(oo(2)-oo(1))
P2=-E8*1.4703*10**4
P3=P+P0+P1+P2

WRITE(6,50)P3
50 FORMAT(F15.8)
30 CONTINUE
WRITE(*,*)P,P0,P1,P2
WRITE(6,601) P3
601 FORMAT(1H , F15.8)
STOP
END

>>102
その質問の仕方じゃ　答えようがない。
もう少し何をやっているのか、根本から家

>>98
出来ました。ありがとうございます！

fortran90です。
いまwindows上で動かしてたプログラムを
Linuxに移そうとしてるんですが

Unable to open MODULE file dflib.mod

とエラーが出てきます。
どうすればいいでしょうか？

Unable to open MODULE file dflib.mod
　　　↓
エキサイト翻訳
　　　↓
MODULEファイルdflib.modを開くことができません。

>>107
ほかの人にやってもらう

>>107
dflib.mod はDECの拡張モジュールだから、そのままではLinuxでは動かない。

DECの拡張ルーチン中で用いているものがPOSIX互換のルーチンだったら
LINUX上に移植できるだろう。

そうでないなら、自分で回避策を考えねばならない。

とりあえず、USE　DFLIB　となっている行を全部つぶすべし。
そうして出てきた未定義ルーチン・エラーをひとつひとつ調べるべし。

>>110
やはりそうですかー。気が遠くなりそうですが、やってみます。
ありがとうがざいました。

>>111
あぁ、もしＬｉｎｕｘ上のＩｎｔｅｌＦｏｒｔｒａｎを使っているなら、ＤＦＯＲＴをＩＦＯＲＴに変えて試してみるといい。
今のＷｉｎ上のＩｎｔｅｌｆｏｒｔｒａｎではＤFORTはＩＦＯＲＴになっている。
Ｉｎｔｅｌが移植用に互換モジュールを用意してくれているかもしれない。

とりあえず一応マニュアルを検索してみるといい。

どっちにしろコマンドライン引数を取るために使っている程度なら、ｇ７７でも移植可能だ。

このプログラムだとわざわざDOループのdaの値を何回も変えてよりP3の値が３に近い値のda
を求めなければいけないのですが、それを実行一回で求めたいのですがどうしたら
よいでしょうか？

>>113
見た感じ非線形だから１回は無理で、はさみうち法とかでイターレションするしかナインでないの？
というかもし１回で求まるなら、はじめっからｄａについての式を逆関数的に定式化してあるのでは？

>>112
なるほど。試してみます。
ありがとうございました！！

下記のプログラムについて、fortran77で行うと実行できないのです＞＜
何がおかしいですかね＞＜
よければ、教えていただけませんか？


dimension ia(100),vw(12),ivw(12),vw2(100)
data am,ix,n/1.0 ,0 ,100/

call poiss(am,ix,ia,n,vw,ivw,vw2)

write(*,*) (ia(i),i=1,n)

stop
end
　
　 subroutine poiss(am,ix,ia,n,vw,ivw,vw2)
dimension ia(*),vw(*),ivw(*),vw2(*)
real*8 x,y
data bm/-1.0/


if (am.ne.bm) then
bm=am
vw(1)=exp(-am)
vw(2)=vw(1)*am
m=2.0*am+10.0

do 10 k=2,m-1
vw(k+1)=vw(k)*am/k
10 continue

do 30 k=1,m-1
if (vw(k).ne.0.0) then
x=vw(k)
do 20 j=k+1,m
y=x+vw(j)
if (y.eq.x) go to 50
x=y
vw(j)=x
20 continue
go to 40
end if
30 continue

k=m
40 j=m

50 do 70 i=j,m
vw(i)=1.0
continue

mm=j-k+1
ai=mm

do 70 i=1,mm
af=(i-1)/ai
do 60 while(af.gt.vw(k))
k=k+1
60 continue
ivw(i)=k
70 continue
end if

call unfm2(ix,vw2,n)

do 90 i=1,n
k=vw2(i)*ai+1.0
if(k.eq.mm+1) k=mm
j=ivw(k)
do 80 while(vw2(i).gt.vw(j))
j=j+1
80 continue
ia(i)=j-1
90 continue

return
end

subroutine unfm2(ix,vw2,n)
implicit real*8 (a-h,o-z)
dimension vw2(*)
real*8 aa/32771.d0/,b/1234567891.d0/,x
real*8 c/2147483648.d0/,ci/4.6566128730773925d-10/

x=ix

do 10 i=1,n
x=dmod(aa*x+b,c)
vw2(i)=x*ci
10 continue
ix=x

return
end

>>116
質問の仕方がなっちょらん。
本文からコンパイルエラーなのか実行時エラーなのかすら読み取れない。
エラーメッセージを出すくらいの知恵も無いのか！ｗｗｗ　

と叱られるよ

>>116
一応DEGITAL Visual Fortranでコンパイルしてみた。
do 70のループが原因でコンパイルが通らない。
後ろの70を170にして、その前のcontinueしかなかったところを70にした。
コンパイルはできるが実行中にArray Bounds Exceededが起こる。

プログラムの中身を追いかけるのは面倒だからやらない。

高校の課題なんですが、妹に教えられなくて兄としての威厳がかかっていますｗ
正直、プログラムはさっぱりなんですが、シスアド初級持っているために安請け合いしてしまったのが間違いですたｗ

以下の問いなんですがよろしくお願いします。


乱数（0〜1までの範囲）を7個発生させなさい。
その発生させた7個の乱数はそれぞれa,b,c,d,e,f,gと定義する。
このa〜gの中で、最大値と最小値をとるものを選別するプログラムを作りなさい。

というものなんですが、例として何か提示していただけませんか？
それを本みながら自分で理解して、妹に説明できるようにしたいんで。

program main
a=rand()
b=rand()
c=rand()
d=rand()
e=rand()
f=rand()
g=rand()
write(6,*) 'MAX ',max(a,b,c,d,e,f,g)
write(6,*) 'MIN ',min(a,b,c,d,e,f,g)
stop
end

>>123

この

write(6,*) 'MAX ',max(a,b,c,d,e,f,g)
write(6,*) 'MIN ',min(a,b,c,d,e,f,g)

の、

maxとminっていうのは最初に宣言しなくても機能するんですか？

>>124
機能する。 最大値／最小値を返す関数の総称名として規格で規定されている。
RAND() の方は使えない可能性がある。

>>125

ありがとうございました。

>>122 Fortran90風には国府田。

PROGRAM chinpoko
IMPLICIT NONE
REAL :: a(7)
CALL RANDOM_NUMBER(a)
WRITE(*, *) 'MAX=', MAXVAL(a)
WRITE(*, *) 'MIN=', MINVAL(a)
STOP
END PROGRAM chinpoko

ちょっと122番と若干かぶる要素があるんだけど、

0から1までの乱数�@〜�Dまで存在すると仮定する。

問1.乱数をそれぞれ、�@*�A、�@*�B、�@＊�C、�@*�D…�C*�Dとすべてのパターンの積を求める問題。

問2.問1で得られた結果に100を掛けて、小数点以下を切り捨てる。

問3.問2で得られた全ての結果に対して50より小さい場合残し、そうでない場合は残さない。残すもの、残さないものを分けろ。

問4.問3で得られた結果（残った数値）を以下のような形でまとめるようにせよ。

yi(y1,y2,y3…yn)

問5.問4で得られた結果からその中で最大値を求めよ。


という5つの問題がある。元は一個の問題なんだが、問題文が無駄に長いので俺が解釈して５つに分割してみた。


質問1.問1を俺なりにやってみたんだが、いまいちプログラムがわからないので自力で全て値を�@*�Aというように入力していったのだが、もっといい方法はないかな？
質問2.50以下のものとそうでないものを分けるとき、使う分は「if」でいいのか？
質問3.yi(y1,y2,y3…yn)とあるけど、これはdimension？使えばいいのか？


わかりにくいけど質問の回答よろしくお願いします。

できたら、問ごとにプログラムの例文作ってくださると理解しやすいのでありがたいです。

>>129
ＦＯＲＴＲＡＮ７７かＦｏｒｔｒａｎ９０か指定はあるか？

特にないです。今使ってるのは、visual fortranってやつなんでたぶん新しい方だと思う。

>>128
PROGRAM q1
IMPLICIT NONE
INTEGER, PARAMETER :: n = 5
INTEGER :: i, j, k
REAL :: x(n), y(n * (n - 1) / 2)
CALL RANDOM_NUMBER(x)
k = 0
DO i = 1, n
DO j = i + 1, n
k = k + 1
y(k) = x(i) * x(j)
END DO
END DO
PRINT *, x
PRINT *
PRINT *, y
STOP
END PROGRAM q1

大文字はＦｏｒｔｒａｎに存在する命令。小文字は自分の定義する変数名など。
RANDOM_NUMBER()は0以上１未満の数を返す。
>>132　インデント

>>128
PROGRAM q2
IMPLICIT NONE
INTEGER, PARAMETER :: n = 5
INTEGER :: i, j, k
REAL :: x(n), w(n * (n - 1) / 2)
REAL, ALLOCATABLE :: y(:)
CALL RANDOM_NUMBER(x)
k = 0
DO i = 1, n
DO j = i + 1, n
k = k + 1
w(k) = x(i) * x(j)
END DO
END DO
w = AINT(w * 100.0)
PRINT *, x
PRINT *
PRINT *, w
STOP
END PROGRAM q2

問２
INT（）は切り捨て整数化。AINTは切り捨て整数化を実数で返す。
>>133

PROGRAM q3
IMPLICIT NONE
INTEGER, PARAMETER :: n = 5
INTEGER :: i, j, k
REAL :: x(n), w(n * (n - 1) / 2), y((n * (n - 1) / 2))
CALL RANDOM_NUMBER(x)
k = 0
DO i = 1, n
DO j = i + 1, n
k = k + 1
w(k) = x(i) * x(j)
END DO
END DO
w = AINT(w * 100.0)
k = 0
DO i = 1, n * (n - 1) / 2
IF (w(i) < 50.0) THEN
k = k + 1
y(k) = w(i)
END IF
END DO
print *, 'no. of less than 50 =', k
PRINT *, y(1:k)
STOP
END PROGRAM q3

すいません。FORTRANについての質問かどうかわからないのですが…。
2次元配列DDDD(512,512)に物理量、1次元配列X(512),Y(512)に座標値が
格納されている場合で、Text column形式で出力するならば
OPEN(10,FILE='TXTCOL.D')
DO L = 1, 512
DO K = 1, 512
WRITE(10,*) X(K),Y(L),DDDD(K,L)
END DO
END DO
というのは理解しているのですが、Binary Columns形式ですとどのように
書けば良いのでしょうか？
Binary Matrix形式ですと
OPEN(20,FILE='BINMAT.D',ACCESS='direct',FORM='unformatted',recl=8*512*512)
WRITE(20,rec=1) ((DDDD(K,L),K=1,512),L=1,512)
とすれば良いと勝手に思っています(※これだと座標値が入ってません)。
そもそもバイナリ書き出しを全然、理解しておりません。
どなたかお教えいただけないでしょうか？

>>135
いまいち何をしようとしているのか理解できないが、単にバイナリー形式で配列を
書き出したいなら書式なし形式で出せばいい。
WRITE(10) DDDD
でおｋ
WRITE(10) X,Y,DDD
なら、Xのすべて、Yのすべて、DDDのすべて
の順で書き出される。

直接アクセス方式は、要素をランダムにアクセスしたいときの形式。
OPEN文のRECLは単一要素の大きさなので５１２＊５１２はいらない。
あと、処理系によってバイトかワードか単位要素か単位が違うのでマニュアルで
確かめる必要がある。
WRITE文でのRECは位置を指定するので、要素ごとに１づつ増やしてゆかねばならない。

直接アクセス方式はしばらく使ってないから細部を忘れたｗ
マニュアルを穴があくまでにらんで読んでくれ。

fortran77を使用して、cos(pi*x)を初期値とした熱方程式の時間発展を計算したいと思います。
範囲は０≦ｘ≦１です。
これで作られたdatファイル（'C:\out.dat'）が(exp((-(pi)**2)*0.05))*cos(pi*x)
とぴったり重なる様な結果を得たいのですが、ずれてしまいます。
改善点をご指摘いただければ幸いです。よろしくお願いします。
program heat equation
implicit none
integer i,j,nstep,n,ndim
parameter (ndim=5001)
double precision flam,f,h,u,c1,c2,fi,t,tmax,dt
double precision a,b,c
dimension a(ndim),b(ndim),c(ndim),u(ndim)
n=10
dt=0.005d0
tmax=0.05d0
h=1.0d0/float(n)
flam=dt/(h**2.0d0)
c1=1.0d0-2.0d0*flam
c2=flam
do 1 i=1,n-1
a(i)=c2
b(i)=c1
c(i)=c2
1 continue
do 2 i=1,n+1
fi=i-1
u(i)=f(fi*h)
2 continue
OPEN(11,FILE='C:\out.dat')
nstep=0.0

3 continue
do 4 i=1,n-1
4 u(i)=b(i)*u(i+1)+a(i)*u(i)+c(i)*u(i+2)
do 5 i=n,1,-1
5 u(i+1)=u(i)
a(1)=0.0d0
b(1)=c1
c(1)=2.0d0*c2
a(n+1)=2.0d0*c2
b(n+1)=c1
c(n+1)=0.0d0
nstep=nstep+1.0
t=dt*float(nstep)
do j=1,n+1
if(mod(nstep,5).eq.0)then
write(11,*) j*h-1.0d0/n,u(j)
end if
end do
write(11,*)' '
if(t+dt.lt.tmax+0.10d0*dt) go to 3
end

function f(x)
implicit none
double precision x,f,pi
pi=3.14159265358979323846264338327950288
f=dcos(pi*x)
end

>>137
ぴったり重なるといっても、どの程度の重なりを求めているのか？
そもそも積分方程式を差分化して数値積分している時点で誤差が溜まってゆくので
相応の誤差はある。見積もり以上の誤差が出ていると言いたいのか？

中身を全く見ていないが、表面だけを見た感じでは
ｆｌｏａｔは単精度を返すので、ＤＢＬＥにしてみそとか、
『fi=i-1』も単精度に変換されているだろうといえる。

＊昔のコンパイラは倍精度に単精度の数を入れると、あまった桁数にはゴミが入るので
結局数値誤差は単精度と変わりないレベルまで汚染されてしまう。

>>137
ご指摘の点、早速改善させて頂きました。
ありがとうございました。
ずれは、x=1付近ではほとんどありませんがx=0では誤差では済まない程あります。
グラフにすると形もcos関数ではなくなってしまっているのですが、
改善する手立てが見つけられないでいます。

>>137
なんでFORTRAN77のプログラムを1カラム目から書いてるの？
整形めんどい。

>>141
s/^/ /

FORTRANﾈﾀじゃないけど、まいっか

どうも。
質問者がこれを暗黙に要求するのはどうかと思いまして。

viを起動して137-138をコピーペーストして、
:%s/^/ /
を実行した後数字のラベルの前の空白を手動で2カラム削除しました。
スレ違いで済みませんが、
「7カラム目が数字の行の先頭を2カラム削除する」というのもsedかviで
自動で処理出来るのでしょうか？

/^......[0-9]/s/^..//

>>144
ありがとうございます。

前述の計算は以下のように書き直したら期待通りのグラフが描けました。
　　　　　　　program heat equation
　　　　　　　implicit none
　　　　　　　integer i,j,nstep,n,ndim
　　　　　　　parameter (ndim=5001)
　　　　　　　double precision flam,f,h,u,c1,c2,fi,t,tmax,dt
　　　　　　　double precision a,b,c
　　　　　　　dimension a(ndim),b(ndim),c(ndim),u(ndim)
　　　　　　　n=10
　　　　　　　dt=0.005d0
　　　　　　　tmax=0.05d0
　　　　　　　h=1.0d0/dble(n)
　　　　　　　flam=dt/(h**2.0d0)
　　　　　　　c1=1.0d0-2.0d0*flam
　　　　　　　c2=flam
　　　　　　　do 1 i=1,n
　　　　　　　a(i)=c2
　　　　　　　b(i)=c1
　　　　　　　c(i)=c2
　　　　　　1 continue
　　　　　　　do 2 i=1,n+2
　　　　　　　fi=i-2.0d0
　　　　　　　u(i)=f(fi*h)
　　　　　　2 continue
　　　　　　　OPEN(11,FILE='C:\out.dat')
　　　　　　　nstep=0.0

　　　　　　3 continue
　　　　　　　do 4 i=1,n
　　　　　　4 u(i)=b(i)*u(i+1)+a(i)*u(i)+c(i)*u(i+2)
　　　　　　　do 5 i=n+1,1,-1
　　　　　　5 u(i+1)=u(i)
　　　　　　　a(1)=0.0d0
　　　　　　　b(1)=c1
　　　　　　　c(1)=2.0d0*c2
　　　　　　　a(n+2)=2.0d0*c2
　　　　　　　b(n+2)=c1
　　　　　　　c(n+2)=0.0d0
　　　　　　　nstep=nstep+1.0
do j=1,n+1
if(mod(nstep,5).eq.0)then
write(11,*) j*h-1.0d0/n,u(j+1)
end if
end do
write(11,*)' '
if(t+dt.lt.tmax+0.10d0*dt) go to 3
end
function f(x)
implicit none
double precision x,f,pi
pi=3.14159265358979323846264338327950288
f=dcos(pi*x)
end
しかしcos(pi*x)としていた初期関数を
if(x.le.0.5) f=x
if(x.gt.0.5) f=1.0d0-x
に変えてみると誤った計算結果が出ます。
改善点等みつかりましたらご指摘頂けると幸いです。

初めてプログラミングをやってみましたが本の通りに入力してもできません（涙）。
こんな事聞いたら馬鹿にされるかもしれませんが、これができなければ前に進めません…
ソースコードをコンパイルできません…何度やっても「javacはコマンドとして認識できません」とエラーが起きます。
Java2SDKの本を買って付属のCDROMよりインストールしたのですが、いったいどこが悪いのかわかりません。
ソースコードにはミスないはずなのにPCが壊れてるのかな？…すいませんわかる方教えてください

>>148
このスレッドはFORTRANという(恐らくキミが生まれるより以前から存在する)言語の
スレッドであって、Javaのスレッドではありません。適切なスレッドに映ってください。

一般的に、コンピュータは同じ事を繰り返すと同じ結果が得られます。
何も変更しないで同じ事を繰り返して、↓こういう結果になるのは極めて当然、当たり前。
> 「javacはコマンドとして認識できません」とエラーが起きます。

登山道がAからFまであるとする。

さいころX,Yを二つ振って、出た目の合計でその登山道を選ぶかどうか決定する。

出た目の合計が7以上であるならば、その登山道を登ることができる、とする。

それぞれの結果を、Aは**、Bは**というように表示し、

Aは8だから行ける、Bは5だから行けないと表示するプログラムを作りなさい。


こんなプログラムなんだが、まずサイコロをどう表現したらいいかわからない点、
次に、でた結果を全てAは**、Bは**というように表示すんのはいいんだけどさ、

行くか行かないかを判別するのって一個一個ifで分岐させてやらないかんの？

追加で書かせてもらうと、

if （ AからFの結果　＞７）　then
write文で表示する
else
データ破棄、


みたいな感じにしたいんだけど…
プログラム全然できなくて日本語はいりまくりですが…
　

おお！こんな奇跡のスレがあるとは…
まじヘルプお願いします。


　　　　　　　　え-------------か
　　　　　　／
　　　　い
　 ／　　　＼
あ　　　　　　　お-------------き
　　＼
　　　　　　／　
　　　　う　-------------------く


上記のような図がある。
あ〜くはそれぞれ独立した乱数である。

あ〜くで、次の経路の距離は（例：あ−い間はあ・いの積を小数点第３位で四捨五入したもの
）とする。

終点までの全経路の距離の合計を表示し、（例：あ+い+え+か、あ+う+く　etc）
もっとも合計距離が小さいものを選ぶプログラムを作りなさい。


なんか上のほうにも似たようなものたくさんあるけど、応用力がなさすぎて作れない…
助けてください…orz

>>152
Fortranで挑戦、途中で挫折・・・
その後、rubyで作成しました。
あいうえおかきく→abcdefghと置き換えて
tree[]にabc...fghの値を入れると答えを見つけてくれます

#------------------------------
a=0; b=1; c=2; d=3; e=4; f=5; g=6; h=7;
null = nil
$ans = []

# a,b,c,d,e,f,g,h
tree=[1,2,3,4,5,6,7,8]

aa = [b,c]; bb = [e,d]
cc = [d,h]; dd = [g,null]
ee = [f,null]; ff = [f,null]
gg = [g,null]; hh = [h,null]

node = [aa,bb,cc,dd,ee,ff,gg,hh]

def len(tree,node,i,n)
tree[i]*tree[node[i][n]]
end

def path(x)
pair = ["a","b","c","d","e","f","g","h"]
pair[x]
end

>>153

def get_ans(node,tree,i,n,sum,_path_)
_next = node[i][n]
if _next == i then
a = [[sum,_path_]]
$ans = $ans + a
elsif _next == nil
nil
else
x = len(tree,node,i,n)
for j in 0..1
get_ans(node,tree,_next,j,sum+x,_path_+path(_next))
end
end
end

# --- main ---
for i in 0..1
get_ans(node,tree,0,i,0,path(0))
end
p $ans.sort
# --- end main ---

#[[27, "ach"], [38, "abdg"], [42, "abef"], [43, "acdg"]]

>>152
real function distance(path)
common /rval/ rval
real rval(10)
integer path(*), i

distance = 1.0
do i=1,1000
if (path(i) .lt. 1) then
distance = int((distance + 0.5) * 100) / 100.0
return
endif
distance = distance * rval(path(i))
enddo
write(*,*) 'internal error'
stop
end

>>155
program f152
common /rval/ rval
real d(6), minimum
integer i, path(10,6)
data path /
+ 1,2,4,6,0,0,0,0,0,0, ! あいえか
+ 1,3,5,2,4,6,0,0,0,0, ! あうおいえか
+ 1,2,5,7,0,0,0,0,0,0, ! あいおき
+ 1,3,5,7,0,0,0,0,0,0, ! あうおき
+ 1,3,8,0,0,0,0,0,0,0, ! あうく
+ 1,2,5,3,8,0,0,0,0,0 ! あいおうく
+ /
real rval(10)
rval(1) = rand(1) ! seed
do i=1,10
rval(i) = rand(0) * 10
enddo
minimum = 3.40282347E+38
do i=1,6
d(i) = distance(path(1,i))
if(minimum .gt. d(i)) minimum = d(i)
write(*,'(i10,f10.3)') i,d(i)
enddo
write(*,'(a10,f10.3)') 'mininum :',minimum
stop
end

質問させてください．
Windowsのgfortranだとコンパイルも実行もできるソースが，
Linux版のgfortranではコンパイルは通るのに実行時に
　 At line 91 of file LHDmake.f90
　 fortran runtime error;Bad adress
となって止まってしまいます．

LHDmake.f90はMODULEで，その91行目は「WRITE(21,*) ''」という文(データを改行するため)
なのですが，なにがおかしいのでしょうか？

ちなみにUNIT=21は
「OPEN(UNIT=21,file='Result/plotLHD1_polygon.dat',status='unknown')」
です．

コンパイラのバージョンは
Win：gcc version 4.2.0 20060401 (experimental)
Linux：gcc version 4.0.4 20060507 (prerelease) (Debian 4.0.3-3)
です

>>157
こんなのがあった
ttp://gcc.gnu.org/ml/gcc-bugs/2006-01/msg00595.html

「0.360e-10」は0.360×(10の-10乗)ですよね。

eがなくて、「0.360+250」のように表示されるのは
どういう意味でしょうか？

指数が3桁になって表示を切り詰めたんじゃないのかな？

>>159
マニュアルでＦＯＲＭＡＴ文のところを嫁

すいません…どなたか150の回答お願いします…

課題で、

-4(V/x)^2*a^3-2*(V/x)^2*3a^2+λv3a^2+a{2-2*(V/x)^2+4λV}+1+λV=0

V=0.1、ｘは0.1-1.0、
未定乗数λ＝1,2,…t適当に入れて決める
求めたいのがaなのですが、
３次方程式なのでまったくプログラムが解りません
お助けねがえないでしょうか？
フォートランのバージョンはFujitsu Fortran V4.0L10です。
知恵をおかしください。　

>>162
program hoge
integer, parameter :: N=6; real :: X(N), Y(N); integer :: Z(N)
character(len=8) :: f(N) = "行けない"
call random_seed(); call random_number(X); call random_number(Y)
Z = int(X*6 + 1) + int(Y*6 + 1); where (Z >= 7) f = "行ける"
print "(1X,A,A,I2,A,A)",(char(ichar("A")+i-1),"は",Z(i),"だから",f(i),i=1,N)
end program hoge

>>163
安直にはニュートン法でとくべし。

１個は必ず実の解がある。残りはたとえば2次方程式として公式で。

fortran77 でのプログラム作成なんですが、

ｘ=１でのf(x)＝exp(x)の一次微係数を２点微分方と４点微分法を用いて、
hを0.1から0.01ずつ掛けながら数値微分するプログラムを作れ。
hは１０＾（ー１５）まで変化させ各hでの値と厳密値とを比較せよ。

f(1)(x）を一次微係数を表す式としますと

２点微分法は

f(1)(x)=(f(x+h)-f(x))/2h

4点微分法は

f(1)(x)=(8f(x+h)-f(x+2h)-8f(x-h)+f(x-2h))/12h

で表されています。

どなたか分かる方お願いします。

double precision f,dx2,dx4,h,x
integer i
f(x) = exp(x)
dx2(h,x) = (f(x+h)-f(x)) / (2*h)
dx4(h,x) = (8*f(x+h)-f(x+2*h)-8*f(x-h)+f(x-2*h)) / (12*h)
x = 1.0
write(*,*) 'exp(1) = ', exp(1.0D0)
do 10 i=1, 15, 2
h = 0.1D0 * 10.0D0**(1-i)
write(*,*) h, dx2(h,x), dx4(h,x)
10 continue
end

>>165
ニュートン法は使ってはいけない
ということなので、
他の方法はないのでしょうか？

>>168
つ ニュートンラプソン法

三次方程式なら普通に代数的解法があるんだからそれを使えばいいんじゃね？

　　　REAL*8 FM(6),X(6),H,T,TI,TF
　　　WRITE(*,201)
　201 FORMAT(' ',5X,'T',9X,'X',11X,'Y',11X,'U',11X,'V')
　　 OPEN(8,FILE='b:alinda.dat')
CALL ICOND(N,TI,TF,H,X)
DO 20 T=TI,TF,H
WRITE(*,'('' '',F8.1,4F12.7)')T,(X(I),I=1,N)
WRITE(8,'(2F10.5)')X(1),X(2)
CALL RUNGE(N,H,T,X,FM)
DO 30 I=1,N
30 X(I) = X(I) + FM(I)*H
20 CONTINUE
CLOSE(8)
END


上記のプログラムで、理解できないところがありまして。
＞WRITE(*,'('' '',F8.1,4F12.7)')T,(X(I),I=1,N)
ここの括弧内なのですが、[4F12.7]は4回実数を繰り返し表示する、という意味であっているのでしょうか?
あと、Tが何のためにあるのかがわかりません。

退職金1000万を源資とし、毎年100万ずつ引き出して使用した場合
預金が無くなるまでの毎年の元利合計を計算するプログラムを作れ
なお、年利率が0.1%、1%、5%の3種類について計算を実行して
預金が無くなるまでの年数を比較しなさい
(預け入れ時点では引き出さず、1年後から引き出すものとする)

という課題が出たのですが、分らなくて困っています。
教えて下さい。

program report11_2
implicit none
integer :: a(3),n,i,j,b
real :: r(3)

read(*,*) a(1)
do i=2,3
a(i)=a(1)
end do
do i=1,3
read(*,*) r(i)
end do
read(*,*) n

b=100.0
do i=2,n
do j=1,3
a(j)=a(j)+a(j)*r(j)/100.0-b
end do
write(*,*) i,(int(a(j)),j=1,3)
end do

stop
end program report11_2

ここまでやってみました。

n原子分子を構成する各原子の質量Miと座標（Xi、Yi、Zi）を読み込み、
質量中心を計算するプログラムを作れ。

っていう問題なんですが、どのようにすればいいでしょうか？
If文とDo文と配列を使わないといけないのですが･･･。
よろしくお願いします。

三つの数値データを与えてその三つの数値の最大値、最小値、中央値、
平均値、標準偏差、不偏分散を求められるプログラムを作成せよ。
という問題が出されてわからなくて大変です。
お願いします

超素人な俺が考えてみた
　PROGRAM DISTS
C
C
C
real x,y,z,max,min,mid,hyoujun,heikin,wa
print *,'enter x,y,z
read x,y,z
wa=x+y+z
***heikinti***
heikin=wa/3
print *,'heikin=wa/3
***max***
if (x .gt. y) .and. (x .gt. z)
print *,'max=x
elseif (y .gt. x) .and. (y .gt. z)
print *,'max=y
elseif (z .gt. y) .and. (z .gt. x)
print *,'max=z
endif
***min***
if (x .gt. z) .and. (y .gt. z)
print *,'min=z
elseif (y .gt. x) .and. (z .gt. x)
print *,'max=x
elseif (z .gt. y) .and. (x .gt. y)
print *,'max=y
endif
***mid***

>>175
program hage
real :: a(3),b,c(1),d(3)=3;read *,a
c=maxloc(a);d(int(c(1)))=0;print *,maxval(a)
c=minloc(a);d(int(c(1)))=0;print *,minval(a)
print *,pack(a,d>1);b=sum(a)/3;print *,b
print *,sqrt(sum((a-b)**2)/3);print *,sum((a-b)**2)/2
end program hage

177を初心者向けのプログラムに出来ませんか？

>>178
十分初心者向け

>>178 >>177は筋が悪い、学んでもかえって馬鹿になるからやめとけ。

fortran90の課題なのですが、分からなくて困っています。
どなたか優しい方、よろしくお願いします。

数字1.2.･･･.nを様々な順序に並べたものを置換という。
例えば、1.2.3の置換は次の6通りである。
（1.2.3）（1.3.2）（2.1.3）（2.3.1）（3.1.2）（3.2.1）
1.2.･･･.nの置換はn!通り存在するが、これらを全て重複なくリストアップする
プログラムを作れ。

というものです。

追加です。

nはデータとして読み込むものとする。

ってのを忘れてました。
よろしくお願いします。

順列 再帰 でググるといいかも

三つの数値データを与えてその三つの数値の最大値、最小値、中央値、
平均値、標準偏差、不偏分散を求められるプログラムを作成せよ。
という問題が出されてわからなくて大変です。
誰か教えてください。

> 数字1.2.･･･.nを様々な順序に並べたものを置換という。
なぜ教育課程でも使われている順列という定義があるのに、独自の定義をするのだろう?

どなたかお願いしますorz

184頼みます

>>184
UZEEE
一般性を考慮しないＦ９０の糞プログラム例だ。
PROGRAM vip
IMPLICIT NONE
REAL :: x(3), xmin, xmax, xmed, xmean, sigma, usigma
x = (/10.0, 1.0, 5.0/)
xmin = MINVAL(x)
xmax = MAXVAL(x)
xmed = SUM(x) - xmin - xmax
xmean = SUM(x) / SIZE(x)
sigma = SUM( (x - xmean)**2 ) / SIZE(x)
usigma = SUM( (x - xmean)**2 ) / ( SIZE(x) - 1 )
WRITE(*, *) xmin, xmax, xmed
WRITE(*, *) xmean, sigma, usigma
STOP
END PROGRAM vip

>>185
有限群論とかでは置換群とか呼ぶから、そっち方面の方言だべ

>>172
subroutine calc(x)
real rate
integer balance, years, interest

balance = 10000000
years = 0
write(*,*) 'rate = ', x, '%'
rate = x / 100.0

10 continue
interest = int(balance * rate + 0.5)
if (interest .ge. 1000000) then
write(*,*) 'rate too much'
stop
endif

balance = balance + interest - 1000000
years = years + 1
if (balance .le. 0) return
write(*,*) years, balance
goto 10
end

program report172
call calc(0.1)
call calc(1.0)
call calc(5.0)
stop
end

>>186
超適当。0.1%, 1%の時は１０年後に残金約５０万円。５％の時は１４年後に残金２０万円。


PROGRAM no_dream
IMPLICIT NONE
REAL, PARAMETER :: rinterest(3) = (/0.1, 1.0, 5.0/), rstart = 1.0e7
REAL :: deposit
INTEGER :: i, j
DO i = 1, 3
WRITE(*, *) 'Case', i, 'Interest =', rinterest(i), '%'
WRITE(*, *) '----------------------------------------------'
deposit = rstart
j = 1
DO
deposit = deposit * (1.0 + rinterest(i) / 100.0) - 1.0e6
WRITE(*, *) 'Year', j, ' Remaining Money ', deposit
j = j + 1
IF (deposit < 1.0e6) EXIT
END DO
WRITE(*, *) '=============================================='
END DO
STOP
END PROGRAM no_dream

x = (/10.0, 1.0, 5.0/)　
この文章の意味を教えてください

n原子分子を構成する各原子の質量Miと座標（Xi、Yi、Zi）を読み込み、
質量中心を計算するプログラムを作れ。

という問題なんですが、どうしていいか分かりません。
If文とDo文と配列を使わないといけないのですが･･･。
よろしくお願いします。

どなたか171お願いします

>>171
>[4F12.7]は4回実数を繰り返し表示する、という意味であっているのでしょうか?
はい、幅12小数部７桁の実数を４回表示する、という意味です。

>あと、Tが何のためにあるのかがわかりません。
何ででしょうね。 時間も表示させたかったんじゃないですかね。

>>194
4F12.7の4は反復数、Tは出力するREAL*8型の変数。

>>192
配列構成子を使って配列に値を代入する。

>>192
x(1)=10.0
x(2)=1.0
x(3)=5.0
と同じこと。

>>193
解法とプログラムのどちらが解らないのかはっきりしたまえ。

>>199
プログラムが分かりません。
よろしくお願いします。

>>200
我輩が知っている解法ではdo文と配列は使用するが、
if文は使用しないのだ。
その、解法を提示しろ。

>>201
i<=2までの場合と2<iの場合で質量による
重心のバランスの出し方って変わりませんかね？
すいません。
ちょっと間違ってるかもしれません。

198ってx(1)が10、x(2)が1,x(3)が5、という定数だという事を示しているんですよね？
定数を指定してしまうと最初から各値が決まってる気がするんですけど・・

>>203
変数x(1)に10を格納する。
変数なので格納される値は書き換える事ができる。

204>>なんとなくわかりました。
もしx(1)に１０を格納しなかったらプログラムとして不成立？

>>202
重心の出し方は n<=2 だろうが 2<n だろうが変わらないはず。
(M1×y1 + M2×y2 + ... + Mk×yk)
-------------------------------- (y座標)
k

x座標もz座標もyをx、yに変更すればO.K.

おそらくif文を使うのは、何個の原子があるかを判定するところで使うのでは？
つまり、プログラムに原子の個数をあらかじめ与えておくのではなく、与えられた初期データを読み込むときに原子が何個あるかをプログラムがカウントするということね。

あるプログラムを作っているのですが、
配列に入れる値をテキストファイルから読み込むことが何度もあります。
その部分をサブルーチン化したいのですが、うまくいきません。
次のようなエラーが出るのですが、どう改善すればようでしょうか？

変数'x'の宣言において,整数式であるべきところが整数式ではありません.
'UNIT'指定子の値または型が正しくありません.
'FILE'指定子の値または型が正しくありません.
'UNIT'指定子の値または型が正しくありません.
'UNIT'指定子の値または型が正しくありません.

subroutine data(x,n,unit,file)

real,intent(inout),dimension(n) :: x
real,intent(in) :: n, unit, file

open(unit,file)
do i = 1,n
read(unit,*) x(i)
end do
close(unit)

end subroutine

xは値を入れたい1次元配列
nは配列の要素数
unitはファイルに付ける番号
fileは読み込むテキストファイルの名前です。

↑の形にこだわらなくてもいいので、
同じことができる方法があれば教えてください。

>>207
fortran90は使ったことないが、
> real,intent(in) :: n, unit, file
realじゃいやーん
といってるんだろ。

>>208
漏れは>>207は釣かとオモタよ。
本気で聞いてるなら、しょん便で顔を洗って出直してくるべき。

y(x)=F(x)*cos(m*x)の形をした関数を数値積分(積分範囲0〜2pi)
をシンプソン則で計算させようと思っているんですが
試しに、F(x)=1、x軸上のサンプリングポイント2000点とし、自然数mを1から1000まで連続的に増やして計算させたところ、
あるm=650近傍で積分値が大幅に振動し、それ以降では振動がある程度収まるという結果になりました。
この振動を抑えるための解決策があればご教授お願いします。

>>211
三角関数の振動の周期と積分のステップ幅の関係を考えねば。
たぶん２π/650〜Δ1/1000くらいなんだろ。
積分範囲が分からんので断言できないが・

>>209
ありがとうございます。
これ書いたときはだいぶ疲れてて気がつきませんでした。

そもそも文字列（ファイル名）を引数にしようとしていました。

>>213
>そもそも文字列（ファイル名）を引数にしようとしていました。
それは可能姉妹

某サイトからの引用ですが

CHARACTER::f*63,q=""""
f="('CHARACTER::f*63,q=',4A/'f=',3A/'PRINT f,q,q,q,q,q,f,q'/'END')"
PRINT f,q,q,q,q,q,f,q
END

を，まわすとプログラムソースと同じ答えが返ってきます．
そこの解説読んでも何故そうなるのかわかりません．
誰か詳しく解説してください．お願いします

CHARACTER(LEN=43) :: f="('hoge',A,A,A,A/'huge',A,A,A/'hage'/'hige')"
WRITE(*,FMT=f) '"', '"', '"', '"', '"', f, '"'
END

>>215
１）ＦＯＲＭＡＴは文字列として与えても良い。　ｆはＦＯＲＭＡＴ文としてと文字列としてと
　　２回違う意味で使用されいる。

２）ＰＲＩＮＴ ｆｍｔ、ｔｅｘｔ ＝　ＷＲＩＴＥ（＊、ｆｍｔ）　ｔｅｘｔ
　　PRINT文の１個目の引数？はフォーマットを表している。

３）文字変数　ｑ　＝　”　（引用符中の２連発の同じ引用符は１個の文字扱いとなる。）

４）PRINT f,q,q,q,q,q,f,q　の意味は、以下と同じ。
１００　FORMAT('CHARACTER::f*63,q=',4A/'f=',3A/'PRINT f,q,q,q,q,q,f,q'/'END')
　　　　WRITE(＊、１００）　ｑ、ｑ、ｑ、ｑ、ｑ、ｆ、ｑ　　

５）FORMAT文の斜線は改行。　
　　文字列は文字列としてそのまま出力。
　　Aは文字変数ｑに対応する文字型指定子。

f2cで変換したソースを読む気がしない

exe にしたファイルを元に戻したいのですが何か良いツールはないでしょうか？
20行くらいの簡単なファイルなのですが数がめちゃくちゃあって一から入力する気になれないです
どなたかそんなツール知っている方いましたらレスお願いします

>>219
つ 「逆アセンブラ」「逆コンパイラ」

まあ、こいつらの出力を見たら一から入力する気になると思うがな。

>>216
>>217
なるほど．．．fortran奥が深いですね

>>220

それをやってくれるツールみたいなのはないのでしょうか？
フリー・シェアのどちらでも構わないのですが・・・

>>222
Googleで「逆アセンブラ」とか「逆コンパイラ」で検索すれば色々見つかるよ。
Fortranの逆コンパイラは知らんが。

但し、実行ファイルから完全に元のソースにもどしてくれるようなものを期待しているのなら
そんなのもは無い。

　来週の月曜までが期限の課題なのですが全く分からなくて困っています。
　内容としては、以下のプログラムでは行列の対角成分に０があると逆行列を
計算できないので、それをできるように改良するというものです。
　行列の対角成分の絶対値が大きくなるように行を入れ替えるピボット選択を
すればいいみたいなのですが…
　２ｃｈに書き込むのは初めてなのでおかしな点があったら申し訳ありません
が、何かヒントをいただければ助かります。

>>224のプログラムです。
　　　PROGRAM INVERSE MATRIX

　　　INTEGER I,J,K
REAL A(3,6)

PRINT *,'Input matrix A'
DO 1000 I=1,3
READ *,(A(I,J),J=1,3)
DO 2000 J=4,6
IF (I+3 .EQ. J) THEN
A(I,J)=1.0
ELSE
A(I,J)=0.0
ENDIF
2000 CONTINUE
1000 CONTINUE

>>225の続きです
　　　DO 3000 K=1,3
DO 4000 J=K+1,K+3
A(K,J)=A(K,J)/A(K,K)
DO 5000 I=1,3
IF (I .NE. K) THEN
A(I,J)=A(I,J)-A(K,J)*A(I,K)
ENDIF
5000 CONTINUE
4000 CONTINUE
3000 CONTINUE

PRINT *,'Result'
DO 6000 I=1,3
PRINT *,(A(I,J),J=4,6)
6000 CONTINUE

END

>>223

ない

ん？ なにが？

文脈から判断すると、fortran逆コンパイラ。

X(O)=1としたとき、K=1,2,3,...の順に
X(K)=COS(X(K-1))
をくりかえして、
|X(K)-X(K-1)|<0.00001
になったら反復を打ち切るプログラムを作れ。

という問題なんですけど、どうすればいいでしょうか？
DO WHILEを使うのは分かるんですが…
どなたかお知恵をお願いします。
FORTRAN90です。

もう少し歯ごたえのある問題をもってこい！

>>230
program test
real :: x(0:1000)
integer :: i,k
k=0
x(k)=1.
do
k=k+1
x(k)=cos(x(k-1))
if(abs(x(k)-x(k-1))<0.00001) exit
enddo
do i=0,k
print *,i,x(i)
enddo
end program test

ご指導ください．
関数副プログラムをF1,F2として，
A＝F1*F2
としたとき，
write(*,*) A,F1*F2
の結果が

NaN,[実数値]

となるのはなぜですか？？？？？
謎です・・・

>>233
ソースがないとなんとも言えないが、
変数の型は大丈夫かな？

>>233
NaNとはNot A Numberの略で、実数として違法な表現が入っている場合などに起きる。

コードをもう少し広く見ないと分からない。
>234の言うように型を間違えている可能性が高い。

>234,235
ご回答ありがとうございます．
そうですか・・・
関数も倍精度，代入している変数も倍精度で定義しています．
なぜ代入してから表示するのと，代入する前の状態で表示するのと
で変わるのでしょうか？

関数は内部関数にしているのですが，なにか関係があるのでしょうか？
コードを貼り付けると非常に長くなりますので，
考えうる可能性としてはどのようなものがあるのでしょうか？

A=F1*F2からwrite(*,*) A,F1*F2までの間で何かAに副作用を及ぼすようなことをしていない？

>>232
ありがとうございます。できました。
DO WHILE 使わなくても出来ますね。全然気がつきませんでした。

>>236
出してる情報が少なすぎるから、エスパー以外答えられんよ。
せめてWrite文の結果とエラーメッセージくらい出せば少しはわかる。
コンパイラの種類もわからないし。


まぁ漏れのESPカードによると、F1*F2は巨大数でオーバーフローして倍精度では表せない
数になった。もしくはアンダーフロー。WRITE文が４倍精度に対応しているのでエラーが出なかった。

最近４倍精度をサポートしている処理系が増えているので、Intel・Lahey（Fujitsu）・PGI・Absoft、ｅｔｃ、
十分ありうる。

とあるプログラムのformat文の中に「1pg」と書いてあったのですが，
これは何を意味しているのでしょうか

>>240
１P　G　のこと。　

それはFとEを自動で切り替えるGフォーマット　に　小数点の位置をずらすP　フォーマット
が組み合わさったもの。

F90時代にはESあたりを使うがよろしい。
プログラム書法では推奨されていたが、もはや時代遅れ。

詳しくはマニャル嫁。

>>241
はい！せんせい

すいませんちょっとお聞きしたいのですが
ちょっと計算機を変えて計算しはじめたのですが
以下のように今まで読めてたopen文が読めなくなってしまいます
これってコンパイラとか計算機によって依存するんでしょうか
なにかヒントいただければと思うのですが・・・
getRegFromUnwindContext: Can't get Gr0 from UnwindContext, using 0
forrtl: severe (59): list-directed I/O syntax error, unit 515, file /・・・・・

すいませんあげさせてください

ついでにエラー部分は
open(515,file='bc2.dat')
cwrite(*,*)'test'

cwrite(*,*)'test2'
read(515,*)iob　　←ここで止まってる
cwrite(*,*)'test3'
read(515,*)(iobb(i),i=1,iob)
close(515)
です

すいません自己解決しました
読み込むファイルのスペースの数が悪いようでした
書き換えます・・

やっぱりめんどくさいです
どうやら改行もだめのようです
こんなことってあるんですか？

Makefileについて質問があります．

test.f90がメインルーチンで，サブルーチンの入ったsub1.f90に依存しているので
次のようなmakefileを作りました．

第2，3行目を追加したら常にclean:が実行されるようになってしまったのですが
何ででしょうか．．
環境はWindows XP SP2 pro，Intel Visual Fortran 9.1で，
MSのNMAKEを使ってます．


FLAGS=%F90FLAGS% %LINK_F90%

.PHONY: all
all: main.exe

main.exe: main.obj sub1.obj
%F90% $(FLAGS) main.obj sub1.obj /o main.exe

main.obj: main.f90
%F90% $(FLAGS) main.f90 /c

sub1.obj: sub1.f90
%F90% $(FLAGS) sub1.f90 /c

.PHONY: clean
clean:
del *.obj

すみません；
× test.f90がメインルーチンで，サブルーチンの入ったsub1.f90に依存しているので
○ main.f90がメインルーチンで，サブルーチンの入ったsub1.f90に依存しているので

あとコマンド行の前にはもちろんタブが入ってます．

>>248
nmakeが .PHONY に対応しておらず、通常のターゲットとみなしてしまうから。

.PHONY: all
を消すか、 all: の後に持って行くかすれば良い。
（nmakeで使う分には残しておいても意味はないが…）

make は結構方言がきついので、油断していると変なところではまってしまうぞ。

>>247
何を言っているのかよく分からんが、一番確実なのは書いたとおりに読み込むことだ。

１行を頭から途中まで読んで、尻のほうを捨てることも出来る。
その場合、次のREAD文では、ファイル内の次の行を頭から読むことになる。

途中まで読んで、その行の残りを別のREAD文で読みたければ、
F90以降で導入されたADVANCE＝”NO”を指定する必要がある。

詳しくはマニャル嫁

>>250
ありがとうございます！
NMAKEの仕様を調べなおします．

ｄ

まず、最初に規定の値（xとここではする）を決める。
そのxに乱数によって導き出された数y（計算方法は下記に記す）を引いていく。
これをz回（zはあらかじめ決めておく）繰り返す。
そのとき、計算過程がマイナス、または0になったとき、「miss」という項目に+1し、
xの値を初期値に戻す。というものです。

たぶんわかりにくいと思うので具体例を挙げて説明します。

（例） 各々の初期値を
x=7
y=int(rand()*10)　※この場合計算結果が3だったとします。
z=10
miss=0
とします。

このとき、x-y（以降、仕事と呼ぶ）をz回繰り返します。

�@7-3=4
�A4-3=1

1-3=miss
miss=0+1
ここで、xの値を7に戻す。

�B7-3=4
�C4-3=1
・
・
という感じで、missになる仕事以外の回数をz回繰り返します。
missになる仕事のとき、xを初期値に戻すためにどういう式を使ったらいいかがわかりません。
よろしくお願いします。

>>254
>missになる仕事のとき、xを初期値に戻すためにどういう式を使ったらいいかがわかりません。
繰り返し処理に入る前に x の初期値を適当な変数に保存しておく。
missに +1するところで保存しておいた値を x に代入する。

>>254
その日本語の解説どおりにFortranで書けばOKだ。
ほとんど完成している。

「miss」か「Z=今までの仕事回数+1」

の判断はif文で分岐させればいいんですよね？

そのためにあるのがIF文です。

副作用のないＬｉｓｐで、変数の型に制約がつけられる、という理解でいいか？

参照透明性とはチューリング機械で言うと内部状態を持たないということか？

今自宅なので、実際にプログラムを起動させて確認を取ることができないので、
構文だけここに記入します。間違いなどを指摘してくだされば幸いです。

※xfirstはxの初期値とする。

program main

integer x,z,xfirst,miss
real y

x=7
xfirst=7
y=int(rand()*10)
z=10
miss=0

do 20 z=1,10

x-y=x
if(x .gt. 0) then
z=z+1

else

miss=miss+1
x=xfirst

end if

20 continue

FORTRAN77です
ビルドしたら<静的変数の合計の大きさまたは共通ブロックの大きさが、許容量を超えました>
というエラーがでたのですが、プログラムを分割するしかないのでしょうか？

>>261
OSの制限なら分割するしかない。

32BitOSだと１配列やコモンブロックあたり２Gバイトまでというのはよくある。
（符号付３２ビット整数の正の最大値＝２G）

そうではなくて、コンパイラのデフォールトの制限でエラーが出ることもある。
その場合はマニュアルを見て、リンカーのオプションを変えればOSの制限値までは増やせるだろう。

>>262
アドヴァイスありがとうございます

>>260
文法のチェックだけでよければ
http://www.lahey.com/check.htm
ここでWEB上でチェックできる。


内容は見ていないが、ざっと見た範囲では
x-y=x
これがおかしい。あと末尾に　END　が必要。

>>264


アドバイスありがとうございます。

x-y=x
というのは、x-yという演算を行い、そのときに出た結果を表現したつもりです。

例文でいうところの、7-3＝4、4-3＝1…と、4や1にあたる表現をしたいのですが、
なかなかうまいこといかないんです…
このような計算を行う場合、たとえばどのように計算してやればいいですかね？

>>265
x = x - y
向きが逆

>>265
266で言われてますが、=は左の変数に
右の計算結果を代入するものです。

>>265
初心者の頃の俺と同じことやってるｗ
= は変数の値を左右が釣り合うようにするワケじゃなくて
= の右側の結果を左側の変数に入れるものだよ。

「←」の代わりだと思えば良いかな。

do i=1,10
うんこ
if (i.eq.1)then
なんちゃら
else if(i.eq.2)then
なんちゃらちゃら
　endif
enddo
なんか一気に何ちゃらちゃらまでまわっちゃうんです＞＜
うんこ→なんちゃら→うんこ→なんちゃらちゃら　の順番になるはずですよね？
というかなんちゃらでエラーが出てます＞＜これが原因でしょうか？
このエラーなおしたらちゃんとまわるようになるんでしょうか？まずはえらー直すことが先だと思うんですけど
構造的にはあってますよね？？

>>264であるサイトいったんですけど、errorってでてダウンロードできなかった…orz
構文はあってますかね？

>>270
ダウンロードするのではなく、自分の書いたソースをアップロードするのだ。

で、>>260をそのままチェックした結果：

Diagnostic messages: program name(main)
1026-S: "SOURCE.F90", line 1: Corresponding END statement missing.
2004-I: "SOURCE.F90", line 8: 'y' is set but never used.
1302-S: "SOURCE.F90", line 14: Not a valid Fortran statement.
1635-S: "SOURCE.F90", line 16, column 7: DO variable 'z' cannot be redefined in range of DO construct.

上から
・ENDがない
・8行目（y=int(rand()*10)）で y に値が代入されているが使われていない（これは次のエラーと関連している）
・14行目（x-y=x）の文が変
・16行目（z=z+1）でzの値を変えようとしているがｚはDO変数なのでDOループ中で変更できない。

構文的におかしいのはこんなところ。
>>254の題意からすると、zの使われ方が少しおかしい。

>>271

zの使い方がおかしいとありますがどのようにおかしいのですか？
z=10とdo文の前で定義しているのがおかしいということですか？

do 20 z=1,10
で、仕事回数zの表現をしたかったんですが…

>>272
もしかして日本語を理解できない人？
> 16行目（z=z+1）でzの値を変えようとしているがｚはDO変数なのでDOループ中で変更できない。
DOループの中でそのDOループのiteratorであるzをz=z+1という文で変更しようとしているが、
それはコンパイラとしては許せないことであると言われている。

では、このような表現ではだめでしょうか？
work(10)は10回までの仕事を表現します。
work(z)というように表現し、do変数zで表現したいと思います。
このプログラムの最終目的は、missの回数をエクセルに出力して、グラフを作成することなんです。
そのために、結果を出力したいとおもっています。

以下に構文を示します。

program main
integer x,z,xfirst,miss,work(10)
real y

x=7
xfirst=7
y=int(rand()*10)
miss=0

do 20 z=1,10
x-y=x

if(x .gt. 0) then
write work(z)　　　→※何回目の仕事かどうか、演算結果を表示したい。
else
miss=miss+1
x=xfirst
end if
end

20 continue

あーダメだ x-y=xを見て諦めた
宿題スレじゃないしこっちは考えるヒントだけ出して自分で完成することで
プログラミングの達成感が味わえれば少しはプログラミングが面白くなるかもとか一瞬思ったが
散々言われてたx-y=xを直す気もなさそうだし相手にするだけ無駄っぽい

>>275

あ、すいません…

間違ったソースコードをコピペしてしまったんで…
いまエラーを一個ずつ直してますが、確かに直っていくと面白いですね。

理解力がなくてすいません。これが最後の質問です。

program main

integer x,z,xfirst,miss,work(10)
real y

x=7
xfirst=7
y=int(rand()*10)
miss=0

do 20 z=1,10

x=x-y
if(x .gt. 0) then

write(*,*) work(z)

else

miss=miss+1
x=xfirst


end if


20 continue


end

上記のような構文を作ってコンパイルしたら、error 0 warning 0となっているのに、
題意のような結果がうまく表示されません。

おそらく、>>271で指摘されている
・8行目で y に値が代入されているが使われていない　という点が原因だと思います。

しかし、ここがイマイチよく理解できません。

とありますが、yは、y=int(rand()*10で値が決定されていますよね？
ここで乱数計算によって、たとえばy=3という結果が出た場合、以後のプログラム内における
yという値は全て3というように決定されるんじゃないんですか？

x=x-yという式によって、ｙという変数はプログラム内で使われているのではないんですか？
この点がイマイチよく理解できないんです。

ドシロートですいません。

>>277
>・8行目（y=int(rand()*10)）で y に値が代入されているが使われていない
というのは
x-y=x とやってた為にコンパイラが

「x-y=x」は不正な文だからとりあえず無視
　→ 他に y を参照している箇所がない
　→ 「yが使われていない」と警告

として出てきたものだ。x-y=xを正しく書き直せば出なくなし、結果がおかしいのとは関係ない。

y=3でテストしたければ「y=int(rand()*10」の後ろに「y=3」と書いて走らせてみたらいい。

ところで、work(z)にはどんな値が入っているかわかるかな？

>>279

仕事回数z、つまりそのときの仕事の結果を表示していくつもりで書いたんです。
ｚが3回目の仕事のときに、3回目の仕事の結果を出す、というようにしたんですが、
間違ってますかね？

workに何も代入してないじゃん。

>>281

つまり、work(1)〜(10)までを定義してやる必要があるということですよね？
(1)〜(10)を「どのように求めるつもりなのか」を明記するってことですか？

この場合だと、x-yでいいのかな？

彼が何をしたいかわからない俺は馬鹿

>>283

バカですんません。
とりあえず、やりたいことというのは、

�@定められた仕事回数の分の仕事（x-y）をdo loopで行う。
�Ax-yを実行する度に、仕事回数ｚの数を増やす

※x-yがマイナスに達した時、xの値を初期設定の値に戻す。仕事回数zに+1する代わりに、
その後、missの項目に+1する。

�B仕事回数が規定に達したらそこでプログラム終了。

何故、こんな意味のわからないプログラムを組むのかというと高校の卒業課題で使うためです。
これによってある計算をしなければならないので…（詳細はここでは書けませんが）

とりあえず、x-yの計算結果と、missの回数すら出力できてないので、僕は非常に混乱しています…
バカでほんとうにすいませんorz

>>284
卒業課題って今年度の？
今の状態で行き当たりばったりにプログラム書くより、教科書なり参考書なりをよく読んで
もう一度基本的な文法を覚えなおした方がいいよ。

それからパソコン持ってるなら家でも予習復習できるように処理系の１つも入れておけ。
Windowsで学校の課題程度ならこれくらいで十分
ttp://www.coastal-env.k.u-tokyo.ac.jp/koibuchi/fortran.htm
他は本スレ参照のこと
http://pc10.2ch.net/test/read.cgi/tech/1163319215/2

構造格子上の3次元の温度データから等数値面をつくり、
その面に沿って面積分を行いたいのですが、うまくいきません。
　そもそも曲面の求め方が悪いのか（勾配から法線ベクトルを出し格子内
で面を張っているだけ）チェック用に完全な球を作り、
面積分を行っても球の表面積と桁が違うような結果になってしまいます。
　曲面の作成と面積分のやり方、あるいはfortranのサンプルプログラム等
に関してご存知でしたら教えてもらえませんか？

倍精度複素数型 cKappa には，
純虚数 か 実数 かの何れかが格納されています．

つまり，cKappaを2乗すると必ず実数になるのですが，
この値を倍精度実数型 dE に代入したいときは，

dE = cKappa * cKappa

で問題ないのでしょうか？

初心者でもぅさっぱり分かりません。
VS2005を入れてSilverfrostのftn95をインストールしてVSでftn95を選択。
一番簡単な
PROGRAM MainProgram
write(*,*) 'TEST'
END
を打ちデバッガ実行するとLINKエラーになります。このLINKエラーって何でしょうか？

288です。失礼
------ Build started: Project: FortranApplication2, Configuration: Debug .NET ------
Building project FortranApplication2...
Updating References...
Linking...
LINK : error - File not found: C:\Documents and Settings\.dbk
こんなエラーです。

>>288
VS2002だが、うまく行っている。
エラーメッセージを見る限り、ファイルそのものが見つかっていないような。
まず単独でビルドしてみたら？


Building project FortranApplication1...
Updating References...
Compiling file: FreeFormat1.f95
Linking...
Assembly file Debug\NET\FortranApplication1.exe for .NET 1.1 produced
Build log saved at: "file://C:\Documents and Settings\All Users\Documents\FORTRAN\FortranApplication1\Debug\NET\buildlog.txt"
FortranApplication1 build succeeded.

>>まず単独でビルドしてみたら？
時間がなく、余り調べていないのですが、自分はVSが初心者です。
コマンドプロンプトからはexeを生成し、実行まで出来るのですが、VSだと･･･orz
上記の「単独でビルド」とは何の作業をさすのでしょうか？

>>291
エラーメッセージから見ると、ソースファイルをプロジェクトに登録していないように見えるんだが。
New で作ったならそういう事は起きないが既存ソースなら登録がうまく言ってないのかも？

s = a**b

の計算が必要なときは

s = exp(b*log(a))

と書いた方が計算が速いと聞いたのですが，本当でしょうか？

>>293
自分でベンチマーク作って確かめろｗ

どっちにしろ一般のべき乗はexp log で計算しているはずなので、
８０８７のようにハードウェア的に実装されたｅｘｐを利用したかしないかの違いだろう。

>>294
うちのPCではどうやら後者の方が2割程度早いようです．
ｺﾝﾊﾟｲﾗはIntel Visual Fortran 9.1です．

それでこれからは全てexp logで計算しようと思ったのですが，
問題が発生しました．

次のような，モジュールで変数を共有する仕様の
プログラムを書いたのですが，コメントしているように，
共有用モジュールの変数を介して値を渡した場合と
直接値を代入した場合で結果が異なってしまいます．

ちなみにb**bで計算した場合はどちらも同じ結果を出力します．

b=0の場合を除いて演算してもいいのですが，
NaNになる原因が分からないので何とも怖くて使えません．

何かアドバイスあったらお願いします．．

module commons
implicit none
real(8) a(1:10)
end module commons

program main
use commons
implicit none
a(1)=0.d0
call sub
end program main

subroutine sub
use commons
implicit none
real(8) b

!こちらを使うと出力は NaN
b=a(1)

!こちらを使うと出力は 1
b=0.d0

write(6,*) exp(b*log(b))
return
end subroutine sub

>>295
組み込み関数 LOG(X) はXが実数の場合、Xは正(X>0)でなければならない。

IVFでは -∞ を返しているが、b * log(b) = 0 * -∞ = NaN となるので、
exp(b * log(b)) = exp(NaN) = NaN となる。

b=0.d0 の時に 1 と出力されるのは、コンパイル時に bの値が既知の定数(=0)なので
最適化によりプログラム中で b*log(b) が計算されず、直接 exp0.0D0) を求めてるため。
最適化なし(-Od)でコンパイルすると両方とも NaN になる。

>>297
おおー！確かにそうなります！

実にすっきりしました．
どうもありがとうございます．

途中で書き込んでしまいました．

あまりに訳が分からないので
さっきVinu Linuxをわざわざインストールしてgfortranでも
同じことをしたんですが，こちらは b=0.d0 と書いても b=a(1) で渡しても
NaNとなりました．

今まで最適化で問題が発生したことが無かったこともあって，
原因がコンパイラの最適化にあったとは全然気がつきませんでした．

Xが0以下の実数の場合、LOG(X) の結果は処理系依存になるので
実行時エラーを出してプログラムを止める処理系もある。

最近の処理系はIEEE準拠になって途中で実行時エラーで止まることが少なくなってありがたい。

せっかくIF文で例外を排除しているのに、コンパイラが最適化して勝手に投機的実行をして、
挙句に０割とかで実行時エラーを起こしてとまるとか、ウンコのようなことが昔はよくあったｗｗｗｗ

学校で今週までの課題で
『ｎ元連立一次方程式をガウス・ザイデル法を利用して解くプログラムを作成し、テストを行え』
という課題がだされたのですが、全くわかりません(;_;)
頑張って本を読んでみたのですが理解できなくて(;_;)
最終手段として、ここに来ました(>_<)
どなたか、ご教授下されば幸いです(;O;)
初めての書き込みなので粗相があったらすみませんm(__)m
ちなみに学校で使っているのはフォートラン77です。

またガウス・ザイデル法かｗ

こんなの教科書にも載りまくりだろう。

誰かFAQにでも入れておけ。

>>302
IMPLICIT REAL*8(A-H, O-Z)
PARAMETER (NMAX=99, LMAX=100, EPS=1.D-10)
REAL*8 A(NMAX,NMAX), B(NMAX), X(NMAX)
OPEN(10, FILE='DATA.TXT', STATUS='OLD')
READ(10,*) N
DO 100 I = 1,N
READ(10,*) (A(I,J), J=1,N), B(I)
100 X(I) = 1.
DO 110 I = 1,N
110 WRITE(*,'(1X,100F12.6)') (A(I,J), J=1,N), B(I)
CLOSE(11)
DO 220 L = 1,LMAX
ERR = 0.
DO 210 I = 1,N
TX=B(I)
DO 200 J = 1,N
200 IF (I .NE. J) TX = TX - A(I,J)*X(J)
TX = TX / A(I,I)
ERR = ERR + ABS(TX - X(I))
210 X(I) = TX
220 IF (ERR .LT. EPS) GOTO 299
STOP
299 WRITE(*,'(1X,F12.6)') (X(I), I=1,N)
STOP
END

まだ77教えてるところなんてあるんだな。センセたちも古い規格しか知らなく
て「古い古い」言ってんだから世話ないな。

古いからこそ実績や信頼性の有るコードが多い。
研究室に入った学生が真っ先にする仕事は古いFORTRANコードのメンテだ。
FORTRAN77を学ぶことは、まだまだ有益なのだ。



決して、ン十年前の講義ノートを使い回しているだけじゃありませ（ｒｙ

古いと、妙な手あかがつきまくって、
難解で珍妙なコードになってることもあるぜ。

FORTRAN は言語的にバグを出しにくくするサポートが未熟で
デバッグしづらくて困る。

>>305
コンピュータ言語の教育という面ではF77は切り捨てたほうがいいのだが、
現実にメンテさせられるコードは７７以前なのだから実用性からはそれもありかなとｗ

まぁクイズとか解いてりゃいい新型言語とは違うしｗ

あ

fcpadでftn77を使っているのですが、あるプログラムをコンパイルして実行
しても結果が一瞬表示されて消えてしまいます。どうしたらいい？

>>310
ｐａｕｓｅ文を入れろ

>>311
これってもしかしosのてシステム上の問題？
xpつかってるからセキュリティのservice pack2が問題かもしれん。
コマンドプロンプトからプログラム名を指定して実行したら成功しますた

>>312
コンソールアプリはそれが普通。
処理を停めなければDOS窓はすぐに消える。
それが嫌ならpauseで停めろ。
それより結果はファイルに出した方が良いぞ。

>>313
ということはつまり結果をテキストファイルみたいにしてだすわけですね？
OPEN文つかうんですよね？

ご指導ありがとうございます

DOS 窓開いてから実行すればいい

ありがっと

OPEN文でファイルからデータを読ませようとしてるんですけど「ファイルがない」って
エラーが出ます。あるのに。


なぜですか

事故解決しました。
UNIT=1にするとエラー
UNIT=2にするとコンパイルはできるが入力データに上書きされてしまう
UNIT=3にするとコンソール画面に結果がでてウマー

ＣＳＱＲＴって使えるの？

>>319
FTN77で。

>>317
一桁目は歴史的な理由から予約されていることが多い。
処理系によっては９０番台も予約されていたりする。

特に理由が無ければ避けるのがよろし。


>>320
たしか、７７の規格には無かったはず。
ＥＸＰを使って解け。

>>320
こんなとき（だけ）役に立つ、岩波FORTRAN辞典〜♪
によるとFORTRAN77でもSQRTの引数複素数時の個別名はCSQRTだ。よかったな。

だが、総称名使え。

>>321
一桁目ってunit=1から9のことｄふぇすか？


複素数のはびっくりした。ああなるほどっておもた。オイラーの公式だね

exp(iθ)=cosθ+i sinθ

>>322
ということは、たとえばD=B**2-4*A*Cという式があって、REAL Dとしているとき、
sqrt(D)を計算しようとしますよね。BとAの与え方によってはD<0となり複素数と
なりますよね。そのときCSQRT(D)とすればいいのですか？

ま、あれですかね。複素数ってプログラムで扱うもんじゃないんですかね

>>323
CSQRTの引数は複素数型。
そしてFORTRANではC言語などの様に引数の型を（関数宣言に合わせて）自動変換したりはしない。

ということは

あらかじめ変数宣言でCOMPLEX Dとしていて
CSQRT(D)ならばありってことですか？

>>326
文法的にはあり

だがまあその場合、SQRT(D)としとけば内部的にCSQRTを呼んでくれるので
わざわざCSQRT(D)とする必要はないけどな。

>まぁクイズとか解いてりゃいい新型言語とは違うし

かっこいいな

ちょっとしつもんです。

txtファイルで
0　0.1　3.2　
6.3
1　2.1　5.6　
5.9
・・・

のようにちょっとずれたテーブルがあるんすけどこれを次のようにグラフ化できるように
整理できるプログラムって作れますか？
0, 0.1 ,3. ,6.3
1, 2.1, 5.6, 5.9
・・・

条件は入力がtxtファイルで出力もtxtファイルかcvsファイルにするものです。

作れるけど、わざわざ FORTRAN で作るもんでもないと思われ。

>>330
おしえてくださいませんか？また、フォトラン使わないならば何を使うのか教えてください。
お願いします先生

awk がちょうどいいんでない？

じゃあawkで

#! /bin/gawk -f
BEGIN { NMAX = 4; n = 0 }

{
for(i=1; i<=NF; i++) { a[n++] = $i; if (n >= NMAX) { putn(a, n); n = 0; } }
}

END { if (n>0) printn(a, n) }

function putn(a, n, i) {
for(i=0; i<n-1; i++) printf("%4.1f, ", a[i])
printf("%4.1f\n", a[i])
}

もっと簡単に書けるかもしれないけど、こんな感じで。

awkってなんですか

awkは語るまい。

awk というのは、加齢臭の漂うおっさんスクリプト言語。

パソコンの処理能力が上がって、データ後処理がパソコンで実用的にできるようになった頃に
流行っていたスクリプト言語がawkだった。あの頃はperlも何も無かった。

大型計算機のFORTRAN出力をパソコンに落として、awkで整形処理、ｘｙプロットソフトで描く。
という流れが８０年代末から９０年代初頭に流行った。

そういうわけでロートルはawkが使える。

じゃぴったりなわけだ。

awkってコマンドプロンプトさえあればできるんですよね？

"フォトラン"ってなんかカワイイ

フォトラン20 平成20年バジョン

>>337
基本的にコマンドラインでスクリプト・ファイルを指定したり、直接コマンドを書き込んだりして使う。

Linuxには大抵、最初から入ってる。
Windows用のもググれば色々見つかる。

にゃるほどシグウィンにするかな。いろいろありがとうごｚぁいましゅ

すげーすげーすげー
オークすげー

ここの人のおかげで先ほどのデータ並べ替えができるようになりました。
ありがとうございました。
オークの勉強もフォトランと一緒にしようとおもいます

>>329
もしFORTRANでやるならどうしますか

先ほどからしつこくてすみません。テキストデータの並べ替えを質問してAWKに感動したものです。
実はフォートランでもプログラムを組んでみました。

REAL A,B,C,D
C 読み込みファイル指定（既存）
OPEN(UNIT=3,FILE='data.TXT',STATUS='OLD')

DO 20 I=1,4
READ(3,*) A,B,C,D

WRITE(*,100) A,B,C,D
C 書き出しのフォーマットを指定
100 FORMAT(F6.3,',',F6.3,',',F6.3,',',F6.3)
20 CONTINUE
C 開けたファイルは閉じましょう
CLOSE(UNIT=3)

END

こんなぐあいでやったらコンパイルはちゃんと通ってデータも並びなおしてました。

でも、ちょっと腑に落ちないところがあります。
DO文のところです。俺は、なんとなくI=1,4として「これで1行から4行まで勝手に読み込んでくれるだろ」
という乗りでやったのですが実際のどげんですか？
このI=1,4はA,B,C,Dを読み込むのを4回やれということになりますか？
それとも1行目から4行目までA,B,C,Dを読むのを繰り返せというような行数を表すのですか？

普通、DO文でI=なんとかとしたら
処理のところにA=A+IとかIの入った式を書いたりするじゃないですか。でもこの問題の場合
Iを処理の回数として勝手に認識してくれているのであれば私のプログラムに間違いはない
ですよね

>>343
AWKは最近流行らないから、本とかも無いだろ。
今更感が無きにしも非ずだがFORTRAN界では大丈夫かｗ

PROGRAM test
IMPLICIT NONE
REAL :: a(4)
INTEGER :: i, ip
DO
READ(9, *, END = 999) a
WRITE(*, '(3(f10.4, ","), f10.4)') a
END DO
999 STOP
END PROGRAM test


>>344
あぁ、もう自作したかｗ

FORTRANでは歴史的な事情により、指定しない限りフォーマット付きI/Oでは、
入力データの空白は０として扱われる。（FORMAT文のNB/NZオプションあたりの説明嫁）

しかし＊での自由書式だと空白が区切り記号として解釈される。
また今の場合のようにデータが足りないと改行して読んでくれる。

大昔は容量を節約するため出力データは空白を開けないで続けて書いたりしていたｗ

質問の解答は『A,B,C,Dを読み込むのを4回やれ』ということになる。
ところがFORMATを指定した場合はこうはならない。
/などで明示的に改行してやらないといかん。

細部は漏れも忘れたので文法書でFortmatのところ嫁。

うわーありがとうございます！

ところでおすすめの書ありますか？

>>345

｡･ﾟ･(ﾉд`)･ﾟ･｡
そりゃあ流行じゃないし、数もそれほどじゃないけど絶版になってないのもあるよ
オライリーのとか。

>>336
うちの研究室では awk を教えられたぞ。

ホムペじゃあまだオークを語ってるところもあるね

FORTRAN よりは若い言語だよな。

>>351
そりゃFORTRANと比べたら……

どなたかご教授下さい！！

FTN95とfcpadの組み合わせで使いたいのですが、
fcpadでコンパイラのパスを「C:\Program Files\Silverfrost\FTN95\ftn95.exe」
として、いざf95のソースファイルをコンパイルすると
[FTN95/Win32 Ver. 5.01.0 Copyright (c) Silverfrost Ltd 1993-2006]
*** Invalid executable file suffix - only .EXE is permitted
と出てコンパイルできません。ファイルの拡張子を.fにしても駄目でした。

どうしたらFTN95+fcpadでコンパイルできるのでしょうか！？
おとなしくplato3を使うしかないのでしょうか・・・

>>353
前にも誰かが聞いていたような・・・・

まぁあのｷﾓｲ顔のプラトン使ってやれよ！

>>353
FTN95のコマンドラインの書式は
> ftn95 <filename> [option [option] ...]
一方、CPad (for FTN77)がコンパイラに渡すコマンドラインの書式は
> ftn95 [option [option] ...] <filename>
書式が違うのでそのままではエラーになる。

解決策としては例えば、↓の様な内容のバッチファイルを

@REM --- ここから ---
@shift
ftn95.exe %1 %2 %3 %4 %5 %6 %7 %8 %9 %0
@REM --- ここまで ---

ftn95.exe と同じディレクトリに作成し、fcpadの「基本設定」の「コンパイラのパス(W)」に
このバッチファイルを指定する。

しかし自由書式のプログラム作成とか考えるとCPad for FTN77を使うよりも
付属のplato3か他のエディタを使うことをお勧めする。

>>353
このスレの57のあたりを見直せば

>>355さん
ありがとうございます！
早速バッチファイルを作ってみましたが、結局うまくいきませんでした。
たしかにfcpadにこだわる必要はなさそうなので、vimを使ってみようと思います。

>>354さん
>>356さん
ご意見ありがとうございます！
前述の通り結局vimを使うことにしました^^;

こんにちは。現在Fortranでプログラムをしているのですが、難しくててこずっています。
そこで、知恵を貸していただきたいと思い、書き込んでいます。

A=x*exp(-u*B)+y*exp(-v*B)+z*B^(-(w+1))
という関数があり、x,y,z,u,v,wは未知数で
AとBはExcelのA列とB列にそれぞれデータが600与えられていて、
test.csvとして保存しています。

このx,y,z,u,v,wの未知数6個を上式の関数とA,Bの配列(?)を用いて,
収束させることによって求めたいと思っています。

ご教授お願いします。

dimension、open、write、read、doなどを使うのはなんとなく分かるのですが、
、難しくてよく分かりません。
どうぞよろしくお願いします。

salford FTN90 とPLATO3のダウンロードできる場所をおしえてくださいお

>>358
まずは問題点を切り分けよう。
ファイルの読み込みと実際の計算部分だ。
どちらならできる？　両方とも難しいかい？

質問させてください
DO i=1,10
　DO j=1,10
　　DO k=1,10
　　　dt_den(i,j,k) = den(i,j,k) * div_vel
　　END DO
　END DO
END DO

というコードを書いたのですが、『ランクが違うので計算できない』という旨のエラーが出ました。
「配列の一つの要素を取り出し、計算して、別の配列の要素へ代入する」
ということをやりたいのですが、どのように書けばいいのですか？

Fortran77のコンパイラの選択に関する質問です。

私はこれからFortran77を学ぼうとしているところです。
利用の目的は学術的な数値の計算です。
有償、無償の様々なコンパイラがあることは調べているうちにわかりました。
しかしなにぶん右も左もわからないもので、どれが良いか決めかねています。
信頼性を重視した場合、やはり有償のものを購入する方がよいのでしょうか。
また、お勧めのコンパイラがありましたら教えていただけますと幸いです。

当方Windowsを使用しており、主な環境は以下の通りです。
OS:WindowsXP Professional
CPU:Pentium D(4?) 2.8GHz
Mem:1GB

なお、購入の場合は5〜10万円程度の予算を考えております。
どうぞよろしくお願い致します。

!これから学ぶのに、いきなり5-10マソも出すなんてすごい！しかも77だなんて。
!どこかの研究者さんでしょうか。メモリ1GBということは、大規模な計算では
!なさそうなので、無料のG95かgfortranでいい気がするのですが、どうなんでしょう？
!学生の私には判断できません・・・。

上に同じ。
安い／無償の環境と高い環境の違いは
・開発のしやすさ
・ライブラリの有無（例：高速な数値演算ライブラリ）
・最適化の性能（＝計算速度）
あたりにあるので、まずはフリーのコンパイラを使ってみたらどうか。

すみません。自分で解決しました。原因は
REAL(8),ALLOCATABLE :: div_vel(:,:,:)
としていたことでした。正しくは
REAL(8) :: div_vel
でした。

>>358
その式を解くとなると、
非線形の最小二乗法（または、自乗法）が
必要になります。
もちろん、FORTRANで計算可能ですが、
一からプログラムを書くとなるととっても大変だよ。

ご返信ありがとうございます。
研究に用いるため信頼性をなによりも重視しており、
予算も高めに覚悟しておりました。
有償無償の間にバグに関して差がないのであれば、
最初は無償のものを使おうと思います。
Salfordのものを使うつもりです。
相談に乗っていただきありがとうございました。

>>367
大規模計算はスパコンで、解析は自分のPCでフリーのコンパイラだけど、
特に問題は出ていない。
少なくとも、今のところ、バグで酷い目にあったことは無いかな。

知り合いは、その昔、スパコンのコンパイラのバグで酷い目にあったそうだから、
有償／無償の差よりは、運／不運のファクターの方が大きいかも。
多くの人が使っているものを選べば、まぁ、大丈夫だと思う。

コンパイラのバグに出会ったことは無いけれど、
大規模な新型スパコンが導入されると、設定ミスが発覚したり、
ジョブ管理のシステムの挙動がおかしかったり、トラブルが続いて
落ち着いて使えるまでに少し時間がかかることもある。
数が多く出ているようなシステムは、こういうことは少ないと思うけど。

>>367
金があるなら、Intel Fortran + MSVisualStudio　の組み合わせでいいんじゃないの？
ソフトは利用者が多い物が勝利するので、寄らば大樹の陰。

VisualStudioはアカデミックなら５,０００円なので無問題。
IntelFortranは、やや高いし、毎年＄２００のお布施でライセンス更新をしないといけないが。
（再インストール出来なくなるだけで、使えなくなるわけではない）
Linux版でよければただで使える。

バックエンドが富士通なら、富士通/Laheyのコンパイラを使うという手もある。

Fortranは未だコンパイラ市場が存在できてい牙城なので、フリーソフトを排撃してこれを守ってゆこうｗ

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

>>370
write(6,'(A)') '1964/1/1 115 89 80 81 104 154 209 256 287 298 285
&260 229 199 176 170 189'
write(6,'(A)') '1964/1/2 155 111 89 81 91 124 171 227 266 293 292
&274 252 221 190 174 181'
write(6,'(A)') '1964/1/3 191 147 112 92 94 113 151 194 241 271 290
& 283 267 241 214 189 176'
write(6,'(A)') '1964/1/4 222 185 147 117 101 107 134 168 208 248 2
&74 284 274 256 233 208 185'
write(6,'(A)') '1964/1/5 238 213 184 156 134 127 139 164 193 227 2
&59 276 274 267 249 226 204'
write(6,'(A)') '1964/1/6 248 234 214 190 165 151 149 161 181 206 2
&34 257 267 266 257 243 222'
write(6,'(A)') '1964/1/7 234 236 229 215 199 184 175 175 184 203 2
&24 241 256 264 259 248 233'
write(6,'(A)') '1964/1/8 216 228 234 229 221 213 203 192 190 197 2
&10 225 239 252 258 254 243'
write(6,'(A)') '1964/1/9 193 212 227 237 238 231 220 209 199 194 1
&95 205 217 232 243 250 246'
write(6,'(A)') '1964/1/10 163 181 203 226 238 243 240 229 217 205
&195 196 202 215 232 245 251'
end

PROGRAM omanko
IMPLICIT NONE
CHARACTER (LEN = 80):: fmt
INTEGER :: iy, im, id, idat(17)
fmt = '(i4, 2("/", i2.0), 17i4 )'
iy = 1964
im = 3
id = 24
idat = (/115, 89, 80, 81, 104, 154, 209, 256, 287, 298, 285, 260, 229, 199, 176, 170, 189/)
WRITE(*, fmt) iy, im, id, idat
STOP
END PROGRAM omanko

>>371
それはDO文でやった方がいいですねｗ

>>372
年月日のところが上手く表示できなかったので
参考にしてみます。

>>371,372
感謝です。

plato3で作成したプログラムのエグゼファイルのアイコンをアイコンエディターでつくった
アイコンに変えるにはどうしたらいいですか？

>>374
プラトン先生はよく知らんが、リソースファイル書いてアイコン指定すれば？

Fortranで何ができるのかも良くわかっていないのですが，
データを出力したときにメモリーの状態はどうなっているのか？
ということを教えて下さい．

その内容は
｢Fortranでエクセル形式データを出力させた場合，データは一度
メモリーに格納され，エクセルへ出力するのだと思っているのですが，
出力した後，メモリーは開放されるのでしょうか？」ということです．

当方の環境は，
コンパイラ　：　fujitsu　fortran77
OS　： WindowsXP Professional


以上，よろしくお願いします．

死にたい

>>376
質問の意味が不明瞭だ。

原則としては、プログラムが終了するなら、メモリーは解放される。　

実行中は、静的に確保した部分は確保されたままだし、動的に確保されたものはDEALLOCATE
すれば開放される。




ただしOSやコンパイラにバグがなければ。　
昔のMS-PowerStationにはメモリーリーク・エラーがあって実行中にI/Oがあると、
徐々に使用メモリーが増えていった。

>>377
死ぬなら中国大使館に爆弾持って特攻だ！！

>>379
FXって言う為替取引で倍率を200倍でやると良いって、死ぬきあるなら博打をしてから詩ねってことだな

　　　　　　 ｒ'::::::::,:::::::::::::::::::::::＼ 　　　　【USBオナホを実現させる会】に入会しませんか＾＾＾
　 　 　 　 　＼::/::;;;:::::::::::::::::::::::::＼
　 　 　 , ─ -ｖ_;;. -─ ─ ── ─-、このオナホが実現すれば、パソコンにUSB接続して
　　　　|　, 　 ｆ-‐ ''_ .二 二二二 二._t 　対応エロゲを立ち上げてペニスを挿入するだけで
　 　 　 Y__,rＹニﾆ -- ‐─ ──‐‐─`ー─‐-- ､画面のキャラクターが喘ぐんだよぉ♪
　　,.．‐Ｋﾉ,ｈ;;::．,-ｒｪ;ェｴｆ７ｆｖ!Ｔl ＴｉＴｉt.､--':::::::::ﾉ　キミの腰の動きに合わせて
　 (:::::::::l::::ヽ.|l.|/,|l.1,!Ｈ､ﾄ|｀ 　H┼I､ｌ|l |l ﾄ.-‐ ''´　わたしたちがアニメーションするの。
　　`,ニｌ:::::::::ﾚｌ/!|lｲｈ_ ｌ|｀｀　,ｲrﾏ,ｲ| |l |ﾚﾉ　　　　　とっても素敵なアイテムでしょっ☆
　 く ／|:::::::::ﾚl.ｨ!ｌl.l. ｰ‐' 　　, kノﾊﾄ,ｲﾚ:::::l
　／ / !:::::::::|1l_l.|l.|ゝ""　(.ｱ　 ,仆!|::´:::::::| 　　ゆっくり動かすと・・・じわじわと感じて、
└ｔ. ｌ　ｌ:::::::::Ｎ'::: Ｎ´ヽ`,,ｰr t１リ::ﾚ::::::::::: !　　　はげしく動かすと・・・失神とかもするんだよ☆
　　`ｰ弋:::::::｀:::::::::ゝ ﾆァ＠'ノ::､:::::::::::::::: ﾉ
　 　　　`ー､::::::: /:::::::,ｲ !ﾄ`ヽ::ヽ:::::, ‐'´　　わたしたちとバーチャルで繋がるから
　　　　　　　 `ー1::::::└Ll_rｰ'::::::〉´　　　　もう、即射精しちゃうよねっ☆

　でも・・・USBオナホールは、まだ実現していないの・・・。
　はやくみんなとエッチなことしたいなぁ・・・。
　USBオナホールを実現させて次元を超えたSEXをしましょうねっ☆
　入会方法はとっても簡単。毎週1回どこかのスレに↑をコピペするだけでいいヨ♪

Intel Visual Fortran9.1を使っています．
Fortran90形式でプログラムを書いており，
マトリックスを表示させたいのですが，これに関して質問があります．

現在，2次元配列に格納されたマトリックスを表示するために，

integer　　　　　　　　　::　　i,j,N
integer,allocatable　　::　　A(:,:)
〜〜〜
N=10
allocate(A(1:N,1:N))
A = 123
do i=1,N,1
　write(6,'(100I3)') (A(i,j),j=1,N,1)
enddo

という風にしています．
つまり，横に並べて表示するformatとして 100I3 という風に反復回数を十分大きくとって
対応しています．
この反復回数100の部分を変数にしたいのですが不可能なのでしょうか？

>>382
文字型変数で書式を作って渡す。
character(len=20) :: fmt
...
write(fmt, "('(', I5, 'I3)')") N
do i=1,N,1
write(6, fmt) (A(i,j), j=1,N,1)
end do

または、Visual Fortran限定で良ければ書式中の <> 内に整数式を書き込む。
do i=1,N,1
write(6, '(<N>I3)') (A(i,j), j=1,N,1)
end do

>>382
これは実はFAQ。このスレでもすでに何回も出てきている。

教科書も少しはこういう事を書けｗ

>>383
ありがとうございます．
助かりました．

式の前に入れる空白の数によってコンパイルエラーが出たりでなかったりするのですが、
これはどのように決まっているのでしょうか？

>>386
ソースやエラーメッセージも貼らずに（ｒｙ

空白を追加したら１行が73桁以上になった、に10ｴｽﾊﾟｰ

73行以上っていつの時代のコンパイラだよｗ

gfortran やIntel Fortran 9.1 でも固定形式ではしっかり73桁目以降を無視してくれますよ。

>>388-389
そうじゃないぞ諸君。
一般的なコンパイラでは拡張子によって、デフォの挙動が変わる。
拡張子が.fや.forの場合は、固定書式になって７２カラムモードとなる。
.f90などの場合は自由書式になる。

あと、ワークステーションなどでは、よく同じコンパイラが、コマンドとしてｆ７７とｆ９０呼び出される設定に
なっていて、デフォの挙動を変えて作動する。この場合も、固定書式と自由書式にそれぞれなっている。


まぁFORTRANを使うなら固定書式の爽やかさを味わってみるのもいいだろう。

n(nx,ny)
という配列をバイナリ形式で
write(8)n
という風に書き出した場合、格納順序はどのようになるのでしょうか？
n(1,1),n(1,2)....
でしょうかそれとも
n(1,1),n(2,1)....
あるいは処理系に依存なのでしょうか？

>>391
n(1,1),n(2,1)....
ただし格納されたデータ形式は処理系依存

>>391
FORTRANでは規格でn(1,1),n(2,1)....の順にメモリー上に入っている。
さすがにこれは教科書にも書いてある。

というか、これを知らないと整合配列がらみでF77などでは実用的なサブルーチンとか書けない。
またライブラリを呼ぶにしても引数の意味が理解できないだろう。

まぁマニュアル嫁。

合同乗積法による乱数生成プログラムを作成しなければならないのですが
どこか参考になるページなどはありませんでしょうか。
かれこれ2、3時間探しているのですが見つかりません。

配列の初期化について質問があります．

以下のプログラム (Fortran90です) のように，
変数宣言部で初期化すると，data文が効かずににMatrixの要素がすべて0になってしまいます．
もちろん変数宣言部を

integer :: Matrix(1:3)

に書き換えるとMatrixの要素はそれぞれ 1, 2, 3 と入ってくれます．

これは一体どういうことなんでしょうか？
最適化を外しても同じです・・
頭が変になりそうです．

!---プログラム-----
program MATRIX_TEST
implicit none
integer :: Matrix(1:3) = 0
data Matrix/1,2,3/
write(6,*) Matrix
end program MATRIX_TEST

色々やった結果，以下のようにしたら期待通りになりました．
変数宣言部で初期化するとsave属性が付くという仕様に気づきませんでした．
data文はsave属性の変数に値を入れることができない(？)．
何がなにやら・・

もうdata文を使うのはやめます．

!---プログラム-----
program MATRIX_TEST
implicit none
integer :: Matrix(1:3) = 0
!data Matrix/1,2,3/ ←やめて，下のようにした．
Matrix = (/1,2,3/)
write(6,*) Matrix
end program MATRIX_TEST

そういえばfortran90ではdata文を使うのはやめましょうっていう記述結構見かけるんだけど、
多次元配列にデータ入れたい場合どうやってます？

たとえば2×2のマトリックスを3個用意する場合。
1次元目がマトリックスを区別し、2次元目が行、3次元目が列を表すものとして、

integer :: A(1:3,1:2,1:2)
data A(1,:,:) /1,2,3,4/
data A(2,:,:) /5,6,7,8/
data A(3,:,:) /9,10,11,12/

に相当するコードをdata文なしでどうやって書きます？
>>396みたく

A(1,:,:) = (/1,2,3,4/)

と書ければいいんだけど、これはコンパイルエラーになる…

変数宣言部で一気に

integer :: A(1:3,1:2,1:2) = (/1,2,3,4,5,6,7,8,9,10,11,12/)

ってやるしかない？？
これだと目的の要素に値入れるのにメモリ上の並びを考えないと
いけないので4次元5次元となると、きついっす。

>>394
合同乗積法による乱数　でググルとトップにFortranのプログラムが出てくるぞ！
しょんべんで顔洗って出直して来い！

>>395-397
DATA文は実行文ではなく初期化文だ。
F90では初期化が直接出来るので、わざわざDATA文を使う理由がない。

宣言文での初期化とDATA文での初期化のどちらが優先されるのかは、規格書を見ないと分からんが
そんな使い方は想定していないので、たぶんメーカーの実装に依存していると思われる。

A(1,:,:) = (/1,2,3,4/)　がうまく行かないのは、次元があっていないから。
RESHAPE関数を使って次元をあわせればおｋ。これはよく教科書にも載っているので調べそ。

二項分布のプログラムについてなのですが（硬貨の表裏）、表を0、裏を1とした場合
それぞれが1/2の確率で現す書き方がわかりません。
超初心者ですみません。

>>398
サンクス！
調べてくるお

すいません、質問させてください。
コンパイル時どうしても
/opt/intel/fc/9.0/lib/for_main.o(.text+0x41): In function `main':
: undefined reference to `MAIN__'
と出てしまうのですが、これはどういう事なのでしょうか。

>>399
０以上１未満を返す乱数のサブルーチンは用意してあるとする。
このとき０以上０．５未満だったら表、０．５以上１未満だったら裏とすればよい。

>>401
確認していないので想像で書くｗ
メインルーチンが無いだけじゃない？プログラム文が無いとか。
昔のHITACのFORTRANだと、サブルーチン名をMAINにすると、
それをスタートアップルーチンと勘違いしていきなりそこから実行が始まるウンコ仕様だったがｗ

なるほど。。
やってみます。ありがとうございました。

NaNを検出するような関数ないですか・・

すみません自己解決しました
isnan(real)でした

>>404-405
IEEEのフラグがらみの関数はFortran2003で導入される予定なので、それはコンパイラ独自の拡張。
コンパイラの種類などを書かないと、誰も答えられないし、何の情報にもならん。
おまんこ。

g77はFortran90の機能も内包しているのですか？

YES
おまんこ。

>>407
ちょっとだけだけどな。

実行時間を計ろうと思って色々検索したら，dtimeというサブルーチンが
用意されていると知り，使ってみたのですが，よくわからないエラーが発生してしまいました．

>>410のプログラムをコンパイル&実行すると，なぜか出力がNaNになってしまいます．

call dtimeを削除すると期待した値が出るのですが，なんでこんなことになるのかわかりません…
どなたか助けて下さい　(´；ω；｀)
ただdtimeがfortran標準でないという記述も見かけたので，使わない方がいいのでしょうか…

コンパイラはIntel Visual Fortran 9.1で，コマンドは f90 /check:all /Od test.f90 でした．

すみません>>410ではなくて>>411です orz

!-- test.f90--
program test
implicit none
real(8) :: x,y
call SUB
x = 11.d0
y = Func(X)
write(6,*) y

contains
subroutine SUB
implicit none
real tarray(2),result
call dtime(tarray,result)
call dtime(tarray,result)
return
end subroutine SUB

complex(8) function Func(x)
implicit none
complex(8),parameter :: IU=(0.d0,1.d0)
real(8),intent(in) :: x
Func = exp(IU*x)/sqrt(2.d0)
return
end function Func
end program test

>>410
Intel Visual Fortranの DTIME は引数が１つの関数だ。サブルーチンじゃない。
SUB の先頭に USE IFPORT をつければちゃんとコンパイル時にチェックしてエラーを出してくれる。

>>412
ｱｯー！たしかにそうです！
ライブラリ・リファレンスに書いてました・・・・・
ありがとうございました！

(1) ２×２の行列A，B の要素を入力し，A＋B とAB およびA−1 計算する部分をそれぞれサブル
ーチンとするプログラムを作成せよ．

A-1ってA-I(もしくはE)のことか？

ああ、逆行列のことか。スマン。

>>414
program main
use la_precision, only: WP => DP
implicit none
integer, parameter :: N = 2
real(WP), dimension(N,N) :: A, B, C, D, E
print *, "input A:"
call sub_read(A)
print *, "input B:"
call sub_read(B)
call sub_add(A, B, C)
print *, "A + B ="
call sub_write(C)
call sub_mul(A, B, D)
print *, "AB ="
call sub_write(D)
call sub_inv(A, E)
print *, "inv(A) ="
call sub_write(E)

contains
subroutine sub_read(X)
real(WP), intent(out) :: X(N,N)
integer :: i, j
do i = 1, n; read *, (X(i,j), j=1,N); end do
end subroutine sub_read
subroutine sub_write(X)
real(WP), intent(in) :: X(N,N)
integer :: i, j
do i = 1, n; print *, (X(i,j), j=1,N); end do
end subroutine sub_write
subroutine sub_add(X, Y, Z)
real(WP), intent(in) :: X(N,N), Y(N,N)
real(WP), intent(out) :: Z(N,N)
Z = X + Y
end subroutine sub_add
subroutine sub_mul(X, Y, Z)
real(WP), intent(in) :: X(N,N), Y(N,N)
real(WP), intent(out) :: Z(N,N)
Z = matmul(X, Y)
end subroutine sub_mul
subroutine sub_inv(X, Y)
use f95_lapack, only: la_getrf, la_getri
real(WP), intent(in) :: X(N,N)
real(WP), intent(inout) :: Y(N,N)
integer :: ipiv(N)
Y = X
call la_getrf(Y, ipiv)
call la_getri(Y, ipiv)
end subroutine sub_inv
end program main

>>417-418
正直なところ、たかがこれしきの行列演算をサブルーチンにする意味が分からん。

おまけに逆行列はLAPACK使ってるし。２＊２限定なら高校で習う公式で十分だろ。

すべてがjokeなのかもしれないが

宿題か何かなんじゃないか？
っていうか超初心者用スレなんだから
サブルーチンにする意味なんか問うな。
Hello worldなんか意味がないとか言っているのと同じだ。

>>421
たぶんＦＯＲＴＲＡＮ７７での宿題なんでね？
Ｆｏｒｔｒａｎ９０では意味を成さない希ガス。

>>421
何その再帰呼び出し的アンカー

>>422
f90 へのあてつけ

最新のgfortranをvistaにインストールし、コンパイルしようとしたら
ld: cannot find -lgfortranbegin
とエラーが出てコンパイルできません。
XPではインストールして問題なくコンパイルできました。
vistaでは使えないのでしょうか？

>>423
f90 は再帰できるじゃん。
F77 へのあてつけなら分かるが。

>>424
状況がよくわからんが Program Files 以外の日本語や半角スペースを含まないフォルダにインストールし直せばいいんじゃね？

>426
それも考えてc直下にインストールし直しても結局同じでした。
vistaだからなのかな･･･誰かvistaで動いてる人は居ないんでしょうか。
g95も同じようにエラーが出て使えませんでした（今手元にvistaマシンが無いのでエラーメッセージ分かりません
FTN95はvistaでも動きました。

>>424
まったく分からんが、ローダーのエラーなのでコンパイルは行ってるんでない？
単にメインプログラムがないだけとか？

>428
XPマシン+gfortranでコンパイル→実行できたソース(非常に簡単なプログラム）をコンパイルしようとしたんですが，
上記のようなエラーが出て何も生成されませんでした。
ソースの指定をしないとno input filesが出るので，インストールもされてるしpathも通ってはいるようですが・・・

とりあえず明日にならないとvistaマシンはさわれないので今はどうしようもないですが，
何か分かる方はアドバイス御願いします。

平方根について質問があります．
複素数の範囲で -1 の平方根をとった場合，
精度の関係で +i が出力される場合と -i が出力される場合があります．
つまり，

Z1 = (-1.d0, 1d-10)
Z2 = (-1.d0, -1d-10)

のような時に，

Q1 = sqrt(Z1)
Q2 = sqrt(Z2)

を計算したような場合，Q1=i，Q2=-i　となります．皆さんはこのような問題を
どうやって排除していますか？
絶対値で見たとき，ある値以下では0.d0で置き換えるという作業も考えられますが，
出る値のオーダがあらかじめ分かっているときしかできません．

それとも初めからルートをとるような計算を
避けるようなアルゴリズムを考えるのでしょうか？

また，平方根をとることによって精度が半分になってしまう場合が多々ありますが，これも諦めるしか
ないのでしょうか．例えば理論上では0となるべきところで，

R = 1d-10

となっているようなときです．

Rと0.d0とを比べたときの誤差は許容範囲内だとしても，sqrt(R)とsqrt(0.d0)とではかなり違ってきており，
sqrtを計算するたびに誤差が馬鹿にならなくなってきます．．．

このようなときに行う定石みたいな処理があればアドバイスお願いします．

i と -i のどちらかでなければダメなのであれば
その問題において何らかの束縛条件があるはずだから、
それに照らし合わせて符号を決めればいい。

どちらでもいい式になってるのなら気にする必要は無い。

どっちかというと、FORTRAN の問題じゃなくて数学の問題かな。
x^2 = -1 の時 x = ±i なわけだけど、
条件によっては +i か -i のどちらかじゃないといけないが、
どちらでもいい、あるいは両方ないといけない場合もある。

>>430
ルートを取ったときに誤差が増えるといっているが、そもそも1＾-10のルートが1^-5なのだから
そうなるのは当然のことだ。

そもそも根本的な発想が間違っている。それゆえに対処法が見当違いで、ど壺にはまっている。
>>431のいうことももっともである。

身近な詳しい人に教えを乞うて、正しい発想をしなさい。

正の整数nを読み込みr=1,2,...,nに対する
nCr=n!/r!(n-r)!
を計算するプログラムを教えてくださいお願いします。
へたれなもんでスマセン。

宿題は自分で（ry

nCr = n(n - 1)・ … ・(n-r+1) / r(r-1)・ … ・1
nC(r-1) = n(n - 1)・ … ・(n-r+2) / (r-1)・ … ・1

nCr = nC(r-1) * ？ / ？
nC0 = 1

nC(n-r) = nCr

これだけ分かってれば解けるはず。

>>433
program hoge
implicit none
integer :: n, r
read *, n
do r = 1, n
print *, n, "C", r, "=", nCr(n,r)
end do
contains
recursive function nCr(n, r) result (c)
integer, intent(in) :: n, r
integer :: c
if (n < 0 .or. r < 0 .or. r > n) then
c = 0
else if (n == 0 .or. r == n) then
c = 1
else
c = nCr(n-1, r-1) + nCr(n-1, r)
end if
end function nCr
end program hoge

パスカルの三角形ってまた重いものをw

>>436
再帰は出来るだけ避けろ！ｗ　　情報工学科じゃあるまいしｗ

再帰が問題じゃないだろー。

gccをfedora core 5、athlon64環境で使っています。
CのfloatとFortranのrealは同じ精度なのでしょうか？

>>440
同じ

質問があります。
今大学の研究室でプログラミングをやっているのですが、
みんなでやっても出来ません。
問題は「fl」というファイルに20個の数字がランダムにあって、
その数字を大きい順に並び替えるという問題です。
ちなみに今私が作ったプログラムは
program ex4
implicit real*8(a-h,o-z)
dimension a(100)
do 10 i=1,19
open (10,file='fl.dat')
read(10,*) a(i)
do 20 j=i+1,20
if(a(k).gt.a(j))k=j
20 continue
b=a(i)
a(i)=a(k)
a(k)=b
10 continue
open(11,file='ans.dat')
write(11,*) a(i)
pause
stop
end
です。どこが間違っているのかわかる方がいましたら教えて下さい。

>>442
院生にきけ

ループの中でファイルを開いてるのが意味分からない。

open文とread文をループの外にやっても正解のファイルが出てきませんでした。
心優しい人の返信待ってます。ほんとに困っているのでよろしくお願いします。

あと今の状態の文は以下のとおりです。
program ex4
implicit real*8(a-h,o-z)
dimension a(100)
open (10,file='fl.dat')
read(10,*) a(i)
do 10 i=1,19
k=i
do 20 j=i+1,20
if(a(k).gt.a(j))k=j
20 continue
b=a(i)
a(i)=a(k)
a(k)=b
10 continue
open(11,file='ans.dat')
write(11,*) a
pause
stop
end

ファイルの中身も貼らないようでは困る。
あと、プログラムは間違いすぎていてどう指摘していいのか分からん。

アドバイスしたいけれど、院生がだした課題でしょ？
これから研究する上でプログラミングが必要なのだから、
できるだけ努力しようよ。
　# 気むずかしい院生に相談するというのも努力の一つ。

一年生なら、初心者救済のために教えたいけど、
四年生なら、今がんばらないと年末に泣くだけだよ。

課題ができなくて院生に怒られるのも、必要な経験のうちだよ。
また、自分達の理解度を正しく伝えるということも、
恥ずかしいけれど必要な作業。

ファイルの中身はこのような感じです。
21.0
54.0
32.0
68.0
64.0
10.0
22.0
98.0
65.0
84.0
81.0
39.0
93.0
74.0
55.0
11.0
71.0
97.0
19.0
50.0
そんなに間違えてるんですか？
よろしかったらプログラムを書いていただいてもよろしいでしょうか？

書くわけないじゃん。

program ex4
i = system('cat fl.dat | sort -n > ans.dat')
end

>>442
とりあえず，

read(10,*) a(i)

の下に

do 999 i=1,20,1
write(6,*) a(i)
999 continue

等と追加し，配列aに入っている値を調べよ．
致命的な間違いに気付くであろう．

ｱﾝｶｰﾐｽった>>446ね．

フィボナッチ数列　a0＝0、a1＝1、a(i+1)＝ai+a(i-1) (i=2,3,…)
を計算するプログラムを作成せよ。なお計算は例はa(i+1)=20まででよい。
条件１：得られた値が奇数であるならば、数値の右に＊をつけ、奇数かつ３の倍数ならば、数値の右に＊＊を付けて出力すること。

どなたかわかる方いませんか？
プログラムお願いします(>_<)

宿題は自分でやれ

研究室の本棚に回答あるぞ。

F77か知らないが、解説書やプログラムの本にも書いてるぞ。
「並び替え」で調べろ。それの英単語で出てくるだろうけど
敢えて書かない。

>>452
ゴミが入って、そのゴミが正解の場合もあるので
read(10,*) a(i)
の前に
a(i)=-1.0
とかを入れた方が良いかな？

>>446
今の段階で効率良いプログラムを書こうとするな。
とにかく動くプログラムを書け。

b=a(i)
a(i)=a(k)
a(k)=b
が、どういう時に作用させたいか、今現在はどのように動いているか
調べろ。

>>456
いや・・・そもそも i に値が入る前にiを使っているのと、
このread文一個では、ファイルからのデータの読み出しがまず失敗している。

てか間違いが多すぎてアドバイスしにくいんだなそもそも。


>>446
まずデータを配列に格納することをちゃんと成功させてから、
ひとつひとつ問題をクリアしていくことをお勧めする。

間違い多すぎるわけではないか。
1か所修正して2行追加するだけでとりあえず答えは出るところまで来てるんだな。

>>467
あれ？
>>442では、ループ回してる中にread文あったけど....

ループ不足は分かってたけど書くと、致命的なことが安直に分かるので
-1.0を入れろってしたんだ。

この人致命的に、理解してないね。
今の文でread(10,*) a(i)の後ろでループ回すの教えると、もっと根本的にミス
するだろうし。

あ、上の
>>467
は>>457です。
orz

fortranでwebの結果を得ることできますか？

get http://www.hoge.com
を行い
結果をファイルに残す

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

>>451
sortのスイッチは-rを追加しないといけないんじゃね？
なんか大きい順って言ってるし．

でもこの3行で「できました〜」って言ったらどう思われるかね．

固有値解析もサブルーチン呼ぶだけだから良いんじゃね？
なわけねーだろ。

>>461
system 関数で wget を呼べば

>>464
なるほど、getは出来ますね。純粋なfortranでは無いですが贅沢は言えません。

戻り値が200(存在する)場合はファイルに書き
304であれば、書かないようにすれば良さそうですが
戻り値を得る方法ありますか？

>>465
wget の出力をファイルにリダイレクトしておいて、それを解析するとか。

>>465
ああ、あと、当然成功すれば 0 が、失敗すれば非 0 が返ってくるはずだから、
それで成功か失敗かは分かる。
エラーの種類は出力を解析する必要あり。

>>466
ああ、そうか。
wcなりで、有るときにファイル作業すれば良いね。
でも、permission や loginの時にどう対処すれば良いのか？
ともかく、一応の道筋は付いた。

ありがとう。

戻り値を得る方法をご存知の方、いらっしゃったら、よろしくです。

問題は fortran で実装する必要性だな。
何らかのスクリプトでいい気がする。

fortranかVBい外のbasicでなら、思うことが自在にできます。
文字列や数値をくっ付けたりとか、あと解析に便利なんです。

つーのが、fortranでする理由です。

そこら辺は perl とか ruby とかの方が得意な気が。
まあ、言語を覚える手間がかかるというのなら仕方がないけど、
bash あたりでもいいという気はする。
ま、いいけどね。

>>461
Windows＋VisualFortranの場合だが、参考例。

オプション：Linker->Input->additinal dependency 　　WinInet.lib

PROGRAM tashiro
USE ifwina
IMPLICIT NONE
INTEGER, PARAMETER :: INTERNET_OPEN_TYPE_PRECONFIG = 0 ! // use registry configuration
INTEGER, PARAMETER :: INTERNET_OPEN_TYPE_DIRECT = 1 ! // direct to net
INTEGER, PARAMETER :: INTERNET_OPEN_TYPE_PROXY = 3 ! // via named proxy
INTEGER, PARAMETER :: INTERNET_OPEN_TYPE_PRECONFIG_WITH_NO_AUTOPROXY = 4 ! // prevent using java/script/INS
!

INTERFACE
INTEGER (HANDLE) FUNCTION InternetOpen(arg1, arg2, arg3, arg4, arg5)
USE ifwinty
!DEC$ ATTRIBUTES DEFAULT, STDCALL, DECORATE, ALIAS:'InternetOpenA' :: InternetOpen
INTEGER (LPLONG) :: arg1
INTEGER (DWORD) :: arg2
INTEGER (LPLONG) :: arg3
INTEGER (LPLONG) :: arg4
INTEGER (DWORD) :: arg5
END FUNCTION
!
INTEGER (HANDLE) FUNCTION InternetOpenUrl(arg1, arg2, arg3, arg4, arg5, arg6)
USE ifwinty
!DEC$ ATTRIBUTES DEFAULT, STDCALL, DECORATE, ALIAS:'InternetOpenUrlA' :: InternetOpenUrl
INTEGER (HANDLE) :: arg1
INTEGER (LPLONG) :: arg2
INTEGER (LPLONG) :: arg3
INTEGER (DWORD) :: arg4
INTEGER (DWORD) :: arg5
INTEGER (DWORD) :: arg6
END FUNCTION
!
INTEGER (BOOL) FUNCTION InternetReadFile(arg1, arg2, arg3, arg4)
USE ifwinty
!DEC$ ATTRIBUTES DEFAULT, STDCALL, DECORATE, ALIAS:'InternetReadFile' :: InternetReadFile
INTEGER (HANDLE) :: arg1
INTEGER (LPVOID) :: arg2
INTEGER (DWORD) :: arg3
INTEGER (LPDWORD) :: arg4
END FUNCTION

!
INTEGER (BOOL) FUNCTION InternetCloseHandle(arg1)
USE ifwinty
!DEC$ ATTRIBUTES DEFAULT, STDCALL, DECORATE, ALIAS:'InternetCloseHandle' :: InternetCloseHandle
INTEGER (HANDLE) :: arg1
END FUNCTION
END INTERFACE

INTEGER (HANDLE) :: hInt, hUrl
INTEGER (BOOL) :: iret
INTEGER (DWORD) :: dwret, dwtot
CHARACTER (LEN = 80) :: agent = 'Mozilla/4.0 (compatible; MSIE 6.0; Windows NT 5.1)'//CHAR(0)
CHARACTER (LEN = 256) :: strurl, buff
CHARACTER (LEN = 256), ALLOCATABLE :: srctext(:)
INTEGER :: i
hInt = InternetOpen(LOC(agent), INTERNET_OPEN_TYPE_PRECONFIG, NULL, NULL, 0)
IF (hInt == 0) PRINT *, 'Open error'
strurl = 'http://www.yahoo.co.jp'C
hUrl = InternetOpenUrl(hInt, LOC(strurl), NULL, -1, 0, 0)
dwtot = 0
ALLOCATE(srctext(1000))
OPEN(9, FORM ='binary')
DO i = 1, 1000
iret = InternetReadFile(hUrl, LOC(buff), LEN(buff), LOC(dwret))
IF (dwret == 0) EXIT
dwtot = dwtot + dwret
WRITE(9) buff(1:dwret)
srctext(i) = buff(1:dwret)
IF (i == 1000) STOP 'buffer overflow '
END DO

print *, srctext(1:i - 1), srctext(i)(1:dwret)
print *, dwtot, 'bytes'
iret = InternetCloseHandle(hUrl)
iret = InternetCloseHandle(hInt)
STOP
END PROGRAM tashiro


面倒なのはWinInet用のインターフェースが用意されていないので、自分でC用のヘッダーファイルを見て、
作ってやらないといけないことか。　

若者を呼び寄せるために
Fortranでエロゲでも作らないか？

>>472
おお、出来るんですね。f77では、ないけど読めます。
ありがとう。
数値解析程度にC/C++やるので、なんとかできると思います。

このプログラム…Ｎを11以上にすると動かなくなってしまいます…
どうしてでしょうか？
ちなみに、コンパイラはsalfordのPlato3を使っています。

program main
integer:: N,k
integer, allocatable:: F(:,:),A(:,:)
real:: eps,epsp1
real, allocatable:: x(:),y(:),b(:,:),r(:,:)
real, parameter:: rad=3.1415926/180
print*, "please entry the N."
read*, N
allocate(F(1,1:N**2),A(1:N**4,1:N**4),x(1:N),y(1:N),b(1:N,1:N),r(1,N**2))
A=0
print*, " "
do i=1,N**2
  A(i,i)=4
  A(i+3,i)=-1
  A(i,i+3)=-1
end do
do i=2,N**2,3
  A(i-1,i)=-1
  A(i,i-1)=-1
  A(i+1,i)=-1
  A(i,i+1)=-1
end do
do i=1,N**2
  print "(10I5)",(A(i,j),j=1,N**2)
end do

!A

do i=1,N
  x(i)=i*(N+1)
  y(i)=i*(N+1)
end do
do i=1,N
  do j=1,N
    b(i,j)=13*(rad**2)*sin(3*rad*x(i))*sin(2*rad*y(j))*((1/N+1)**2)
  end do
end do
print*, " "
!b
do k=1,N
  do i=k+1,N
    do j=k+1,N
      A(i,j)=A(i,j)-A(k,j)*A(k,j)/A(k,k)
    end do
    b(i,i)=b(i,i)-b(k,k)*A(i,k)/A(k,k)
  end do
end do
do k=N, 1, -1
  do j=k+1, N
    b(k,k)=b(k,k)-A(k,j)*b(j,j)
  end do
  b(k,k)=b(k,k)/A(k,k)
end do
print*, b
print*, " "
!gauss

エラーメッセージは

Access Violation
the instruction at address 036238f0 attempted to write to location 00000000

 036238b4 __ALLOCATE_U [+003c]
 036ef495 get_non_checked_memory(int)#8F [+001a]
 036ef4db get_heap_block(int)#8F [+0056]
 036ef6ea __PALLOCATE2 [+008b]
 main -  in file freeformat1.f95 at line 11 [+0299]

です。

>>478
メモリーがアロケートできなかったんじゃないかな？

11**4〜10^4
A(1:N**4,1:N**4)だから、必要とされるメモリーは、4*10^4*10^4＜10^9Byte〜1Gbyte

メモリが足りないと言うことは、私のＰＣが原因なんですね・・・
なんとか改善策はありませんでしょうか？

>>476
正直もう Fortran は消えて欲しいのでやめてください。

N=2において、
b(1)=x(1)y(1)
b(2)=x(2)y(1)
b(3)=x(1)y(2)
b(4)=x(2)y(2)

N=3において、
b(1)=x(1)y(1)
b(2)=x(2)y(1)
b(3)=x(3)y(1)
b(4)=x(1)y(2)
b(5)=x(2)y(2)
b(6)=x(3)y(2)
b(7)=x(1)y(3)
b(8)=x(2)y(3)
b(9)=x(3)y(3)

ってな感じに配列を格納したいのですが、下のプログラムでは
例えばN=2の場合、
b(1)、b(2)、b(3)、b(4)全ての配列にx(2)y(2)のときの値が格納されてしまいました…。
どうしてでしょうか？

do k=1,N**2
 do i=1,N
  do j=1,N
    b(k)=13*rad**2*sin(3*rad*x(i))*sin(2*rad*y(j))*(1/N+1)**2
  end do
 end do
end do

>>482
よく見てないので、配列ミスないならメモリを節約するアルゴリズムの
本でも見なさい。

コンパイルオプションで配列ミスチェックできるなら付けてみる。

unix系でunlimitedしてみたら？

>>484
「どうしてでしょうか？」なので原因だけ書く
君がそう書いてるから、そうなる。
上の442か？

k=1の時、iもjも2まで回るから、そうなる。k=2の時も同様。
以下同じ

>>476
まずはヌける絵を描くんだ。
話はそれからだ。

>>476
DO「ああんっ……CONTINUE君っ……」
細い腰を震わせ、DOは歓喜の牝声をあげた。
奥深く突き刺さった行番号を、DOとIが招き入れるように吸い付いている。
「いいんっ……もっと……もっと下さい……」
好きなCONTINUEの行番号をDOとIの間に入れられることがこんなにも気持ちいいなんて。
CONTINUEの熱い肉体をもっと感じたくて、ついいやらしく激しい突きをせがんでしまう。

こういう感じ？

だがEND DOに寝取られる

FORTRANのイメージならホモゲーだよな

アッー！

>>478
A をアロケートしようとしているところが 4GB を超えてしまっているのが原因。
４GB以下になるように分割して扱う必要がある。

>>485>>491
アドバイスありがとうございます！
なんとか、N=11以上のとき動かなくなるエラーは亡くすことができたのですが、
>>484が未だに解決できないので、こちらの方のアドバイスも
大変恐縮ですが、アドバイスをよろしくお願いします。

>>492
紙と鉛筆を用意して、N=2の時 >>484のプログラムがどんな動きをするのか
1ステップずつ k, i, j, b の値を紙の上に書き出して確認してみましょう。

>>492
i と　j を k から作り出すように書き換えるのがたぶん簡単。

FORTRANを盛り上げていこうぜ！

お断りします

そう言われると、かえって>>496の菊穴にCOMMON文を無理やり押し込んでやりたいｗ　うふふ

最近フリーのGfortranを使い始めたんですが
サブルーチンの使い方がわかりません・・・．
LAPACKを使いたいんですけど，プログラムの始めに何か書くのでしょうか？
いきなりcall文で呼び出せないでしょうし．

ちなみに，windowsのコマンドプロンプトでコンパイルしています．
fortran初心者なものですいません・・・

>>498
本を読まずに1から全部聞いて済ませようとしているのか？

>>498
ヘッダファイル（インクルードファイル）がないのかというようなことを聞いているのでしょうが
ないんじゃないかな。
リンクするのを忘れなければいきなりcall できると思う。
f90(f95)なら　use lapack_f95 みたいなモジュールの読み込みだけはいると思う。

Fortran IVの方でも聞いていたのですが、返事がなかったのでこちらでも
質問させてください。

C(Not C++ / MSVC)で作られたライブラリの関数が使用したいのです。
そいつは STDCALL で作られていて、大文字小文字が混在しているものもあります。

それをg95またはgFortran もしくは g77 で呼び出したいと思っています。
_cdecl で 小文字なら呼び出せるのですが 呼び出し規約と関数名の両方について
対処方があるならば知りたいと思っています。

ラッパ関数を C で作るとか。

素早いお返事ありがとうございます。
やはりそれが普通でしょうか。
関数がわりかし多めな+ IVF から呼べるインターフェース類はあるので
再利用ができるとうれしいなという感じなのですが。

ラッパ関数生成スクリプトでも書くとか。

>>504
ですね。
あきらめて、書き始めておりました。

すいません。学校でFORTRANを習っているのですが、
解けない問題があるのでよろしくお願いします。


問:三角形の底辺A、高さHを入力し、
面積を求めて印刷するプログラムを作りなさい。
と、いう問題なのですが、TeraPadで
C toi
read(5,*)A
read(5,*)H
S=A*H/2
write(6,*)'TATE=',A
write(6,*)'YOKO=',H
write(6,*)'MENSEKI=',S
STOP
END

と、入力しコマンドプロンプトで
cd fortran
g77 -o ファイル名 ファイル名.forとして、
エンターを押すと、エラーメッセージは発生せず、
その後、ファイル名→エンターをすると、
いつも出ている FORTRAN〜が消えます。
計算は出来ませんでした。アドバイスお願いします。

>>506
g77　云々　以降の画面に出ている文字を全部コピペして書け。
これだけでは判断のしようも無い。
プログラム自体は間違ってないように見える。

DO文を用いて整数NからM（N<M）までに存在する3の倍数の和を求める
プログラムを作成したいのですが、うまい方法がわかりません。よければ実行例を教えてください。

自分なりに考えたのですが、この数列の初項と個数が分かれば解けそうなのですが、
Nに最も近くてN以上の3の倍数が初項というのを、
IF（MOD(N+1,3).eq.0）THEN
初項＝N+1
のようにするとNとN+2の時の場合分けをしなければならないし、個数に至っては、
末項/3-初項/3+1で出そうとすると、かなりの数のIFを使わなければならないので
頭がこんがらがってうまく作れませんでした。

>>506
実行するとコマンドプロンプトが修了してしまうという意味でしょうか？
だとすると、たぶんマシンが変だと思う。
違うマシンで試せるなら試してみて。

SUM=0
DO I=N,M
IF (MOD(I,3) .EQ. 0) THEN
SUM=SUM+I
END IF
END DO
じゃだめ？

>>510
上手くいきました。
DO文の中にIF文入れて簡単に表せるんですね。とても勉強になりました。
ご教授ありがとうございました。

FORTRAN77です。
一様乱数を１０個ずつ発生させ、その平均を計算する。
これをｎ回繰り返す。
計算した平均値をｎが１千回、１万回、１０万回、１００万回についてそれぞれ分布を求めよ。
という問題が分かりません。
乱数の発生は調べたら[０，１）区間の一様乱数が
　　　real function unifrd(i)
　　　integer l,c,t30
　　　real mu
　　　parameter (l=843314861,c=453816693,t30=2**30,mu=2.0**31)
　　　i=l*i+c
　　　if(i.lt.0) i=(i+t30)+t30
　　　unifrd=real(i)/mu
　　　end
で出るらしく、これを応用すれば出せそうですが、１０個ずつ発生させ、その平均を計算する
プログラムの書き方がよく分かりません。
幼稚な質問ですが、どなたか教えて下さい。
お願いします。

C　S仮
C　SUM合計

S=0
SUM=0
N=10
DO I=1,N
　乱数発生ルーチン(ここで乱数はRとする)
　S=S+R
CONTINUE

SUM=S/N

でNを返れ

用すませて思い出したけど
変数は
REAL*8
とした方が良いんじゃね？

あ、これで動くは動くが、変数名とか恥ずかしい

C　A平均
C　SM合計
REAL*8　A,SUM,R
A=0
SUM=0
N=10
DO I=1,N
　乱数発生ルーチン(ここで乱数はRとする)
　SUM=SUM+R
CONTINUE

A=SUM/N

わからないので助けてください；；

整変数I,J,Kをキーボードより与え、下記を実行するプログラムを実行せよ。
x=i^2/4 + j^3/27
y=x/(1+k^2)
z=(yの3乗根)klogk

条件　ｋあるいはｙの値が0以下の場合はｚを求める計算ができない、その場合z=(絶対値yの3乗根)(絶対値k)log(絶対値k)
を計算するように条件わけをすること。


IMPLICIT REAL(A-H,O-Z)
DATA I,J,K,X,Y,Z/5,6,1,2,1,6/

x=i**2/4+j**3/27
y=x/(1+k**2)
z=y**(1/3)*k*LOG(k)

WRITE(6,'(4I10)')I,J,K
WRITE(6,'(4F10)')X,Y,Z
WRITE(6,'(2I10,2F10)7)M,L,E,D
STOP
END

この辺までは組んでみたのですがこの後どうすればいいのでしょうか？
すみません、よろしくおねがいします><

>>516
>整変数I,J,Kをキーボードより与え、
READ文

>ｋあるいはｙの値が0以下の場合はｚを求める計算ができない、その場合z=(絶対値yの3乗根)(絶対値k)log(絶対値k)
>を計算するように条件わけをすること。
IF文

>>516
IF文でy=0 または z=0 のところを別にして
それ以外は　ｙとkをabs でくくる。
という風にやるのが期待されているんだろうけど、

y<=0 でも計算できる気がするぞ。問題間違ってない？

>>507さん レスありがとうございます。遅くなって申し訳ないですが、
もう一度やってみました。
>>506後
W:\fortran>g77 -o ファイル名 ファイル名.for

W:\fortran>ファイル名、と打ちエンターを押すと、
何も書かれてない下の行にカーソルが移動し、
エンターを押しても同じことで、もう一度ファイル名→エンターでは、
ファイル名
invalid number: Incomprehensible list input
apparent state: unit 5 (unnamed)
last format: list io
lately reading direct formatted external IO

abnormal program termination

とエラー?が出ます。
すいませんがよろしくお願いします。

>>519
２回数値を入力してEnterを押すと答えが表示される。
READ文の復習をしましょう。

あと、細かいことだが 'TATE', 'YOKO' と底辺A, 高さHの対応はそれでいいのか？

>>520さん、レスありがとうございます。
しかし、>>「２回数値を入力して」と、書かれていますが、
数値とはどこの数値をどこに入力すればよいのでしょうか?
何度もすいません…。

>>521
君はいったい何をするプログラムを作っているのかね？
もう一回問題をよく読みなさい。

>>522さん
すいません。ようやく分かりました。ありがとうございます。

恥ずかしながら、出力結果が S=AH/2 になるのかと思ってました。

迷惑おかけしました。すいませんでした。

すみません、使いたいソフトのソースコードがFortran90で書かれていますが、
コンパイラがありません。
しかもUNIXじゃなくてWindowsXPなんです。Cygwinはあるんですが。
というわけで、Fortran90のWindows用フリーコンパイラなどをご存じでしたら
教えてください。

g95 -minGW で決まり。

３×３の行列の積を解くプログラムがわかりません＞＜
ちなみにTextファイルから行列を読み込むプログラムです。
どなたかご教授願います。

宿題頑張ってね

fortranでフィボナッチ数列A0 = 0 , A1 = 1 , Ai+1 = Ai + Ai-1 (i =2,3...)を計算するプログラムを作っています。
得られた値が奇数の時には、数値の右横に＊をつけ、奇数かつ３の倍数なら＊＊をつけます。以下のようにプログラムしましたが、＊印が偶数時にもついてしまい、また＊＊の付け方にも困っています。どこが、間違っているのか、アドバイス手本よろしくお願いします。
　　　IMPLICIT REAL (A-H,O-Z)
A0 = 0
A1 = 1
A2 = 1
WRITE(6,*)'A 0 = 0'
WRITE(6,*)'A 1 = 1'
DO 20 L = 2,20,1
A2 = A1 + A0
I = A2-INT(A2/2)*2
J = A2-INT(A2/3)*3
IF(J.NE.O.AND.I.NE.O) THEN
WRITE(6,*) 'A',L,'=',A2,'*'
ELSE IF(J.NE.O.AND.I.EQ.O) THEN
WRITE(6,*) 'A',L,'=',A2,'*'
ELSE
WRITE(6,*) 'A',L,'=',A2
END IF
A0 = A1
A1 = A2
20 CONTINUE
STOP
END
　　結果
A 0 = 0　　　A 1 = 1　 A 2= 1.*　　 A 3= 2.*
A 4= 3.*　　　A 5= 5.*　A 6= 8.*　　　A 7= 13.*
A 8= 21.*　　　A 9= 34.*　A 10= 55.*　　　A 11= 89.*
A 12= 144.*　　　A 13= 233.*　A 14= 377.*　　　A 15= 610.*
A 16= 987.*　　　A 17= 1597.*　A 18= 2584.*　　　A 19= 4181.*
A 20= 6765.*

>>528
・FORTRANにはMODという便利な関数がある。
・「奇数」、「奇数かつ3の倍数」の判定条件が変。
・せっかくIF/ELSE IFで条件分けしているのにWRITE文がまったく同じ。
・プログラムの頭に「IMPLICIT NONE」をつけてみよう。

>>528
実数で計算する必要があるのか？

>>528
Ai+1 = Ai + Ai-1
を計算するなら、次のように配列を使う方が簡単だと思います。
A(i+1) = A(i) + A(i-1)

>>530
絶対にない


>>526
絶対に3x3なら式を９個並べろ
規則性が分かったらループ回せ
さすればkxn ・nxjでも簡単に計算できる

>>526
行列の積の定義を線形代数の本で見よ．
定義式をそのまま計算すればよろしい．

3x3が保障されてるなら532の言うように，式を9本かけばよい．

テキストファイルから行列を読み込む方法が分らないんじゃないかな。

open(2) "mat.txt"

C　ループ回すなら工夫して回せ

read(2,*) A(I,J)
read(2,*) B(K,L)

C calc statment

どーのこーの

close(2)
stop
end


I,JとK,Lで変えたけど分かってるならB(I,J)でも良い
open文でstatus=oldかunknownが必要なら書け。

つーか、そしたら質問者は何も理解してないな。

>>526
c = matmul(a, b)
じゃだめなの？

>>536
多分７７の課題なんでね？それでいいとすると、あまりにあまり（正笑

!(Φ＿Φ+){？？？}

Ai+1 = Ai + Ai-1
　 ↓↓↓
「全てに一致する」 i『option』 + 1 ≫『代入する』≪ 「全てに一致する」 i『option』 + 「全てに一致する」 i『option』 -1

A(i+1) = A(i) + A(i-1)
　 ↓↓↓
「全てに一致する」 (i『option』 + 1) ≫『代入する』≪ 「全てに一致する」 (i『option』) + 「全てに一致する」 (i『option』- 1)

!(Φ＿Φ+){？？？}

『「全てに一致する」(command) i{option}』 == 『文字列、配列順に一致を優先する + i{option}』

　　　　 〆{この様に？}

　　　Ж　　√ζ
　　　　　 　 {Ж}

do i=1,10
end do
の場合、このループを抜けた後iの値は11になるのですか？

yes

大変初歩的なこととは思いますが質問です

x : 1 　2　 3 　4　 5 ...
y : -0.1 -0.02 0.03 0.1 0.7 ...

のように，無数のxに対応するyの値を記録させるアルゴリズムがすでにあります．

ここで，yの符号が変わったときの，対応するxを記録する（ここでいうx=3）
命令を組みたいのですが・・・

どうもわかりません．
わかる方いれば教えていただけないでしょうか？

>>543
データ量が少ないなら、頭から１こづつ見ていっても対して変わらん。

それが嫌なら、単純なのは二分法。ぐぐればよろし。

×づつ
○ずつ

>>545
歴史的仮名遣いでは『づつ』の方が正しい。
Fortran使いは、古いほうを好むべし。

x-1とxのyを掛けて負ならば記憶させたら良いだけじゃ？

グローバル変数みたいなことしたいなら、commonじゃなくモジュール使えって言われますが
commonだとどの変がまずいんですか？正確に複写しないとダメ、という他に
速度などの面でデメリットとかあるんでしょうか？

>>548
COMMONは色々厄介なのよ。
メモリーの先頭番地を合わせるだけなので、まったく別の名前で別の並びでも成立する。
サイズが合わなくてもいい。（処理系に依存するが最初に出た大きさになることが多いと思う）

昔はメモリー節約のために、今で言う動的割付のようにプログラマーがCOMMONを
自己責任で様々に再利用した。プログラムの前半と後半で、同じCOMMONブロックを
まったく別の割付で用いるとか。COMMONの前半分は共通で、後ろ半分が次々変化するとか。

そういうプログラムは非常にデバッグが難しい。読まされるほうは地獄。

そういう過去の怨念が溜まっているので、COMMONはとても忌み嫌われているのだ。

単なるGlobal変数として使う分には、その恐ろしさが分からないだろうが、COOMONは
ほとんど無制限のメモリー領域の共有だと理解したほうが適切。

速度的にはむしろ早いんではないかと思ったりもする。

>>549
＞プログラムの前半と後半で、同じCOMMONブロックを
＞まったく別の割付で用いるとか。COMMONの前半分は共通で、後ろ半分が次々変化するとか。
昔はそんなことをしてたんですか・・・壮絶ですね

でも逆に考えると、割と単純な科学技術計算とかで、複雑なことはしないが
速度はひたすらに要求されるような場合は、あえてcommonを使ったりするんでしょうか？

今は普通に引数で渡してますが、これ変えたら速くなるかなぁ

あらゆるところで変更されうるから、
common に起因するバグがあった場合、
バグのあるとろを探すのが大変、
というのはよく言われるところ。

ただ、common を避ける事で引数がもの凄く多くなるとなると、
それはそれで問題になる事もあったり。

>>550
最適化の問題などもあるので、今となってはＣＯＭＭＯＮにする利益は無いのではないかと思う。
MODULEで渡しても変わらないのではないかと推測する。

FＯＲＴＲＡＮは参照渡し、すなわち変数のアドレスしか渡さないので、
変数のコピーを渡す値渡しと違って、それほど引数渡しのオーバーヘッドの問題は少ないと思う。

ただFortran90になってから、配列を受け渡すときに、サイズや次元の情報まで渡すようになったが・・
（FORTRAN７７までは配列の先頭要素の番地しか渡さない）。




Fortran90の配列引渡しのとき、x(:)とかx(5:10)などと : 指定子を使って部分配列を渡そうとすると
いまのFortranの実装では、呼び出し前に暗黙のダミー配列にコピーをつくり、それを渡して
戻ってきた結果を再び元の配列にコピーしている。したがってオーバーヘッドが大きくなる。

意味的には同じだが、xを渡す場合とx(:)を渡す場合だと、後者のほうが少し遅くなる。

巨大な配列は部分配列で受け渡さないほうがいい。

詳しい解説ありがとうございます。勉強になります。

いつまでもFORTRAN77ではさすがに時代遅れになるかも・・・という漠然とした不安から、
上のバージョンへの移行を考えていたのですが、よく検討したほうが良さそうですね。
Fortran95とか新しいやつはオブジェクト指向とかも取り入れてるんでしたっけ。
速度を求めてる人にとっては論外なんでしょうね。。

>>553
Fortran95はFORTRAN77を含んでいるので、中途半端に移行してもいいんで内科医。
部分配列の引数受け渡しが遅くなるといっても、よほどのことが無い限り差はでないだろう。

オブジェクト指向はFortran2000からで、Fortran95ならPASCALが半分混じったような感じなので、
８０年代とかに教育を受けた人には親しみやすいところがあると思う。

FORTRAN７７を長く使った後で、Fortran90/95を使うと、痒いところに手が届くというか、
こんな機能があったらな、と思っていたような機能が色々入っているので、学んで損は無いと思う。

IntelのコンパイラがVer10になった？

全日空の旧システムはＦＯＲＴＲＡＮで書かれていた！
ttp://itpro.nikkeibp.co.jp/article/JIREI/20070530/273051/

うｎ
レジストレーション・センターに置いてあった。

IntelやXLsoftからはまだアナウンスがないけど、Windows版はなにやら
Microsoft Visual Studio 2005 Premier Partner Edition Included
だそうで。

>>557
>Microsoft Visual Studio 2005 Premier Partner Edition Included

kwsk

>>558
ttp://www.intel.com/cd/software/products/asmo-na/eng/compilers/fwin/278834.htm
VSの簡易版(?)が付いて単体で使えるようになった。

あとPro版にMKLが付くようになった。

Product name Full product Support renewal
IntelR Visual Fortran Compiler, Standard Edition, for Windows* $599 $2801
IntelR Visual Fortran Compiler, Professional Edition, for Windows*
Includes IntelR Math Kernel Library $699 $280
IntelR Visual Fortran Compiler, Professional Edition, for Windows* with IMSL*
Includes IntelR Math Kernel Library and IMSL Fortran Library $1,599 $749

1Support renewal qualifies you for upgrade to the Professional Edition

ttp://www.intel.com/cd/software/products/asmo-na/eng/compilers/279613.htm

よく分からんが、いままでのスタンダード版の更新料が＄２８０に値上げになるが、
ＭＫＬも付くようになるという事でいいのかな？

fortran90って外部コマンド実行できますか？
perlで言うsystem()みたいな

ちなみにlinux環境です

何でとりあえずやってみようとしないの？

fortranで外部コマンドなんて使えるの？いったいどういう仕組みで・・・
perlとかなら分かるけど

なんで使えないと思うのかが理解できない。

Fortran77の宿題なんですが・・、下のような行列を絶対値の小さいものから1行目に順に
並べるプログラムをつくれというものなんですが、分かりません。どなたかよろしくお願いします。

‖2.0 1.0 -4.5‖
‖0.0 -0.5 1.1‖
　　　↓
‖0.0 -0.5 1.0‖
‖1.1 2.0 -4.5‖

>>560
スタンダード版は新規のみでサポート更新がなくなったからプロ版（のサポート更新）にアップグレードしろ
ということかと。

アップグレードで得られる権利：

各 プロフェッショナル・エディション に含まれる製品に対する使用権を所有します。

サポート期間が所有製品の現在の終了日より 1 年間更新されます。

例 ) インテル C++ コンパイラー 9.1 Windows 版 を 2006 年 10 月に購入した場合
( この場合、サポートサービスが 2007 年 10 月まで有効 ) 、インテル C++ コンパイラー
Windows 版プロフェッショナル・エディション SSR を購入することで、インテル C++
コンパイラー Windows 版プロフェッショナル・エディション に含まれるインテル C++
コンパイラーおよびすべてのライブラリー (IPP 、 MKL 、 TBB) のサポートサービスが
2008 年 10 月まで有効になります。

ttp://www.xlsoft.com/jp/products/intel/compilers/comp10_faq.html

よく和歌欄が、今すぐ更新すれば従来の更新期間１年延長で、ＭＫＬ付きの版に変わるという事か？

>>565
２行３列に見えるのは、俺の目の錯覚？

「ソート」で調べて絶対値を被せたら良いだけじゃ？

古いコンパイラで　-Nl50　というオプションの効果がわかる方がいれば教えてください。

>>569
コンパイラの種類は？バージョンは？

バージョンは分かりません。種類はFortranのコンパイラです。

>>571
お前はアホか！　エスパーだって答えられないぞ！
しょんべんで顔洗って出直せ！


コンパイラのヘルプを見れなり、バナーを出すなりして、それをこのスレにコピペしろ。
それも出来ないなら、せめてコンピュータ名やOSくらいは書け。

脳みそがたりないようだから、何も考えるな、言われたままやれ。

>571
まあ、とりあえず、そのコンパイラでコンパイルするときに打つコマンドは何？
frt とか g77 とか、そんなんでわかるばあいもあるし。

windows xp のコマンドプロンプトで
f90のコンパイルして実行するとstack overflowしちゃいます．
visual studio だとスタックサイズ増やせたんですが，
コマンドプロンプト上でやりたいんです．
どうしたらスタックサイズ増やせますか．

>>574
>>572-573

>>575
俺の頭が悪いのか理解できません
噛み砕くと？

事故解決しました．おさわがせしました．

>>574
解決結果を書いておくと、備忘録になり、他の人の助けにもなる。

まぁ一般的に言えば、大概のコンパイラでは /help とか /? で、オプション一覧が出るので
それを見れば解決するのだが・・・

じゃぁ。
自分ができた方法を・・・

普通にコンパイル
　　　　↓
editbin /stack:十分な値（デフォルトは１Mだから、それ以上の数字） program.exe
　　　　↓
program.exe

で、できました。もっと効率のよい方法もあるのかも、、、

>>579
コンパイラは何？
Visual Studio上で動くFortranコンパイラにしても色々あるわけだが

　　　　　　　　　　　　　 ry､
　　　　　　　　　　　　 / / }
　　　　　　　　　　　_/ノ.. /、
　　　　　　　　　　 /　 < 　 }
　 　 　 ry､ 　　　 ｛ｋ＿　_/`;,　　ﾉﾉ　ﾊﾟﾝﾊﾟﾝ
　　　 / / } 　　　　　;'　　　　 `i､　
　　 _/ノ../、　　　＿/ 入/ /　　 ｀ヽ,　ﾉﾉ
　 /　r;ｧ 　}''i"￣.　 ￣r'_ノ"'ヽ.i 　　）　―☆
　｛ｋ＿　_/,,.'　　;.　 :.　　　　　 l､ 　ﾉ　　
　　　 ＼　｀ 　､　　,i. 　　　.:, :, ' / /　＼
　　　　 ,;ゝr;,;_二∠r;,_ｪ=-ｰ'" r,_,/　　　☆


【ラッキーレス】
このレスを見た人はコピペでもいいので
10分以内に3つのスレへ貼り付けてください。
そうすれば１４日後好きな人から告白されるわ宝くじは当たるわ
出世しまくるわ体の悪い所全部治るわでえらい事です

>>580
まぁ、前に書いたんですが・・・・
f90です

( ﾟдﾟ)

>582
俺が不勉強なのかもしれんがf90というコンパイラは聞いたことがない。
そのソフト買ったとき、箱に f90 って書いてあった？
なんていうメーカーの製品？

その製品を知ってれば答えられるが知らないなら答えられないというなら、
その製品が何かを聞き返すこと自体、意味がない気がしないでもない。

f90 は汎用コマンド名であって、コンパイラ名ではないな。

任意の文字列C、整数Nを入力して、これらを連結してファイル名に使いたいのですが、
どうやればできるでしょうか？

>>587
WRITE文を使えば簡単ですよ

PROGRAM test
  IMPLICIT NONE
  CHARACTER(20) :: str1 = "test", str2
  INTEGER :: i
  DO i=1,10
    WRITE(str2,'(A,I3.3,".dat")') trim(str1), i
    WRITE(*,'(A)') trim(str2)
  END DO
END PROGRAM test

出力結果

test001.dat
test002.dat
test003.dat
test004.dat
test005.dat
test006.dat
test007.dat
test008.dat
test009.dat
test010.dat

ありがとうございます。trim関数ですか
ただFORTRAN77なんですよね・・・

>>589
PROGRAM HOGE
IMPLICIT NONE
CHARACTER*10 C
CHARACTER*20 STR
INTEGER N, I
READ(*,*) C
READ(*,*) N
I = INDEX(C,' ') - 1
IF (I .LE. 0) I = LEN(C)
WRITE(STR, '(A, I3.3, ''.dat'')') C(1:I), N
WRITE(*,*) STR
STOP
END

>>587
またこの質問かｗｗｗｗ
次からテンプレに入れておけ！（正笑

>>590
知りませんでした・・・＞INDEX
ありがとうございました

//でくっ付けろよ。

>>593
まぁ、内部ファイルでやると１行で出来るからな。

N行N列の逆行列を求めるプログラムの作り方を教えてください・・・
ちなみにfortran90です
LAPACK使えって言うのは無しで・・・

>>595
線形代数の行列の所で、掃き出し法によって逆行列を求める方法をやったでしょ。
それをそのままプログラムにすればOK、習って無くても、どの線形代数の本にもに載ってる
と思うからから、調べてみて。数値計算だとガウス・ジョルダン法と呼ばれてる方法かな？

LU分解
掃き出し法
スカイライン

色々あるよ。つーか数値計算の本読め。

返事が遅くなって申し訳ない
助言ありがとうございます
とりあえず、自力でまずはプログラム化してみます

教えて下さい。
H(l)を以下のように計算したいのですが、走らせても答えが「NaN」と表示されます。
どこがおかしいか教えて下さい。お願いします。

program FILTER
c
implicit real*8(a-h,o-z)
dimension H(100000)
open(115,file='H.dat')
c
do 10 l=1,n/2+1
f=df*(l-1)
e=f/10.0
H(l)=((1.0/f)**0.5)*((1.0+(0.694*(e**2))+(0.241*(e**4))
&+(0.0557*(e**6))+(0.009664*(e**8))+(0.00134*(e**10))
&+(0.000155*(e**12)))**(-0.5))*((1.0-(dexp(-(f/0.5)**3)))**0.5)
write(115,*) H(l)
10 continue
c
pause
stop
end

（注）＆は数式が長くなったため、改行している事を意味しているようです。

>>599
n の値が不定
df の値も不定

＞＞600
ありがとうございました。

>>599
暇なんで書き換えてみた。
i=1 のときは、f がゼロだから 1/f が計算できないからパス。

      program filter
      implicit integer(i-n)
      implicit real*8(a-h,o-z)
      real*8 ek(1:6)
      ek = (/ 0.69400d0, 0.24100d0, 0.055700d0,
     &        0.09664d0, 0.00134d0, 0.000155d0 /)
      df = 0.01d0
      n  = 20
      do 1 i = 2, n / 2 + 1
        f = df * (i-1)
        ee = f * f * 0.01d0
        ex = 1.d0
        ss = 1.d0
        do 2 j = 1, 6
          ex = ex * ee
          ss = ss + ex * ek(j)
2       continue
        sx = 1.d0 - dexp(-(f+f)**3)
        H = dsqrt( sx / (f*ss) )
        write(*,*) i, H
1     continue
      stop
      end

たぶんファイルからNを読み。
H.datの中のN個分のデータからH(I)を作りたいじゃ？

もとのプログラムではHを配列にする必要ないけど。

dexp(-(f/0.5)**3)のf/0.5は　2*fの方が見やすいと思うけど。
元の定義式がこうなってるのかな？

複素数
C=(1.d0/3.d0,0.d0)のように分数などを打つとエラーが出るのはなんでですか？
ちなみにC=(dlog(2.d0),0.d0)もエラーがでました。

>>604
cmplx()にしたらコンパイルできる。
次に質問するときはエラーメッセージくらいは
書いたほうが良いぞ。

倍精度にしたそうだから dcmplx じゃないかな。
それともどっちでも自動的に補正してくれたっけ？

>>604
ダメな理由はしらんが、

C = cmplx(1.d0/3.d0, 0.d0)
でＯＫ。

>>605
>>606
>>607
ありがとうございました。
はじめにcomplex*16::c(100,100)と宣言しているので、そのままいけると思ってました。

要するに (***, ***) で１つのリテラルだから、
その中に演算子とか（符号以外）入れれないんだろう。

質問させていただきます。
以下のプログラムで、y1とy3は答えとともにオープンしてくれるのですが、y2のオープンしたファイルが何も表示されません。
サブルーチン文はy1y2y3ともに引数は表示してくれてます。何か問題がありましたら教えて下さい。
お願いします。
program IFFT
implicit real*8(a-h,o-z)
dimension y1(1000000),y2(1000000),y3(1000000)
complex*8 w1(1000000),w2(1000000),w3(1000000)
dimension H(1000000),f(1000000)
c
open(114,file='d.dat')
open(115,file='e.dat')
open(116,file='f.dat')
c
m=15.d0
dt=0.01
n=2.0**m
tp=n*dt
df=1.0/tp
c
call SPECFILT(w1,w2,w3)
do 21 j=n/2+2,n
w1(j)=conjg(w1(n-j+2))
w2(j)=conjg(w2(n-j+2))
w3(j)=conjg(w3(n-j+2))
21 continue
c

call nlogn(m,w1,1.0)
do 24 k=1,n
y1(k)=real(w1(k))/dt
t=dt*(k-1)
write(114,*) y1(k)
24 continue
call nlogn(m,w2,1.0)
do 25 k=1,n
y2(k)=real(w2(k))/dt
t=dt*(k-1)
write(115,*) y2(k)
25 continue
call nlogn(m,w3,1.0)
do 26 k=1,n
y3(k)=real(w3(k))/dt
t=dt*(k-1)
write(116,*) y3(k)
26 continue
c
stop
end

>>610-611
入力ファイルが空なんだろう。

だいたい３つのファイルが同じ処理をしているなら、一つ用のサブルーチンを作って
それを３回呼ぶようにすれば、このような疑問は生じない。

どうせ同じことだと横着して、同じ処理を３回書くから、どこかで間違ったかもと余計な時間を食う。

反省しろ！　Fortranの神様に謝罪と賠償しろ！

612>
それが入力ファイルもオープンしたらちゃんと値が表示されていたんですよ。

>>610
質問と関係ないけど、

m=15.d0
dt=0.01
n=2.0**m
tp=n*dt
df=1.0/tp

これ気持ち悪い。整数と倍精度実数はっきり区別して書く。

      m  = 15
      n  = 2**m
      dt = 1.0d-2
      tp = dble(n) * dt
      df = 1.0d0 / tp

f77も複素数型も使うことないので、あやふやだけど、
complex*8 は、実部と虚部の合計で８バイトで、実部と虚部それぞれは４バイト（単精度実数）、
complex*16は、実部と虚部の合計で１６バイトで、実部と虚部それぞれは８バイト（倍精度実数）じゃなかったっけ？
y1-y3 は real*8 なんだから、complex*16 のが良い気がする。

ちなみに、F90 以降の complex(8) みたいな宣言の方法だと、
実部・虚部それぞれが８バイトだったような気がする。

Microsoft PowerstationをWindowsXpにインストールしてFortranを使おう
と思ったのですが、コンパイル時に
Error executing fl32.exe.
Text1.obj - 1 error(s), 0 warning(s)
とエラーが出て先に進めません。
どうすれば解決するのでしょうか
よろしくお願いします

>>615
どの文でエラーが起こったか分からないと
答えようがないんで…
文法間違いでないならサンプルで試してみたら
良いのでは？

>>611

do 25 k=1,n
の直前にnの値が何か書くように汁。

>>615
PowerStationはバグがあるから、新しいコンパイラ買え。

DO 100 L = 1 , 5
READ(IOT)
& (((FTR(I,J,K,L),I=1,A),J=1,B),K=1,C)
100 CONTINUE

というプログラムで、値は全てファイルから読み込ませます。
予めAとBとCは読み込みをさせています。
それでこのプログラムを実行させようとすると

input statement requires too much data, unit200, file "ファイル名"

と表示されて実行出来ません。どうやらデータの読みすぎのようなのですが、
心当たりがありません。
考えられる原因は何があるのでしょうか。

A,B,Cは200×100×50くらいの値です。

>>619

ttp://www.xlsoft.com/jp/products/intel/cvf/docs/vf-html/em/em01_03.htm

>evere (67): Input statement requires too much data
>FOR$IOS_INPSTAREQ。書式なし READ 文または PAD='NO' で開かれたファイルからの書式付き順番 READ 文で既存の記録以上のデータを読み取ろうとしました。

>>619
A*B*C〜1,000,000〜１M
単精度なら４倍して４M

データファイルの大きさがこれくらいあるか？

まぁファイルの読み取りが、ここだけとしての話だがｗ

ありがとうございました

ハードウェアで４倍精度演算をサポートしてないと劇遅だよな

end file premeture
で終わるタイプと違う？ただのエラー表示の違い？

>>623
当然だろｗ
コプロが無かった時代のパソコンでの計算を思い出せｗ

初歩的な質問でしたら申し訳ないのですが、調べてもわからないので質問させて下さい。
fortranで作ったプログラムを実行ファイル(EXE形式)にするにはどうすればいいのでしょうか。
（プログラムを実行するたびにfortranを起動して[コンパイル]→[ビルド]→[実行]するのは面倒なので、
　.exeをダブルクリックするだけにしたいのです）

>>626
一度、ビルドしてたら、実行ファイルできてないですか？
フォルダが作成されてたらその中とかに。

>>627
「Debug」フォルダの中に実行ファイルできてました。
どうやらアホな質問をしてしまったようで申し訳ないです。
どうもありがとうございました。

>>628
練習用の短いプログラムならいいけど、本番用や他のマシンに持って行くプログラムは
ちゃんと「Release」でビルドして使え。

質問させてください。
OPEN文の中にDO文を適用するのは可能でしょうか？
例えば、
OPEN(10,FILE='TEST1.TXT',STATUS='UNKNOWN')という文章の中の
↑ ↑
矢印で指した数値にDO文を適用して、たくさんのデータを
開きたいと考えています。
例えば
DO 1 I=1,N
OPEN(I,FILE='TESTI.TXT',STATUS='UNKNOWN')
1 CONTINUE
というような感じです。エロい人、宜しくお願いします。

>>630
>>587-590

>>630
もう千回くらい既出。

エクセルから数値を読み込んで処理させようとしたのですが、
うまく読み込んでくれません。ちなみに2次元ベクトルです。
ファイル形式は.csvでいいのでしょうか？

open(10,file='〜.csv')
do i=1,10
read(10,*) (A(i,j),j=1,10)
end do
close(10)

という感じでしました。
よろしくお願いします。

>>633
どう「うまく読み込んでくれ」ないのか詳しく。
全く読めないのか、特定のデータが来ると読めないのか、ランタイムエラーで止まるのか・・・

>>634
write文のように数字が表示され、止まってしまうんです。
同じ形式であと3つほど書いたので、たぶん1つ目と同じで止まると思います。

おまえの頭では分かっていても他人に説明できてないのを理解しろ

A(i,j)は２次元ベクトルでないと言うのはOKか？

>>633では数値データは100個分だけだよな？
i,jの成分で200個有るとかじゃないよな？

エラーメッセージ書け

今から学校に行ってまた送ります。

>>635
EXCELのCVS形式の読み書きの話はここレスがあったから、まず探して鯉！

>>635
EXCELのCVS形式の読み書きの話はここレスがあったから、まず探して鯉！