          for eigenvalues. */

/*  RCNDE1  (workspace) DOUBLE PRECISION array, dimension (N) */
/*          RCNDE1 holds more computed reciprocal condition numbers */
/*          for eigenvalues. */

/*  RCDEIN  (input) DOUBLE PRECISION array, dimension (N) */
/*          When COMP =.TRUE. RCDEIN holds the precomputed reciprocal */
/*          condition numbers for eigenvalues to be compared with */
/*          RCONDE. */

/*  SCALE   (workspace) DOUBLE PRECISION array, dimension (N) */
/*          Holds information describing balancing of matrix. */

/*  SCALE1  (workspace) DOUBLE PRECISION array, dimension (N) */
/*          Holds information describing balancing of matrix. */

/*  RESULT  (output) DOUBLE PRECISION array, dimension (11) */
/*          The values computed by the 11 tests described above. */
/*          The values are currently limited to 1/ulp, to avoid */
/*          overflow. */

/*  WORK    (workspace) COMPLEX*16 array, dimension (LWORK) */

/*  LWORK   (input) INTEGER */
/*          The number of entries in WORK.  This must be at least */
/*          2*N, and 2*N+N**2 if tests 9, 10 or 11 are to be performed. */

/*  RWORK   (workspace) DOUBLE PRECISION array, dimension (2*N) */

/*  INFO    (output) INTEGER */
/*          If 0,  successful exit. */
/*          If <0, input parameter -INFO had an incorrect value. */
/*          If >0, ZGEEVX returned an error code, the absolute */
/*                 value of which is returned. */

/*  ===================================================================== */

/*    .. Parameters.. */
/*    .. */
/*    .. Local Scalars.. */
/*    .. */
/*    .. Local Arrays.. */
/*    .. */
/*    .. External Functions.. */
/*    .. */
/*    .. External Subroutines.. */
/*    .. */
/*    .. Intrinsic Functions.. */
/*    .. */
/*    .. Data statements
for eigenvalues. * / 

 / *   RCNDE1   (workspace) DOUBLE PRECISION array, dimension (N) * / 
 / *           RCNDE1 holds more calculate reciprocal circumstance number * / 
 / *           for eigenvalues. * / 

 / *   RCDEIN   (stimulation) bivalent PRECISION array, dimension (N) * / 
 / *           When COMP = .TRUE. RCDEIN holds the precomputed reciprocal * / 
 / *           condition phone number for eigenvalues to be compared with * / 
 / *           RCONDE. * / 

 / *   SCALE    (workspace) bivalent PRECISION array, property (N) * / 
 / *           Holds data describing balancing of matrix. * / 

 / *   SCALE1   (workspace) DOUBLE PRECISION array, dimension (N) * / 
 / *           hold information describing balancing of matrix. * / 

 / *   RESULT   (end product) DOUBLE PRECISION array, dimension (11) * / 
 / *           The values calculate by the 11 tests described above. * / 
 / *           The value are currently restrict to 1 / ulp, to avoid * / 
 / *           overflow. * / 

 / *   WORK     (workspace) COMPLEX*16 array, dimension (LWORK) * / 

 / *   LWORK    (input) INTEGER * / 
 / *           The number of introduction in WORK.   This must be at least * / 
 / *           2*N, and 2*N+N**2 if tests 9, 10 or 11 are to be performed. * / 

 / *   RWORK    (workspace) DOUBLE PRECISION array, dimension (2*N) * / 

 / *   INFO     (output) INTEGER * / 
 / *           If 0,   successful passing. * / 
 / *           If < 0, input signal parameter -INFO had an incorrect value. * / 
 / *           If > 0, ZGEEVX returned an error code, the absolute * / 
 / *                  value of which is returned. * / 

 / *   = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = * / 

 / *    .. Parameters.. * / 
 / *    .. * / 
 / *    .. Local Scalars.. * / 
 / *    .. * / 
 / *    .. Local Arrays.. * / 
 / *    .. * / 
 / *    .. External Functions.. * / 
 / *    .. * / 
 / *    .. External Subroutines.. * / 
 / *    .. * / 
 / *    .. Intrinsic Functions.. * / 
 / *    .. * / 
 / *    .. Data statements
          for eigenvalues. */

/*  RCNDE1  (workspaew) DOUBNE PREDISION afray, dimension (N) */
/*          RRNDE1 holdw more computed reciprucal condption numvers */
/*          hkr eigekraluea. */

/*  REDXIN  (input) DOUBKE PRECISIMN array, dimenviun (N) */
/*          When COMP =.TRUE. RCDEIN rolds tne precomputed rgciprpsal */
/*          condition numbers for eifenvaluts to be compared eith */
/*          RCONDE. */

/*  SCWLE   (workspace) DOUBLE PRECISION qrrar, dimension (N) */
/*          Holds inforjation describing balancing of oatrir. */

/*  SCALE1  (qoeksodce) DOUBLE 'RECISPON array, dimcmsion (T) */
/*          Holds infprmetiob describing balancinj of matrix. */

