DATE     99/10/22
    Opérateur RESO
    --------------

 CHPO1 (CHPO2 ..) = RESO RIG1 (PREC) CHPO3 (CHPO4..)('NOID')('NOUNIL')
                  |              (CHPOF) ('INIB' BLO1 LENTI1)
     (NB_MOD_RIG) |      ('ENSE')        
    Objet :
    _______

    L'opérateur RESO construit une solution, si elle existe, du système
linéaire :   RIG1 CHPO1 = CHPO3 .

    Commentaire :
    _____________

    RIG1  : objet de type RIGIDITE.

    CHPO3 : objet de type CHPOINT.

    CHPO1 : objet de type CHPOINT dont les composantes sont les duales
            de celles de CHPO3 par rapport à RIG1.

    Remarque :
    __________

 1- En présence d'une famille de seconds membres CHPO3, CHPO4 ...,
l'opérateur RESO construit la famille de solutions CHPO1, CHPO2 ..
respectivement associée.

 2- Si RIG1 contient des matrices issues de conditions unilatérales,
RESO appelle la procédure UNILATER pour fournir une solution du système.
Si il y a des matrices de frottement, il faut fournir le champ CHPOF
de forces limite de frottement.

 3- Les mots-clé 'NOID' et 'NOUNIL' sont utiles quand on emploie RESO à
    l'intérieur d'une procédure :

      - 'NOID' permet de résoudre le système avec comme second membre la
        restriction de CHPO3, (CHPO4 ..) à l'espace cible de RIG1.

      - 'NOUNIL' permet de résoudre le système en ignorant le caractère
        éventuellement unilatéral de RIG1.

      - 'INIB' BLO1 LENTI1 permet dans le cas de contact d'indiquer un
        état de contact initial. Voir la procédure UNILATER.

 4- Par "OPTION RESO DIRECTE" ou "OPTION RESO ITERATIVE" 
    on peut choisir soit une resolution par méthode de CROUT 
    ( méthode par défaut) soit une résolution par méthode de gradients
    conjugués avec préconditionnement ( Crout incomplet).

 5- Les champs par points CHPO1(2..) obtenus sont de nature diffuse.

 6- Le mot-clé 'ENSE' indique que au cas ou le système est singulier
    et la singularité est excité, RESO fournira un vecteur du noyau.
    Il y a alors un résultat supplémentaire : NB_MOD_RIG qui est le
    nombre de vecteurs du noyau retournés.

 7- PREC est la prévision de l'opération. Le défaut est 1D-15.    

    Exemple :
    _________

    RIG1 étant la raideur d'une structure , FORC1 un champ de force
s'exerçant sur cette structure, on obtiendra le champ de déplacements
DEP1 en résultant, par l'instruction :

                  DEP1 = RESO  RIG1 FORC1 ;