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

passes mexext to k_order and reads and unpacks data passed back 

git-svn-id: https://www.dynare.org/svn/dynare/trunk@2446 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
george 2009-03-03 11:14:36 +00:00
parent 3945918e47
commit 68535ae212
1 changed files with 8 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
matlab/kordpert/dr1_k_order.m
View File

@ -133,13 +133,14 @@ function [dr,info,M_,options_,oo_] = dr1(dr,task,M_,options_,oo_)
% z = z(iyr0) ; % z = z(iyr0) ;
% oo_.dyn_ys=z; % extended ys % oo_.dyn_ys=z; % extended ys
try try
[ysteady, gx, gu]=k_order_perturbation(dr,task,M_,options_, oo_ ); [ysteady, ghx_u]=k_order_perturbation(dr,task,M_,options_, oo_ , ['.' mexext]);
load(M_.fname); % load(M_.fname);
ghxu = eval([M_.fname '_g_1']); % dr.ys=eval([M_.fname '_ss']);
sss= size(ghxu,2); dr.ys=ysteady;
dr.ghx= ghxu(:,1:sss-M_.exo_nbr); % ghxu = eval([M_.fname '_g_1']);
dr.ghu= ghxu(:,sss-M_.exo_nbr+1:end); sss= size(ghx_u,2);
dr.ys=eval([M_.fname '_ss']); dr.ghx= ghx_u(:,1:sss-M_.exo_nbr);
dr.ghu= ghx_u(:,sss-M_.exo_nbr+1:end);
catch catch
disp('*************************************************************************************'); disp('*************************************************************************************');
% disp('Problem with using k_order perturbation solver - Using Dynare solver instead'); % disp('Problem with using k_order perturbation solver - Using Dynare solver instead');