diff --git a/src/eems2.cpp b/src/eems2.cpp index e059fea..78c1568 100644 --- a/src/eems2.cpp +++ b/src/eems2.cpp @@ -138,7 +138,7 @@ void EEMS2::initialize_state(const MCMC &mcmc) { nowmrateS = params.max_omegam; - nowqrateS = params.min_omegaq; + nowqrateS = params.max_omegaq; int niters = mcmc.num_iters_to_save(); mRates = MatrixXd::Zero(niters, d); diff --git a/src/util.cpp b/src/util.cpp index 41b47f6..4f97866 100644 --- a/src/util.cpp +++ b/src/util.cpp @@ -60,7 +60,7 @@ Params::Params(const string ¶ms_file, const long seed_from_command_line) { min_omegaq = -4; min_omegam = -4; - max_omegaq = log10(0); + max_omegaq = log10(1); max_omegam = log10(1.5); } ostream& operator<<(ostream& out, const Params& params) {