From 9ee2a23a5dbcce9c476b1ed05950096644cef9e6 Mon Sep 17 00:00:00 2001 From: george Date: Thu, 12 Feb 2009 17:21:38 +0000 Subject: [PATCH] Updated for passing qz_criterium and unit root models git-svn-id: https://www.dynare.org/svn/dynare/trunk@2415 ac1d8469-bf42-47a9-8791-bf33cf982152 --- mex/sources/korderpert/src/k_order_perturbation.cpp | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/mex/sources/korderpert/src/k_order_perturbation.cpp b/mex/sources/korderpert/src/k_order_perturbation.cpp index 1ba6823d5..e46d89e0c 100644 --- a/mex/sources/korderpert/src/k_order_perturbation.cpp +++ b/mex/sources/korderpert/src/k_order_perturbation.cpp @@ -114,14 +114,14 @@ extern "C" { double qz_criterium = 1+1e-6; mxFldp = mxGetField(options_, 0,"qz_criterium" ); if (mxIsNumeric(mxFldp)) - qz_criterium = (int)mxGetScalar(mxFldp); + qz_criterium = (double)mxGetScalar(mxFldp); mxFldp = mxGetField(M_, 0,"params" ); double * dparams = (double *) mxGetData(mxFldp); int npar = (int)mxGetM(mxFldp); Vector * modParams = new Vector(dparams, npar); #ifdef DEBUG - mexPrintf("k_ord_perturbation: nParams=%d .\n",npar); + mexPrintf("k_ord_perturbation: qz_criterium=%g, nParams=%d .\n",qz_criterium,npar); for (int i = 0; i < npar; i++) { mexPrintf("k_ord_perturbation: Params[%d]= %g.\n", i, (*modParams)[i]); } #endif @@ -271,7 +271,7 @@ extern "C" { // DynamicFn * pDynamicFn = loadModelDynamicDLL (fname); DynamicModelDLL dynamicDLL(fName,nEndo, jcols, nMax_lag, nExog); #ifdef DEBUG - mexPrintf("k_order_perturbation: Calling dynare constructor.\n"); + mexPrintf("k_order_perturbation: Calling dynare constructor .\n"); #endif // make KordpDynare object KordpDynare dynare(endoNamesMX, nEndo, exoNamesMX, nExog, nPar, // paramNames, @@ -299,9 +299,9 @@ extern "C" { // 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 with qz_criterium=%f \n", qz_criterium); #endif - Approximation app(dynare, journal, nSteps); + Approximation app(dynare, journal, nSteps, false, qz_criterium); // run stochastic steady #ifdef DEBUG mexPrintf("k_order_perturbation: Calling walkStochSteady.\n");