質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
86.12%
FORTRAN

FORTRAN(フォートラン)は科学時術計算に向いた手続き型プログラミング言語です。 並列計算の最適化が行いやすい特性上、数値予報および気候モデルなどの大規模な計算を行う分野のスーパーコンピュータで使われています。

解決済

fortranで正規乱数を発生させるプログラムを実行するとエラーが出る

ochikonbu
ochikonbu

総合スコア2

FORTRAN

FORTRAN(フォートラン)は科学時術計算に向いた手続き型プログラミング言語です。 並列計算の最適化が行いやすい特性上、数値予報および気候モデルなどの大規模な計算を行う分野のスーパーコンピュータで使われています。

1回答

0リアクション

0クリップ

395閲覧

投稿2022/07/18 17:31

前提

fortranでボックス=ミューラー法による正規乱数を発生させたいのですが、ネットに載っていた以下のコードを私の環境で実行するとエラーが発生してしまいました。
エラーの内容は「データ型の誤り」なのですが、何をどう改善すれば良いのかわからなかったため、質問させていただきました。
初歩的な質問ですみませんがご回答よろしくお願いします。

 

発生している問題・エラーメッセージ

エラー: boxm.f90(38): 組込みRANDOM_SEEDへの引数SIZEで誤ったデータ型が指定されています エラー: boxm.f90(41): 組込みRANDOM_SEEDへの引数PUTで誤ったデータ型が指定されています

該当のソースコード

!**************************************************** ! 正規乱数生成(by ボックス=ミューラー法) ! : 0 〜 20 の正規乱数をボックス=ミューラー法により ! 計算し、ヒストグラムを出力する。 ! ! date name version ! 2018.12.03 mk-mode.com 1.00 新規作成 ! ! Copyright(C) 2018 mk-mode.com All Rights Reserved. !**************************************************** ! module const implicit none ! SP: 単精度(4), DP: 倍精度(8) integer, parameter :: SP = kind(1.0) integer(SP), parameter :: DP = selected_real_kind(2 * precision(1.0_SP)) integer(SP), parameter :: M = 10 ! 平均 real(DP), parameter :: S = 2.5_DP ! 標準偏差 integer(SP), parameter :: N = 10000 ! 発生させる乱数の個数 real(DP), parameter :: PI = 4.0_DP * atan(1.0_DP) ! 円周率 real(DP), parameter :: SC = N / 100.0_DP ! ヒストグラム用スケール end module const module box_muller use const, only : SP, DP, S, M, PI implicit none private public :: rnd_seed, rnd contains ! 乱数の種の設定 subroutine rnd_seed implicit none integer(SP) :: seed_size, clock integer(SP), allocatable :: seed(:) call system_clock(clock) call random_seed(size=seed_size) allocate(seed(seed_size)) seed = clock call random_seed(put=seed) deallocate(seed) end subroutine rnd_seed ! 正規乱数生成 ! ! :param(out) integer(4) r(2) subroutine rnd(r) integer(SP), intent(out) :: r(2) real(DP) :: r_u(2) ! [0, 1) の一様乱数 2 個 call random_number(r_u(1)) call random_number(r_u(2)) r(1) = int(S * sqrt(-2 * log(r_u(1))) * cos(2 * PI * r_u(2)) + M) r(2) = int(S * sqrt(-2 * log(r_u(1))) * sin(2 * PI * r_u(2)) + M) end subroutine rnd end module box_muller program rndnum_bm use const, only : SP, N, M, SC use box_muller implicit none integer(SP) :: i, j, r(2), hist(0:M * 2) ! 乱数の種の設定 ! (正規乱数生成に使用する一様乱数の種) call rnd_seed ! 正規乱数の生成 hist(:) = 0 do i = 1 , N call rnd(r) hist(r) = hist(r) + 1 end do ! 結果出力 do i = 0, M * 2 write (*, '(I3, ":", I4, " | ", A)') & & i, hist(i), repeat("*", int(hist(i) / SC)) end do end program rndnum_bm

以下のような質問にはリアクションをつけましょう

  • 質問内容が明確
  • 自分も答えを知りたい
  • 質問者以外のユーザにも役立つ

リアクションが多い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

気になる質問をクリップする

クリップした質問は、後からいつでもマイページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

下記のような質問は推奨されていません。

  • 間違っている
  • 質問になっていない投稿
  • スパムや攻撃的な表現を用いた投稿

適切な質問に修正を依頼しましょう。

2022/07/18 23:18

こちらの質問が他のユーザーから「過去に投稿した質問と同じ内容の質問」という指摘を受けました。

2022/07/19 01:14

こちらの質問が他のユーザーから「やってほしいことだけを記載した丸投げの質問」という指摘を受けました。

まだ回答がついていません

会員登録して回答してみよう

アカウントをお持ちの方は

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
86.12%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問

同じタグがついた質問を見る

FORTRAN

FORTRAN(フォートラン)は科学時術計算に向いた手続き型プログラミング言語です。 並列計算の最適化が行いやすい特性上、数値予報および気候モデルなどの大規模な計算を行う分野のスーパーコンピュータで使われています。