/*--------------------------------------------------------------------------
**  name this file PsZtgSen.tst
**
**
**
**
**
**
**
**
**
**
**
**
**  Paul Soderlind, 18 April 1997
----------------------------------------------------------------------------*/



output file = PsZtgSen.out reset;
/*-----------------------------------------------------------------------*/


"\l\l------------------Example 1------------------------";
format 8,5;
pii_sii = 1;
GtIt = 1;
"\lpii_sii and GtIt: ";;pii_sii~GtIt;

A =
  {  2.1385   -0.5344    0.8260,
   -0.5344    0.1372   -0.1815,
    0.8260   -0.1815    0.4875 };
"\lA ";A;
A = complex( A, zeros(rows(A),rows(A)) );
B =
  {  0.6834    0.8023   -1.0064,
    0.8023    0.9419   -1.1732,
   -1.0064   -1.1732    1.9998 };
"\lB: ";B;
B = complex( B, zeros(rows(B),rows(B)) );

dlibrary -a PsZtgSen;
{ S,P,Q,Z } = ZtgSen( A,B,1.01,pii_sii,GtIt );

"\lS: ";;S;
"\lP: ";;P;
"\lQ: ";;Q;
"\lZ: ";;Z;

"\lreal( Q*S*Z' ) "; real( Q*S*Z' );
"\lreal( Q*P*Z' )";  real( Q*P*Z' );
"\lreal( Q'A*Z );"; real( Q'A*Z );
"\lreal( Q'B*Z );"; real( Q'B*Z );

"\l\l------------------Example 2------------------------";
pii_sii = 1;
GtIt = 1;
"\lpii_sii and GtIt: ";;pii_sii~GtIt;

A = eye(3);
"\lA ";A;
A = complex( A, zeros(rows(A),rows(A)) );
B =
  {  0.6834    0.8023   -1.0064,
    0.8023    0.9419   -1.1732,
   -1.0064   -1.1732    1.9998 };
"\lB: ";B;
B = complex( B, zeros(rows(B),rows(B)) );

{ S,P,Q,Z } = ZtgSen( A,B,1.01,pii_sii,GtIt );

"\lSii/Pii and Pii/Sii";
lambda = (diag(S)./diag(P));
lambda~(1./lambda);

"\lS: ";;S;
"\lP: ";;P;
"\lQ: ";;Q;
"\lZ: ";;Z;

"\lreal( Q*S*Z' ) "; real( Q*S*Z' );
"\lreal( Q*P*Z' )";  real( Q*P*Z' );
"\lreal( Q'A*Z );"; real( Q'A*Z );
"\lreal( Q'B*Z );"; real( Q'B*Z );


"\lEigenvalue decomposition";
{ lambda, Eigenvectors } = eigv( B );
lambda';
Eigenvectors;
/*-----------------------------------------------------------------------*/

end;