From cf933608d0aed705a89b40fd125410f09e8bfb32 Mon Sep 17 00:00:00 2001 From: michel Date: Mon, 7 Dec 2009 12:25:24 +0000 Subject: [PATCH] 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-bf33cf982152 --- matlab/qz/mjdgges.m | 9 +++------ 1 file changed, 3 insertions(+), 6 deletions(-) diff --git a/matlab/qz/mjdgges.m b/matlab/qz/mjdgges.m index 577255af7..8d188c999 100644 --- a/matlab/qz/mjdgges.m +++ b/matlab/qz/mjdgges.m @@ -39,11 +39,8 @@ function [ss,tt,w,sdim,eigval,info] = mjdgges(e,d,qz_criterium) % along with Dynare. If not, see . if exist('OCTAVE_VERSION') - % Octave can only compute a real QZ with the function qz() - % while by default MATLAB computes a complex QZ - % 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.') + [ss,tt,w,eigval]=qz(e,d,'S'); + sdim =sum(abs(eigval) < qz_criterium); end % Chek number of inputs and outputs. @@ -82,4 +79,4 @@ try sdim = sum( abs(eigval) < qz_criterium ); catch info = 1;% Not as precise as lapack's info! -end \ No newline at end of file +end