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

v4 dr1: corrected typo in 'computer == 'PCWIN' 

git-svn-id: https://www.dynare.org/svn/dynare/dynare_v4@1479 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
michel 2007-12-07 13:20:33 +00:00
parent b5669458e6
commit 1e226a9783
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
matlab/dr1.m
View File

@ -454,7 +454,7 @@ zu=[zu; eye(M_.exo_nbr);zeros(M_.exo_det_nbr,M_.exo_nbr)];
[nrzx,nczx] = size(zx);
if ~exist('sparse_hessian_times_B_kronecker_C') | ((strcmp(version('-release'),'2007b') ...
~= 1) & compute == 'PCWIN')
~= 1) & computer == 'PCWIN')
if nrzx*nrzx*nczx*nczx > 1e7
rhs = zeros(M_.endo_nbr,nczx*nczx);
k1 = 1;
@ -534,7 +534,7 @@ hu = dr.ghu(nstatic+1:nstatic+npred,:);
%kk = kk(1:npred,1:npred);
%rhs = -hessian*kron(zx,zu)-f1*dr.ghxx(end-nyf+1:end,kk(:))*kron(hx(1:npred,:),hu(1:npred,:));
if ~exist('sparse_hessian_times_B_kronecker_C') | ((strcmp(version('-release'),'2007b') ...
~= 1) & compute == 'PCWIN')
~= 1) & computer == 'PCWIN')
if nrzx*nrzx*nczx*M_.exo_nbr > 1e7
rhs = zeros(M_.endo_nbr,nczx*M_.exo_nbr);
k1 = 1;
@ -556,7 +556,7 @@ hu1 = [hu;zeros(np-npred,M_.exo_nbr)];
[nrhx,nchx] = size(hx);
[nrhu1,nchu1] = size(hu1);
if ~exist('A_times_B_kronecker_C') | ((strcmp(version('-release'),'2007b') ...
~= 1) & compute == 'PCWIN')
~= 1) & computer == 'PCWIN')
if nrhx*nrhu1*nchx*nchu1 > 1e7
B1 = zeros(size(dr.ghxx,1),nchx*nchu1);
k1 = 1;
@ -584,7 +584,7 @@ dr.ghxu = A\rhs;
kk = reshape([1:np*np],np,np);
kk = kk(1:npred,1:npred);
if ~exist('sparse_hessian_times_B_kronecker_C') | ((strcmp(version('-release'),'2007b') ...
~= 1) & compute == 'PCWIN')
~= 1) & computer == 'PCWIN')
if nrzx*nrzx*M_.exo_nbr*M_.exo_nbr > 1e7
rhs = zeros(M_.endo_nbr,M_.exo_nbr*M_.exo_nbr);
k1 = 1;
@ -601,7 +601,7 @@ else
rhs = sparse_hessian_times_B_kronecker_C(hessian,zu);
end
if ~exist('A_times_B_kronecker_C') | ((strcmp(version('-release'),'2007b') ...
~= 1) & compute == 'PCWIN')
~= 1) & computer == 'PCWIN')
if nrhu1*nrhu1*nchu1*nchu1 > 1e7
B1 = zeros(size(dr.ghxx,1),nchu1*nchu1);
k1 = 1;
@ -662,7 +662,7 @@ for i=1:M_.maximum_endo_lead
find(M_.lead_lag_incidenceordered(M_.maximum_endo_lag+j+1,:));
nk3a = length(k3a);
if ~exist('sparse_hessian_times_B_kronecker_C') | ((strcmp(version('-release'),'2007b') ...
~= 1) & compute == 'PCWIN')
~= 1) & computer == 'PCWIN')
if nk3a*nk3a*M_.exo_nbr*M_.exo_nbr > 1e7
B1 = zeros(M_.endo_nbr,M_.exo_nbr*M_.exo_nbr);
k1 = 1;
@ -695,7 +695,7 @@ for i=1:M_.maximum_endo_lead
gu = dr.ghx*Gu;
[nrGu,ncGu] = size(Gu);
if ~exist('A_times_B_kronecker_C') | ((strcmp(version('-release'),'2007b') ...
~= 1) & compute == 'PCWIN')
~= 1) & computer == 'PCWIN')
if nrGu*nrGu*ncGu*ncGu > 1e7
G1 = zeros(M_.endo_nbr,ncGu*ncGu);
G2 = zeros(size(hxx,1),ncGu*ncGu);