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

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

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

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

Q&A

解決済

1回答

1067閲覧

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

ochikonbu

総合スコア2

FORTRAN

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

0グッド

0クリップ

投稿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

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

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

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

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

fortran

1 ! SP: 単精度(4), DP: 倍精度(8) 2 integer, parameter :: SP = kind(1.0)

SP は単精度実数の種別となっていますが integer(SP) のように interger の種別に設定しているところが意味的に正しくないです。コメント内ですでに指摘されているように単に integer としましょう。

元のプログラムの誤りだと思います(環境によっては動作します)。それはともかく、プログラムの出典は質問内でちゃんと明示しましょう。

投稿2022/07/21 02:38

chirimen

総合スコア189

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

ただいまの回答率
85.54%

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

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

質問する

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

FORTRAN

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