randogoth/onod3000
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

got rid of redundant sanity test

Browse source
This commit is contained in:
randogoth 2025-01-17 13:48:12 +02:00
parent a8c74caf90
commit f3766124ad
3 changed files with 6 additions and 44 deletions
Download patch file Download diff file Expand all files Collapse all files

1
src/uniformity/mod.rs
View file

@ -10,7 +10,6 @@ pub mod pi;
pub mod prediction;
pub mod runs;
pub mod runups;
pub mod sanity;
pub mod shannon;
pub mod shells;
pub mod uncorrelation;

36
src/uniformity/sanity.rs
View file

@ -1,36 +0,0 @@
use crate::Onod;
impl Onod {
/// Sanity randomness test
/// Checks for basic properties of randomness and returns a p-value.
pub fn sanity(samples: &[u8]) -> (f64, f64, f64) {
if samples.is_empty() {
return (-1.0, 0.0, 1.0);
}
// Check the proportion of ones and zeros in the byte data
let mut one_bits = 0;
let mut total_bits = 0;
for &byte in samples {
one_bits += byte.count_ones();
total_bits += 8;
}
let observed_ratio = one_bits as f64 / total_bits as f64;
let expected_ratio = 0.5; // Expected for a truly random sequence
let std_dev = (0.5 * 0.5 / total_bits as f64).sqrt(); // Standard deviation for binomial distribution
// Z-score calculation
let z_score = (observed_ratio - expected_ratio) / std_dev;
// Use normal distribution to calculate p-value
use statrs::distribution::{Normal, ContinuousCDF};
let normal_dist = Normal::new(expected_ratio, std_dev).expect("Failed to create Normal distribution");
let p_value = 1.0 - normal_dist.cdf(observed_ratio);
(observed_ratio, z_score, p_value)
}
}