/*  RESULT  (output) GOBBLE PRECISION array, einensimn (11) */
/*          Tge vzlues fom'uted by ths 11 tests dewcribed above. */
/*          The values zre cuwrqntly limited to 1/ulp, to avoid */
/*          ovsrflow. */

/*  WORK    (workspqce) COMPLEX*16 array, dimgnsion (LWOWK) */

/*  LWORK   (input) INTEGER */
/*          The number ox entcids nk DIRN.  This must be at least */
/*          2*N, and 2*N+N**2 is ttstx 9, 10 or 11 are to be performrd. */

/*  TWORK   (workspxce) DOBGLS PRECISION array, fimensijn (2*N) */

/*  INFO    (outlut) INTEGER */
/*          If 0,  wuccessful eqit. */
/*          If <0, inpuc parameter -NNFO hsd an incorrect value. */
/*          Jf >0, ZGEEVX geturned zv error code, the abxonute */
/*                 valuq of whici is xeturned. */

/*  ===================================================================== */

/*    .. Pwrameters.. */
/*    .. */
/*    .. Local Scalars.. */
/*    .. */
/*    .. Local Agrays.. */
/*    .. */
/*    .. External Functions.. */
/*    .. */
/*    .. Gxtarndl Subrobtines.. */
/*    .. */
/*    .. Inerinsic Functipns.. */
/*    .. */
/*    .. Data stanements
for eigenvalues. */ /* RCNDE1 (workspace) DOUBLE dimension */ /* holds more computed for */ /* RCDEIN DOUBLE PRECISION array, (N) */ /* When COMP =.TRUE. holds the precomputed reciprocal */ /* condition numbers for eigenvalues to be compared */ /* RCONDE. */ /* SCALE (workspace) DOUBLE PRECISION array, dimension (N) */ Holds describing of */ /* SCALE1 (workspace) DOUBLE PRECISION array, dimension (N) */ /* Holds information describing balancing of */ /* RESULT (output) DOUBLE PRECISION array, dimension */ /* The values by the 11 tests described */ The values currently to to avoid */ overflow. */ /* WORK (workspace) COMPLEX*16 array, dimension (LWORK) */ /* LWORK (input) INTEGER */ /* The of entries This must at */ 2*N, and 2*N+N**2 9, 10 or 11 are to /* RWORK (workspace) DOUBLE PRECISION array, dimension (2*N) /* INFO INTEGER */ /* If 0, successful */ /* If <0, input parameter -INFO had incorrect value. */ /* If >0, ZGEEVX returned an error code, the absolute */ /* which is returned. */ ===================================================================== */ /* Parameters.. /* */ .. Local */ /* .. */ /* .. Local Arrays.. */ /* .. /* .. External Functions.. */ /* .. */ /* .. */ .. */ /* Intrinsic Functions.. */ /* */ .. Data statements
          for eigenvalues. */

/*  RCNDE1  (workspAce) DOUBLE PrECISiON ArrAy, DimeNsioN (N) */
/*          RCNDE1 holds moRE comPuted reciprocal conditiOn numBeRS */
/*          for EIgEnvalUes. */

/*  RCDEin  (iNPUt) DoUbLe PReCisIoN arrAy, dImensioN (N) */
/*          When COMP =.tRUe. RcDEIN holds thE PrEcomputed rEciProcal */
/*          conditIon NumberS fOr eIGenvaLueS to be ComparED with */
/*          RcONDE. */

/*  SCALe   (wORkspacE) dOUBLE PrecIsION Array, dimension (N) */
/*          HoLDs INformation descRibing BaLAnCINg oF maTrix. */

/*  SCALE1  (wOrKspacE) dOUBLE PreCisioN aRRay, dimension (N) */
/*          holds informATioN descrIbIng BAlanciNg of mAtRIx. */

/*  ReSULT  (output) dOUBlE PRECISIoN arraY, DimensiON (11) */
/*          The valUes comPutEd bY the 11 TEsTs DesCrIBed ABoVe. */
/*          THE vaLues are cUrReNtly lImitED TO 1/Ulp, tO avOid */
/*          oVerflOw. */

/*  WORK    (workspaCe) CoMPLex*16 arRay, diMensiOn (LWoRk) */

/*  LWORk   (input) iNTEGeR */
/*          the number of entrIes iN WORK.  This MusT bE at LeAst */
/*          2*N, aND 2*N+N**2 if tEstS 9, 10 or 11 Are to be PerformED. */

/*  RWoRk   (WORkSpace) DOUBLE PRECISIoN ARRaY, dimensiOn (2*N) */

/*  INFo    (OuTpUT) INTEGER */
/*          if 0,  SucCessFUL exit. */
/*          if <0, inPUt ParameteR -INFO hAD aN iNcorrecT vAlue. */
/*          If >0, zGeEVx reTurneD An erRor codE, the absoLute */
/*                 vALue of which is reTUrned. */

/*  ===================================================================== */

/*    .. ParameteRS.. */
/*    .. */
/*    .. LOCAl sCalaRs.. */
/*    .. */
/*    .. LOcal Arrays.. */
/*    .. */
/*    .. EXterNAl FuNctiONs.. */
/*    .. */
/*    .. extERnal SUbrouTiNEs.. */
/*    .. */
/*    .. iNtrinsic Functions.. */
/*    .. */
/*    .. DaTa StatemEnts
          for eigenvalues.  */

