diff --git a/CHANGELOG.md b/CHANGELOG.md index 04fa905c..6d11df23 100644 --- a/CHANGELOG.md +++ b/CHANGELOG.md @@ -1,6 +1,12 @@ # Changelog The semantic versioning is kind of random. +## 4.3.3 +- create a new type for the probabilities +### BugFix +- saving timepoints when the simulations run for only 1 year +- remove the factor 2 for the fit variants introduced by `v4.2.0` (not clear why we shouldn't divide by 2) + ## 4.3.2 ### BugFix - Do not multiply by tau when converting the rate per year into rate per division diff --git a/Cargo.lock b/Cargo.lock index 9559b283..048a0b2b 100644 --- a/Cargo.lock +++ b/Cargo.lock @@ -311,7 +311,7 @@ checksum = "95505c38b4572b2d910cecb0281560f54b440a19336cbbcb27bf6ce6adc6f5a8" [[package]] name = "hsc" -version = "4.3.2" +version = "4.3.3" dependencies = [ "anyhow", "chrono", diff --git a/Cargo.toml b/Cargo.toml index d2ede160..55cc27c0 100644 --- a/Cargo.toml +++ b/Cargo.toml @@ -1,6 +1,6 @@ [package] name = "hsc" -version = "4.3.2" +version = "4.3.3" edition = "2021" [dependencies] diff --git a/src/lib.rs b/src/lib.rs index 9c50ae81..ebfc6d08 100644 --- a/src/lib.rs +++ b/src/lib.rs @@ -87,7 +87,7 @@ impl Probs { mu_division, mu, }; - let mut u = mu / (cells as f32); + let u = mu / (cells as f32); let probs = if (asymmetric - 0.).abs() > f32::EPSILON { Probs::Asymmetric { u, @@ -95,13 +95,11 @@ impl Probs { asymmetric, } } else { - u *= 0.5; Probs::Symmetric { u, probs_per_year } }; if verbosity > 0 { println!("probs {:#?}", probs); } - dbg!(u); assert!((0f32..1.).contains(&u), "Invalid u: u>=0 and u<1"); assert!( (0f32..=1.).contains(&asymmetric),