4.1: solution for QZ in Octave (to be tested)
git-svn-id: https://www.dynare.org/svn/dynare/trunk@3202 ac1d8469-bf42-47a9-8791-bf33cf982152time-shift
parent
c402bd8969
commit
cf933608d0
|
@ -39,11 +39,8 @@ function [ss,tt,w,sdim,eigval,info] = mjdgges(e,d,qz_criterium)
|
||||||
% along with Dynare. If not, see <http://www.gnu.org/licenses/>.
|
% along with Dynare. If not, see <http://www.gnu.org/licenses/>.
|
||||||
|
|
||||||
if exist('OCTAVE_VERSION')
|
if exist('OCTAVE_VERSION')
|
||||||
% Octave can only compute a real QZ with the function qz()
|
[ss,tt,w,eigval]=qz(e,d,'S');
|
||||||
% while by default MATLAB computes a complex QZ
|
sdim =sum(abs(eigval) < qz_criterium);
|
||||||
% Here we need a complex QZ (so that ss and tt are uppeer
|
|
||||||
% triangular), so fail under Octave
|
|
||||||
error('You must compile the MEX file "mjddges", I can''t solve the model without that.')
|
|
||||||
end
|
end
|
||||||
|
|
||||||
% Chek number of inputs and outputs.
|
% Chek number of inputs and outputs.
|
||||||
|
@ -82,4 +79,4 @@ try
|
||||||
sdim = sum( abs(eigval) < qz_criterium );
|
sdim = sum( abs(eigval) < qz_criterium );
|
||||||
catch
|
catch
|
||||||
info = 1;% Not as precise as lapack's info!
|
info = 1;% Not as precise as lapack's info!
|
||||||
end
|
end
|
||||||
|
|
Loading…
Reference in New Issue