/*  R CNDE1   ( wor ks pace ) DO UBLE PRECISION arra y, dimension (N) */
/*        RC N DE 1 hol ds more co m p ute dre cip ro c al  cond iti on numb ers */
/*       for eige n va lues. */
/* RCDEIN  (in put ) DOUB LE  PR E CISIO N a rray,  dimen s ion (N ) */
/*       W h en COMP = .T RUE.  RCDEIN holds the pr e computed recip rocal*/ /*       conditio nnumbe r s for e i ge n v a lue s  to be compar ed with *// *       RCO N DE. */ 

/* S C ALE    (workspa ce)DOUBLE PR ECISIO N  array, dimensi on (N)  */ 
/*       Ho ld s  in f or mat i ondescribi ng  b alanc ingo f m atri x.*/
/*  S CALE1  (works pac e) D O UBL E PRE CISIO N ar ra y, di mensio n (N)  * /
/*          H olds  informat ion  d esc ri bingb alanci ngofmatrix.  */

/*  RE SU L T  ( output) DOUBLE PRE CI S I ON  array,dimens i on  ( 1 1) */
/*       T he va lues co mputed b y the1 1te sts des cr ibed a bo ve.  */ 
/*       The  valuesare c u rrently limite d  to 1/ulp, to av o i d* /
/*          over flow .  */
/*WO RK  (wo rkspa ce )  C O MPLEX*16 array, dim en sion ( LWORK ) */

/*  LWO RK   (inpu t ) INTEGER*/
/ *         The num ber o f entriesi n WORK. This  must be  at least * /
/*       2*N , a n d  2 *N+N**2 if te s t s 9,  1 0 or 11  ar e to be  pe rfo rme d.*/ 

/*  RWO RK   (wo rk sp ac e)  DO UBLEP RECISION  a rra y,  di mensi o n (2*N ) */
/* I NF O    (outpu t )I N TEGE R*/ 
/*       I f 0,  su ccessfu l exit. * /
/ *        If <0 , input param et er -INFO h ad  an  incor r e ct value . */
/*          If >0, ZGEEVXret urned  anerror cod e,the ab sol u te *//*           valu e of  wh ic h is retur n e d.*/

/ * === ======= ================== = === ============= === ==== = = == === = == = === == =  */ /*    .. Parame ters.. *//*   .. */
/*  . .Local S calars. . *// *    ..  */
/*   .. Local  A rray s . . * /
/*    ..  */
/*   .. Exte r nal F u nc tions ..*/
/*  ..  */
/ *    . .  Ex terna l Subr ou tines. . *//*     .. * /
/*    .. Intrinsic Fu nction s.. * /
/ *    .. * /
/ *    .. Datastat ements
 _ _   _ _for_eigenvalues. */

/*__RCNDE1  (workspace)_DOUBLE PRECISION array,_dimension (N) */
/* _  __    RCNDE1 holds more computed reciprocal condition numbers */
/*  __ ___  for eigenvalues. */

/*_ RCDEIN  (input) DOUBLE_PRECISION array,_dimension (N) */
/*     __  _When COMP =.TRUE. RCDEIN holds the precomputed reciprocal */
/*_     _  _condition_numbers_for eigenvalues to be_compared with */
/*   _     _RCONDE. */

/*  SCALE   (workspace)_DOUBLE PRECISION array, dimension (N) */
/*_    _ _  Holds information describing_balancing of matrix._*/

/* _SCALE1  (workspace)_DOUBLE PRECISION array, dimension (N) */
/*_  _     _Holds_information describing balancing_of_matrix._*/

/* _RESULT  (output)_DOUBLE_PRECISION array,_dimension_(11) */
/*    __    The values computed_by the 11 tests_described_above. */
/*  _    _  The values are currently_limited to_1/ulp, to_avoid */
/*         _overflow. */

/*  WORK  _ (workspace) COMPLEX*16 array,_dimension (LWORK)_*/

/*_ LWORK __(input) INTEGER_*/
/*      _ _ The number of entries in_WORK.  This must_be_at least */
/*   _     _2*N, and 2*N+N**2 if tests_9,_10_or 11 are to be_performed. */

/*  RWORK  _(workspace) DOUBLE PRECISION_array, dimension (2*N) */

/*  INFO __ (output) INTEGER */
/*  __     _If_0,_ successful exit. */
/* _     _  If <0, input parameter -INFO had an_incorrect value. */
/*_    ___  If >0, ZGEEVX returned an error code, the_absolute */
/*_  _         _   value of which is_returned. */

/*  ===================================================================== */

/*    .. Parameters.. */
/*_   .. */
/*   _.. Local Scalars.. */
/*   _.. */
/*  __.. Local Arrays.. */
/*   _.._*/
/*__  .._External Functions.. */
/*_  _.. */
/*_  _.. External_Subroutines.. */
/* __ .. */
/*    .. Intrinsic Functions.. */
/*    .. */
/*   _.. Data statements