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

Merge remote-tracking branch 'ratto/master'

time-shift
Sébastien Villemot 2012-09-18 11:13:50 +02:00
parent fd7386b593 fce6c5cc31
commit c41ce4ae69
2 changed files with 7 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
matlab/disp_identification.m
View File

@ -285,7 +285,7 @@ if SampleSize==1 && advanced,
pax(i,dumpindx)=idemoments.cosnJ(i,j);
end
end
fprintf('%-15s [%s] %10.3f\n',name{i},namx,idemoments.cosnJ(i,j))
fprintf('%-15s [%s] %14.7f\n',name{i},namx,idemoments.cosnJ(i,j))
end
end
end

21
matlab/getH.m
View File

@ -388,21 +388,12 @@ if kronflag==1, % kronecker products
H = -DfDtau\DfDthet;
x = reshape(H(1:m*m,:),m,m,param_nbr);
y = reshape(Dm*H(m*m+1:end,:),m,m,param_nbr);
if nauxe,
x = x(nauxe+1:end,nauxe+1:end,:);
y = y(nauxe+1:end,nauxe+1:end,:);
dA = x;
dOm = y;
m = size(y,1);
x = reshape(x,m*m,param_nbr);
Dm = duplication(m);
DmPl = inv(Dm'*Dm)*Dm';
y = DmPl*reshape(y,m*m,param_nbr);
H = [x;y];
else
dA = x;
dOm = y;
end
dA = x;
dOm = y;
% convert to dyn_vech
tmpH = Dm*H(m*m+1:end,:);
Index = find(triu(ones(m)));
H(m*m+1:end,:) = tmpH(Index,:);
Hx = [];
if ~isempty(indexo),