From cf898e0d7ce2b36ad87a3992d58c500a1e8dd9ce Mon Sep 17 00:00:00 2001
From: halasadi <halasadi89@gmail.com>
Date: Sun, 11 Feb 2018 10:38:18 -0600
Subject: [PATCH] oops fixed really stupid bug I just introduced

---
 src/eems2.cpp | 2 +-
 src/util.cpp  | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

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 &params_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) {