CALL SEQUENCE...
       meofvec  A [B] AVAL AVEC INPUT OUTPUT \ 
                      [ -kind VAL1 [VAL2]  -a VAL3 \ 
                          -b  VAL4 -c VAL5 -d VAL6 ]
 
 
      AUTHOR - B.Dugas, mars 17/94.
      LAST REVISION: $Header: /home/cvs/RMNLIB/utils/rdiag/lspgm/meofvec.ptn,v 4.13 2014/09/25 20:20:51 dugas Exp $
 
 
      DESCRIPTION...
       MEOFVEC calcule les valeurs propres et optionnellement, les 
       vecteurs propres droits de la matrice/fichier A ou bien, dans
       le cas symmetrique generalise (KIND=3), du produit matriciel B*A.
       B est alors definie positive et A est symmetrique. Dans les deux
       autres cas, la matrice A est la seule a etre lue.
 
       Les enregistrements des matrices correspondents aux colonnes 
       des matrices. Il doit donc y avoir autant d'echantillons 
       que de variables.
 
 
      FICHIERS DE SORTIE...
       AVAL = Valeurs propres de a en ordre descendant d'importance.
              Dans le cas general (KIND=1), ces valeurs propres sont
              sous formes de paires de conjuguees complexes ou purement
              reelles. Les parties reeles et imaginaires sont alors
              sauvees separement.
       AVEC = Vecteurs propres de a correspondants a aval
 
       NOTE: SI a est un matrice generale, les valeurs propres sont
             complexes. Les parties reeles et imaginaires sont alors
             sauvees separement dans aval (noms = VALR et VALI).
 
 
      PARAMETRES...
        KIND = >=0 , Le programme ne calculera que les vecteurs propres;
         (1)   < 0 , Le programme calculera les vecteurs et valeurs popres;
               0,1 , Si la matrice A est generale (DEFAUT);
               2   ,    "     "     "    symmetrique;
               3   ,    "     "     "    symmetrique
                     et la matrice B est definie positive.
        KIND  != 0 , Les vecteurs propres sont multitplies par la
         (2)         racine carree de leurs valeurs propres a la sortie.
                     AVEC contient alors la matrice des facteurs ("factor
                     loadings") qui est utilisee par EOFROT
                     (SSI KIND(1) = 2).
 
        Les trois parametres suivants n'ont de sens que dans le cas
        SYMMETRIQUE. De plus, notez que la paire A,B a preseance sur
        le parametre C.
 
        A,B  = Limites inferieures et superieures des valeurs
               propres recherchees.
        C    = Fraction de la variance totale exclue.
        D    = 0,     trier les valeurs propres reelles (lorsque
                      KIND=2,3) de la plus grande a la plus petite
               sinon, trier ces valeurs propres selon leur module
 
        Si aucun des parametres A,B ou C n'est specifie, le programme
        tentera de trouver les modes expliquant 99% de la variance.
 
        Dans ce dernier cas, de meme que lorsqu'on utilise le parametre C,
        la matrice doit en plus etre definie positive puisque le programme
        utilise la trace comme mesure de la variance totale.
               
        READ(5,5000) KIND,A,B,C,D
 5000   FORMAT(10X,I5,3E10.0,I5)
 
 
      EXEMPLE DE CARTE DE DIRECTIVE...
 
  MEOFVEC    -1
 
         0123456789012345678901234567890123456789012345678901234567890
         1         2         3         4         5         6         7
 
 
     EXIT CONDITIONS...
        0 NO PROBLEMS
        1 UNABLE TO READ FIRST SET IN FILE A
        2 UNABLE TO READ FIRST SET IN FILE B
       12 ORDER OF A AND B DIFFER
        3 NOT ENOUGH MEMORY
        4 UNABLE TO READ A COLUMN FROM A
       14 UNABLE TO READ A COLUMN FROM B
        5 MORE THAN ONE SIZE IN FILE A
       15 MORE THAN ONE SIZE IN FILE B
        6 UNABLE TO PERFORM TRANSFORMATION
        7 UNABLE TO PERFORM FACTORIZATION
        8 UNABLE TO PERFORM BACK-TRANSFORMATION
        9 BAD EIGENVALUE RANGE SPECIFICATION
       10 ERROR READING DIRECTIVES