Merge pull request #464 from JohannesPfeifer/steady_display
Bug fixes for display of problems with steady statetime-shift
commit
3c8daf1cd7
|
@ -65,8 +65,15 @@ end
|
||||||
[dr.ys,params,check1]=evaluate_steady_state(oo.steady_state,M,options,oo,1);
|
[dr.ys,params,check1]=evaluate_steady_state(oo.steady_state,M,options,oo,1);
|
||||||
|
|
||||||
% testing for problem
|
% testing for problem
|
||||||
if check1
|
if check1(1)
|
||||||
disp('model diagnostic can''t obtain the steady state')
|
disp('model diagnostic can''t obtain the steady state')
|
||||||
|
if any(isnan(dr.ys))
|
||||||
|
disp(['model diagnostic obtains a steady state with NaNs'])
|
||||||
|
end
|
||||||
|
if any(isinf(dr.ys))
|
||||||
|
disp(['model diagnostic obtains a steady state with Inf'])
|
||||||
|
end
|
||||||
|
return;
|
||||||
end
|
end
|
||||||
|
|
||||||
if ~isreal(dr.ys)
|
if ~isreal(dr.ys)
|
||||||
|
|
|
@ -79,9 +79,9 @@ if info(1)
|
||||||
end
|
end
|
||||||
|
|
||||||
[steady_state,M_.params,info] = steady_(M_,options_,oo_);
|
[steady_state,M_.params,info] = steady_(M_,options_,oo_);
|
||||||
|
oo_.steady_state = steady_state;
|
||||||
|
|
||||||
if info(1) == 0
|
if info(1) == 0
|
||||||
oo_.steady_state = steady_state;
|
|
||||||
if options_.noprint == 0
|
if options_.noprint == 0
|
||||||
disp_steady_state(M_,oo_);
|
disp_steady_state(M_,oo_);
|
||||||
end
|
end
|
||||||
|
|
Loading…
Reference in New Issue