From ebe81eb6c90837914205ff003f5b2fde7a35b86b Mon Sep 17 00:00:00 2001 From: Johannes Pfeifer Date: Mon, 11 Oct 2021 17:45:03 +0200 Subject: [PATCH] display_problematic_vars_Jacobian.m: fix display when auxiliary equations outside of Ramsey are present --- matlab/display_problematic_vars_Jacobian.m | 36 +++++++++++----------- 1 file changed, 18 insertions(+), 18 deletions(-) diff --git a/matlab/display_problematic_vars_Jacobian.m b/matlab/display_problematic_vars_Jacobian.m index 53a6edc35..388f5d020 100644 --- a/matlab/display_problematic_vars_Jacobian.m +++ b/matlab/display_problematic_vars_Jacobian.m @@ -1,5 +1,5 @@ function []=display_problematic_vars_Jacobian(problemrow,problemcol,M_,x,type,caller_string) -% []=display_problematic_vars_Jacobian(problemrow,problemcol,M_,ys,caller_string) +% []=display_problematic_vars_Jacobian(problemrow,problemcol,M_,x,caller_string) % print the equation numbers and variables associated with problematic entries % of the Jacobian % @@ -16,7 +16,7 @@ function []=display_problematic_vars_Jacobian(problemrow,problemcol,M_,x,type,ca % none. % -% Copyright (C) 2014-2018 Dynare Team +% Copyright (C) 2014-2021 Dynare Team % % This file is part of Dynare. % @@ -37,7 +37,7 @@ skipline(); if nargin<6 caller_string=''; end -aux_eq_nbr=M_.eq_nbr-M_.orig_eq_nbr; +initial_aux_eq_nbr=M_.ramsey_eq_nbr; if strcmp(type,'dynamic') for ii=1:length(problemrow) if problemcol(ii)>max(M_.lead_lag_incidence) @@ -54,46 +54,46 @@ if strcmp(type,'dynamic') type_string='lead of'; end if problemcol(ii)<=max(max(M_.lead_lag_incidence)) && var_index<=M_.orig_endo_nbr - if problemrow(ii)<=aux_eq_nbr + if problemrow(ii)<=initial_aux_eq_nbr eq_nbr = problemrow(ii); fprintf('Derivative of Auxiliary Equation %d with respect to %s Variable %s (initial value of %s: %g) \n', ... eq_nbr, type_string, M_.endo_names{var_index}, M_.endo_names{var_index}, x(var_index)); else - eq_nbr = problemrow(ii)-aux_eq_nbr; + eq_nbr = problemrow(ii)-initial_aux_eq_nbr; fprintf('Derivative of Equation %d with respect to %s Variable %s (initial value of %s: %g) \n', ... eq_nbr, type_string, M_.endo_names{var_index}, M_.endo_names{var_index}, x(var_index)); end elseif problemcol(ii)<=max(max(M_.lead_lag_incidence)) && var_index>M_.orig_endo_nbr % auxiliary vars if M_.aux_vars(1,problemcol(ii)-M_.orig_endo_nbr).type==6 %Ramsey Lagrange Multiplier - if problemrow(ii)<=aux_eq_nbr + if problemrow(ii)<=initial_aux_eq_nbr eq_nbr = problemrow(ii); fprintf('Derivative of Auxiliary Equation %d with respect to %s of Langrange multiplier of equation %s (initial value: %g) \n', ... eq_nbr, type_string, M_.aux_vars(1,problemcol(ii)-M_.orig_endo_nbr).eq_nbr, x(problemcol(ii))); else - eq_nbr = problemrow(ii)-aux_eq_nbr; + eq_nbr = problemrow(ii)-initial_aux_eq_nbr; fprintf('Derivative of Equation %d with respect to %s of Langrange multiplier of equation %s (initial value: %g) \n', ... eq_nbr, type_string, M_.aux_vars(1,problemcol(ii)-M_.orig_endo_nbr).eq_nbr, x(problemcol(ii))); end else - if problemrow(ii)<=aux_eq_nbr + if problemrow(ii)<=initial_aux_eq_nbr eq_nbr = problemrow(ii); orig_var_index = M_.aux_vars(1,var_index-M_.orig_endo_nbr).orig_index; fprintf('Derivative of Auxiliary Equation %d with respect to %s Variable %s (initial value of %s: %g) \n', ... eq_nbr, type_string, M_.endo_names{orig_var_index}, M_.endo_names{orig_var_index}, x(orig_var_index)); else - eq_nbr = problemrow(ii)-aux_eq_nbr; + eq_nbr = problemrow(ii)-initial_aux_eq_nbr; orig_var_index = M_.aux_vars(1,var_index-M_.orig_endo_nbr).orig_index; fprintf('Derivative of Equation %d with respect to %s Variable %s (initial value of %s: %g) \n', ... eq_nbr, type_string, M_.endo_names{orig_var_index}, M_.endo_names{orig_var_index}, x(orig_var_index)); end end elseif problemcol(ii)>max(max(M_.lead_lag_incidence)) && var_index<=M_.exo_nbr - if problemrow(ii)<=aux_eq_nbr + if problemrow(ii)<=initial_aux_eq_nbr eq_nbr = problemrow(ii); fprintf('Derivative of Auxiliary Equation %d with respect to %s shock %s \n', ... eq_nbr, type_string, M_.exo_names{var_index}); else - eq_nbr = problemrow(ii)-aux_eq_nbr; + eq_nbr = problemrow(ii)-initial_aux_eq_nbr; fprintf('Derivative of Equation %d with respect to %s shock %s \n', ... eq_nbr, type_string, M_.exo_names{var_index}); end @@ -102,40 +102,40 @@ if strcmp(type,'dynamic') end end fprintf('\n%s The problem most often occurs, because a variable with\n', caller_string) - fprintf('%s exponent smaller than 1 has been initialized to 0. Taking the derivative\n', caller_string) + fprintf('%s exponent smaller than 0 has been initialized to 0. Taking the derivative\n', caller_string) fprintf('%s and evaluating it at the steady state then results in a division by 0.\n', caller_string) fprintf('%s If you are using model-local variables (# operator), check their values as well.\n', caller_string) elseif strcmp(type, 'static') for ii=1:length(problemrow) if problemcol(ii)<=M_.orig_endo_nbr - if problemrow(ii)<=aux_eq_nbr + if problemrow(ii)<=initial_aux_eq_nbr eq_nbr = problemrow(ii); fprintf('Derivative of Auxiliary Equation %d with respect to Variable %s (initial value of %s: %g) \n', ... eq_nbr, M_.endo_names{problemcol(ii)}, M_.endo_names{problemcol(ii)}, x(problemcol(ii))); else - eq_nbr = problemrow(ii)-aux_eq_nbr; + eq_nbr = problemrow(ii)-initial_aux_eq_nbr; fprintf('Derivative of Equation %d with respect to Variable %s (initial value of %s: %g) \n', ... eq_nbr, M_.endo_names{problemcol(ii)}, M_.endo_names{problemcol(ii)}, x(problemcol(ii))); end else %auxiliary vars if M_.aux_vars(1,problemcol(ii)-M_.orig_endo_nbr).type ==6 %Ramsey Lagrange Multiplier - if problemrow(ii)<=aux_eq_nbr + if problemrow(ii)<=initial_aux_eq_nbr eq_nbr = problemrow(ii); fprintf('Derivative of Auxiliary Equation %d with respect to Lagrange multiplier of equation %d (initial value: %g) \n', ... eq_nbr, M_.aux_vars(1,problemcol(ii)-M_.orig_endo_nbr).eq_nbr, x(problemcol(ii))); else - eq_nbr = problemrow(ii)-aux_eq_nbr; + eq_nbr = problemrow(ii)-initial_aux_eq_nbr; fprintf('Derivative of Equation %d with respect to Lagrange multiplier of equation %d (initial value: %g) \n', ... eq_nbr, M_.aux_vars(1,problemcol(ii)-M_.orig_endo_nbr).eq_nbr, x(problemcol(ii))); end else - if problemrow(ii)<=aux_eq_nbr + if problemrow(ii)<=initial_aux_eq_nbr eq_nbr = problemrow(ii); orig_var_index = M_.aux_vars(1,problemcol(ii)-M_.orig_endo_nbr).orig_index; fprintf('Derivative of Auxiliary Equation %d with respect to Variable %s (initial value of %s: %g) \n', ... eq_nbr, M_.endo_names{orig_var_index}, M_.endo_names{orig_var_index}, x(problemcol(ii))); else - eq_nbr = problemrow(ii)-aux_eq_nbr; + eq_nbr = problemrow(ii)-initial_aux_eq_nbr; orig_var_index = M_.aux_vars(1,problemcol(ii)-M_.orig_endo_nbr).orig_index; fprintf('Derivative of Equation %d with respect to Variable %s (initial value of %s: %g) \n', ... eq_nbr, M_.endo_names{orig_var_index}, M_.endo_names{orig_var_index}, x(problemcol(ii)));