From 19b2619d0613e7d1628d8a9a64f301a6daf3312d Mon Sep 17 00:00:00 2001
From: Johannes Pfeifer <jpfeifer@gmx.de>
Date: Sat, 2 Dec 2023 11:08:22 +0100
Subject: [PATCH] Remove unused gsa/gsa_speed.m
---
license.txt | 1 -
matlab/gsa/gsa_speed.m | 58 ------------------------------------------
2 files changed, 59 deletions(-)
delete mode 100644 matlab/gsa/gsa_speed.m
diff --git a/license.txt b/license.txt
index c39165710b..71083f2ab4 100644
--- a/license.txt
+++ b/license.txt
@@ -173,7 +173,6 @@ Files: matlab/gsa/cumplot.m
matlab/gsa/filt_mc_.m
matlab/gsa/gsa_plotmatrix.m
matlab/gsa/gsa_skewness.m
- matlab/gsa/gsa_speed.m
matlab/gsa/log_trans_.m
matlab/gsa/map_calibration.m
matlab/gsa/map_ident_.m
diff --git a/matlab/gsa/gsa_speed.m b/matlab/gsa/gsa_speed.m
deleted file mode 100644
index 33b14695dd..0000000000
--- a/matlab/gsa/gsa_speed.m
+++ /dev/null
@@ -1,58 +0,0 @@
-function [tadj, iff] = gsa_speed(A,B,mf,p)
-% [tadj, iff] = gsa_speed(A,B,mf,p),
-%
-% Written by Marco Ratto
-% Joint Research Centre, The European Commission,
-% marco.ratto@ec.europa.eu
-
-% Copyright © 2012 European Commission
-% Copyright © 2012-2017 Dynare Team
-%
-% This file is part of Dynare.
-%
-% Dynare is free software: you can redistribute it and/or modify
-% it under the terms of the GNU General Public License as published by
-% the Free Software Foundation, either version 3 of the License, or
-% (at your option) any later version.
-%
-% Dynare is distributed in the hope that it will be useful,
-% but WITHOUT ANY WARRANTY; without even the implied warranty of
-% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
-% GNU General Public License for more details.
-%
-% You should have received a copy of the GNU General Public License
-% along with Dynare. If not, see <https://www.gnu.org/licenses/>.
-
-nvar=length(mf);
-nstate= size(A,1);
-nshock = size(B,2);
-nrun=size(B,3);
-
-iff=zeros(nvar,nshock,nrun);
-tadj=iff;
-disp('Computing speed of adjustement ...')
-h = dyn_waitbar(0,'Speed of adjustement...');
-
-for i=1:nrun
- irf=zeros(nvar,nshock);
- a=squeeze(A(:,:,i));
- b=squeeze(B(:,:,i));
- IFF=inv(eye(nstate)-a)*b;
- iff(:,:,i)=IFF(mf,:);
- IF=IFF-b;
-
- t=0;
- while any(any(irf<0.5))
- t=t+1;
- IFT=((eye(nstate)-a^(t+1))*inv(eye(nstate)-a))*b-b;
- irf=IFT(mf,:)./(IF(mf,:)+eps);
- irf = irf.*(abs(IF(mf,:))>1.e-7)+(abs(IF(mf,:))<=1.e-7);
- %irf=ft(mf,:);
- tt=(irf>0.5).*t;
- tadj(:,:,i)=((tt-tadj(:,:,i))==tt).*tt+tadj(:,:,i);
- end
- dyn_waitbar(i/nrun,h)
-end
-skipline()
-disp('.. done !')
-dyn_waitbar_close(h)
--
GitLab