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

Added routines related to auxiliary variables.

time-shift
Stéphane Adjemian(Charybdis) 2018-05-29 15:44:41 +02:00
parent d5862f93bf
commit 474556e0ca
4 changed files with 173 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

48
matlab/get_aux_variable_id.m Normal file
View File

@ -0,0 +1,48 @@
function ida = get_aux_variable_id(var)
% Returns the index of an auxiliary variable in M_.aux_vars
%
% INPUTS
% - var [string, integer] Variable name or index in M_.endo_names.
%
% OUTPUTS
% - ida [integer] Corresponding index in M_.aux_vars.
% Copyright (C) 2018 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 <http://www.gnu.org/licenses/>.
global M_
if ischar(var)
id = strmatch(var, M_.endo_names, 'exact');
if isempty(var)
error('Variable %s is unkown.', var)
else
var = id;
end
else
if ~isint(var) || var>M_.endo_nbr || var<1
error('Input must be the name of an endogenous variable or an integer between 1 and %s!', num2str(M_.endo_nbr))
end
end
if var<=M_.orig_endo_nbr
% var is not an auxiliary variable.
ida = 0;
else
ida = find([M_.aux_vars.endo_index]==var);
end

56
matlab/get_difference_order.m Normal file
View File

@ -0,0 +1,56 @@
function n = get_difference_order(var)
% Returns true iff endogenous variable `var` is a variable in difference.
%
% INPUTS
% - var [string, integer] Variable name or index in M_.endo_names.
%
% OUTPUTS
% - boo [logical] true/false.
% Copyright (C) 2018 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 <http://www.gnu.org/licenses/>.
global M_
% Get index of the auxiliary
id = get_aux_variable_id(var);
% Set default difference order
n = 0;
% If var is not an auxiliary variable it cannot be a first difference.
if ~id, return, end
while id
if M_.aux_vars(id).type==8
% diff auxiliary variable.
n = n+1;
v = M_.aux_vars(id).orig_index;
if v<=M_.orig_endo_nbr
id = 0;
else
id = get_aux_variable_id(v);
end
elseif M_.aux_vars(id).type==9
% lagged diff auxiliary variable
v = M_.aux_vars(id).orig_index;
id = get_aux_variable_id(v);
else
break
end
end

32
matlab/get_equation_number_by_tag.m Normal file
View File

@ -0,0 +1,32 @@
function eqnumber = get_equation_number_by_tag(eqname)
% Translates an equation name into an equation number.
%
% INPUTS
% - eqname [string] Name of the equation.
%
% OUTPUTS
% - eqnumber [integer] Equation number.
% Copyright (C) 2018 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 <http://www.gnu.org/licenses/>.
global M_
eqnumber = strmatch(eqname, M_.equations_tags(strmatch('name', M_.equations_tags(:,2), 'exact'), 3), 'exact');
if isempty(eqnumber), eqnumber = 0; end

37
matlab/isdiff.m Normal file
View File

@ -0,0 +1,37 @@
function boo = isdiff(var)
% Returns true iff endogenous variable `var` is a v ariable in difference.
%
% INPUTS
% - var [string, integer] Variable name or index in M_.endo_names.
%
% OUTPUTS
% - boo [logical] true/false.
% Copyright (C) 2018 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 <http://www.gnu.org/licenses/>.
global M_
ida = get_aux_variable_id(var);
boo = false;
if ~ida, return, end
if ismember(M_.aux_vars(ida).type, [8, 9])
boo = true;
end