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

fixing bug in test case of standalone c++

time-shift
Michel Juillard 2014-10-14 10:25:07 +02:00
parent 8244862bd1
commit c8ac7007f0
1 changed files with 11 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
others/cpp/tests/test1.cc
View File

@ -32,20 +32,29 @@ main(int argc, char **argv)
int nback = zeta_back.size();
int nmixed = zeta_mixed.size();
int nstatic = zeta_static.size();
int sdyn = nfwrd + nback + 2*nmixed;
int jacob_cols = zeta_back.size()+zeta_fwrd.size()+zeta_mixed.size()+exo_nbr;
int jacob_cols = sdyn + endo_nbr + exo_nbr;
double *exo_steady_state = new double[exo_nbr];
double *jacob_data = new double[endo_nbr*jacob_cols];
std::cout << endo_nbr << " " << jacob_cols << endl;
FirstDerivatives(steady_state, exo_steady_state, 0, params, steady_state, 1, NULL, jacob_data, NULL, NULL);
std::cout << "g1[44] = " << jacob_data[44] << endl;
DecisionRules dr(endo_nbr, exo_nbr, zeta_fwrd, zeta_back, zeta_mixed,
zeta_static, qz_criterium);
int sdyn = nfwrd + nback + 2*nmixed;
std::cout << "g1[44] = " << jacob_data[44] << endl;
int jacobian_col_nbr = sdyn + endo_nbr + exo_nbr;
std::cout << "g1[44] = " << jacob_data[44] << endl;
MatrixView jacob_tmp(jacob_data, endo_nbr, jacobian_col_nbr, endo_nbr);
std::cout << "g1[44] = " << jacob_data[44] << endl;
std::cout << "jacob_tmp(2,7) = "<< jacob_tmp(2,7) << endl;
Matrix jacobian(endo_nbr, jacobian_col_nbr), g_y(endo_nbr, nback+nmixed), g_u(endo_nbr, exo_nbr);
jacobian = jacob_tmp;