diff --git a/license.txt b/license.txt
index 51b37da97b0b9067955540c6c0ec037bc28f0454..120c04eff6f7eb890a82abf0e2732e953a02de66 100644
--- a/license.txt
+++ b/license.txt
@@ -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
  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/+bvar/toolbox.m matlab/partial_information/PI_gensys.m matlab/partial_information/qzswitch.m
- matlab/partial_information/qzdiv.m
+ matlab/+bvar/toolbox.m matlab/partial_information/PI_gensys.m matlab/partial_information/PI_qzswitch.m
+ matlab/partial_information/PI_qzdiv.m
 Copyright: 1993-2009 Christopher Sims
-           2006-2023 Dynare Team
+           2006-2024 Dynare Team
 License: GPL-3+
 
 Files: matlab/optimization/cmaes.m
diff --git a/matlab/partial_information/PI_gensys.m b/matlab/partial_information/PI_gensys.m
index 01639bf7d7076d6eeaa1da1e986b335d587802ce..361ba717e8e999f59672391109433c084e321c5b 100644
--- a/matlab/partial_information/PI_gensys.m
+++ b/matlab/partial_information/PI_gensys.m
@@ -195,7 +195,7 @@ for i=1:nn
 end
 
 if ~zxz
-    [a, b, ~, z]=qzdiv(div,a,b,q,z);
+    [a, b, ~, z]=PI_qzdiv(div,a,b,q,z);
 end
 
 gev=[diag(a) diag(b)];
diff --git a/matlab/partial_information/qzdiv.m b/matlab/partial_information/PI_qzdiv.m
similarity index 89%
rename from matlab/partial_information/qzdiv.m
rename to matlab/partial_information/PI_qzdiv.m
index 4c802eb8e50b87942a3aa05e5a31ea35afb9782c..710df76d23f29886e6c16ce30492529043a92005 100644
--- a/matlab/partial_information/qzdiv.m
+++ b/matlab/partial_information/PI_qzdiv.m
@@ -1,5 +1,5 @@
-function [A,B,Q,Z] = 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)
+%function [A,B,Q,Z] = PI_qzdiv(stake,A,B,Q,Z)
 %
 % 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
@@ -10,7 +10,7 @@ function [A,B,Q,Z] = qzdiv(stake,A,B,Q,Z)
 % http://sims.princeton.edu/yftp/gensys/mfiles/qzdiv.m
 
 % Copyright © 1993-2007 Christopher Sims
-% Copyright © 2008-2017 Dynare Team
+% Copyright © 2008-2024 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -43,7 +43,7 @@ for i = n:-1:1
         return
     end
     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);
         root(k,2) = root(k+1,2);
         root(k+1,2) = tmp;
diff --git a/matlab/partial_information/qzswitch.m b/matlab/partial_information/PI_qzswitch.m
similarity index 94%
rename from matlab/partial_information/qzswitch.m
rename to matlab/partial_information/PI_qzswitch.m
index 60869de79de2775bcd435e37545db0ea4d2df800..74132abb9441096a02b39238149786c199169618 100644
--- a/matlab/partial_information/qzswitch.m
+++ b/matlab/partial_information/PI_qzswitch.m
@@ -1,5 +1,5 @@
-function [A,B,Q,Z] = 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)
+%function [A,B,Q,Z] = PI_qzswitch(i,A,B,Q,Z)
 %
 % Takes U.T. matrices A, B, orthonormal matrices Q,Z, interchanges
 % 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)
 % http://sims.princeton.edu/yftp/gensys/mfiles/qzswitch.m
 
 % Copyright © 1993-2007 Christopher Sims
-% Copyright © 2008-2011 Dynare Team
+% Copyright © 2008-2024 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -79,4 +79,4 @@ B(i:i+1,:) = xy*B(i:i+1,:);
 A(:,i:i+1) = A(:,i:i+1)*wz;
 B(:,i:i+1) = B(:,i:i+1)*wz;
 Z(:,i:i+1) = Z(:,i:i+1)*wz;
-Q(i:i+1,:) = xy*Q(i:i+1,:);
\ No newline at end of file
+Q(i:i+1,:) = xy*Q(i:i+1,:);