read文を用いて input datから読み込む（外積計算fortran)

Question

下のプログラムに外積の計算のプログラムのコードがあります。
それで下のプログラムで２点の座標値を
read文を使ってファイルのinput.datから読み込みなさいと言われました。

program gaisekikeisann
	implicit none
	 integer, parameter :: vsize=3                             
	 real(8) vector1( vsize ), vector2( vsize )
	 real(8) gaiseki( vsize )
	 
	read(*, *) vector1( 1 ), vector1( 2 ), vector1( 3 )							!3つのベクトルの読み込み
	read(*, *) vector2( 1 ), vector2( 2 ), vector2( 3 )
	
	 !gaiseki(外積）の計算
	 gaiseki( 1 ) = vector1( 2 ) * vector2( 3 ) - vector1( 3 ) * vector2( 2 )
	 gaiseki( 2 ) = vector1( 3 ) * vector2( 1 ) - vector1( 1 ) * vector2( 3 )
	 gaiseki( 3 ) = vector1( 1 ) * vector2( 2 ) - vector1( 2 ) * vector2( 1 )
	 write(*,*) gaiseki(1), gaiseki(2), gaiseki(3)

end program gaisekikeisann

それで

    read(*, *) vector1( 1 ), vector1( 2 ), vector1( 3 )                            !3つのベクトルの読み込み
    read(*, *) vector2( 1 ), vector2( 2 ), vector2( 3 )

これを何かしら変えると思うのですが、
input.datを読み込めというのが全く分かりません。
どのように手順でプログラムコードや何かしらの操作をして、
read文を使ってファイルのinput.datから読み込むのですか。
回答お願いします。

Answer 1

科学計算を行う際、Fortranを使用することが非常に多いですが、
その時に使用する値の大半が、端末からの直接入力ではなく、
外部ファイル（.csvや.dat,.txtファイルなど）から値を読み込むものです。

外部ファイルから読み込みを行うために、まずはファイルをオープンする必要があります。

現状のコードですと、

read(*, *) vector1( 1 ), vector1( 2 ), vector1( 3 )
```となっているため、端末からの標準入力のみとなっています。

そこで、このように書き換えると値を読み込むことができます。
```fortran
integer, parameter :: points=2, vsize=3                             
real(8) vector( points, vsize )
open(10, file=input.dat)
  do i = 1, 2
    read(10, *) vector((i,j), j=1,3)) 
  end do
close(10)

とすれば、２つのベクトルを１つに２次元配列にできますので、計算も楽になりますよ^^

Comments

[mercurian-teto]

メモ帳で何も書いていないものを

input.datで保存したほうがいいですか。

[TakuyaHidaka]

メモ帳には
1,2,3
4,5,6
と書いておくだけで大丈夫です。

Answer 2

単純にこれだけでは？

program gaisekikeisann
    implicit none
     integer, parameter :: vsize=3                             
     real(8) vector1( vsize ), vector2( vsize )
     real(8) gaiseki( vsize )
     open(10, file='input.dat', status='old')

    read(10, *) vector1( 1 ), vector1( 2 ), vector1( 3 )                            !3つのベクトルの読み込み
    read(10, *) vector2( 1 ), vector2( 2 ), vector2( 3 )
    close(10)

     !gaiseki(外積）の計算
     gaiseki( 1 ) = vector1( 2 ) * vector2( 3 ) - vector1( 3 ) * vector2( 2 )
     gaiseki( 2 ) = vector1( 3 ) * vector2( 1 ) - vector1( 1 ) * vector2( 3 )
     gaiseki( 3 ) = vector1( 1 ) * vector2( 2 ) - vector1( 2 ) * vector2( 1 )
     write(*,*) gaiseki(1), gaiseki(2), gaiseki(3)

end program gaisekikeisann

Comments

[mercurian-teto]

input.datというファイルは作るのでしょうか。

[_Victorique__]

プログラムがあるフォルダにinput.datファイルを作成してください。

形式は以下です。

1 2 3
4 5 6

空白区切で3つ、改行で2行
その他余計なデータは書き込まないでください。