4.2: corrected bugs for terminal_condition == 2
git-svn-id: https://www.dynare.org/svn/dynare/trunk@3298 ac1d8469-bf42-47a9-8791-bf33cf982152time-shift
parent
2a75491196
commit
9b3718b50f
|
@ -78,7 +78,8 @@ elseif isstruct(compute_linear_solution)
|
||||||
end
|
end
|
||||||
|
|
||||||
if compute_linear_solution
|
if compute_linear_solution
|
||||||
ghx = dr.ghx(end-dr.nfwrd+1:end,:);
|
ghx(dr.order_var,:) = dr.ghx;
|
||||||
|
ghx = ghx(iyf,:);
|
||||||
end
|
end
|
||||||
|
|
||||||
periods = options_.periods;
|
periods = options_.periods;
|
||||||
|
@ -134,7 +135,7 @@ for iter = 1:options_.maxit_
|
||||||
ic = ic + ny;
|
ic = ic + ny;
|
||||||
icp = icp + ny;
|
icp = icp + ny;
|
||||||
s = jacobian(:,is);
|
s = jacobian(:,is);
|
||||||
s(:,iyp-nyp) = s(:,iyp-nyp)+jacobian(:,isf)*ghx;
|
s(:,iyp) = s(:,iyp)+jacobian(:,isf)*ghx;
|
||||||
c (ic,:) = s\jacobian(:,isf1);
|
c (ic,:) = s\jacobian(:,isf1);
|
||||||
end
|
end
|
||||||
c = bksup0(c,ny,nrc,iyf,icf,periods);
|
c = bksup0(c,ny,nrc,iyf,icf,periods);
|
||||||
|
|
Loading…
Reference in New Issue