2011-10-12 21:46:50 +02:00
function [ys,params,info] = evaluate_steady_state ( ys_init,M,options,oo,steadystate_check_flag)
2012-04-21 10:12:00 +02:00
% function [ys,params,info] = evaluate_steady_state(ys_init,M,options,oo,steadystate_check_flag)
2012-03-06 19:05:59 +01:00
% Computes the steady state
%
2011-10-12 21:46:50 +02:00
% INPUTS
% ys_init vector initial values used to compute the steady
% state
% M struct model structure
% options struct options
% oo struct output results
% steadystate_check_flag boolean if true, check that the
% steadystate verifies the
2012-03-06 19:05:59 +01:00
% static model
%
2011-10-12 21:46:50 +02:00
% OUTPUTS
% ys vector steady state
% params vector model parameters possibly
% modified by user steadystate
% function
% info 2x1 vector error codes
%
% SPECIAL REQUIREMENTS
% none
2013-04-11 17:07:39 +02:00
% Copyright (C) 2001-2013 Dynare Team
2011-10-12 21:46:50 +02:00
%
% 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/>.
info = 0 ;
check = 0 ;
2012-03-06 19:05:59 +01:00
2011-10-12 21:46:50 +02:00
steadystate_flag = options . steadystate_flag ;
params = M . params ;
exo_ss = [ oo . exo_steady_state ; oo . exo_det_steady_state ] ;
2012-03-06 19:05:59 +01:00
2011-10-12 21:46:50 +02:00
if length ( M . aux_vars ) > 0
2012-03-06 19:05:59 +01:00
h_set_auxiliary_variables = str2func ( [ M . fname ' _set_auxiliary_variables' ] ) ;
2011-10-12 21:46:50 +02:00
ys_init = h_set_auxiliary_variables ( ys_init , exo_ss , M . params ) ;
end
2012-03-06 19:05:59 +01:00
2011-10-12 21:46:50 +02:00
if options . ramsey_policy
[ ys , params ] = dyn_ramsey_static ( ys_init , M , options , oo ) ;
elseif steadystate_flag
% explicit steady state file
2012-05-10 20:48:44 +02:00
[ ys , params , info ] = evaluate_steady_state_file ( ys_init , exo_ss , M , ...
2012-04-21 11:39:21 +02:00
options ) ;
if info ( 1 )
return ;
end
2011-10-12 21:46:50 +02:00
elseif ( options . bytecode == 0 && options . block == 0 )
if options . linear == 0
% non linear model
2011-10-13 17:14:41 +02:00
[ ys , check ] = dynare_solve ( [ M . fname ' _static' ] , ...
2011-10-12 21:46:50 +02:00
ys_init , ...
2012-03-06 19:05:59 +01:00
options . jacobian_flag , ...
2011-10-12 21:46:50 +02:00
exo_ss , params ) ;
else
% linear model
2011-10-13 17:14:41 +02:00
fh_static = str2func ( [ M . fname ' _static' ] ) ;
[ fvec , jacob ] = fh_static ( ys_init , exo_ss , ...
2012-04-21 10:12:00 +02:00
params ) ;
2013-06-01 16:32:40 +02:00
ii = find ( ~ isfinite ( fvec ) ) ;
if ~ isempty ( ii )
ys = fvec ;
check = 1 ;
disp ( [ ' STEADY: numerical initial values or parameters incompatible with the following' ...
' equations' ] )
disp ( ii ' )
disp ( ' Check whether your model in truly linear' )
elseif isempty ( ii ) && max ( abs ( fvec ) ) > 1e-12
2011-10-12 21:46:50 +02:00
ys = ys_init - jacob \ fvec ;
else
ys = ys_init ;
end
2013-07-28 00:09:11 +02:00
if options . debug
if any ( any ( isinf ( jacob ) | isnan ( jacob ) ) )
[ infrow , infcol ] = find ( isinf ( jacob ) | isnan ( jacob ) ) ;
fprintf ( ' \nSTEADY: The Jacobian contains Inf or NaN. The problem arises from: \n\n' )
for ii = 1 : length ( infrow )
if infcol ( ii ) < = M . orig_endo_nbr
fprintf ( ' STEADY: Derivative of Equation %d with respect to Variable %s (initial value of %s: %g) \n' , infrow ( ii ) , deblank ( M . endo_names ( infcol ( ii ) , : ) ) , deblank ( M . endo_names ( infcol ( ii ) , : ) ) , ys_init ( infcol ( ii ) ) )
else %auxiliary vars
orig_var_index = M . aux_vars ( 1 , infcol ( ii ) - M . orig_endo_nbr ) . orig_index ;
fprintf ( ' STEADY: Derivative of Equation %d with respect to Variable %s (initial value of %s: %g) \n' , infrow ( ii ) , deblank ( M . endo_names ( orig_var_index , : ) ) , deblank ( M . endo_names ( orig_var_index , : ) ) , ys_init ( infcol ( ii ) ) )
end
end
disp ( ' STEADY: Check whether your model in truly linear\n' )
end
end
2011-10-12 21:46:50 +02:00
end
else
% block or bytecode
2012-04-23 16:57:30 +02:00
[ ys , check ] = dynare_solve_block_or_bytecode ( ys_init , exo_ss , params , ...
options , M ) ;
2011-10-12 21:46:50 +02:00
end
if check
if options . steadystate_flag
info ( 1 ) = 19 ;
resid = check1 ;
else
info ( 1 ) = 20 ;
resid = evaluate_static_model ( ys_init , exo_ss , params , M , options ) ;
end
info ( 2 ) = resid ' * resid ;
2013-06-01 16:32:40 +02:00
if isnan ( info ( 2 ) )
info ( 1 ) = 22 ;
end
2011-10-12 21:46:50 +02:00
return
end
2013-04-11 17:07:39 +02:00
% If some equations are tagged [static] or [dynamic], verify consistency
if M . static_and_dynamic_models_differ
% Evaluate residual of *dynamic* model using the steady state
% computed on the *static* one
z = repmat ( ys , 1 , M . maximum_lead + M . maximum_lag + 1 ) ;
zx = repmat ( [ oo . exo_simul oo . exo_det_simul ] , M . maximum_lead + M . maximum_lag + 1 , 1 ) ;
if options . bytecode
[ chck , r , junk ] = bytecode ( ' dynamic' , ' evaluate' , z , zx , M . params , ys , 1 ) ;
mexErrCheck ( ' bytecode' , chck ) ;
elseif options . block
[ r , data ] = feval ( [ M . fname ' _dynamic' ] , z ' , zx , M . params , ys , M . maximum_lag + 1 , data ) ;
else
iyv = M . lead_lag_incidence ' ;
iyr0 = find ( iyv ( : ) ) ;
xys = z ( iyr0 ) ;
r = feval ( [ M . fname ' _dynamic' ] , z ( iyr0 ) , zx , M . params , ys , M . maximum_lag + 1 ) ;
end
% Fail if residual greater than tolerance
if max ( abs ( r ) ) > options . solve_tolf
info ( 1 ) = 25 ;
return
end
end
2011-10-12 21:46:50 +02:00
if ~ isreal ( ys )
info ( 1 ) = 21 ;
info ( 2 ) = sum ( imag ( ys ) .^ 2 ) ;
ys = real ( ys ) ;
return
end
if ~ isempty ( find ( isnan ( ys ) ) )
info ( 1 ) = 22 ;
info ( 2 ) = NaN ;
return
end