大きなN. Common Lispの結果が不正なコード

Question

以下のコードは間違った答えです。これは、0と1の間の多くの乱数を持つ配列の平均である約0.5を与えるはずです。問題は、Nが "大きい"か、おそらく生成された乱数の精度ですか？このコードは、N（10^7,10^6など）の小さな値に対してうまく機能します。いくつかのアドバイスが参考になります。

ありがとうございます。

(defun randvec(n) 
    (let ((arr (make-array n))) 
     (dotimes (i n) 
      (setf (aref arr i) (random 1.0)) 
     ) 
     arr 
    ) 
) 

(defparameter N (expt 10 8)) 
(setf *random-state* (make-random-state t)) 
(defparameter vector1 (randvec N)) 
(format t "~a~%" (/ (reduce #'+ vector1) (length vector1))) 

Answer 1

あなたは単精度浮動小数点数で計算されている浮動小数点数

の精度。すべての乱数を足し合わせることで、1つの浮動小数点数が得られます。追加する数値が多いほど、浮動小数点数は大きくなります。結果的に精度が足りなくなります。

1.0d0のようなダブルフロートは、1.0s0のようなシングルフロートよりも高い精度を持っています。デフォルトでは、1.0はシングルフロートとして読み込まれます。 (RANDOM 1.0d0)はdouble floatを計算します。

(defun randvec (n) 
    (let ((v (make-array n))) 
    (dotimes (i n v) 
     (setf (aref v i) (random 1.0d0))))) ; create a double float random number 

(defun test (&optional (n 10)) 
    (setf *random-state* (make-random-state t)) 
    (let ((v (randvec n))) 
    (/ (reduce #'+ v) (length v)))) 

例：

CL-USER 58 > (test (expt 10 8)) 
0.4999874882753848D0 

スタイル

のCommon Lispでプログラミングする際の一般的なLispのプログラミングスタイルを使用してください：

• 、グローバル変数をそうでない場合は使用しないでください必要。代わりにローカル変数で関数を記述してください。
• defparameterでグローバル変数を定義する場合は、nではなく、*n*という名前を付けます。
• 形式でコードを正しくインデントしてください。インデントはエディタのヘルプで行う必要があります。
• 括弧を使用しないでください。

上記の私の例を参照してください。