diff --git a/src/observation.rs b/src/observation.rs index a0e0fb5..11b3bda 100644 --- a/src/observation.rs +++ b/src/observation.rs @@ -32,7 +32,7 @@ pub(crate) struct Core { pub indices: Vec, pub ns_cav: Vec, pub xs_cav: Vec, - pub t: f64, + pub _t: f64, pub logpart: f64, pub exp_ll: f64, } @@ -49,7 +49,7 @@ impl Core { .collect(), ns_cav: (0..elems.len()).map(|_| 0.0).collect(), xs_cav: (0..elems.len()).map(|_| 0.0).collect(), - t, + _t: t, logpart: 0.0, exp_ll: 0.0, } diff --git a/src/storage.rs b/src/storage.rs index 09c7347..526848b 100644 --- a/src/storage.rs +++ b/src/storage.rs @@ -22,6 +22,7 @@ pub trait Storage { } */ +#[derive(Default)] pub struct Storage { keys: HashMap, items: Vec, @@ -67,12 +68,3 @@ impl Storage { self.keys.keys() } } - -impl Default for Storage { - fn default() -> Self { - Self { - keys: HashMap::new(), - items: Vec::new(), - } - } -} diff --git a/src/utils.rs b/src/utils.rs index a7a4834..3c675a9 100644 --- a/src/utils.rs +++ b/src/utils.rs @@ -88,6 +88,7 @@ pub fn logphi(z: f64) -> (f64, f64) { } } +/* pub fn logsumexp2(xs: &[f64], bs: &[f64]) -> f64 { let a = xs.iter().max_by(|a, b| a.partial_cmp(b).unwrap()).unwrap(); @@ -98,6 +99,7 @@ pub fn logsumexp2(xs: &[f64], bs: &[f64]) -> f64 { .sum::() .ln() } +*/ pub const ROOTS_HERMITENORM_XS: [f64; 30] = [ -9.706236,