From 5e9eea55ad74871162c07b0f50b36d16f5259c97 Mon Sep 17 00:00:00 2001 From: hyh Date: Fri, 28 Mar 2025 15:20:10 +0800 Subject: [PATCH 1/2] update sp800_22_non_overlapping_template_matching_test.py --- sp800_22_non_overlapping_template_matching_test.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/sp800_22_non_overlapping_template_matching_test.py b/sp800_22_non_overlapping_template_matching_test.py index 76ee3ed..7666d5c 100644 --- a/sp800_22_non_overlapping_template_matching_test.py +++ b/sp800_22_non_overlapping_template_matching_test.py @@ -105,9 +105,9 @@ def non_overlapping_template_matching_test(bits): chisq = 0.0 # Compute Chi-Square for j in range(N): - chisq += ((W[j] - mu)**2)/(sigma**2) + chisq += ((W[j] - mu)**2)/(sigma) - p = gammaincc(N/2.0, chisq/2.0) # Compute P value + p = gammaincc((N-1)/2.0, chisq/2.0) # Compute P value success = ( p >= 0.01) return (success,p,None) From 2b2f5572034c9c98d41e272f7fad0d488cff6e52 Mon Sep 17 00:00:00 2001 From: foreverstream <90609287+goonstream@users.noreply.github.com> Date: Sun, 30 Mar 2025 14:19:23 +0800 Subject: [PATCH 2/2] Update sp800_22_non_overlapping_template_matching_test.py --- sp800_22_non_overlapping_template_matching_test.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/sp800_22_non_overlapping_template_matching_test.py b/sp800_22_non_overlapping_template_matching_test.py index 7666d5c..9cc085a 100644 --- a/sp800_22_non_overlapping_template_matching_test.py +++ b/sp800_22_non_overlapping_template_matching_test.py @@ -101,13 +101,13 @@ def non_overlapping_template_matching_test(bits): W.append(count) mu = float(M-m+1)/float(2**m) # Compute mu and sigma - sigma = M * ((1.0/float(2**m))-(float((2*m)-1)/float(2**(2*m)))) + variance = M * ((1.0/float(2**m))-(float((2*m)-1)/float(2**(2*m)))) chisq = 0.0 # Compute Chi-Square for j in range(N): - chisq += ((W[j] - mu)**2)/(sigma) + chisq += ((W[j] - mu)**2)/(variance) - p = gammaincc((N-1)/2.0, chisq/2.0) # Compute P value + p = gammaincc(N/2.0, chisq/2.0) # Compute P value success = ( p >= 0.01) return (success,p,None)