From aec140697c69614c7f2344a623f348c2df5774ad Mon Sep 17 00:00:00 2001
From: edmundcrawley <edmundcrawley@gmail.com>
Date: Wed, 17 Jun 2020 12:19:59 -0400
Subject: [PATCH 1/2] Bug fix to drawDiscrete

These two lines can be run outside the loop, avoiding an error when N<events.size and also saving time.
The error occurs when j=0 and event_draws is empty (X[event_draws] throws an error).
---
 HARK/distribution.py | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/HARK/distribution.py b/HARK/distribution.py
index 314591630..5ef12d872 100644
--- a/HARK/distribution.py
+++ b/HARK/distribution.py
@@ -514,9 +514,9 @@ def drawDiscrete(self, N,X=None,exact_match=False,seed=0):
                 bot = top
                 top = cutoffs[j]
                 event_list += (top-bot)*[events[j]]
-                # Randomly permute the event indices and store the corresponding results
-                event_draws = RNG.permutation(event_list)
-                draws = X[event_draws]
+            # Randomly permute the event indices and store the corresponding results
+            event_draws = RNG.permutation(event_list)
+            draws = X[event_draws]
         else:
             indices = self.draw_events(N, seed=seed)
             if J > 1:

From 3f69148cfd3311b26da891f1cb12ba34722a0f7c Mon Sep 17 00:00:00 2001
From: edmundcrawley <edmundcrawley@gmail.com>
Date: Wed, 17 Jun 2020 13:54:55 -0400
Subject: [PATCH 2/2] Updates testAggShockMarkovConsumerType due to random
 number change

---
 HARK/ConsumptionSaving/tests/test_ConsAggShockModel.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/HARK/ConsumptionSaving/tests/test_ConsAggShockModel.py b/HARK/ConsumptionSaving/tests/test_ConsAggShockModel.py
index 4aa2c06f3..8cf28c183 100644
--- a/HARK/ConsumptionSaving/tests/test_ConsAggShockModel.py
+++ b/HARK/ConsumptionSaving/tests/test_ConsAggShockModel.py
@@ -79,5 +79,5 @@ def test_economy(self):
         
         self.economy.AFunc = self.economy.dynamics.AFunc
         self.assertAlmostEqual(self.economy.AFunc[0].slope,
-                         1.0921217053006234)
+                         1.0845554708377696)