From ed88a8e528d9993a725ac8f73eb9ae0eea6cdb9f Mon Sep 17 00:00:00 2001
From: Houtan Bastani <houtanb@gmail.com>
Date: Wed, 9 Jun 2010 15:09:55 +0200
Subject: [PATCH] Bug fix: catch the case when a parameter follows a beta
 distribution with mean = std err = 0.5

---
 ComputingTasks.cc | 9 +++++++++
 1 file changed, 9 insertions(+)

diff --git a/ComputingTasks.cc b/ComputingTasks.cc
index ef3bec57..836d6db8 100644
--- a/ComputingTasks.cc
+++ b/ComputingTasks.cc
@@ -317,6 +317,15 @@ EstimatedParamsStatement::EstimatedParamsStatement(const vector<EstimationParams
   estim_params_list(estim_params_list_arg),
   symbol_table(symbol_table_arg)
 {
+  for (vector<EstimationParams>::const_iterator it = estim_params_list.begin();
+       it != estim_params_list.end(); it++)
+    if (it->prior == "1") //BETA_PDF is associated with "1" in DynareBison.yy
+      if (dynamic_cast<NumConstNode *>(it->mean)->isNumConstNodeEqualTo(0.5) &&
+          dynamic_cast<NumConstNode *>(it->std)->isNumConstNodeEqualTo(0.5))
+        {
+          cerr << "The prior density is not defined for the beta distribution when the mean = standard deviation = 0.5." << endl;
+          exit(EXIT_FAILURE);
+        }
 }
 
 void