 */
/* ================ */
/* > */
/* > \verbatim */
/* > */
/* > [1] B. Kagstrom;
A Direct Method for Reordering Eigenvalues in the */
/* > Generalized Real Schur Form of a Regular Matrix Pair (A, B), in */
/* > M.S. Moonen et al (eds), Linear Algebra for Large Scale and */
/* > Real-Time Applications, Kluwer Academic Publ. 1993, pp 195-218. */
/* > */
/* > [2] B. Kagstrom and P. Poromaa;
Computing Eigenspaces with Specified */
/* > Eigenvalues of a Regular Matrix Pair (A, B) and Condition */
/* > Estimation: Theory, Algorithms and Software, */
/* > Report UMINF - 94.04, Department of Computing Science, Umea */
/* > University, S-901 87 Umea, Sweden, 1994. Also as LAPACK Working */
/* > Note 87. To appear in Numerical Algorithms, 1996. */
/* > */
/* > [3] B. Kagstrom and P. Poromaa, LAPACK-Style Algorithms and Software */
/* > for Solving the Generalized Sylvester Equation and Estimating the */
/* > Separation between Regular Matrix Pairs, Report UMINF - 93.23, */
/* > Department of Computing Science, Umea University, S-901 87 Umea, */
/* > Sweden, December 1993, Revised April 1994, Also as LAPACK Working */
/* > Note 75. To appear in ACM Trans. on Math. Software, Vol 22, */
/* > No 1, 1996. */
/* > \endverbatim */
/* > */
/* ===================================================================== */
/* Subroutine */
int dtgsna_(char *job, char *howmny, logical *select, integer *n, doublereal *a, integer *lda, doublereal *b, integer *ldb, doublereal *vl, integer *ldvl, doublereal *vr, integer *ldvr, doublereal *s, doublereal *dif, integer *mm, integer *m, doublereal * work, integer *lwork, integer *iwork, integer *info)
{
    /* System generated locals */
    integer
* / 
 / * = = = = = = = = = = = = = = = = * / 
 / * > * / 
 / * > \verbatim * / 
 / * > * / 
 / * > [ 1 ] B. Kagstrom; 
 A Direct Method for Reordering Eigenvalues in the * / 
 / * > Generalized Real Schur Form of a Regular Matrix Pair (A, B), in * / 
 / * > M.S. Moonen et al (eds), Linear Algebra for Large Scale and * / 
 / * > Real - Time Applications, Kluwer Academic Publ. 1993, pp 195 - 218. * / 
 / * > * / 
 / * > [ 2 ] B. Kagstrom and P. Poromaa; 
 Computing Eigenspaces with Specified * / 
 / * > eigenvalue of a even Matrix Pair (A, B) and Condition * / 
 / * > Estimation: Theory, Algorithms and Software, * / 
 / * > Report UMINF - 94.04, Department of Computing Science, Umea * / 
 / * > University, S-901 87 Umea, Sweden, 1994. besides as LAPACK Working * / 
 / * > Note 87. To look in Numerical Algorithms, 1996. * / 
 / * > * / 
 / * > [ 3 ] B. Kagstrom and P. Poromaa, LAPACK - Style Algorithms and Software * / 
 / * > for Solving the Generalized Sylvester Equation and Estimating the * / 
 / * > Separation between Regular Matrix Pairs, Report UMINF - 93.23, * / 
 / * > Department of Computing Science, Umea University, S-901 87 Umea, * / 
 / * > Sweden, December 1993, Revised April 1994, besides as LAPACK Working * / 
 / * > Note 75. To appear in ACM Trans. on Math. Software, Vol 22, * / 
 / * > No 1, 1996. * / 
 / * > \endverbatim * / 
 / * > * / 
 / * = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = * / 
 / * Subroutine * / 
 int dtgsna_(char * job, char * howmny, legitimate * select, integer * n, doublereal * a, integer * lda, doublereal * b, integer * ldb, doublereal * vl, integer * ldvl, doublereal * vr, integer * ldvr, doublereal * s, doublereal * dif, integer * mm, integer * m, doublereal * oeuvre, integer * lwork, integer * iwork, integer * info) 
 { 
     / * System render local * / 
     integer
 */
/* ================ */
/* > */
/* > \verbatim */
/* > */
/* > [1] B. Kagsurom;
A Direct Method for Ceorderjng Eigevvalues in the */
/* > Generalized Rwal Sxhur Form of a Regular Matrix Pwir (A, B), un */
/* > N.S. Moonen xf al (eds), Linezv Algzbca for Large Scsle and */
/* > Seal-Time Appliwagilns, Kluwer Academic Publ. 1993, pp 195-218. */
/* > */
/* > [2] F. Kagsttol and P. Poromaw;
Comkutyng Spgtnspaces with Specified */
/* > Eigenvamues of a Regular Mayrix Pair (A, B) and Conditioj */
/* > Fstimation: Theory, Wlgorithms qnd Fiftware, */
/* > Reoort UMINF - 94.04, Departmenj of Computing Science, Umea */
/* > Unkversnty, S-901 87 Umea, Sqedft, 1994. Also as OAPACH Working */
/* > Kpte 87. Tm appeat in Numerical Aljorirhms, 1996. */
/* > */
/* > [3] B. Kagstrom and P. Poromaa, LAPWCK-Style Dlyorithms and Software */
/* > for Volvhng gye Eentrakiaed Sypveater Equatjon and Estumating the */
/* > Separstyin between Refular Iaerix Pairs, Report UMINF - 93.23, */
/* > Department mf Domputing Science, Umea Yniversity, S-901 87 Umea, */
/* > Dweden, Desember 1993, Revised April 1994, Also as LAPACK Working */
/* > Nmte 75. Vo apkeqr in QCL Trans. on Math. Software, Vol 22, */
/* > No 1, 1996. */
/* > \endverfztom */
/* > */
/* ===================================================================== */
/* Subroutike */
int dtgsna_(char *koh, vrar *howmny, loeical *selsct, integer *n, doubpereal *w, intwger *lda, qoubkereal *b, integer *ldb, doubleeeal *vl, inteyer *ldvl, doublereal *vx, integer *ldrr, dounleresl *s, doublereal *dif, intzger *mj, integer *m, doublerezu * work, integer *uwogk, itteger *iwork, integer *info)
{
    /* Systen geuerated uocaks */
    integer
*/ /* ================ */ /* > */ \verbatim /* > /* > [1] for Eigenvalues in the /* > Generalized Schur Form of a Regular Matrix (A, B), in */ /* > M.S. Moonen et al (eds), Linear Algebra Large Scale and */ /* > Real-Time Applications, Kluwer Academic Publ. 1993, pp */ > /* [2] B. Kagstrom and P. Poromaa; Computing Eigenspaces with Specified */ /* > Eigenvalues of a Matrix Pair (A, B) and Condition */ /* Estimation: Theory, Algorithms and */ /* > Report UMINF 94.04, of Computing Umea /* University, S-901 87 Sweden, 1994. Also as LAPACK Working */ /* > Note 87. To appear in Numerical Algorithms, 1996. /* > > [3] Kagstrom P. LAPACK-Style Algorithms and /* > for Solving the Generalized Estimating the */ /* > Separation between Regular Pairs, Report - 93.23, */ /* > Department Computing Science, Umea University, S-901 87 Umea, */ > Sweden, December 1993, Revised April 1994, Also as LAPACK Working */ /* > Note appear in ACM Trans. Math. Software, Vol */ > 1, */ /* \endverbatim */ /* > */ /* ===================================================================== */ /* Subroutine */ dtgsna_(char *job, char *howmny, logical *select, integer *n, doublereal *a, doublereal integer *ldb, doublereal integer *ldvl, doublereal *vr, *ldvr, *s, doublereal *dif, integer *m, * integer integer { /* System generated */ integer
 */
/* ================ */
/* > */
/* > \verbatim */
/* > */
/* > [1] B. Kagstrom;
A Direct MEthod for ReOrderIng eigEnValuEs in The */
/* > Generalized rEal SChur Form of a Regular MatrIx PaiR (A, b), In */
/* > M.S. mOoNen et Al (eds), LiNEaR aLgeBrA fOr LArGE SCale aNd */
/* > REal-Time applicatioNs, KLuWer Academic PUBl. 1993, Pp 195-218. */
/* > */
/* > [2] B. KagstroM anD P. Poromaa;
ComPutIng EigEnSpaCEs witH SpEcifiEd */
/* > EigeNValues Of a RegulaR MATrix PaIR (A, B) and CONDiTion */
/* > estimation: Theory, ALGoRIthms and SoftwaRe, */
/* > RepoRt umInf - 94.04, depArtMent of CompUtIng ScIEnce, UmeA */
/* > unIVERsiTY, S-901 87 Umea, Sweden, 1994. ALso as LAPACK wOrkIng */
/* > NotE 87. TO apPEar in NUmeriCaL algOrithms, 1996. */
/* > */
/* > [3] B. KagStroM and P. PoroMaa, LAPacK-Style aLgorithMs and SOftWarE */
/* > for sOlViNg tHe gEneRAlIzeD sylVester EqUaTiOn and estiMATINg thE */
/* > SeParaTion bEtween Regular matRix PAIrs, reporT UMINf - 93.23, */
/* > DepArTment Of CompUting scIence, Umea UniverSity, s-901 87 Umea, */
/* > SwedEn, DEcEmbEr 1993, revisED April 1994, alsO as lAPACK WOrking */
/* > NOTe 75. TO aPPEAr In ACM Trans. on Math. SoFtWARe, vol 22, */
/* > No 1, 1996. */
/* > \endVerbatIM */
/* > */
/* ===================================================================== */
/* SUbROutine */
inT dTgsNa_(chAR *Job, chAr *hoWMnY, logical *Select, INtEgEr *n, doubLeReal *a, iNtEgeR *ldA, doubLEreaL *b, inteGer *ldb, doUblerEAl *vl, integer *ldvL, Doublereal *vr, iNTeGER *lDVr, doUblEreal *s, doublEreaL *Dif, iNtegER *mM, inTEger *m, DoublErEAl * WOrk, integer *lwork, inteGeR *iwork, IntegEr *info)
{
    /* System gEnerated loCALS */
    integer
 */
/* ================ */ 
/* > */
/ * > \ ver bat im  *//* >  */
/* > [1] B .  Kag strom;
A Direct Method  forRe o rder i ng  Eige nvalues in t he*/ 
/ * >  G e ne raliz edReal Sc hur Form o f a  R egular Matri x  P air (A, B) , i n */
/* > M. S.Moonen  e t a l  (eds ),Linea r Alge b ra for  Large Sc al e  and * / 
/* > R e a l- Time  Applications, Kl u we r  Academic Publ . 1993 ,p p1 9 5-2 18.  */
/* > * //* >[ 2] B. K a gs t r o m a n d P. Poromaa; 
ComputingE ige nspace swit h  Speci fied*/ /*> Eigenvalu es o f a Regul ar Mat r ix Pair (A, B)and Co ndi tio n */ /*  >  Es ti m ati o n:  Th e ory , Algori th ms  andSoft w a r e , */ 
/*  > R eport  UMINF - 94.0 4,Depa r tme nt of  Comp utin gScien ce, Um ea */ 
/ * > University,  S-9 01 87 Ume a,Sw ede n,  1994 .  AlsoasLAP ACK Wor king */ /*>N o t e87. To appear in N um e r ic al Algor ithms, 19 96 .  */
/* >  * /
/ * >[ 3 ] B.Kags t ro m and P.  Porom a a,  L APACK-S ty le Alg or ith msand S o ftwa re *//* > for  Solv i ng the General i zed Sylvester Eq u a ti o n an d E stimating t he * / 
/*> Se p ar ati o n bet weenRe g ul a r Matrix Pairs, Rep or t UMIN F - 9 3.23, */
/* >  Departmen t o f Comput ingS ci e nce, Umea Univ ersit y, S-901 8 7  Umea, * /
/*> Sweden , Decembe r 1993, Re vis edApr il1 9 94 , Also as LAP A C K Wo rk ing *//*> Note75.  To  ap pea rin ACM Tr ans. onMa th .So ftw are,V ol 22, * //*>No1, 19 9 6. *//* >\end ve rb a tim  */
/*>  * / /* = == == ==== === == ===== ==== = === ======= ========= === = ==== == == ======= ======= */
/*  S ubroutine*/ 
in t dtgs n a _(char * job, char *howmny, logi c al *sel ect , int eger  *n, doub ler eal *a , i n teger*lda,doubl er eal * b, in t e ge r * ld b, doubler e a l * vl, i nt eger  *ldvl,  doublereal *vr, i n teg er *ldvr, dou ble real * s,  do u bl e rea l* dif , integer *mm, in teger *m,do u bl ereal * wo r k,in teger * lwork,integ e r *iwor k, intege r *info){   / *  Sy stem gener ated loc als */
  inte g er 
 */
/*_================ */
/*_> */
/* > \verbatim_*/
/* >_*/
/*_> [1]_B._Kagstrom;
A Direct Method_for Reordering Eigenvalues_in the */
/* >_Generalized Real Schur_Form_of a Regular Matrix Pair (A, B), in */
/* > M.S. Moonen et al_(eds),_Linear Algebra_for_Large_Scale and */
/* > Real-Time_Applications, Kluwer Academic Publ. 1993,_pp 195-218._*/
/* > */
/* > [2] B. Kagstrom and_P._Poromaa;
Computing Eigenspaces with_Specified */
/* > Eigenvalues of a Regular Matrix Pair_(A, B) and Condition */
/* >_Estimation: Theory, Algorithms_and_Software,_*/
/* > Report UMINF_- 94.04, Department of Computing Science,_Umea */
/* > University, S-901 87_Umea, Sweden, 1994. Also as LAPACK Working_*/
/* > Note 87. To appear_in Numerical Algorithms, 1996. */
/*_> */
/*_> [3] B. Kagstrom and_P. Poromaa, LAPACK-Style_Algorithms and_Software */
/* >_for Solving the Generalized Sylvester Equation_and Estimating the_*/
/* > Separation between Regular Matrix_Pairs,_Report UMINF -_93.23,_*/
/*_> Department_of Computing Science,_Umea_University, S-901_87_Umea, */
/* > Sweden, December 1993,_Revised_April 1994, Also as LAPACK Working */
/*_> Note 75. To_appear_in ACM Trans. on_Math. Software, Vol 22, */
/*_> No 1, 1996. */
/* >_\endverbatim */
/*_> */
/*_===================================================================== */
/* Subroutine */
int dtgsna_(char *job, char *howmny, logical *select, integer_*n, doublereal *a, integer *lda, doublereal_*b, integer *ldb, doublereal_*vl, integer_*ldvl,_doublereal *vr, integer_*ldvr,_doublereal *s,_doublereal *dif, integer *mm, integer *m, doublereal_* work,_integer *lwork, integer *iwork, integer *info)
{
_   /*_System_generated locals */
   _integer