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

Correct bugs in dr_block: 

- the index of lead_lag_incidence for forward looking blocks
   - the block type index
time-shift
Ferhat Mihoubi 2011-09-11 14:37:06 +02:00
parent 0c72764932
commit e51d8f7a56
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
matlab/dr_block.m
View File

@ -216,8 +216,8 @@ for i = 1:Size;
%% ------------------------------------------------------------------
%Evaluate Backward
if maximum_lead > 0 && n_fwrd > 0
indx_r = find(M_.block_structure.block(i).lead_lag_incidence(2,:));
indx_c = M_.block_structure.block(i).lead_lag_incidence(2,indx_r);
indx_r = find(M_.block_structure.block(i).lead_lag_incidence(3,:));
indx_c = M_.block_structure.block(i).lead_lag_incidence(3,indx_r);
data(i).eigval = 1 ./ diag(jacob(indx_r, indx_c));
data(i).rank = sum(abs(data(i).eigval) > 0);
else
@ -315,7 +315,7 @@ for i = 1:Size;
end;
dr.rank = dr.rank + data(i).rank;
dr.eigval = [dr.eigval ; data(i).eigval];
case 5
case 6
%% ------------------------------------------------------------------
%Solve Forward complete
if maximum_lag > 0 && n_pred > 0
@ -382,7 +382,7 @@ for i = 1:Size;
end
end
end
case 6
case 7
%% ------------------------------------------------------------------
%Solve Backward complete
if maximum_lead > 0 && n_fwrd > 0
@ -395,7 +395,7 @@ for i = 1:Size;
end;
dr.rank = dr.rank + data(i).rank;
dr.eigval = [dr.eigval ; data(i).eigval];
case 8
case {5,8}
%% ------------------------------------------------------------------
%The lead_lag_incidence contains columns in the following order:
% static variables, backward variable, mixed variables and forward variables