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

v4: added test exist('gensylv') 

git-svn-id: https://www.dynare.org/svn/dynare/dynare_v4@479 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
michel 2005-10-03 09:45:26 +00:00
parent a6ec6bd2ec
commit 9550d6f798
1 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
matlab/dr1.m
View File

@ -440,12 +440,16 @@ end
[junk,k1,k2] = find(M_.lead_lag_incidence(M_.maximum_lag+M_.maximum_lead+1,order_var));
A(1:M_.endo_nbr,nstatic+1:nstatic+npred)=...
A(1:M_.endo_nbr,nstatic+[1:npred])+jacobia_(:,k2)*gx1(k1,1:npred);
%C = kron(hx,hx);
C = hx;
D = [rhs; zeros(n-M_.endo_nbr,size(rhs,2))];
%x0 = sylvester3(A,B,C,D);
%dr.ghxx = sylvester3a(x0,A,B,C,D);
dr.ghxx = gensylv(2,A,B,C,D);
if exist('gensylv')
dr.ghxx = gensylv(2,A,B,C,D);
else
C = kron(hx,hx);
x0 = sylvester3(A,B,C,D);
dr.ghxx = sylvester3a(x0,A,B,C,D);
end
%ghxu
%rhs