From b5aeefe526af0b024b2a52d09cb819d831822497 Mon Sep 17 00:00:00 2001
From: Johannes Pfeifer <jpfeifer@gmx.de>
Date: Fri, 24 Jul 2015 16:23:40 +0200
Subject: [PATCH] Make prior_draw.m check for uniform distribution violation of
 bounds

Follows the logic for all other distributions
---
 matlab/prior_draw.m | 8 +++++++-
 1 file changed, 7 insertions(+), 1 deletion(-)

diff --git a/matlab/prior_draw.m b/matlab/prior_draw.m
index 9a22dc39c..937730748 100644
--- a/matlab/prior_draw.m
+++ b/matlab/prior_draw.m
@@ -1,4 +1,4 @@
-function pdraw = prior_draw(init,uniform)
+function pdraw = prior_draw(init,uniform) % --*-- Unitary tests --*--
 % This function generate one draw from the joint prior distribution.
 % 
 % INPUTS 
@@ -93,6 +93,11 @@ end
 
 if uniform_draws
     pdraw(uniform_index) = rand(length(uniform_index),1).*(p4(uniform_index)-p3(uniform_index)) + p3(uniform_index);  
+    out_of_bound = find( (pdraw(uniform_index)'>ub(uniform_index)) | (pdraw(uniform_index)'<lb(uniform_index)));
+    while ~isempty(out_of_bound),
+        pdraw(uniform_index) = rand(length(uniform_index),1).*(p4(uniform_index)-p3(uniform_index)) + p3(uniform_index);
+        out_of_bound = find( (pdraw(uniform_index)'>ub(uniform_index)) | (pdraw(uniform_index)'<lb(uniform_index)));
+    end
 end
 
 if gaussian_draws
@@ -143,3 +148,4 @@ if inverse_gamma_2_draws
         out_of_bound = find( (pdraw(inverse_gamma_2_index)'>ub(inverse_gamma_2_index)) | (pdraw(inverse_gamma_2_index)'<lb(inverse_gamma_2_index)));
     end
 end
+