From 75440028ebe26a31ec8e34dd657fc7990529ceec Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?S=C3=A9bastien=20Villemot?= <sebastien.villemot@ens.fr> Date: Thu, 16 Sep 2010 15:37:09 +0200 Subject: [PATCH] Preprocessor: fix crash when mean or stderr of an estimated parameter is specified as an expression instead of a plain value (bug introduced in 28d30d9a626a8d7be0866038924bfc52a935ee27) --- ComputingTasks.cc | 16 ++++++++++++---- 1 file changed, 12 insertions(+), 4 deletions(-) diff --git a/ComputingTasks.cc b/ComputingTasks.cc index 8030fe59..483b3d89 100644 --- a/ComputingTasks.cc +++ b/ComputingTasks.cc @@ -362,12 +362,20 @@ EstimatedParamsStatement::checkPass(ModFileStructure &mod_file_struct) if (it->name == "dsge_prior_weight") mod_file_struct.dsge_prior_weight_in_estimated_params = true; + // Handle case of degenerate beta prior 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)) + try { - cerr << "ERROR: The prior density is not defined for the beta distribution when the mean = standard deviation = 0.5." << endl; - exit(EXIT_FAILURE); + if (it->mean->eval(eval_context_type()) == 0.5 && + it->std->eval(eval_context_type()) == 0.5) + { + cerr << "ERROR: The prior density is not defined for the beta distribution when the mean = standard deviation = 0.5." << endl; + exit(EXIT_FAILURE); + } + } + catch (ExprNode::EvalException &e) + { + // We don't have enough information to compute the numerical value, skip the test } } } -- GitLab