Skip to content
Snippets Groups Projects
Verified Commit 60ab11ea authored by Sébastien Villemot's avatar Sébastien Villemot
Browse files

Partial information: rename modified copies of Sims’ qzdiv and qzswitch

This prevents a name clash with qzdiv and qzswitch shipped with MS-SBVAR code.

(cherry picked from commit 2aad980f)
parent 8953f1ab
Branches
Tags
No related merge requests found
Pipeline #10829 failed
...@@ -89,10 +89,10 @@ License: public-domain-aim ...@@ -89,10 +89,10 @@ License: public-domain-aim
Files: matlab/optimization/bfgsi1.m matlab/optimization/csolve.m matlab/optimization/csminit1.m matlab/optimization/numgrad2.m Files: matlab/optimization/bfgsi1.m matlab/optimization/csolve.m matlab/optimization/csminit1.m matlab/optimization/numgrad2.m
matlab/optimization/numgrad3.m matlab/optimization/numgrad3_.m matlab/optimization/numgrad5.m matlab/optimization/numgrad3.m matlab/optimization/numgrad3_.m matlab/optimization/numgrad5.m
matlab/optimization/numgrad5_.m matlab/optimization/csminwel1.m matlab/+bvar/density.m matlab/optimization/numgrad5_.m matlab/optimization/csminwel1.m matlab/+bvar/density.m
matlab/+bvar/toolbox.m matlab/partial_information/PI_gensys.m matlab/partial_information/qzswitch.m matlab/+bvar/toolbox.m matlab/partial_information/PI_gensys.m matlab/partial_information/PI_qzswitch.m
matlab/partial_information/qzdiv.m matlab/partial_information/PI_qzdiv.m
Copyright: 1993-2009 Christopher Sims Copyright: 1993-2009 Christopher Sims
2006-2023 Dynare Team 2006-2024 Dynare Team
License: GPL-3+ License: GPL-3+
Files: matlab/optimization/cmaes.m Files: matlab/optimization/cmaes.m
......
...@@ -239,7 +239,7 @@ for i=1:nn ...@@ -239,7 +239,7 @@ for i=1:nn
end end
div ; div ;
if ~zxz if ~zxz
[a, b, ~, z]=qzdiv(div,a,b,q,z); [a, b, ~, z]=PI_qzdiv(div,a,b,q,z);
end end
gev=[diag(a) diag(b)]; gev=[diag(a) diag(b)];
......
function [A,B,Q,Z] = qzdiv(stake,A,B,Q,Z) function [A,B,Q,Z] = PI_qzdiv(stake,A,B,Q,Z)
%function [A,B,Q,Z] = qzdiv(stake,A,B,Q,Z) %function [A,B,Q,Z] = PI_qzdiv(stake,A,B,Q,Z)
% %
% Takes U.T. matrices A, B, orthonormal matrices Q,Z, rearranges them % Takes U.T. matrices A, B, orthonormal matrices Q,Z, rearranges them
% so that all cases of abs(B(i,i)/A(i,i))>stake are in lower right % so that all cases of abs(B(i,i)/A(i,i))>stake are in lower right
...@@ -10,7 +10,7 @@ function [A,B,Q,Z] = qzdiv(stake,A,B,Q,Z) ...@@ -10,7 +10,7 @@ function [A,B,Q,Z] = qzdiv(stake,A,B,Q,Z)
% http://sims.princeton.edu/yftp/gensys/mfiles/qzdiv.m % http://sims.princeton.edu/yftp/gensys/mfiles/qzdiv.m
% Copyright © 1993-2007 Christopher Sims % Copyright © 1993-2007 Christopher Sims
% Copyright © 2008-2017 Dynare Team % Copyright © 2008-2024 Dynare Team
% %
% This file is part of Dynare. % This file is part of Dynare.
% %
...@@ -43,7 +43,7 @@ for i = n:-1:1 ...@@ -43,7 +43,7 @@ for i = n:-1:1
return return
end end
for k=m:1:i-1 for k=m:1:i-1
[A, B, Q, Z] = qzswitch(k,A,B,Q,Z); [A, B, Q, Z] = PI_qzswitch(k,A,B,Q,Z);
tmp = root(k,2); tmp = root(k,2);
root(k,2) = root(k+1,2); root(k,2) = root(k+1,2);
root(k+1,2) = tmp; root(k+1,2) = tmp;
......
function [A,B,Q,Z] = qzswitch(i,A,B,Q,Z) function [A,B,Q,Z] = PI_qzswitch(i,A,B,Q,Z)
%function [A,B,Q,Z] = qzswitch(i,A,B,Q,Z) %function [A,B,Q,Z] = PI_qzswitch(i,A,B,Q,Z)
% %
% Takes U.T. matrices A, B, orthonormal matrices Q,Z, interchanges % Takes U.T. matrices A, B, orthonormal matrices Q,Z, interchanges
% diagonal elements i and i+1 of both A and B, while maintaining % diagonal elements i and i+1 of both A and B, while maintaining
...@@ -14,7 +14,7 @@ function [A,B,Q,Z] = qzswitch(i,A,B,Q,Z) ...@@ -14,7 +14,7 @@ function [A,B,Q,Z] = qzswitch(i,A,B,Q,Z)
% http://sims.princeton.edu/yftp/gensys/mfiles/qzswitch.m % http://sims.princeton.edu/yftp/gensys/mfiles/qzswitch.m
% Copyright © 1993-2007 Christopher Sims % Copyright © 1993-2007 Christopher Sims
% Copyright © 2008-2011 Dynare Team % Copyright © 2008-2024 Dynare Team
% %
% This file is part of Dynare. % This file is part of Dynare.
% %
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment