stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

trunk preprocessor: fixed bug with estimated_params_init and estimated_params_bounds in conjunction with correlation of endogenous variables 

git-svn-id: https://www.dynare.org/svn/dynare/trunk@2366 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
sebastien 2009-01-22 10:10:39 +00:00
parent 3f36ec5b07
commit 50985b004e
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
preprocessor/ComputingTasks.cc
View File

@ -1,5 +1,5 @@
/* /*
* Copyright (C) 2003-2008 Dynare Team * Copyright (C) 2003-2009 Dynare Team
* *
* This file is part of Dynare. * This file is part of Dynare.
* *
@ -495,7 +495,7 @@ EstimatedParamsInitStatement::writeOutput(ostream &output, const string &basenam
else if (symbol_table.getType(it->name) == eEndogenous) else if (symbol_table.getType(it->name) == eEndogenous)
{ {
output << "tmp1 = find((estim_params_.corrn(:,1)==" << symbol_table.getID(it->name)+1 << ")) & (estim_params_.corrn(:,2)==" << symbol_table.getID(it->name2)+1 << ";\n"; output << "tmp1 = find((estim_params_.corrn(:,1)==" << symbol_table.getID(it->name)+1 << ")) & (estim_params_.corrn(:,2)==" << symbol_table.getID(it->name2)+1 << ";\n";
output << "estim_params_.corrx(tmp1,3) = " << it->init_val << ";\n"; output << "estim_params_.corrn(tmp1,3) = " << it->init_val << ";\n";
} }
} }
} }
@ -547,8 +547,8 @@ EstimatedParamsBoundsStatement::writeOutput(ostream &output, const string &basen
else if (symbol_table.getType(it->name) == eEndogenous) else if (symbol_table.getType(it->name) == eEndogenous)
{ {
output << "tmp1 = find((estim_params_.corrn(:,1)==" << symbol_table.getID(it->name)+1 << ")) & (estim_params_.corrn(:,2)==" << symbol_table.getID(it->name2)+1 << ";\n"; output << "tmp1 = find((estim_params_.corrn(:,1)==" << symbol_table.getID(it->name)+1 << ")) & (estim_params_.corrn(:,2)==" << symbol_table.getID(it->name2)+1 << ";\n";
output << "estim_params_.corrx(tmp1,4) = " << it->low_bound << ";\n"; output << "estim_params_.corrn(tmp1,4) = " << it->low_bound << ";\n";
output << "estim_params_.corrx(tmp1,5) = " << it->up_bound << ";\n"; output << "estim_params_.corrn(tmp1,5) = " << it->up_bound << ";\n";
} }
} }
} }