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

Update: test main which writes MAT file too 

git-svn-id: https://www.dynare.org/svn/dynare/trunk@2379 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
george 2009-01-26 19:48:53 +00:00
parent ad1608c9af
commit 88afe87f0c
1 changed files with 31 additions and 26 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

57
mex/sources/korderpert/src/k_order_test_main.cpp
View File

@ -197,39 +197,44 @@ int main(int argc, char* argv[])
, ySteady, vCov, modParams, nStat, nPred, nForw, nBoth
, jcols, nSteps, kOrder, journal, dynamicDLL, sstol, var_order
, llincidence );
try {
// intiate tensor library
// intiate tensor library
#ifdef DEBUG
mexPrintf("k_order_perturbation: Call tls init\n");
mexPrintf("k_order_perturbation: Call tls init\n");
#endif
tls.init(dynare.order(),
dynare.nstat()+2*dynare.npred()+3*dynare.nboth()+
2*dynare.nforw()+dynare.nexog());
// construct main K-order approximation class
// FistOrderApproximation app(dynare, journal, nSteps);
tls.init(dynare.order(),
dynare.nstat()+2*dynare.npred()+3*dynare.nboth()+
2*dynare.nforw()+dynare.nexog());
// construct main K-order approximation class
// FistOrderApproximation app(dynare, journal, nSteps);
#ifdef DEBUG
mexPrintf("k_order_perturbation: Call Approximation constructor \n");
mexPrintf("k_order_perturbation: Call Approximation constructor \n");
#endif
Approximation app(dynare, journal, nSteps);
// run stochastic steady
Approximation app(dynare, journal, nSteps);
// run stochastic steady
#ifdef DEBUG
mexPrintf("k_order_perturbation: Calling walkStochSteady.\n");
mexPrintf("k_order_perturbation: Calling walkStochSteady.\n");
#endif
app.walkStochSteady();
} catch (const KordException& e) {
// tell about the exception and continue
printf("Caught (not yet fatal) Kord exception: ");
e.print();
JournalRecord rec(journal);
rec << "Solution routine not finished (" << e.get_message()
<< "), see what happens" << endrec;
} catch (const TLException& e) {
mexErrMsgTxt("Caugth TL exception.");
} catch (SylvException& e) {
mexErrMsgTxt("Caught Sylv exception.");
}
app.walkStochSteady();
// open mat file
std::string matfile(fName);//(params.basename);
matfile += ".mat";
FILE* matfd = NULL;
if (NULL == (matfd=fopen(matfile.c_str(), "wb"))) {
fprintf(stderr, "Couldn't open %s for writing.\n", matfile.c_str());
exit(1);
}
std::string ss_matrix_name(fName);//params.prefix);
ss_matrix_name += "_steady_states";
ConstTwoDMatrix(app.getSS()).writeMat4(matfd, ss_matrix_name.c_str());
// write the folded decision rule to the Mat-4 file
app.getFoldDecisionRule().writeMat4(matfd, fName);//params.prefix);
fclose(matfd);
// get latest ysteady
double * dYsteady = (dynare.getSteady().base());
ySteady = (Vector*)(&dynare.getSteady());