gschur() was declared obsolete in Scilab 4, now removed
[scilab.git] / scilab / modules / compatibility_functions / macros / mtlb_qz.sci
1 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
2 // Copyright (C) ???? - INRIA - Scilab
3 //
4 // Copyright (C) 2012 - 2016 - Scilab Enterprises
5 //
6 // This file is hereby licensed under the terms of the GNU GPL v2.0,
7 // pursuant to article 5.3.4 of the CeCILL v.2.1.
8 // This file was originally licensed under the terms of the CeCILL v2.1,
9 // and continues to be available under such terms.
10 // For more information, see the COPYING file which you should have received
11 // along with this program.
12
13 function [AA, BB, Q, Z, V] = mtlb_qz(A,B)
14     [lhs,rhs] = argn(0)
15     [AA, BB, Q, Z] = schur(A,B)
16     if lhs==5 then
17         n  = size(A,1)
18         LA = diag(AA)
19         LB = diag(BB)
20         V  = zeros(n,n)
21         for k = 1:n
22             K = kernel(A*LB(k)-B*LA(k))
23             if size(K,2)<>1 then
24                 error(msprintf(gettext("%s: pencil is not diagonalizable.\n"),"mtlb_qz"))
25             end
26             V(:,k)=K
27         end
28     end
29 endfunction