2013-06-01から1ヶ月間の記事一覧

配列構成子の書き方

なじみの薄い配列構成子ですが、うまく使えばコードが綺麗になりそうだったので、サンプルコードを書いてみた。 program main implicit none integer, pointer :: arr_p(:) integer :: arr(3) integer, parameter :: n = 3 allocate(arr_p(n)) arr_p = (/ 1,…

Fortran特有の行列演算の書き方

Fortran90で行列を扱ってみます。配列ではそれぞれの配列要素の演算を簡単な記述で表すことができます。main.f90 program main integer(kind=8), pointer :: matrix(:,:) integer(kind=8), parameter :: nele = 3 allocate(matrix(nele, nele)) call set_mat…